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

Timestamp:

Mar 27, 2003, 2:44:17 PM (22 years ago)

Author:

lmdzadmin

Message:

Un peu de menage pour avoir moins de print dans les sorties
LF

File:

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

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

@@ -206 +206 @@
   integer, save        :: error
   integer              :: ii, index
-  logical,save              :: check = .true.
+  logical,save              :: check = .false.
   real, dimension(klon):: cal, beta, dif_grnd, capsol
 !!$PB  real, parameter      :: calice=1.0/(5.1444e+06*0.15), tau_gl=86400.*5.
@@ -218 +218 @@
   real, DIMENSION(klon):: zfra
   logical              :: cumul = .false.
-!IM cf LF
-  integer :: i
+  INTEGER,dimension(1) :: iloc
+  INTEGER                 :: isize
+  real, dimension(klon):: fder_prev
 
   if (check) write(*,*) 'Entree ', modname
@@ -407 +408 @@
       cumul = .false.
 
+      iloc = maxloc(fder(1:klon))
+      if (check) then
+        if (fder(iloc(1))> 0.) then
+          WRITE(*,*)'**** Debug fder ****'
+          WRITE(*,*)'max fder(',iloc(1),') = ',fder(iloc(1))
+        endif
+      endif
+!!$
+!!$      where(fder.gt.0.) 
+!!$        fder = 0.
+!!$      endwhere
+
       call interfoce(itime, dtime, cumul, &
       & klon, iim, jjm, nisurf, pctsrf, knon, knindex, rlon, rlat, &
@@ -437 +450 @@
      &   petAcoef, peqAcoef, petBcoef, peqBcoef, &
      &   tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
-    
-    fder = fder + dflux_s + dflux_l
+
+    fder_prev = fder    
+    fder = fder_prev + dflux_s + dflux_l
+
+      iloc = maxloc(fder(1:klon))
+        if (check.and.fder(iloc(1))> 0.) then
+          WRITE(*,*)'**** Debug fder****'
+          WRITE(*,*)'max fder(',iloc(1),') = ',fder(iloc(1))
+          WRITE(*,*)'fder_prev, dflux_s, dflux_l',fder_prev(iloc(1)), &
+     &                        dflux_s(iloc(1)), dflux_l(iloc(1))
+        endif
+!!$
+!!$      where(fder.gt.0.) 
+!!$        fder = 0.
+!!$      endwhere
 
 !
@@ -488 +514 @@
 
       cumul =.false.
+
+      iloc = maxloc(fder(1:klon))
+      if (check.and.fder(iloc(1))> 0.) then
+        WRITE(*,*)'**** Debug fder ****'
+        WRITE(*,*)'max fder(',iloc(1),') = ',fder(iloc(1))
+      endif
+!!$
+!!$      where(fder.gt.0.) 
+!!$        fder = 0.
+!!$      endwhere
 
       call interfoce(itime, dtime, cumul, &
@@ -573 +609 @@
     ENDIF
 
-    fder = fder + dflux_s + dflux_l
+    fder_prev = fder    
+    fder = fder_prev + dflux_s + dflux_l
+
+      iloc = maxloc(fder(1:klon))
+      if (check.and.fder(iloc(1))> 0.) then
+        WRITE(*,*)'**** Debug fder ****'
+        WRITE(*,*)'max fder(',iloc(1),') = ',fder(iloc(1))
+        WRITE(*,*)'fder_prev, dflux_s, dflux_l',fder_prev(iloc(1)), &
+     &                        dflux_s(iloc(1)), dflux_l(iloc(1))
+      endif
+!!$      where(fder.gt.0.) 
+!!$        fder = 0.
+!!$      endwhere
 
 !
@@ -2000 +2048 @@
   character (len = 20),save :: fich ='limit.nc'
   logical,save     :: newlmt = .false.
-  logical,save     :: check = .true.
+  logical,save     :: check = .false.
 ! Champs lus dans le fichier de CL
   real, allocatable , save, dimension(:) :: alb_lu, rug_lu
@@ -2182 +2230 @@
   REAL, parameter :: chasno = 3.334E+05/(2.3867E+06*0.15)
 !
-  logical, save         :: check = .true.
+  logical, save         :: check = .false.
   character (len = 20)  :: modname = 'calcul_fluxs'
   logical, save         :: fonte_neige = .false.