Changeset 4575

Timestamp:

Jun 14, 2023, 5:06:15 PM (11 months ago)

Author:

Laurent Fairhead

Message:

Sortir les profils de flux sensible et d'eau dus à la diffusion turbulente (déja calculés auparavant)
FC

Location:

LMDZ6/trunk

Files:

• DefLists/field_def_lmdz.xml (modified) (1 diff)
• libf/phylmd/phys_local_var_mod.F90 (modified) (3 diffs, 1 prop)
• libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
• libf/phylmd/phys_output_write_mod.F90 (modified) (3 diffs)
• libf/phylmd/physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/DefLists/field_def_lmdz.xml

@@ -773 +773 @@
         <field id="qcontr2"     long_name="Contrail qcontr2"     unit="kg/kg" />
         <field id="Tcontr"     long_name="Contrail Tcontr"     unit="K" />
+        <field id="fluxt"     long_name="flux h"     unit="W/m2" />
+        <field id="fluxq"     long_name="flux q"     unit="-" />
     </field_group>
 

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -280 +280 @@
 !$OMP THREADPRIVATE(sens, flwp, fiwp)
 !!
+!FC 
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zxfluxt, zxfluxq
+!$OMP THREADPRIVATE(zxfluxt, zxfluxq)
+!FC
 !!         Wake variables
       REAL,ALLOCATABLE,SAVE,DIMENSION(:)            :: alp_wake
@@ -855 +859 @@
       ALLOCATE(wfevap(klon, nbsrf), wfrain(klon,nbsrf), wfsnow(klon, nbsrf))
       ALLOCATE(evap_pot(klon, nbsrf))
+! FC 
+      ALLOCATE(zxfluxq(klon,klev),zxfluxt(klon,klev))
 !
 !  Deep convective variables used in phytrac
@@ -1150 +1156 @@
       DEALLOCATE(alp_bl_stat, n2, s2)
       DEALLOCATE(proba_notrig, random_notrig)
+!FC
+      DEALLOCATE(zxfluxq,zxfluxt)
 
       DEALLOCATE(dnwd0)

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1551 +1551 @@
   TYPE(ctrl_out), SAVE :: o_stratomask = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'stratomask', 'Stratospheric fraction', '1', (/ ('', i=1, 10) /))
+!FC
+  TYPE(ctrl_out), SAVE :: o_zxfluxt = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'fluxt', 'flux h ', 'W/m2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_zxfluxq = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'fluxq', 'flux q ', 'kg/(s*m2)', (/ ('', i=1, 10) /))
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -223 +223 @@
          o_flx_co2_ff, o_flx_co2_bb, & 
          o_delta_sst, o_delta_sal, o_ds_ns, o_dt_ns, o_dter, o_dser, o_tkt, &
-         o_tks, o_taur, o_sss
+         o_tks, o_taur, o_sss, &
+!FC
+         o_zxfluxt,o_zxfluxq
 
 #ifdef CPP_StratAer
@@ -357 +359 @@
 #endif
          ep, epmax_diag, &  ! epmax_cape
-         p_tropopause, t_tropopause, z_tropopause
+         p_tropopause, t_tropopause, z_tropopause, &
+         zxfluxt,zxfluxq
+!FC
 
 #ifdef CPP_StratAer 
@@ -1800 +1804 @@
        CALL histwrite_phy(o_ovap, q_seri)
        CALL histwrite_phy(o_oliq, ql_seri)
+       !FC
+       CALL histwrite_phy(o_zxfluxt, zxfluxt)
+       CALL histwrite_phy(o_zxfluxq, zx_tmp_fi3d)
+       !FC
 
        IF (vars_defined) zx_tmp_fi3d = ql_seri+qs_seri

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -323 +323 @@
        beta_prec,  &
        rneb,  &
-       zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic
+       zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic, &
+       zxfluxt,zxfluxq 
        !
       USE output_physiqex_mod, ONLY: output_physiqex
@@ -844 +845 @@
     REAL fluxqbs(klon,klev, nbsrf)   ! flux turbulent de neige soufflee
     !
-    REAL zxfluxt(klon, klev)
-    REAL zxfluxq(klon, klev)
+    !FC    REAL zxfluxt(klon, klev)
+    !FC    REAL zxfluxq(klon, klev)
     REAL zxfluxqbs(klon,klev)
     REAL zxfluxu(klon, klev)