Changeset 385

Timestamp:

05/25/16 12:07:23 (8 years ago)

Author:

ymipsl

Message:

Some fixes in dcmip2016 physics

YM

File:

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

codes/icosagcm/trunk/src/physics_dcmip2016.f90

@@ -151 +151 @@
     IF (testcase==moist_baroclinic .OR. testcase==cyclone .OR. testcase==supercell ) THEN
        DO ij=1,ngrid
-          lastz=0 ;
+          lastz=0 
           DO l=1,llm
            ll=llm+1-l
-           rho(l) = gas_constant*Temp(ij,ll) / pmid(ij,ll)
+           rho(l) = gas_constant*Temp(ij,l) / pmid(ij,ll)
            z(l)=lastz+ (p(ij,l)-p(ij,l+1) /g) / rho(l)
            lastz=z(l)
-           theta(l)= Tfi(ij,l)*(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
           ENDDO
 
@@ -165 +165 @@
           DO l=1,llm
            ll=llm+1-l
-           Tfi(ij,l) = 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)) * cpp / pk(ij,l)
           ENDDO
        ENDDO