- 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_apbp.f90
r221 r234 76 76 WRITE(*,*) "bp()=",bp 77 77 WRITE(*,*) "Approximative mid-layer pressure, assuming a surface pressure preff=",preff," Pa" 78 WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scale height/1000," (km)"78 WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scale_height/1000," (km)" 79 79 DO l=1,llm 80 WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),' Z ~ ',log(preff/presnivs(l))*scale height/1000, &81 ' DZ ~ ',scale height/1000*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10))80 WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),' Z ~ ',log(preff/presnivs(l))*scale_height/1000, & 81 ' DZ ~ ',scale_height/1000*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10)) 82 82 ENDDO 83 83 !$OMP END MASTER
Note: See TracChangeset
for help on using the changeset viewer.