Changeset 545
- Timestamp:
- 08/10/17 18:13:14 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/disvert_strato.f90
r467 r545 12 12 SUBROUTINE init_disvert 13 13 USE icosa 14 USE getin_mod 14 15 IMPLICIT NONE 15 16 REAL(rstd) :: dsigmin 16 17 ALLOCATE(ap(llm+1)) 17 18 ALLOCATE(bp(llm+1)) 18 19 ALLOCATE(presnivs(llm)) 19 20 20 CALL disvert(ap,bp,presnivs) 21 dsigmin=1.0 ! Should be 0.3 for CMIP5 22 CALL getin('disvert_dsigmin', dsigmin) 23 CALL disvert(dsigmin,ap,bp,presnivs) 21 24 22 25 END SUBROUTINE init_disvert 23 26 24 27 25 SUBROUTINE disvert( ap,bp,presnivs)28 SUBROUTINE disvert(dsigmin,ap,bp,presnivs) 26 29 USE icosa 27 30 USE mpipara 28 31 USE earth_const 29 32 IMPLICIT NONE 33 REAL(rstd) :: dsigmin 30 34 REAL(rstd),INTENT(OUT) :: ap(:) 31 35 REAL(rstd),INTENT(OUT) :: bp(:) … … 41 45 DO l = 1, llm 42 46 x = Pi*(l-0.5)/(llm+1) 43 dsig(l) = ( 1.0+ 7.0 * SIN(x)**2) * (0.5* (1.-tanh( (x-Pi/2)/(Pi/2) ) ))**247 dsig(l) = (dsigmin + 7.0 * SIN(x)**2) * (0.5* (1.-tanh( (x-Pi/2)/(Pi/2) ) ))**2 44 48 snorm = snorm + dsig(l) 45 49 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.