Changeset 2047 for trunk/LMDZ.VENUS/libf/phyvenus/gwstress.F

Timestamp:

Nov 29, 2018, 4:47:07 PM (6 years ago)

Author:

slebonnois

Message:

SL: VENUS, ajout des modifs apportees par Thomas Navarro pour la parametrisation des ondes de gravite orographiques

File:

• trunk/LMDZ.VENUS/libf/phyvenus/gwstress.F (modified) (6 diffs)

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

@@ -5 +5 @@
      *         , prho  , pstab , pvph  , pstd, psig
      *         , pmea , ppic , pval  , ptfr  , ptau  
-     *         , pgeom1 , pgamma , pd1  , pd2   , pdmod , pnu )
+     *         , pgeom1 , pgamma , pd1  , pd2   , pdmod , pnu 
+     *         , zeff )
 c
 c**** *gwstress*
@@ -74 +75 @@
       real pmea(nlon),ppic(nlon),pval(nlon)
       real pdmod(nlon)
+      real zeff(nlon) ! effective height seen by the flow when there is blocking
 c
 c-----------------------------------------------------------------------
@@ -79 +81 @@
 c*       0.2   local arrays
 c              ------------
-c  zeff--real: effective height seen by the flow when there is blocking
 
       integer jl
-      real zeff  
 c
 c-----------------------------------------------------------------------
@@ -101 +101 @@
 c
 c
+      zeff = 0.
       do 301 jl=kidia,kfdia
       if(ktest(jl).eq.1) then
@@ -106 +107 @@
 c  effective mountain height above the blocked flow
   
-         zeff=ppic(jl)-pval(jl)
+         zeff(jl)=ppic(jl)-pval(jl)
          if(kkenvh(jl).lt.klev)then
-         zeff=amin1(GFRCRIT*pvph(jl,klev+1)/sqrt(pstab(jl,klev+1))
-     c              ,zeff)
+         zeff(jl)=amin1(GFRCRIT*pvph(jl,klev+1)/sqrt(pstab(jl,klev+1))
+     c              ,zeff(jl))
          endif
 
@@ -116 +117 @@
      *     *psig(jl)*pdmod(jl)/4./pstd(jl)
      *     *pvph(jl,klev+1)*sqrt(pstab(jl,klev+1))
-     *     *zeff**2
+     *     *zeff(jl)**2