Changeset 271

Timestamp:

09/09/14 11:28:29 (10 years ago)

Author:

dubos

Message:

Fixed DCMIP1 initial condition

File:

• codes/icosagcm/trunk/src/etat0_dcmip1.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/etat0_dcmip1.f90

@@ -281 +281 @@
 
     SUBROUTINE hadleyq(hx) 
-    IMPLICIT NONE 
-    REAL(rstd)::hx(iim*jjm,llm) 
-    REAL(rstd),PARAMETER:: zz1=3500.,zz2=6500.,zz0=0.5*(zz1+zz2)
-   
+      IMPLICIT NONE 
+      REAL(rstd)::hx(iim*jjm,llm) 
+      REAL(rstd),PARAMETER:: zz1=2000.,zz2=5000.,zz0=0.5*(zz1+zz2)
+      
       DO l=1,llm
-        IF ( ( zz1 .LT. zrl(l) ) .and. ( zrl(l) .LT. zz2 ) )  THEN 
-          hx(:,l) = 0.5*(1. + cos(0.002*pi*(zrl(l)-zz0)/3.))  
-        ELSE
-          hx(:,l) = 0.0 
-        END IF
-      END DO
-    END SUBROUTINE hadleyq 
+         IF ( ( zz1 .LT. zrl(l) ) .and. ( zrl(l) .LT. zz2 ) )  THEN 
+            hx(:,l) = 0.5*(1. + cos(2*pi*(zrl(l)-zz0)/(zz2-zz1)))  
+         ELSE
+            hx(:,l) = 0.0 
+         END IF
+      END DO
+    END SUBROUTINE hadleyq
 
   END SUBROUTINE compute_etat0_ncar
-  
 
 END MODULE etat0_dcmip1_mod