Changeset 2124 for LMDZ5/trunk/libf

Timestamp:

Sep 19, 2014, 7:51:02 PM (10 years ago)

Author:

lguez

Message:

In guide_main, there is the line:

f_out=((MOD(itau,iguide_sav).EQ.0).AND.guide_sav)

The Fortran standard says that the result of mod(a, p) with p = 0 is
processor dependent. With ifort (on Ada at IDRIS), it produces a
run-time error. This is now fixed. The user must still choose
iguide_sav=0 in order to write only once.

Location:

LMDZ5/trunk/libf

Files:

• dyn3d/guide_mod.F90 (modified) (1 diff)
• dyn3dmem/guide_loc_mod.F90 (modified) (1 diff)
• dyn3dpar/guide_p_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/dyn3d/guide_mod.F90

@@ -109 +109 @@
     ! frequences f>0: fx/jour; f<0: tous les f jours; f=0: 1 seule fois.
     IF (iguide_sav.GT.0) THEN
-        iguide_sav=day_step/iguide_sav
-    ELSE
-        iguide_sav=day_step*iguide_sav
+       iguide_sav=day_step/iguide_sav
+    ELSE if (iguide_sav == 0) then
+       iguide_sav = huge(0)
+    else
+       iguide_sav=day_step*iguide_sav
     ENDIF
 

LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90

@@ -114 +114 @@
     ! frequences f>0: fx/jour; f<0: tous les f jours; f=0: 1 seule fois.
     IF (iguide_sav.GT.0) THEN
-        iguide_sav=day_step/iguide_sav
+       iguide_sav=day_step/iguide_sav
+    ELSE if (iguide_sav == 0) then
+       iguide_sav = huge(0)
     ELSE
-        iguide_sav=day_step*iguide_sav
+       iguide_sav=day_step*iguide_sav
     ENDIF
 

LMDZ5/trunk/libf/dyn3dpar/guide_p_mod.F90

@@ -113 +113 @@
     ! frequences f>0: fx/jour; f<0: tous les f jours; f=0: 1 seule fois.
     IF (iguide_sav.GT.0) THEN
-        iguide_sav=day_step/iguide_sav
+       iguide_sav=day_step/iguide_sav
+    ELSE if (iguide_sav == 0) then
+       iguide_sav = huge(0)
     ELSE
-        iguide_sav=day_step*iguide_sav
+       iguide_sav=day_step*iguide_sav
     ENDIF