Changeset 1284

Timestamp:

Dec 17, 2009, 10:37:48 AM (15 years ago)

Author:

idelkadi

Message:

-Rajout de sorties dans les fichiers histoires :

re : Cloud droplet effective radius
fl : Denominator of re

Location:

LMDZ4/trunk/libf/phylmd

Files:

• phys_output_mod.F90 (modified) (2 diffs)
• phys_output_write.h (modified) (1 diff)

LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

@@ -398 +398 @@
   type(ctrl_out),save :: o_pr_lsc_l     = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_lsc_l')
   type(ctrl_out),save :: o_pr_lsc_i     = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_lsc_i')
+  type(ctrl_out),save :: o_re           =ctrl_out((/ 5, 10, 10, 10, 10 /),'re')
+  type(ctrl_out),save :: o_fl           =ctrl_out((/ 5, 10, 10, 10, 10 /),'fl')
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
@@ -1006 +1008 @@
  CALL histdef3d(iff,o_pr_lsc_l%flag,o_pr_lsc_l%name, "Large scale precipitation lic", " ")
  CALL histdef3d(iff,o_pr_lsc_i%flag,o_pr_lsc_i%name, "Large scale precipitation ice", " ")
+!Cloud droplet effective radius
+ CALL histdef3d(iff,o_re%flag,o_re%name, "Cloud droplet effective radius","um")
+ CALL histdef3d(iff,o_fl%flag,o_fl%name, "Denominator of Cloud droplet effective radius"," ")
 !FH Sorties pour la couche limite
      if (iflag_pbl>1) then

LMDZ4/trunk/libf/phylmd/phys_output_write.h

@@ -762 +762 @@
      s        o_pr_lsc_i%name,itau_w,psfl(:,1:klev))
       ENDIF
+
+      IF (o_re%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),o_re%name,itau_w,re)
+      ENDIF
+
+      IF (o_fl%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),o_fl%name,itau_w,fl)
+      ENDIF
+