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

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/physiq.F (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -722 +722 @@
 
       REAL bils(klon) ! bilan de chaleur au sol
+
       REAL wfbilo(klon,nbsrf) ! bilan d'eau, pour chaque
 C                             ! type de sous-surface et pondere par la fraction
@@ -751 +752 @@
       SAVE lmt_pas                ! frequence de mise a jour
 c$OMP THREADPRIVATE(lmt_pas) 
-      real zmasse(klon, llm) 
+      real zmasse(klon, llm),exner(klon, llm) 
 C     (column-density of mass of air in a cell, in kg m-2)
       real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
@@ -3748 +3749 @@
 c Accumuler les variables a stocker dans les fichiers histoire:
 c
-c+jld ec_conser
-      DO k = 1, klev
-      DO i = 1, klon
-        ZRCPD = RCPD*(1.0+RVTMP2*q_seri(i,k))
-        d_t_ec(i,k)=0.5/ZRCPD
-     $      *(u(i,k)**2+v(i,k)**2-u_seri(i,k)**2-v_seri(i,k)**2)
-      ENDDO
-      ENDDO
-
-      DO k = 1, klev
-      DO i = 1, klon
-        t_seri(i,k)=t_seri(i,k)+d_t_ec(i,k)
-        d_t_ec(i,k) = d_t_ec(i,k)/dtime
-       END DO 
-      END DO 
-c-jld ec_conser
+
+!================================================================
+! Conversion of kinetic and potential energy into heat, for
+! parameterisation of subgrid-scale motions
+!================================================================
+
+      d_t_ec(:,:)=0.
+      forall (k=1: llm) exner(:, k) = (pplay(:, k)/paprs(:,1))**RKAPPA
+      CALL ener_conserv(klon,klev,pdtphys,u,v,t,qx(:,:,ivap),
+     s        u_seri,v_seri,t_seri,q_seri,
+     s        zmasse,exner,d_t_ec)
+      t_seri(:,:)=t_seri(:,:)+d_t_ec(:,:)
+
 cIM
       IF (ip_ebil_phy.ge.1) THEN 
@@ -3833 +3831 @@
       END IF
 
-c=============================================================
+
 c
 c Convertir les incrementations en tendances