Changeset 208 for codes/icosagcm/trunk
- Timestamp:
- 07/10/14 11:28:31 (10 years ago)
- Location:
- codes/icosagcm/trunk/src
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/disvert_apbp.f90
r207 r208 33 33 REAL(rstd),INTENT(OUT) :: presnivs(:) 34 34 35 INTEGER,PARAMETER :: unit=4236 35 CHARACTER(len=255) :: filename 37 36 INTEGER :: l,ok … … 69 68 ! tell the world about it 70 69 IF (is_mpi_root) THEN 70 !$OMP MASTER 71 71 WRITE(*,*) "ap()=",ap 72 72 WRITE(*,*) "bp()=",bp 73 73 WRITE(*,*) "Approximative mid-layer pressure, assuming a surface pressure preff=",preff," Pa" 74 WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scaleheight ," (km)"74 WRITE(*,*) "and approximative mid-layer height, assuming an atmospheric scale height of ",scaleheight/1000," (km)" 75 75 DO l=1,llm 76 WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),' Z ~ ',log(preff/presnivs(l))*scaleheight , &77 ' DZ ~ ',scaleheight *log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10))76 WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),' Z ~ ',log(preff/presnivs(l))*scaleheight/1000, & 77 ' DZ ~ ',scaleheight/1000*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10)) 78 78 ENDDO 79 !$OMP END MASTER 79 80 ENDIF 80 81 -
codes/icosagcm/trunk/src/disvert_std.f90
r186 r208 26 26 USE icosa 27 27 USE mpipara 28 USE earth_const 28 29 IMPLICIT NONE 29 30 REAL(rstd),INTENT(OUT) :: ap(:) … … 59 60 ap(1)=0. 60 61 ap(llm+1) = pa * ( sig(llm+1) - bp(llm+1) ) 61 62 IF (is_mpi_root) PRINT*,'ap',ap63 IF (is_mpi_root) PRINT*,'bp',bp64 65 IF (is_mpi_root) PRINT*, 'Niveaux de pressions approximatifs aux centres des'66 IF (is_mpi_root) PRINT*, 'couches calcules pour une pression de surface =', preff67 IF (is_mpi_root) PRINT*, 'et altitudes equivalentes pour une hauteur d echelle de'68 IF (is_mpi_root) PRINT*, '8km'69 70 62 DO l = 1, llm 71 63 presnivs(l) = 0.5 *( ap(l)+bp(l)*preff + ap(l+1)+bp(l+1)*preff ) 72 73 IF (is_mpi_root) PRINT*, 'PRESNIVS(',l,')=',presnivs(l),' Z ~ ',log(preff/presnivs(l))*8., &74 ' DZ ~ ',8.*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10))75 64 ENDDO 65 66 ! tell the world about it 67 IF (is_mpi_root) THEN 68 !$OMP MASTER 69 WRITE(*,*) "ap()=",ap 70 WRITE(*,*) "bp()=",bp 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 ",scaleheight/1000," (km)" 73 DO l=1,llm 74 WRITE(*,*) 'PRESNIVS(',l,')=',presnivs(l),' Z ~ ',log(preff/presnivs(l))*scaleheight/1000, & 75 ' DZ ~ ',scaleheight/1000*log((ap(l)+bp(l)*preff)/ max(ap(l+1)+bp(l+1)*preff,1.e-10)) 76 ENDDO 77 !$OMP END MASTER 78 ENDIF 76 79 77 80 END SUBROUTINE disvert -
codes/icosagcm/trunk/src/earth_const.f90
r207 r208 11 11 REAL(rstd),SAVE :: preff=101325. 12 12 REAL(rstd),SAVE :: pa=50000. 13 REAL(rstd),SAVE :: scaleheight=8 . ! atmospheric scale height (km)13 REAL(rstd),SAVE :: scaleheight=8000. ! atmospheric scale height (m) 14 14 REAL(rstd),SAVE :: scale_factor=1. 15 15
Note: See TracChangeset
for help on using the changeset viewer.