Changeset 2046 for LMDZ5/trunk/libf/phylmd/1DUTILS.h

Timestamp:

May 18, 2014, 3:18:19 PM (10 years ago)

Author:

fhourdin

Message:

Nouvelles options des thermiques en route vers les stratocumulus

+ sorties de la discretisation verticale dans le vieux disvert_old

New options in thermals for stratocumulus

File:

• LMDZ5/trunk/libf/phylmd/1DUTILS.h (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/1DUTILS.h

@@ -822 +822 @@
       INTEGER np,ierr
       REAL pi,x
+      REAL :: scaleheight=8.
  
 !-----------------------------------------------------------------------
@@ -924 +925 @@
       PRINT *,  ap
  
-      DO l = 1, llm
-       dpres(l) = bp(l) - bp(l+1)
-       presnivs(l) = 0.5 *( ap(l)+bp(l)*preff + ap(l+1)+bp(l+1)*preff )
-      ENDDO
+   print*,'scaleheight=',scaleheight
+  DO l = 1, llm
+     dpres(l) = bp(l) - bp(l+1)
+     presnivs(l) = 0.5 *( ap(l)+bp(l)*preff + ap(l+1)+bp(l+1)*preff )
+     write(*, *)'PRESNIVS(', l, ')=', presnivs(l), ' Z ~ ', &
+          log(preff/presnivs(l))*scaleheight &
+          , ' DZ ~ ', scaleheight*log((ap(l)+bp(l)*preff)/ &
+          max(ap(l+1)+bp(l+1)*preff, 1.e-10))
+  ENDDO
+
  
       PRINT *,' PRESNIVS '