Changeset 180 for LMDZ.3.3

Timestamp:

Mar 27, 2001, 11:09:52 AM (23 years ago)

Author:

lmdzadmin

Message:

Interversion dans les sorties netcdf des champs venant du coupleur
MAFO/LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90 (modified) (7 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

@@ -1205 +1205 @@
       CALL histwrite(nidcs,cl_read(1),itime,read_sst,iim*(jjm+1),ndexcs)
       CALL histwrite(nidcs,cl_read(2),itime,read_sic,iim*(jjm+1),ndexcs)
-      CALL histwrite(nidcs,cl_read(3),itime,read_sit,iim*(jjm+1),ndexcs)
-      CALL histwrite(nidcs,cl_read(4),itime,read_alb_sic,iim*(jjm+1),ndexcs)
+      CALL histwrite(nidcs,cl_read(3),itime,read_alb_sic,iim*(jjm+1),ndexcs)
+      CALL histwrite(nidcs,cl_read(4),itime,read_sit,iim*(jjm+1),ndexcs)
       CALL histsync(nidcs)
 ! pas utile      IF (npas-itime.LT.nexca )CALL histclo(nidcs)
@@ -1271 +1271 @@
 !
 ! allocation memoire
-    if (nisurf == is_oce) then
+    if (nisurf == is_oce .and. (.not. cumul) ) then
       sum_error = 0
       allocate(tmp_sols(iim,jjm+1,2), stat=error); sum_error = sum_error + error
@@ -1947 +1947 @@
 ! Traitement neige et humidite du sol
 !
-!    if (first .and. nisurf == is_ter) then
-!      do i = 1, knon
-!        write(67,*)i, snow(i), precip_snow(i), evap(i)
-!      enddo
-!    endif
-!    if (second .and. nisurf == is_ter) then
-!      do i = 1, knon
-!        write(77,*)i, snow(i), precip_snow(i), evap(i)
-!      enddo
-!    endif
     if (nisurf == is_oce) then
       snow = 0.
@@ -2004 +1994 @@
      & * p1lay(i)/(RD*t1lay(i))
 
-!    if (first .and. nisurf == is_ter) then
-!      write(43,*) &
-!     &i,zx_coef(i),coef1lay(i),u1lay(i),v1lay(i),p1lay(i),t1lay(i)
-!    endif
-!    if (second .and. nisurf == is_ter) then
-!      write(53,*) &
-!     &i,zx_coef(i),coef1lay(i),u1lay(i),v1lay(i),p1lay(i),t1lay(i)
-!    endif
-    
   ENDDO
 
@@ -2042 +2023 @@
     zx_nh(i) = - (zx_k1(i) * RCPD * zx_pkh(i))/ zx_oh(i)
 
-!    if (first .and. nisurf == is_ter) then
-!      write(41,*) &
-!     &            i,zx_k1(i),petAcoef(i),petBcoef(i),&
-!     & zx_oh(i),zx_mh(i),zx_nh(i)
-!    endif
-
-!    if (second .and. nisurf == is_ter) then
-!      write(51,*) &
-!     &            i,zx_k1(i),petAcoef(i),petBcoef(i),&
-!     & zx_oh(i),zx_mh(i),zx_nh(i)
-!    endif
 
 ! Tsurface
@@ -2081 +2051 @@
     dflux_l(i) = (zx_sl(i) * zx_nq(i))
 
-!    if (first .and. nisurf == is_ter) then
-!      write(42,*) &
-!     &            i,tsurf_new(i),d_ts(i),zx_h_ts(i),zx_q_0(i),&
-!     & evap(i),fluxsens(i)
-!    endif
-
-!    if (second .and. nisurf == is_ter) then
-!      write(52,*) &
-!     &            i,tsurf_new(i),d_ts(i),zx_h_ts(i),zx_q_0(i),&            
-!     & evap(i),fluxsens(i) 
-!    endif
 !
 ! en cas de fonte de neige
@@ -2108 +2067 @@
     qsol(i) = min(qsol(i), max_eau_sol) 
   ENDDO
-
-!  if (nisurf == is_oce .and. second) then
-!    second =.f.
-!  endif
-
-!  if (nisurf == is_oce .and. first) then
-!    first = .f.
-!    second = .t.
-!  endif
 
   END SUBROUTINE calcul_fluxs