Changeset 397 for codes/icosagcm/trunk/src/physics_dcmip2016.f90

Timestamp:

06/06/16 20:41:52 (8 years ago)

Author:

ymipsl

Message:

Prepare DCMIP2016 output by XIOS2

YM

File:

• codes/icosagcm/trunk/src/physics_dcmip2016.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/physics_dcmip2016.f90

@@ -120 +120 @@
     REAL(rstd) :: qfi(ngrid,llm,5)
 
-    REAL(rstd)  :: rho(llm), z(llm), theta(llm)
+    REAL(rstd)  :: rho(llm), z(llm), theta(llm), qv(llm),qc(llm),qr(llm)
     REAL(rstd)  :: lastz
     REAL(rstd)  :: dcl1,dcl2
-    INTEGER :: l,ll,ij
+     INTEGER :: l,ll,ij
     REAL(rstd) :: dt_phys, inv_dt
-
+    INTEGER :: simple_physic_testcase
+    
     ! prepare input fields and mirror vertical index      
     ps(:) = p(:,1) ! surface pressure
@@ -147 +148 @@
     ENDDO
     
+    
     IF (testcase==moist_baroclinic .OR. testcase==cyclone  ) THEN
-      CALL simple_physics(ngrid, llm, dt_phys, lat, tfi, qfi(:,:,1) , ufi, vfi, pmid, pint, pdel, 1./pdel, ps, precl, 1, .FALSE., .FALSE.)
+      IF (testcase==moist_baroclinic) simple_physic_testcase=1
+      IF (testcase==cyclone) simple_physic_testcase=0
+      CALL simple_physics(ngrid, llm, dt_phys, lat, tfi, qfi(:,:,1) , ufi, vfi, pmid, pint, pdel, 1./pdel, ps, precl, &
+                          simple_physic_testcase, .TRUE., .FALSE.)
     ENDIF
 
@@ -164 +169 @@
           ENDDO
 
-          CALL KESSLER(theta(:), qfi(ij,llm:1:-1,1), qfi(ij,llm:1:-1,2), qfi(ij,llm:1:-1,3), rho(:),  &
-                       pk(ij,:), dt_phys, z(:), llm, precl(ij)) 
+          qv(:)=qfi(ij,llm:1:-1,1)
+          qc(:)=qfi(ij,llm:1:-1,2)
+          qr(:)=qfi(ij,llm:1:-1,3)
+          
+!          CALL KESSLER(theta(:), qv, qc, qr, rho(:),  &
+!                       pk(ij,:), dt_phys, z(:), llm, precl(ij)) 
+          
           
           DO l=1,llm
@@ -172 +182 @@
            Tfi(ij,ll) = theta(l)  * ( pk(ij,l) / cpp)
           ENDDO
+
+          qfi(ij,:,1)=qv(llm:1:-1)
+          qfi(ij,:,2)=qc(llm:1:-1)
+          qfi(ij,:,3)=qr(llm:1:-1)
+
        ENDDO
     ENDIF