Changeset 2042

Timestamp:

May 14, 2014, 11:18:03 AM (10 years ago)

Author:

fhourdin

Message:

Ajout d'un diagnostique pour la conservation de l'energie.
Added diagnostic for energie conservation.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• ener_conserv.F90 (modified) (6 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_var_mod.F90 (modified) (4 diffs)
• phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/ener_conserv.F90

@@ -20 +20 @@
 ! From module
 USE phys_local_var_mod, ONLY : d_u_vdf,d_v_vdf,d_t_vdf,d_u_ajs,d_v_ajs,d_t_ajs,d_u_con,d_v_con,d_t_con,d_t_diss
-USE phys_output_var_mod, ONLY : bils_ec,bils_tke,bils_kinetic,bils_enthalp,bils_latent,bils_diss
+USE phys_local_var_mod, ONLY : d_t_eva,d_t_lsc,d_q_eva,d_q_lsc
+USE phys_output_var_mod, ONLY : bils_ec,bils_ech,bils_tke,bils_kinetic,bils_enthalp,bils_latent,bils_diss
 
 IMPLICIT none
@@ -41 +42 @@
 REAL, DIMENSION(klon,klev+1) :: fluxu,fluxv,fluxt
 REAL, DIMENSION(klon,klev+1) :: dddu,dddv,dddt
-REAL, DIMENSION(klon,klev) :: d_u,d_v,d_t,zv,zu
+REAL, DIMENSION(klon,klev) :: d_u,d_v,d_t,zv,zu,d_t_ech
 REAL ZRCPD
 
@@ -131 +132 @@
 
    do k=1,klev
-      d_t_ec(:,k)=-(dddu(:,k)+dddu(:,k+1)+dddv(:,k)+dddv(:,k+1) &
-   &  +rcpd*(dddt(:,k)+dddt(:,k+1)))/(2.*rcpd*masse(:,k))
+      d_t_ech(:,k)=-(rcpd*(dddt(:,k)+dddt(:,k+1)))/(2.*rcpd*masse(:,k))
+      d_t_ec(:,k)=-(dddu(:,k)+dddu(:,k+1)+dddv(:,k)+dddv(:,k+1))/(2.*rcpd*masse(:,k))+d_t_ech(:,k)
    enddo
-! d_t_ec=0.
 
 ENDIF
@@ -141 +141 @@
 !  Computation of integrated enthalpie and kinetic energy variation
 !  FH (hourdin@lmd.jussieu.fr), 2013/04/25
+!  bils_ec : energie conservation term
+!  bils_ech : part of this term linked to temperature
+!  bils_tke : change of TKE
+!  bils_diss : dissipation of TKE (when activated)
+!  bils_kinetic : change of kinetic energie of the column
+!  bils_enthalp : change of enthalpie
+!  bils_latent  : change of latent heat. Computed between
+!          after reevaporation (at the beginning of the physics)
+!          and before large scale condensation (fisrtilp)
 !================================================================
 
       bils_ec(:)=0.
+      bils_ech(:)=0.
       bils_tke(:)=0.
       bils_diss(:)=0.
@@ -151 +161 @@
       DO k=1,klev
         bils_ec(:)=bils_ec(:)-d_t_ec(:,k)*masse(:,k)
+        bils_ech(:)=bils_ech(:)-d_t_ech(:,k)*masse(:,k)
         bils_tke(:)=bils_tke(:)+0.5*(dtke(:,k)+dtke(:,k+1))*masse(:,k)
         bils_diss(:)=bils_diss(:)-d_t_diss(:,k)*masse(:,k)
@@ -157 +168 @@
      &            -puo(:,k)*puo(:,k)-pvo(:,k)*pvo(:,k))
         bils_enthalp(:)= &
-     &  bils_enthalp(:)+masse(:,k)*(ptn(:,k)-pto(:,k)+d_t_ec(:,k))
+     &  bils_enthalp(:)+masse(:,k)*(ptn(:,k)-pto(:,k)+d_t_ec(:,k)-d_t_eva(:,k)-d_t_lsc(:,k))
+!    &  bils_enthalp(:)+masse(:,k)*(ptn(:,k)-pto(:,k)+d_t_ec(:,k))
         bils_latent(:)=bils_latent(:)+masse(:,k)* &
