Changeset 1591 for trunk/LMDZ.VENUS/libf/phyvenus/clmain.F

Timestamp:

Aug 31, 2016, 4:31:29 PM (8 years ago)

Author:

slebonnois

Message:

SL: Mise a jour de la haute atmosphere, du transfert radiatif (solaire=Rainer Haus; IR: reglages de juin 2016), et implementation des variations de temperature potentielle dans la basse atmosphere (variation de la masse molaire moyenne)

File:

• trunk/LMDZ.VENUS/libf/phyvenus/clmain.F (modified) (8 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/clmain.F

@@ -291 +291 @@
 c        call coefkzmin(knon,ypaprs,ypplay,yu,yv,yt,ycoefm
 c    .   ,ycoefm0,ycoefh0)
-c      call dump2d(iim,jjm-1,ycoefm(2:klon-1,2), 'KZ         ')
-c      call dump2d(iim,jjm-1,ycoefm0(2:klon-1,2),'KZMIN      ')
+c      call dump2d(nbp_lon,nbp_lat-2,ycoefm(2:klon-1,2), 'KZ         ')
+c      call dump2d(nbp_lon,nbp_lat-2,ycoefm0(2:klon-1,2),'KZMIN      ')
  
          ycoefm0 = 1.e-3
@@ -349 +349 @@
         ENDIF
 
-c   iflag_pbl peut etre utilise comme longuer de melange
+c   iflag_pbl peut etre utilise comme longueur de melange
 
          if (iflag_pbl.ge.11) then
@@ -558 +558 @@
       PARAMETER (check=.false.)
 C
+c----------------------
+c ATMOSPHERE PROFONDE
+      real p_lim,p_ref
+      real rsmu_mod(klon,klev),zrsmu_mod(klon,klev)
+      real ratio_mod(klon,klev) 
+! theta_mod = theta * ratio_mod => ratio_mod = mmm0/mmm(p)
+
+      p_lim = 6e6
+      p_ref = 8.9e6
+
+      DO k = 1, klev
+      DO i = 1, klon
+        ratio_mod(i,k) = 1.
+        rsmu_mod(i,k) = RD
+       if (pplay(i,k).gt.p_lim) then
+          ratio_mod(i,k) = RMD /
+     &  (RMD+0.56*(log(pplay(i,k))-log(p_lim))/(log(p_ref)-log(p_lim)))
+           ratio_mod(i,k) = max(ratio_mod(i,k),RMD/(RMD+0.56))
+           rsmu_mod(i,k) = RD*ratio_mod(i,k)
+       endif
+      ENDDO
+      ENDDO
+c----------------------
+
       if (iflag_pbl.eq.1) then
         do k = 3, klev
@@ -584 +608 @@
       DO i = 1, knon
             zt(i,k)    = (t(i,k)+t(i,k-1)) * 0.5
+c----------------------
+c ATMOSPHERE PROFONDE
+            zrsmu_mod(i,k)    = (rsmu_mod(i,k)+rsmu_mod(i,k-1)) * 0.5
+c----------------------
       ENDDO
       ENDDO
@@ -612 +640 @@
       DO i = 1, knon
          zx_coef(i,k) = coef(i,k)*RG/(pplay(i,k-1)-pplay(i,k))
-     .                  *(paprs(i,k)/zt(i,k)/RD)**2
+c----------------------
+c ATMOSPHERE PROFONDE
+     .                  *(paprs(i,k)/zt(i,k)/zrsmu_mod(i,k))**2
+c----------------------
          zx_coef(i,k) = zx_coef(i,k) * dtime*RG
       ENDDO
       ENDDO
 c
-c Preparer les flux lies aux contre-gardients
+c Preparer les flux lies aux contre-gardients  OBSOLETE
 c
       DO k = 2, klev
@@ -984 +1015 @@
       DATA appel1er /.TRUE./
 c
+c----------------------
+c ATMOSPHERE PROFONDE
+      real p_lim,p_ref
+      real modgam(klon,klev)
+
+      p_lim = 6e6
+      p_ref = 8.9e6
+
+      DO k = 1, klev
+      DO i = 1, klon
+        modgam(i,k) = 0.
+       if (pplay(i,k).gt.p_lim) then
+           modgam(i,k) = 0.56/(log(p_ref)-log(p_lim))/R
+       endif
+      ENDDO
+      ENDDO
+c----------------------
+
       isommet=klev
 
@@ -1045 +1094 @@
             zmgeom(i,k)= zgeop(i,k)-zgeop(i,k-1)
             zdphi      = zmgeom(i,k)/2.
-            ztvd(i,k)  = (t(i,k)   + zdphi/cpdet(zt(i,k)))
-            ztvu(i,k)  = (t(i,k-1) - zdphi/cpdet(zt(i,k)))
+c----------------------
+c ATMOSPHERE PROFONDE
+            ztvd(i,k)  = t(i,k)   + zdphi*
+     &                      (1./cpdet(zt(i,k))+modgam(i,k))
+            ztvu(i,k)  = t(i,k-1) - zdphi*
+     &                      (1./cpdet(zt(i,k))+modgam(i,k))
+c----------------------
             zpk(i,k)   = ppk(i,k)*(paprs(i,k)/pplay(i,k))**RKAPPA
       ENDDO
@@ -1054 +1108 @@
         zt(i,1)   = ts(i)
         ztvu(i,1) = ts(i) 
-        ztvd(i,1) = t(i,1)+zgeop(i,1)/cpdet(zt(i,1)) 
+c----------------------
+c ATMOSPHERE PROFONDE
+        ztvd(i,1) = t(i,1)+zgeop(i,1)*
+     &                      (1./cpdet(zt(i,1))+modgam(i,1))
+c----------------------
         zpk(i,1)  = ppk(i,1)*(paprs(i,1)/pplay(i,1))**RKAPPA
       ENDDO