Changeset 208 for codes/icosagcm/trunk/src/disvert_apbp.f90

Timestamp:

07/10/14 11:28:31 (10 years ago)

Author:

ymipsl

Message:

Atmosphere scaleheight is not hard coded in disvert_std, and given in meters instead km
YM

File:

• codes/icosagcm/trunk/src/disvert_apbp.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/disvert_apbp.f90

@@ -33 +33 @@
   REAL(rstd),INTENT(OUT) :: presnivs(:)
   
-  INTEGER,PARAMETER :: unit=42
   CHARACTER(len=255) :: filename
   INTEGER :: l,ok
@@ -69 +68 @@
     ! tell the world about it
     IF (is_mpi_root) THEN
+!$OMP MASTER
       WRITE(*,*) "ap()=",ap
       WRITE(*,*) "bp()=",bp
       WRITE(*,*) "Approximative mid-layer pressure, assuming a surface pressure preff=",preff," Pa"
-      WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scaleheight," (km)"
+      WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scaleheight/1000," (km)"
       DO l=1,llm
-        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))
+        WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),'  Z ~ ',log(preff/presnivs(l))*scaleheight/1000,       &
+                   ' DZ ~ ',scaleheight/1000*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10))
       ENDDO
+!$OMP END MASTER
     ENDIF