Changeset 111 for codes/icosagcm/trunk
- Timestamp:
- 08/08/12 02:21:52 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/guided_ncar_mod.f90
r73 r111 6 6 7 7 REAL(rstd), PARAMETER :: alpha=0.0 ! tilt of solid-body rotation 8 REAL(rstd), PARAMETER :: tau = 12*daysec ! 12 days ! see p. 168 REAL(rstd), PARAMETER :: tau_hadley=daysec, tau = 12*daysec ! 12 days ! see p. 16 9 9 REAL(rstd), PARAMETER :: w0_deform = 23000*pi/tau, b=0.2, ptop=25494.4 ! see p. 16 10 10 REAL(rstd), PARAMETER :: u0_hadley=40.,w0_hadley=0.15 ,ztop= 12000. … … 76 76 zr = -kappa*cpp*ncar_T0/g*log(pr/ncar_p0) ! reciprocal of (1) p. 13, isothermal atmosphere 77 77 u1 = w0_deform*radius/b/ptop*cos(2*pitbytau)*(exp((ptop-pr)/b/ptop)-exp((pr-ncar_p0)/b/ptop)) 78 v0 = -radius*w0_hadley*pi/(5.0*ztop)*(ncar_p0/pr)*cos(pi*zr/ztop)*cos(pi tbytau) ! for Hadley cell78 v0 = -radius*w0_hadley*pi/(5.0*ztop)*(ncar_p0/pr)*cos(pi*zr/ztop)*cos(pi*tt/tau_hadley) ! for Hadley cell 79 79 80 80 DO j=jj_begin-1,jj_end+1
Note: See TracChangeset
for help on using the changeset viewer.