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

Timestamp:

Jun 19, 2008, 12:25:57 PM (17 years ago)

Author:

lmdzadmin

Message:

Initialisations : concvl, cv3_routines, cva_driver, physiq
Correction bug i0 + ajout tests : cv3p1_closure
Ajout sorties : ale, alp, cin, wape
Ajout variables wake : phyetat0, phyredem
IM

File:

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

LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

@@ -68 +68 @@
   integer, dimension(nfiles),save   :: flag_ndayrain     = (/ 1, 10, 10, 10 /)
   integer, dimension(nfiles),save  :: flag_precip(nfiles)= (/ 1, 1, 1, 1 /)
-  integer,  dimension(nfiles), save :: flag_plul(nfiles) = (/ 1, 1, 10, 1 /)
-
-  integer, dimension(nfiles) , save :: flag_pluc         = (/ 1, 1, 10, 1 /)
+  integer,  dimension(nfiles), save :: flag_plul(nfiles) = (/ 1, 1, 1, 1 /)
+
+  integer, dimension(nfiles) , save :: flag_pluc         = (/ 1, 1, 1, 1 /)
   integer, dimension(nfiles) , save :: flag_snow         = (/ 1, 1, 10, 1 /) 
   integer, dimension(nfiles) , save :: flag_evap         = (/ 1, 1, 10, 1 /)
@@ -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 /)
+  integer, dimension(nfiles) , save :: flag_ale_bl    = (/ 1, 1, 1, 1 /)
+  integer, dimension(nfiles) , save :: flag_alp_bl    = (/ 1, 1, 1, 1 /)
+  integer, dimension(nfiles) , save :: flag_ale_wk    = (/ 1, 1, 1, 1 /)
+  integer, dimension(nfiles) , save :: flag_alp_wk    = (/ 1, 1, 1, 1 /)
+
+  integer, dimension(nfiles) , save :: flag_ale       = (/ 1, 1, 1, 1 /)
+  integer, dimension(nfiles) , save :: flag_alp       = (/ 1, 1, 1, 1 /)
+  integer, dimension(nfiles) , save :: flag_cin       = (/ 1, 1, 1, 1 /)
+  integer, dimension(nfiles) , save :: flag_wape       = (/ 1, 1, 1, 1 /)
+
 
 ! Champs interpolles sur des niveaux de pression ??? a faire correctement
@@ -227 +233 @@
   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_ovap         = (/ 2, 3, 4, 1 /)
   integer, dimension(nfiles) , save :: flag_wvapp        = (/ 2, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_geop         = (/ 2, 3, 10, 1 /)
@@ -604 +610 @@
  CALL histdef2d(iff,flag_alp_wk,"alp_wk","ALP WK","m2/s2") 
 
+ CALL histdef2d(iff,flag_ale,"ale","ALE","m2/s2")
+ CALL histdef2d(iff,flag_alp,"alp","ALP","W/m2")
+ CALL histdef2d(iff,flag_cin,"cin","Convective INhibition","m2/s2")
+ CALL histdef2d(iff,flag_wape,"WAPE","WAPE","m2/s2")
+
  CALL histdef2d(iff,flag_weakinv, "weakinv","Weak inversion", "-")
  CALL histdef2d(iff,flag_dthmin,"dthmin","dTheta mini", "K/m")
@@ -645 +656 @@
  CALL histdef3d(iff,flag_cldtau,"cldtau","Cloud optical thickness","1")
  CALL histdef3d(iff,flag_cldemi,"cldemi","Cloud optical emissivity","1")
- CALL histdef3d(iff,flag_pr_con_l,"pmflxr","Convective precipitation lic"," ")
- CALL histdef3d(iff,flag_pr_con_i,"pmflxs","Convective precipitation ice"," ")
- CALL histdef3d(iff,flag_pr_lsc_l,"prfl","Large scale precipitation lic"," ")
- CALL histdef3d(iff,flag_pr_lsc_i,"psfl","Large scale precipitation ice"," ")
+!IM: bug ?? dimensionnement variables (klon,klev+1) pmflxr, pmflxs, prfl, psfl
+! CALL histdef3d(iff,flag_pr_con_l,"pmflxr","Convective precipitation lic"," ")
+! CALL histdef3d(iff,flag_pr_con_i,"pmflxs","Convective precipitation ice"," ")
+! CALL histdef3d(iff,flag_pr_lsc_l,"prfl","Large scale precipitation lic"," ")
+! CALL histdef3d(iff,flag_pr_lsc_i,"psfl","Large scale precipitation ice"," ")
 
 !FH Sorties pour la couche limite
@@ -778 +790 @@
       end subroutine histdef3d
 
-
 END MODULE phys_output_mod