Changeset 1417 for trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

Timestamp:

Apr 17, 2015, 1:02:55 PM (10 years ago)

Author:

sglmd

Message:

correction on plevrad at the top of the model: option ptopzero= .true. to activate it

File:

• trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

@@ -20 +20 @@
       use comcstfi_mod, only: pi, mugaz, cpp
       use callkeys_mod, only: varactive,diurnal,tracer,water,nosurf,varfixed,satval,    &
-                kastprof,strictboundcorrk,specOLR,CLFvarying
+                kastprof,strictboundcorrk,specOLR,CLFvarying,ptopzero
 
       implicit none
@@ -607 +607 @@
       
       plevrad(1) = 0.
-      plevrad(2) = max(pgasmin,0.0001*plevrad(3))
+
+      if (ptopzero) then
+         plevrad(2) = 0. 
+      else
+         plevrad(2) = max(pgasmin,0.0001*plevrad(3))
+      endif
 
       tlevrad(1) = tlevrad(2)
@@ -614 +619 @@
       tmid(1) = tlevrad(2)
       tmid(2) = tlevrad(2)
+      tmid(3) = tlevrad(2)
+      tmid(4) = tlevrad(2)
       pmid(1) = plevrad(2)
       pmid(2) = plevrad(2)
-      
-      DO l=1,L_NLAYRAD-1
+      pmid(3) = plevrad(2)
+      pmid(4) = plevrad(2)
+     
+      DO l=2,L_NLAYRAD-1
          tmid(2*l+1) = tlevrad(2*l)
          tmid(2*l+2) = tlevrad(2*l)