Changeset 386 for codes

Timestamp:

05/25/16 17:58:59 (8 years ago)

Author:

ymipsl

Message:

Some fixes in dcmip2016 physics

YM

Location:

codes/icosagcm/trunk/src

Files:

• etat0_dcmip2016_baroclinic_wave.f90 (modified) (1 diff)
• physics_dcmip2016.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/etat0_dcmip2016_baroclinic_wave.f90

@@ -33 +33 @@
        STOP
     END IF
-    CALL set_advection_scheme(3,advect_none)
+!    CALL set_advection_scheme(3,advect_none)
     
   END SUBROUTINE getin_etat0

codes/icosagcm/trunk/src/physics_dcmip2016.f90

@@ -65 +65 @@
     ! Physics-specific data
     ALLOCATE(precl_packed(ngrid))
+    precl_packed(:)=0.
     CALL allocate_field(f_precl, field_t,type_real)
 
@@ -84 +85 @@
     CALL unpack_field(f_precl, precl_packed)
     CALL output_field("precl",f_precl)
+    precl_packed(:)=0.
+    
   END SUBROUTINE write_physics
 
@@ -154 +157 @@
           DO l=1,llm
            ll=llm+1-l
-           rho(l) = gas_constant*Temp(ij,l) / pmid(ij,ll)
-           z(l)=lastz+ (p(ij,l)-p(ij,l+1) /g) / rho(l)
+           rho(l) = pmid(ij,ll)/(287*Temp(ij,l))
+           z(l)=lastz+ (p(ij,l)-p(ij,l+1)) /g / rho(l)
            lastz=z(l)
-           theta(l)= Tfi(ij,ll)*(1+0.608*qfi(ij,ll,1)) * pk(ij,l) / cpp
+!           theta(l)= Tfi(ij,ll)*(1+0.608*qfi(ij,ll,1)) / ( pk(ij,l) / cpp)
+           theta(l)= Tfi(ij,ll) / ( pk(ij,l) / cpp)
           ENDDO
 
@@ -165 +169 @@
           DO l=1,llm
            ll=llm+1-l
-           Tfi(ij,ll) = theta(l) /(1+0.608*qfi(ij,ll,1)) * cpp / pk(ij,l)
+!           Tfi(ij,ll) = theta(l) /(1+0.608*qfi(ij,ll,1)) * ( pk(ij,l) / cpp)
+           Tfi(ij,ll) = theta(l)  * ( pk(ij,l) / cpp)
           ENDDO
        ENDDO