Changeset 1813 for LMDZ5/trunk/libf/phylmd/traclmdz_mod.F90

Timestamp:

Jul 24, 2013, 1:49:19 PM (11 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/traclmdz_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/traclmdz_mod.F90

@@ -337 +337 @@
   SUBROUTINE traclmdz(nstep, julien, gmtime, pdtphys, t_seri, paprs, pplay, &
        cdragh,  coefh, yu1, yv1, ftsol, pctsrf, xlat, xlon, couchelimite, sh, &
-       rh, pphi, ustar, zu10m, zv10m, &
+       rh, pphi, ustar, wstar, ale_bl, ale_wake,  zu10m, zv10m, &
 !!          tr_seri, source, solsym, d_tr_cl, zmasse)                      !RomP
           tr_seri, source, solsym, d_tr_cl,d_tr_dec, zmasse)               !RomP
@@ -385 +385 @@
     REAL,DIMENSION(klon,klev),INTENT(IN) :: pphi    ! geopotentie
     REAL,DIMENSION(klon),INTENT(IN)      :: ustar   ! ustar (m/s)
+    REAL,DIMENSION(klon),INTENT(IN)      :: wstar,ale_bl,ale_wake   ! wstar (m/s) and Avail. Lifti. Energ.
     REAL,DIMENSION(klon),INTENT(IN)      :: zu10m   ! vent zonal 10m (m/s)
     REAL,DIMENSION(klon),INTENT(IN)      :: zv10m   ! vent zonal 10m (m/s)