Changeset 1813 for LMDZ5/trunk/libf/phylmd/phys_output_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_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -1 +1 @@
 ! $Id$
 !
+
+MODULE phys_output_mod 
+  USE indice_sol_mod
+  USE phys_output_var_mod
+  USE aero_mod, only : naero_spc,name_aero
+  USE phys_output_write_mod, ONLY : phys_output_write
+
 ! Abderrahmane 12 2007
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -9 +16 @@
 ! histins.nc : valeurs instantanees
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
-MODULE phys_output_mod 
-  USE indice_sol_mod
-  USE phys_output_var_mod
-  USE aero_mod, only : naero_spc,name_aero
-
-  IMPLICIT NONE
 
 CONTAINS
@@ -41 +41 @@
     USE mod_phys_lmdz_para
     USE aero_mod, only : naero_spc,name_aero
-    USE phys_output_write_mod
     USE phys_output_ctrlout_mod
 
@@ -124 +123 @@
     IF (.NOT. ALLOCATED(o_trac)) ALLOCATE(o_trac(nqtot))
     IF (.NOT. ALLOCATED(o_trac_cum)) ALLOCATE(o_trac_cum(nqtot))
+    ALLOCATE(o_dtr_the(nqtot),o_dtr_con(nqtot),o_dtr_lessi_impa(nqtot))
+    ALLOCATE(o_dtr_lessi_nucl(nqtot),o_dtr_insc(nqtot),o_dtr_bcscav(nqtot))
+    ALLOCATE(o_dtr_evapls(nqtot),o_dtr_ls(nqtot),o_dtr_trsp(nqtot))
+    ALLOCATE(o_dtr_sscav(nqtot),o_dtr_sat(nqtot),o_dtr_uscav(nqtot))
+    ALLOCATE(o_dtr_dry(nqtot),o_dtr_vdf(nqtot))
+
 
     levmax = (/ klev, klev, klev, klev, klev, klev /)
@@ -330 +335 @@
             o_trac(iq-2) = ctrl_out((/ 4, 5, 1, 1, 1, 10 /),tname(iiq),'Tracer '//ttext(iiq), "-",&
                   (/ '', '', '', '', '', '' /))
+
+            o_dtr_vdf(iq-2) = ctrl_out((/ 5, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_vdf' &
+               ,'Tendance tracer '//ttext(iiq), "-" , (/ '', '', '', '', '', '' /))
+
+            o_dtr_the(iq-2) = ctrl_out((/ 5, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_the' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_con(iq-2) = ctrl_out((/ 5, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_con' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_lessi_impa(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_lessi_impa' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_lessi_nucl(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_lessi_nucl' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_insc(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_insc' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_bcscav(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_bcscav' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_evapls(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_evapls' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_ls(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_ls' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_trsp(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_trsp' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_sscav(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_sscav' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_sat(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_sat' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_uscav(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'d'//trim(tname(iq))//'_uscav' &
+               ,'Tendance tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
+            o_dtr_dry(iq-2) = ctrl_out((/ 7, 7, 7, 7, 10, 10 /),'cum'//'d'//trim(tname(iq))//'_dry' &
+               ,'tracer tendency dry deposition'//ttext(iiq), "-", (/ '', '', '', '', '', '' /) )
+
             o_trac_cum(iq-2) = ctrl_out((/ 3, 4, 10, 10, 10, 10 /),'cum'//tname(iiq),&
                   'Cumulated tracer '//ttext(iiq), "-", (/ '', '', '', '', '', '' /))