Changeset 4843 for LMDZ6/trunk/libf/phylmd/lmdz_thermcell_height.F90

Timestamp:

Mar 4, 2024, 6:58:03 PM (3 months ago)

Author:

crio

Message:

Nouvelle formulation du strig et correction thermiques montent trop haut

File:

• LMDZ6/trunk/libf/phylmd/lmdz_thermcell_height.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_thermcell_height.F90

@@ -2 +2 @@
 CONTAINS
 
-      SUBROUTINE thermcell_height(ngrid,nlay,lalim,lmin,linter,lmix,  &
-     &           zw2,zlev,lmax,zmax,zmax0,zmix,wmax)
+      SUBROUTINE thermcell_height(ngrid,nlay,lalim,lmin,linter,lcong,lintercong,lmix,  &
+     &           zw2,zlev,lmax,zmax,zmax0,zmix,wmax,zcong)
       IMPLICIT NONE
 
@@ -14 +14 @@
 ! Entree
       integer, intent(in) :: ngrid,nlay
-      real, intent(in), dimension(ngrid) :: linter
+      real, intent(in), dimension(ngrid) :: linter,lintercong
       real, intent(in), dimension(ngrid,nlay+1) :: zlev
 ! Sortie
-      real, intent(out), dimension(ngrid) :: wmax,zmax,zmax0,zmix
-      integer, intent(out), dimension(ngrid) :: lmax
+      real, intent(out), dimension(ngrid) :: wmax,zmax,zmax0,zmix,zcong
+      integer, intent(out), dimension(ngrid) :: lmax,lcong
 ! Les deux
      integer, intent(inout), dimension(ngrid) :: lmix,lalim,lmin
@@ -24 +24 @@
 
 ! local
-     real, dimension(ngrid) :: num,denom,zlevinter
+     real, dimension(ngrid) :: num,denom,zlevinter,zlevintercong
      integer ig,l
 
@@ -113 +113 @@
       enddo
 
+!CR:calcul du zcong
+      do  ig=1,ngrid
+! calcul de zlevintercong
+         zlevintercong(ig)=(zlev(ig,lcong(ig)+1)-zlev(ig,lcong(ig)))*  &
+     &    lintercong(ig)+zlev(ig,lcong(ig))-lcong(ig)*(zlev(ig,lcong(ig)+1)  &
+     &    -zlev(ig,lcong(ig)))
+         zcong(ig)=zlevintercong(ig)-zlev(ig,1)
+!         print*,"calcul zcong",lcong(ig),lintercong(ig),zlevintercong(ig),zcong(ig)
+      enddo
 
       endif