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

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

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callcorrk.F90 (modified) (3 diffs)
• callkeys_mod.F90 (modified) (1 diff)
• inifis.F (modified) (1 diff)

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)

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

@@ -7 +7 @@
      &   , callstats,calleofdump                                        &
      &   , callgasvis,continuum,H2Ocont_simple,graybody                 &
-     &   , strictboundcorrk                                        
+     &   , strictboundcorrk,ptopzero                                        
 
       logical enertest

trunk/LMDZ.GENERIC/libf/phystd/inifis.F

@@ -228 +228 @@
          call getin_p("strictboundcorrk",strictboundcorrk)
          write(*,*) "strictboundcorrk = ",strictboundcorrk
+
+         write(*,*) "correction on plevrad in the top levels",
+     &              " (callcorrk)?"
+         ptopzero=.false. ! default value
+         call getin_p("ptopzero",ptopzero)
+         write(*,*) "ptopzero = ",ptopzero
 
          write(*,*) "call gaseous absorption in the visible bands?",