Changeset 110 for codes/icosagcm/trunk/src/physics_dcmip.f90

Timestamp:

08/07/12 19:10:05 (12 years ago)

Author:

ymipsl

Message:

Correction for dcmip moist physics
Temperature is ouput instead of virtual temperature.

YM

File:

• codes/icosagcm/trunk/src/physics_dcmip.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/physics_dcmip.f90

@@ -172 +172 @@
     
 !    q=0
-    out_i=T
+!    out_i=T
     
     CALL simple_physics(iim*jjm, llm, dt, lat, tfi, qfi , ufi, vfi, pmid, pint, pdel, 1/pdel, ps, precl, testcase) 
@@ -198 +198 @@
     ENDDO       
 
-    out_i=q
+!    out_i=q
         
     utemp=utemp-u
@@ -496 +496 @@
          do i=1,pcols
             qsat = epsilo*e0/pmid(i,k)*exp(-latvap/rh2o*((1._r8/t(i,k))-1._r8/T0))  ! saturation specific humidity
+            out_i(i,llm+1-k)=q(i,k)-qsat
             if (q(i,k) > qsat) then                                                 ! saturated?
                tmp  = 1._r8/dtime*(q(i,k)-qsat)/(1._r8+(latvap/cpair)*(epsilo*latvap*qsat/(rair*t(i,k)**2)))