Changeset 1753 for LMDZ5/trunk/libf/phylmd/phys_output_write.h

Timestamp:

May 13, 2013, 2:53:28 PM (11 years ago)

Author:

idelkadi

Message:

Concerns energy conservation.

---

The source terms of the TKE prognostic equation are diagnosed
from the tendencies (du, dv, dT) associated with subrgrid scale
motions and treated as an additional heat source.
Controled by a new key, iflag_ener_conserv :

0 no conservation

-1 old adhoc correction for kinetic E only (used for CMIP5)

1 conservation
101 conversion from kinetic to heat only
110 conversion from potential to heat only

iflag_ener_conserv=-1 kept as default value for a test period
iflag_ener_conserv=1 is the advisable value
Concerns clesphys.h, and conf_phys.F90
New routine : ener_conserv.F90, called by physic.
New outputs :
bils_ec, contribution to the energy budget of the column of the

additional heat source (in W/m2)

bils_kinetic : change kinetic energy of the column in physics (W/m2)
bils_enthalp : idem for the total column enthalphy
bils_latent : idem for latent heat
Modified files : clesphys.h, conf_phys_m.F90, physiq.F,

phys_output_mod.F90, phys_output_var_mod.F90, phys_output_write.h,
ener_conserv.F90

File:

• LMDZ5/trunk/libf/phylmd/phys_output_write.h (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_write.h

@@ -358 +358 @@
       CALL histwrite_phy(nid_files(iff),clef_stations(iff),
      $o_bils%name,itau_w,bils)
+       ENDIF
+
+       IF (o_bils_ec%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),clef_stations(iff),
+     $o_bils_ec%name,itau_w,bils_ec)
+       ENDIF
+
+       IF (o_bils_kinetic%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),clef_stations(iff),
+     $o_bils_kinetic%name,itau_w,bils_kinetic)
+       ENDIF
+
+       IF (o_bils_latent%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),clef_stations(iff),
+     $o_bils_latent%name,itau_w,bils_latent)
+       ENDIF
+
+       IF (o_bils_enthalp%flag(iff)<=lev_files(iff)) THEN
+      CALL histwrite_phy(nid_files(iff),clef_stations(iff),
+     $o_bils_enthalp%name,itau_w,bils_enthalp)
        ENDIF
 
@@ -2021 +2041 @@
 
        IF (o_dtec%flag(iff)<=lev_files(iff)) THEN
-      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
+      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)/pdtphys
       CALL histwrite_phy(nid_files(iff),clef_stations(iff),
      $o_dtec%name,itau_w,zx_tmp_fi3d)