Changeset 35

Timestamp:

07/27/12 16:10:12 (12 years ago)

Author:

ymipsl

Message:

add compute_temperature2theta_rhodz functionnlity

YM

File:

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

codes/icosagcm/trunk/src/theta_rhodz.f90

@@ -171 +171 @@
     
   END SUBROUTINE compute_theta_rhodz2temperature
-  
+
+  SUBROUTINE compute_temperature2theta_rhodz(ps,temp,theta_rhodz,offset)
+  USE icosa
+  USE pression_mod
+  USE exner_mod
+  IMPLICIT NONE
+    REAL(rstd),INTENT(IN) :: ps(iim*jjm)
+    REAL(rstd),INTENT(OUT) :: theta_rhodz(iim*jjm,llm)
+    REAL(rstd),INTENT(IN) :: temp(iim*jjm,llm)
+    INTEGER,INTENT(IN) :: offset
+    INTEGER :: i,j,ij,l
+    REAL(rstd) :: p(iim*jjm,llm+1)
+    REAL(rstd) :: pk(iim*jjm,llm)
+    REAL(rstd) :: pks(iim*jjm)
+
+    CALL compute_pression(ps,p,offset)
+    CALL compute_exner(ps,p,pks,pk,offset)
+        
+    DO    l    = 1, llm
+      DO j=jj_begin-offset,jj_end+offset
+        DO i=ii_begin-offset,ii_end+offset
+          ij=(j-1)*iim+i
+!          temp(ij,l) = ( theta_rhodz(ij,l) / ((p(ij,l)-p(ij,l+1))/g) ) * pk(ij,l) / cpp
+          theta_rhodz(ij,l) = temp(ij,l) * ((p(ij,l)-p(ij,l+1))/g) / (pk(ij,l) / cpp )
+        ENDDO
+      ENDDO
+    ENDDO
+    
+  END SUBROUTINE compute_temperature2theta_rhodz
+   
 END MODULE theta2theta_rhodz_mod