Changeset 1604 for LMDZ5/trunk/libf/phylmd/coef_diff_turb_mod.F90

Timestamp:

Dec 22, 2011, 12:26:02 PM (12 years ago)

Author:

lguez

Message:

Removed two "fi" with no corresponding "if" in "makdim".

In procedure "coefkz", when calculating cloud fraction "zfr", "zq" can
be zero at high altitude. Avoid division by zero, set "zfr" to 0 when
"zq" is 0.

In "physiq", allocatable arrays "tabijgcm", "longcm", "latgcm",
"igcm", "jgcm" cannot be arguments of "phys_output_open" if they have
not been allocated. Allocate them with zero size when
'npCFMIP_param.data' cannot be opened.

File:

• LMDZ5/trunk/libf/phylmd/coef_diff_turb_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/coef_diff_turb_mod.F90

@@ -389 +389 @@
 !           calculer la fraction nuageuse (processus humide):
 !
-          zfr = (zq+ratqs*zq-zqs) / (2.0*ratqs*zq)
+          if (zq /= 0.) then
+             zfr = (zq+ratqs*zq-zqs) / (2.0*ratqs*zq)
+          else
+             zfr = 0.
+          end if
           zfr = MAX(0.0,MIN(1.0,zfr))
           IF (.NOT.richum) zfr = 0.0