Changeset 1813 for LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_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/phys_output_ctrlout_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1049 +1049 @@
   TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac(:)
   TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac_cum(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_vdf(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_the(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_con(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_lessi_impa(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_lessi_nucl(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_insc(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_bcscav(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_evapls(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_ls(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_trsp(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_sscav(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_sat(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_uscav(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_dry(:)
+
   TYPE(ctrl_out), SAVE :: o_rsu = ctrl_out((/ 4, 10, 10, 10, 10, 10 /), &
         'rsu', 'SW upward radiation', 'W m-2', (/ ('', i=1, 6) /))