- Timestamp:
- 07/17/14 16:41:08 (10 years ago)
- Location:
- codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src
- Files:
-
- 3 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 -
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 -
codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/earth_const.f90
r221 r234 11 11 REAL(rstd),SAVE :: preff=101325. 12 12 REAL(rstd),SAVE :: pa=50000. 13 REAL(rstd),SAVE :: scale height=8000. ! atmospheric scale height (m)13 REAL(rstd),SAVE :: scale_height=8000. ! atmospheric scale height (m) 14 14 REAL(rstd),SAVE :: scale_factor=1. 15 15 … … 30 30 CALL getin("cpp",cpp) 31 31 CALL getin("preff",preff) 32 CALL getin("scale height",scaleheight)32 CALL getin("scale_height",scale_height) 33 33 34 34 boussinesq=.FALSE.
Note: See TracChangeset
for help on using the changeset viewer.