Changeset 2891 for trunk/LMDZ.GENERIC/libf

Timestamp:

Feb 7, 2023, 4:13:22 PM (23 months ago)

Author:

jleconte

Message:

potential energy in rain_generic added

correction : generic tracer is not always CH4

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• physiq_mod.F90 (modified) (2 diffs)
• rain_generic.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -1671 +1671 @@
             zdqssnow_generic(1:ngrid,1:nq)    = 0.0
 
-            call rain_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,pt,pdt,pq,pdq,            &
+            call rain_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,pphi,pt,pdt,pq,pdq,            &
                         zdtrain_generic,dq_rain_generic_vap,dq_rain_generic_cld,              &
                         zdqsrain_generic,zdqssnow_generic,reevap_precip_generic,cloudfrac)
@@ -2692 +2692 @@
          CALL send_xios_field('h2o_ice',zq(:,:,igcm_h2o_ice))
       ENDIF
-      
+
       CALL send_xios_field("ISR",fluxtop_dn)
       CALL send_xios_field("OLR",fluxtop_lw)

trunk/LMDZ.GENERIC/libf/phystd/rain_generic.F90

@@ -1 @@
-subroutine rain_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,t,pdt,pq,pdq,d_t,dq_rain_generic_vap,dq_rain_generic_cld,dqsrain_generic,dqssnow_generic,reevap_precip,rneb)
+subroutine rain_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,pphi,t,pdt,pq,pdq,d_t,dq_rain_generic_vap,dq_rain_generic_cld,dqsrain_generic,dqssnow_generic,reevap_precip,rneb)
 
 
@@ -33 +33 @@
    real,intent(in) :: pplev(ngrid,nlayer+1) ! inter-layer pressure (Pa)
    real,intent(in) :: pplay(ngrid,nlayer)   ! mid-layer pressure (Pa)
+   real,intent(in) :: pphi(ngrid,nlayer)   ! mid-layer geopotential
    real,intent(in) :: t(ngrid,nlayer) ! input temperature (K)
    real,intent(in) :: pdt(ngrid,nlayer) ! input tendency on temperature (K/s)      
@@ -283 +284 @@
                         precip_rate(i)  = precip_rate_tmp(i)
                      end if
+#ifdef MESOSCALE 
+                     d_t(i,k) = d_t(i,k)+(pphi(i,k+1)-pphi(i,k))*precip_rate(i)*ptimestep/(RCPD*dmass(i,k))
+                     ! JL22. Accounts for gravitational energy of falling precipitations (probably not to be used in the GCM
+                     !   where the counterpart is not included in the dynamics.)
+#endif
                   endif ! of if (precip_rate(i) .GT.0.)
                enddo