Changeset 114

Timestamp:

08/08/12 17:46:28 (12 years ago)

Author:

ymipsl

Message:

correction for initial temperature on dcmip testcase 5

YM

File:

• codes/icosagcm/trunk/src/etat0_dcmip5.f90 (modified) (5 diffs)

codes/icosagcm/trunk/src/etat0_dcmip5.f90

@@ -64 +64 @@
   USE wind_mod
   USE theta2theta_rhodz_mod
+  USE exner_mod
   IMPLICIT NONE  
   REAL(rstd),INTENT(OUT) :: ps(iim*jjm)
@@ -84 +85 @@
   REAL(rstd) :: ulat(3*iim*jjm,llm)
   REAL(rstd) :: lon,lat
+  REAL(rstd) :: pks(iim*jjm) 
+  REAL(rstd) :: pk(iim*jjm,llm) 
 
 
@@ -149 +152 @@
           ENDIF
           T(ij,l)=Tave+Tp 
+          out_i(ij,l)=Tave+Tp
         ENDDO
       ENDDO
@@ -230 +234 @@
 !          ulon(ij+u_ldown,l)=vt*d1/d
 !          ulat(ij+u_ldown,l)=vt*d2/d
-          out_i(ij,l)=sqrt((vt*d1/d)**2+(vt*d2/d)**2)
+!          out_i(ij,l)=sqrt((vt*d1/d)**2+(vt*d2/d)**2)
         ENDDO
       ENDDO
@@ -236 +240 @@
     
     CALL compute_wind_perp_from_lonlat_compound(ulon, ulat, ue)
-    CALL compute_temperature2theta_rhodz(T,ps,theta_rhodz,0)
-      
+    CALL compute_temperature2theta_rhodz(ps,T,theta_rhodz,0)
+
   END SUBROUTINE compute_etat0_dcmip5