Changeset 550 for codes/icosagcm/devel
- Timestamp:
- 09/18/17 13:58:05 (7 years ago)
- Location:
- codes/icosagcm/devel/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dissip/dissip_gcm.f90
r533 r550 87 87 88 88 89 INTEGER :: i,j,ij,ind,it,iter 89 INTEGER :: i,j,ij,ind,it,iter,M 90 90 91 91 rayleigh_friction_key='none' … … 159 159 u=f_u(ind) 160 160 161 ! set random seed to get reproductibility when using a different number of process 162 CALL RANDOM_SEED(size=M) 163 CALL RANDOM_SEED(put=(/(i,i=1,M)/)) 164 161 165 DO j=jj_begin,jj_end 162 166 DO i=ii_begin,ii_end … … 247 251 u=f_u(ind) 248 252 249 DO j=jj_begin,jj_end 253 ! set random seed to get reproductibility when using a different number of process 254 CALL RANDOM_SEED(size=M) 255 CALL RANDOM_SEED(put=(/(i,i=1,M)/)) 256 257 DO j=jj_begin,jj_end 250 258 DO i=ii_begin,ii_end 251 259 ij=(j-1)*iim+i … … 335 343 CALL swap_geometry(ind) 336 344 theta=f_theta(ind) 337 345 346 ! set random seed to get reproductibility when using a different number of process 347 CALL RANDOM_SEED(size=M) 348 CALL RANDOM_SEED(put=(/(i,i=1,M)/)) 349 338 350 DO j=jj_begin,jj_end 339 351 DO i=ii_begin,ii_end -
codes/icosagcm/devel/src/vertical/disvert_strato.f90
r531 r550 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.