Changeset 2703 for LMDZ5/trunk/libf/phylmd

Timestamp:

Nov 17, 2016, 11:25:51 AM (8 years ago)

Author:

musat

Message:

Ajout iflag_oldbug_fisrtilp pour tester l'ancien bug sur fisrtilp
Par defaut iflag_oldbug_fisrtilp=0 (sans bug)

File:

• LMDZ5/trunk/libf/phylmd/fisrtilp.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/fisrtilp.F90

@@ -16 +16 @@
   USE print_control_mod, ONLY: prt_level, lunout
   USE cloudth_mod 
+  USE ioipsl_getin_p_mod, ONLY : getin_p
   IMPLICIT none
   !======================================================================
@@ -146 +147 @@
   !$OMP THREADPRIVATE(appel1er)
   !
+! iflag_oldbug_fisrtilp=0 enleve le BUG par JYG : tglace_min -> tglace_max
+! iflag_oldbug_fisrtilp=1 ajoute le BUG
+  INTEGER,SAVE :: iflag_oldbug_fisrtilp=0 !=0 sans bug
+  !$OMP THREADPRIVATE(iflag_oldbug_fisrtilp)
   !---------------------------------------------------------------
   !
@@ -189 +194 @@
   if (prt_level>9)write(lunout,*)'NUAGES4 A. JAM'
   IF (appel1er) THEN
+     CALL getin_p('iflag_oldbug_fisrtilp',iflag_oldbug_fisrtilp)
+     write(lunout,*)' iflag_oldbug_fisrtilp =',iflag_oldbug_fisrtilp 
      !
      WRITE(lunout,*) 'fisrtilp, ninter:', ninter
@@ -631 +638 @@
                     zdelta = MAX(0.,SIGN(1.,t_glace_min_old-Tbef(i)))
                     else if (iflag_t_glace.ge.1) then
-                    zdelta = MAX(0.,SIGN(1.,t_glace_max-Tbef(i)))
-! BUG corrige par JYG   zdelta = MAX(0.,SIGN(1.,t_glace_min-Tbef(i)))
+                       if (iflag_oldbug_fisrtilp.EQ.0) then
+                          zdelta = MAX(0.,SIGN(1.,t_glace_max-Tbef(i)))
+                       else
+!avec bug : zdelta = MAX(0.,SIGN(1.,t_glace_min-Tbef(i)))
+                          zdelta = MAX(0.,SIGN(1.,t_glace_min-Tbef(i)))
+                       endif
                     endif
                  endif