Changeset 458 for LMDZ.3.3/branches

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

Location:

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

Files:

• clmain.F (modified) (3 diffs)
• interface_surf.F90 (modified) (8 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F

@@ -236 +236 @@
       character (len = 20) :: modname = 'clmain'
       LOGICAL check
-      PARAMETER (check=.true.)
+      PARAMETER (check=.false.)
 C
       if (check) THEN
@@ -872 +872 @@
       character (len = 20) :: modname = 'Debut clqh'
       LOGICAL check
-      PARAMETER (check=.true.)
+      PARAMETER (check=.false.)
 C
       if (check) THEN
@@ -1388 +1388 @@
       ENDDO
 
-      PRINT*,' isommet=',isommet,' knon=',knon
+c$$$      PRINT*,' isommet=',isommet,' knon=',knon
 
       DO k = 2, isommet

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.