Changeset 412 for codes

Timestamp:

06/09/16 19:24:04 (8 years ago)

Author:

ymipsl

Message:

Positivate traceur value before entering to dcmip2016 kessler physics.

YM

File:

• codes/icosagcm/trunk/src/physics_dcmip2016.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/physics_dcmip2016.f90

@@ -170 +170 @@
            theta(l)= Tfi(ij,ll) / ( pk(ij,l) / cpp)
           ENDDO
-
-          qv(:)=qfi(ij,llm:1:-1,1)
-          qc(:)=qfi(ij,llm:1:-1,2)
-          qr(:)=qfi(ij,llm:1:-1,3)
+          
+          DO l=1,llm-1
+           z(l)= 0.5*(z(l)+z(l+1))
+          ENDDO
+          z(llm)=z(llm)+(z(llm)-z(llm-1))
+          
+          qv(:)=max(qfi(ij,llm:1:-1,1),0.)
+          qc(:)=max(qfi(ij,llm:1:-1,2),0.)
+          qr(:)=max(qfi(ij,llm:1:-1,3),0.)
           
           CALL KESSLER(theta(:), qv, qc, qr, rho(:),  &