-     &             (pqn(:,k)-pqo(:,k))
+!    &             (pqn(:,k)-pqo(:,k))
+     &             (pqn(:,k)-pqo(:,k)-d_q_eva(:,k)-d_q_lsc(:,k))
       ENDDO
       bils_ec(:)=rcpd*bils_ec(:)/pdtphys
+      bils_ech(:)=rcpd*bils_ech(:)/pdtphys
       bils_tke(:)=bils_tke(:)/pdtphys
       bils_diss(:)=rcpd*bils_diss(:)/pdtphys

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -237 +237 @@
     'bils_diss', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
   TYPE(ctrl_out), SAVE :: o_bils_ec = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
-    'bils_ec', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))
+    'bils_ec', 'Surf. total heat flux correction', 'W/m2', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_bils_ech = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
+    'bils_ech', 'Surf. total heat flux correction', 'W/m2', (/ ('', i=1, 9) /))
   TYPE(ctrl_out), SAVE :: o_bils_kinetic = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11 /), &
     'bils_kinetic', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 9) /))

LMDZ5/trunk/libf/phylmd/phys_output_var_mod.F90

@@ -17 +17 @@
   !$OMP THREADPRIVATE(itau_con)
   REAL, ALLOCATABLE :: bils_ec(:) ! Contribution of energy conservation
+  REAL, ALLOCATABLE :: bils_ech(:) ! Contribution of energy conservation
   REAL, ALLOCATABLE :: bils_tke(:) ! Contribution of energy conservation
   REAL, ALLOCATABLE :: bils_diss(:) ! Contribution of energy conservation
@@ -22 +23 @@
   REAL, ALLOCATABLE :: bils_enthalp(:) ! bilan de chaleur au sol
   REAL, ALLOCATABLE :: bils_latent(:) ! bilan de chaleur au sol
-  !$OMP THREADPRIVATE(bils_ec,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
+  !$OMP THREADPRIVATE(bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
 
 
@@ -83 +84 @@
     allocate(snow_o(klon), zfra_o(klon))
     allocate(itau_con(klon))
-    allocate (bils_ec(klon),bils_tke(klon),bils_diss(klon),bils_kinetic(klon),bils_enthalp(klon),bils_latent(klon))
+    allocate (bils_ec(klon),bils_ech(klon),bils_tke(klon),bils_diss(klon),bils_kinetic(klon),bils_enthalp(klon),bils_latent(klon))
 
     IF (ok_gwd_rando) allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon))
@@ -95 +96 @@
 
     deallocate(snow_o,zfra_o,itau_con)
-    deallocate (bils_ec,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
+    deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
 
   END SUBROUTINE phys_output_var_end

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -45 +45 @@
          o_LWupSFC, o_LWdnSFC, o_LWupSFCclr, &
          o_LWdnSFCclr, o_bils, o_bils_diss, &
-         o_bils_ec, o_bils_tke, o_bils_kinetic, &
+         o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &
          o_bils_latent, o_bils_enthalp, o_sens, &
          o_fder, o_ffonte, o_fqcalving, o_fqfonte, &
@@ -217 +217 @@
 
     USE phys_output_var_mod, only: vars_defined, snow_o, zfra_o, bils_diss, &
-         bils_ec, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
+         bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
          itau_con, nfiles, clef_files, nid_files, zvstr_gwd_rando
     USE indice_sol_mod, only: nbsrf
@@ -467 +467 @@
        CALL histwrite_phy(o_bils_diss, bils_diss)
        CALL histwrite_phy(o_bils_ec, bils_ec)
+       CALL histwrite_phy(o_bils_ech, bils_ech)
        CALL histwrite_phy(o_bils_tke, bils_tke)
        CALL histwrite_phy(o_bils_kinetic, bils_kinetic)