Changeset 945 for LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

Timestamp:

Apr 11, 2008, 3:53:10 PM (16 years ago)

Author:

lmdzadmin

Message:

Ajout variables thermiques/wake/autres FH
IM

File:

• LMDZ4/trunk/libf/phylmd/phys_output_mod.F90 (modified) (8 diffs)

LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

@@ -162 +162 @@
   integer, dimension(nfiles) , save :: flag_slab_bils    = (/ 1, 1, 10, 10 /)
 
+  integer, dimension(nfiles) , save :: flag_ale_bl    = (/ 1, 1, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_alp_bl    = (/ 1, 1, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_ale_wk    = (/ 1, 1, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_alp_wk    = (/ 1, 1, 10, 10 /)
+
 ! Champs interpolles sur des niveaux de pression ??? a faire correctement
 ! if=1 on ecrit u v w phi sur 850 700 500 200 au niv 1
@@ -221 +226 @@
   integer, dimension(nfiles) , save :: flag_iwcon        = (/ 2, 5, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_temp         = (/ 2, 3, 4, 1 /)
+  integer, dimension(nfiles) , save :: flag_theta         = (/ 2, 3, 4, 1 /)
   integer, dimension(nfiles) , save :: flag_ovap         = (/ 2, 3, 4, 10 /)
   integer, dimension(nfiles) , save :: flag_wvapp        = (/ 2, 10, 10, 10 /)
@@ -257 +263 @@
   integer, dimension(nfiles) , save :: flag_ducon        = (/ 4, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dqcon        = (/ 4, 5, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_dtwak        = (/ 4, 5, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_dqwak        = (/ 4, 5, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dtlsc        = (/ 4, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dtlschr      = (/ 4, 10, 10, 10 /)
@@ -267 +275 @@
   integer, dimension(nfiles) , save :: flag_ratqs        = (/ 4, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dtthe        = (/ 4, 10, 10, 10 /)
+  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_d_th        = (/ 4, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_f0_th        = (/ 4, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_zmax_th        = (/ 4, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dqthe        = (/ 4, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dtajs        = (/ 4, 10, 10, 10 /)
@@ -585 +598 @@
      ENDIF
 
+ CALL histdef2d(iff,flag_ale_bl,"ale_bl","ALE BL","m2/s2") 
+ CALL histdef2d(iff,flag_alp_bl,"alp_bl","ALP BL","m2/s2") 
+ CALL histdef2d(iff,flag_ale_wk,"ale_wk","ALE WK","m2/s2") 
+ CALL histdef2d(iff,flag_alp_wk,"alp_wk","ALP WK","m2/s2") 
+
  CALL histdef2d(iff,flag_weakinv, "weakinv","Weak inversion", "-")
  CALL histdef2d(iff,flag_dthmin,"dthmin","dTheta mini", "K/m")
@@ -608 +626 @@
  CALL histdef3d(iff,flag_iwcon,"iwcon","Cloud ice water content","kg/kg")
  CALL histdef3d(iff,flag_temp,"temp","Air temperature","K" )
+ CALL histdef3d(iff,flag_theta,"theta","Potential air temperature","K" )
  CALL histdef3d(iff,flag_ovap,"ovap","Specific humidity","kg/kg" )
  CALL histdef3d(iff,flag_geop,"geop","Geopotential height","m2/s2")
@@ -645 +664 @@
  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")
  CALL histdef3d(iff,flag_dtlsc,"dtlsc","Condensation dT", "K/s")
  CALL histdef3d(iff,flag_dtlschr,"dtlschr","Large-scale condensational heating rate","K/s")
@@ -655 +676 @@
  CALL histdef3d(iff,flag_ratqs,"ratqs", "RATQS"," ")
  CALL histdef3d(iff,flag_dtthe,"dtthe","Dry adjust. dT", "K/s")
+ 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_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 histdef3d(iff,flag_dqthe,"dqthe","Dry adjust. dQ","(kg/kg)/s")
  CALL histdef3d(iff,flag_dtajs,"dtajs","Dry adjust. dT", "K/s")