- Timestamp:
- 09/09/14 11:28:29 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/etat0_dcmip1.f90
r186 r271 281 281 282 282 SUBROUTINE hadleyq(hx) 283 IMPLICIT NONE284 REAL(rstd)::hx(iim*jjm,llm)285 REAL(rstd),PARAMETER:: zz1=3500.,zz2=6500.,zz0=0.5*(zz1+zz2)286 283 IMPLICIT NONE 284 REAL(rstd)::hx(iim*jjm,llm) 285 REAL(rstd),PARAMETER:: zz1=2000.,zz2=5000.,zz0=0.5*(zz1+zz2) 286 287 287 DO l=1,llm 288 IF ( ( zz1 .LT. zrl(l) ) .and. ( zrl(l) .LT. zz2 ) ) THEN289 hx(:,l) = 0.5*(1. + cos(0.002*pi*(zrl(l)-zz0)/3.))290 ELSE291 292 END IF293 END DO 294 END SUBROUTINE hadleyq 288 IF ( ( zz1 .LT. zrl(l) ) .and. ( zrl(l) .LT. zz2 ) ) THEN 289 hx(:,l) = 0.5*(1. + cos(2*pi*(zrl(l)-zz0)/(zz2-zz1))) 290 ELSE 291 hx(:,l) = 0.0 292 END IF 293 END DO 294 END SUBROUTINE hadleyq 295 295 296 296 END SUBROUTINE compute_etat0_ncar 297 298 297 299 298 END MODULE etat0_dcmip1_mod
Note: See TracChangeset
for help on using the changeset viewer.