Changeset 964 for LMDZ4/trunk/libf

Timestamp:

Jun 9, 2008, 11:55:44 AM (16 years ago)

Author:

Laurent Fairhead

Message:

Sorties de certaines variables que quand elles sont bien activees
LF

Location:

LMDZ4/trunk/libf/phylmd

Files:

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

LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

@@ -319 +319 @@
   include "advtrac.h"
   include "clesphys.h"
+  include "thermcell.h"
 
   integer                               :: jjmp1, nqmax
@@ -605 +606 @@
  CALL histdef2d(iff,flag_weakinv, "weakinv","Weak inversion", "-")
  CALL histdef2d(iff,flag_dthmin,"dthmin","dTheta mini", "K/m")
- CALL histdef2d(iff,flag_meantaucld,"meantaucld","ISCCP mean cloud optical thickness","1")
+ if (ok_isccp) then
+   CALL histdef2d(iff,flag_meantaucld,"meantaucld","ISCCP mean cloud optical thickness","1")
+ endif
  CALL histdef2d(iff,flag_rh2m,"rh2m","Relative humidity at 2m", "%" )
  CALL histdef2d(iff,flag_qsat2m,"qsat2m","Saturant humidity at 2m", "%")
@@ -664 +667 @@
  CALL histdef3d(iff,flag_ducon,"ducon","Convection du","m/s2")
  CALL histdef3d(iff,flag_dqcon,"dqcon","Convection dQ", "(kg/kg)/s")
- CALL histdef3d(iff,flag_dtwak,"dtwak","Wake dT","K/s")
- CALL histdef3d(iff,flag_dqwak,"dqwak","Wake dQ","(kg/kg)/s")
+ IF (iflag_wake == 1) THEN
+   CALL histdef3d(iff,flag_dtwak,"dtwak","Wake dT","K/s")
+   CALL histdef3d(iff,flag_dqwak,"dqwak","Wake dQ","(kg/kg)/s")
+ endif
  CALL histdef3d(iff,flag_dtlsc,"dtlsc","Condensation dT", "K/s")
  CALL histdef3d(iff,flag_dtlschr,"dtlschr","Large-scale condensational heating rate","K/s")

LMDZ4/trunk/libf/phylmd/phys_output_write.h

@@ -966 +966 @@
        ENDIF
 
-       IF (iflag_wake.ge.1) THEN
-        IF (flag_dtwak(iff)<=lev_files(iff)) THEN
-        zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)/pdtphys
-        CALL histwrite_phy(nid_files(iff),"dtwak",itau_w,zx_tmp_fi3d)
-        ENDIF
+       IF (iflag_wake == 1) THEN
+         IF (flag_dtwak(iff)<=lev_files(iff)) THEN
+           zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)           &
+     &                                        /pdtphys
+           CALL histwrite_phy(nid_files(iff),"dtwak",itau_w,zx_tmp_fi3d)
+         ENDIF
 
         IF (flag_dqwak(iff)<=lev_files(iff)) THEN