Changeset 1813 for LMDZ5/trunk/libf/phylmd/physiq.F

Timestamp:

Jul 24, 2013, 1:49:19 PM (12 years ago)

Author:

idelkadi

Message:

1. transform phytrac into a module, in order to pass some variables

(tracer tendencies) to the standard physiq ouput codes.

2. Correct a (big) bug in the call to phytrac.
3. Include w*, and ALEs in the call to phytrac and traclmdz.

physiq.F

• Bug correction in the call of phytrac from the physics u10m,v10m, ustar -> zu10m, zv10m, zustar

phytrac.F90 -> phytrac_mod.F90

• Tranformation of routine phytrac into a module phytrac_mod, in order to tranfer the tracer tendencies from phytrac to

phys_output...

• Inclusion of w*, Ale bl/wake in the call to phytrac and traclmdz

(to be used for dust emmission)

by respectively, vertical diffusion, thermal plumes and convection

• desactivation of ini_histrac.h and write_histrac.h
• USE phys_output_mod removed since it was creating a circular

dependency

between phytrac_mod and phys_output_mod.
So the automatic computation of ecrit_tra is desactivated

ini_histrac.h and write_histrac.h

Descactivated in phytrac but kept for backard compatibility
couchelimite -> iflag_vdf_trac>0

phys_output_ctrlout_mod.F90

New variables : o_dtr_vdf, o_dtr_the ... for output of tracer tendencies

phys_output_mod.F90

Default definition for these new output variables.

phys_output_write_F90.h

disapears, included directly in phys_output_write_mod.F90

phys_output_write_mod.F90

writing of the tracer tendencies

phys_state_var_mod.F90

New declaration (wstar)

traclmdz_mod.F90

• Inclusion of w*, Ale bl/wake in the call to phytrac and traclmdz

(to be used for dust emmission)

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -50 +50 @@
 #endif
       USE indice_sol_mod
+      USE phytrac_mod, ONLY : phytrac
 
 !IM stations CFMIP
@@ -1523 +1524 @@
 c$OMP BARRIER
 
+#undef histISCCP
 #ifdef histISCCP
 #include "ini_histISCCP.h"
 #endif
 
+#undef histNMC
 #ifdef histNMC
 #include "ini_histhfNMC.h"
@@ -3637 +3640 @@
        END IF
 
+      wstar=0. ! FH 2013/07/22 a corriger des que dans pbl_surface_mod
       call phytrac (
      I     itap,     days_elapsed+1,    jH_cur,   debut,
@@ -3642 +3646 @@
      I     paprs,    pplay,     pmfu,     pmfd,
      I     pen_u,    pde_u,     pen_d,    pde_d,
-     I     cdragh,   coefh(:,:,is_ave),     fm_therm, entr_therm,
+     I     cdragh,   coefh(:,:,is_ave),   fm_therm, entr_therm,
      I     u1,       v1,        ftsol,    pctsrf,
-     I     ustar,     u10m,      v10m,
+     I     zustar,   zu10m,     zv10m,
+     I     wstar(:,is_ave),    ale_bl_stat,         ale_wake,
      I     rlat,     rlon,
      I     frac_impa,frac_nucl, beta_prec_fisrt,beta_prec,