Changeset 4843 for LMDZ6/trunk/libf/phylmd/lmdz_thermcell_env.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_env.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_thermcell_env.F90

@@ -3 +3 @@
 
    SUBROUTINE thermcell_env(ngrid,nlay,po,pt,pu,pv,pplay,  &
-     &           pplev,zo,zh,zl,ztv,zthl,zu,zv,zpspsk,pqsat,lev_out)
+     &           pplev,zo,zh,zl,ztv,zthl,zu,zv,zpspsk,pqsat,lcong,lintercong,lev_out)
 
 !--------------------------------------------------------------
@@ -22 +22 @@
    real,intent(out), dimension(ngrid,nlay) :: zo,zl,zh,ztv,zthl
    real,intent(out), dimension(ngrid,nlay) :: zpspsk,zu,zv,pqsat
-   
+   real, intent(out), dimension(ngrid) :: lintercong
+   integer, intent(out), dimension(ngrid) :: lcong   
 ! Local
 
@@ -78 +79 @@
       enddo
    enddo
+!CR: Calcul du niveau de congelation
+   do ig=1,ngrid
+      lcong(ig)=1
+   enddo
+   do ig=1,ngrid
+      do ll=1,nlay-1
+         if ((zh(ig,ll).gt.273.15).and.(zh(ig,ll+1).le.273.15)) then
+            lcong(ig)=ll+1
+            lintercong(ig)=(ll*(zh(ig,ll+1)-zh(ig,ll))  &
+     &               -zh(ig,ll)+273.15)/(zh(ig,ll+1)-zh(ig,ll))
+         endif
+      enddo
+   enddo 
  
  RETURN