Changeset 4146 for trunk/LMDZ.GENERIC/libf/phygeneric/vdif_kc.F

Timestamp:

Mar 19, 2026, 2:35:46 PM (11 days ago)

Author:

gmilcareck

Message:

Thermodynamics update on LMDZ.GENERIC

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/vdif_kc.F (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/phygeneric/vdif_kc.F

@@ -1 @@
-      SUBROUTINE vdif_kc(ngrid,nlay,nq,dt,g,zlev,zlay,u,v,
-     &                       zq,teta,cd,q2,km,kn)
+      SUBROUTINE vdif_kc(ngrid,nlay,nq,dt,g,zlev,zlay,pplev,pplay,u,v,
+     &                       zt,zq,teta,cd,q2,km,kn)
       use generic_cloud_common_h, only: epsi_generic
       use generic_tracer_index_mod, only: generic_tracer_index
       use callkeys_mod, only: generic_condensation, 
-     &                        virtual_theta_correction
+     &                        virtual_theta_correction,thermo_phy
+      use thermo_mod, only: rcp
+      use tracer_h
       IMPLICIT NONE
 c.......................................................................
@@ -35 +37 @@
       REAL,INTENT(IN) :: zlev(ngrid,nlay+1)
       REAL,INTENT(IN) :: zlay(ngrid,nlay)
+      REAL,INTENT(IN) :: pplev(ngrid,nlay+1)
+      REAL,INTENT(IN) :: pplay(ngrid,nlay)
       REAL,INTENT(IN) :: u(ngrid,nlay)
       REAL,INTENT(IN) :: v(ngrid,nlay)
+      REAL,INTENT(IN) :: zt(ngrid,nlay)
       REAL,INTENT(IN) :: zq(ngrid,nlay,nq)
       REAL,INTENT(IN) :: teta(ngrid,nlay)
@@ -59 +64 @@
       REAL unsdzdec(ngrid,nlay+1)
       REAL q(ngrid,nlay+1)
+      REAL ztv(ngrid,nlay)
       REAL tetav(ngrid,nlay)
       integer iq,igcm_generic_vap, igcm_generic_ice
@@ -263 +269 @@
 c
 c.......................................................................
+c  Brunt-Vaisala frequency
+c.......................................................................
+c
+
+      SELECT CASE (TRIM(thermo_phy))
+      
+      CASE ('thermo_uni_ideal')
+
+
 c  Virtual theta correction
-c.......................................................................
-c
 
       if((generic_condensation) .and. (virtual_theta_correction)) THEN
@@ -283 +296 @@
       endif 
 
-c
-c-----------------------------------------------------------------------
       DO ilev=2,nlev-1
        DO igrid=1,ngrid
-c-----------------------------------------------------------------------
-c
+
       if((generic_condensation) .and. (virtual_theta_correction)) THEN
         n2(igrid,ilev)=g*unsdzdec(igrid,ilev)
@@ -317 +327 @@
         m(igrid,ilev)=sqrt(m2(igrid,ilev))
         mpre(igrid,ilev)=m(igrid,ilev)
-c
-c-----------------------------------------------------------------------
+
        ENDDO
       ENDDO
-c-----------------------------------------------------------------------
-c
+      
+      END SELECT
+c-----------------------------------------------------------------------
+c
+
+
       DO igrid=1,ngrid
         m2(igrid,nlev)=m2(igrid,nlev-1)