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

Timestamp:

08/05/14 15:56:49 (10 years ago)

Author:

ymipsl

Message:

Synchronize trunk and Saturn branch.
Merge modification from Saturn branch to trunk

YM

File:

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

codes/icosagcm/trunk/src/disvert_apbp.f90

@@ -57 +57 @@
       ENDIF
     ENDDO
+    
+    CLOSE(unit)
 !$OMP END MASTER
     IF (omp_in_parallel()) THEN
@@ -74 +76 @@
       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/1000," (km)"
+      WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scale_height/1000," (km)"
       DO l=1,llm
-        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))
+        WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),'  Z ~ ',log(preff/presnivs(l))*scale_height/1000,       &
+                   ' DZ ~ ',scale_height/1000*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10))
       ENDDO
 !$OMP END MASTER