- Timestamp:
- 07/17/14 16:41:08 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/disvert_std.f90
r221 r234 70 70 WRITE(*,*) "bp()=",bp 71 71 WRITE(*,*) "Approximative mid-layer pressure, assuming a surface pressure preff=",preff," Pa" 72 WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scale height/1000," (km)"72 WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scale_height/1000," (km)" 73 73 DO l=1,llm 74 WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),' Z ~ ',log(preff/presnivs(l))*scale height/1000, &75 ' DZ ~ ',scale height/1000*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10))74 WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),' Z ~ ',log(preff/presnivs(l))*scale_height/1000, & 75 ' DZ ~ ',scale_height/1000*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10)) 76 76 ENDDO 77 77 !$OMP END MASTER
Note: See TracChangeset
for help on using the changeset viewer.