Changeset 1028 for LMDZ4/trunk/libf/phylmd

Timestamp:

Oct 24, 2008, 5:15:32 PM (16 years ago)

Author:

lmdzadmin

Message:

Sorties en plus
FH/IM

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

@@ -276 +276 @@
   integer, dimension(nfiles) , save :: flag_f_th        = (/ 4, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_e_th        = (/ 4, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_w_th        = (/ 4, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_lambda_th        = (/ 4, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_q_th        = (/ 4, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_a_th        = (/ 4, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_d_th        = (/ 4, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_f0_th        = (/ 4, 10, 10, 10 /)
@@ -678 +682 @@
  CALL histdef3d(iff,flag_f_th,"f_th","Thermal plume mass flux", "K/s")
  CALL histdef3d(iff,flag_e_th,"e_th","Thermal plume entrainment", "K/s")
+ CALL histdef3d(iff,flag_w_th,"w_th","Thermal plume vertical velocity", "m/s")
+ CALL histdef3d(iff,flag_lambda_th,"lambda_th","Thermal plume vertical velocity", "m/s")
+ CALL histdef3d(iff,flag_q_th,"q_th","Thermal plume total humidity", "kg/kg")
+ CALL histdef3d(iff,flag_a_th,"a_th","Thermal plume fraction", "")
  CALL histdef3d(iff,flag_d_th,"d_th","Thermal plume detrainment", "K/s")
- CALL histdef2d(iff,flag_d_th,"f0_th","Thermal closure mass flux", "K/s")
- CALL histdef2d(iff,flag_d_th,"zmax_th","Thermal plume height", "K/s")
+ CALL histdef2d(iff,flag_f0_th,"f0_th","Thermal closure mass flux", "K/s")
+ CALL histdef2d(iff,flag_zmax_th,"zmax_th","Thermal plume height", "K/s")
  CALL histdef3d(iff,flag_dqthe,"dqthe","Dry adjust. dQ","(kg/kg)/s")
  CALL histdef3d(iff,flag_dtajs,"dtajs","Dry adjust. dT", "K/s")

LMDZ4/trunk/libf/phylmd/phys_output_write.h

@@ -1003 +1003 @@
         ENDIF
 
-!       IF (flag_d_th(iff)<=lev_files(iff)) THEN
-!       CALL histwrite_phy(nid_files(iff),"d_th",itau_w,detr_therm)
-!       ENDIF
+        IF (flag_w_th(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),"w_th",itau_w,zw2)
+        ENDIF
+
+        IF (flag_q_th(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),"q_th",itau_w,zqasc)
+        ENDIF
+
+        IF (flag_lambda_th(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),"lambda_th",itau_w,lambda_th)
+        ENDIF
+
+        IF (flag_a_th(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),"a_th",itau_w,fraca)
+        ENDIF
+
+       IF (flag_d_th(iff)<=lev_files(iff)) THEN
+       CALL histwrite_phy(nid_files(iff),"d_th",itau_w,detr_therm)
+       ENDIF
+
        ENDIF !iflag_thermals