Changeset 111

Timestamp:

08/08/12 02:21:52 (12 years ago)

Author:

dubos

Message:

Fixed DCMIP1.2

File:

• codes/icosagcm/trunk/src/guided_ncar_mod.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/guided_ncar_mod.f90

@@ -6 +6 @@
   
   REAL(rstd), PARAMETER :: alpha=0.0 ! tilt of solid-body rotation
-  REAL(rstd), PARAMETER :: tau = 12*daysec ! 12 days               ! see p. 16
+  REAL(rstd), PARAMETER :: tau_hadley=daysec, tau = 12*daysec ! 12 days               ! see p. 16
   REAL(rstd), PARAMETER :: w0_deform = 23000*pi/tau, b=0.2, ptop=25494.4  ! see p. 16 
   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
        u1 = w0_deform*radius/b/ptop*cos(2*pitbytau)*(exp((ptop-pr)/b/ptop)-exp((pr-ncar_p0)/b/ptop))
-       v0 = -radius*w0_hadley*pi/(5.0*ztop)*(ncar_p0/pr)*cos(pi*zr/ztop)*cos(pitbytau) ! for Hadley cell
+       v0 = -radius*w0_hadley*pi/(5.0*ztop)*(ncar_p0/pr)*cos(pi*zr/ztop)*cos(pi*tt/tau_hadley) ! for Hadley cell
 
        DO j=jj_begin-1,jj_end+1