Changeset 441

Timestamp:

06/17/16 22:00:08 (8 years ago)

Author:

ymipsl

Message:

Fix virtual temperature for supercell initial state;

YM

File:

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

codes/icosagcm/trunk/src/etat0_dcmip2016_supercell.f90

@@ -28 +28 @@
     USE omp_para
     USE dcmip2016_supercell_mod, ONLY : supercell_init, supercell_test, supercell_z
+    USE terminator, ONLY: initial_value_Terminator
     INTEGER, INTENT(IN) :: ngrid
     REAL(rstd),INTENT(IN) :: lon(ngrid)
@@ -50 +51 @@
          p=0.5*(ap(l)+ap(l+1) + (bp(l)+bp(l+1)) * ps(ij))
          CALL supercell_test(lon(ij),lat(ij),p,z,0,ulon(ij,l),ulat(ij,l),temp(ij,l),thetav, ps(ij),rho, q(ij,l,1),perturbation) 
+        
+        IF (physics_thermo==thermo_fake_moist) temp(ij,l)=Temp(ij,l)*(1+0.608*q(ij,l,1))
+        q(ij,l,2)=0.
+        q(ij,l,3)=0.
+        CALL initial_value_Terminator(lat(ij),lon(ij),q(ij,l,4),q(ij,l,5))
+
        END DO
     END DO
+ 
     
   END SUBROUTINE compute_etat0