Changeset 360 for trunk/LMDZ.MARS/util/zrecast.F90

Timestamp:

Nov 8, 2011, 3:14:39 PM (13 years ago)

Author:

acolaitis

Message:

Added possibility for Zrecast to work if there is no temperature, but potential temperature.

Temperature is recomputed from potential temperature and pressure field, using constant r/cp = 0.256793.
R can be taken constant for this conversion from teta to temp because this is what is done in the model in the output section of physiq.F

File:

• trunk/LMDZ.MARS/util/zrecast.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/util/zrecast.F90

@@ -101 +101 @@
 real,dimension(:,:,:,:),allocatable :: press ! GCM atmospheric pressure 
 real,dimension(:,:,:,:),allocatable :: temp ! GCM atmospheric temperature
+real,dimension(:,:,:,:),allocatable :: teta ! GCM atmospheric potential temperature
 real,dimension(:,:,:,:),allocatable :: rho ! GCM atmospheric density
 real,dimension(:,:,:,:),allocatable :: za_gcm ! GCM above areoid levels (m)
@@ -773 +774 @@
 endif
 
+allocate(teta(lonlength,latlength,altlength,timelength),stat=ierr)
+if (ierr.ne.0) then
+  write(*,*) "Error: Failed to allocate teta(lonlength,latlength,altlength,timelength)"
+  write(*,*) "       lonlength=",lonlength," latlength=",latlength
+  write(*,*) "       altlength=",altlength," timelength=",timelength
+  stop
+endif
+
+
 ierr=NF_INQ_VARID(infid,"temp",tmpvarid)
 if (ierr.ne.NF_NOERR) then
@@ -779 +789 @@
   ierr=NF_INQ_VARID(infid,"t",tmpvarid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get t ID"
-  else
-    ierr=NF_GET_VAR_REAL(infid,tmpvarid,temp)
-    if (ierr.ne.NF_NOERR) then
-      stop "Error: Failed reading atmospheric temperature"
-    endif
+     ierr=NF_INQ_VARID(infid,"teta",tmpvarid)
+     if (ierr.ne.NF_NOERR) then
+       stop "Error: Failed to get t or teta ID"
+     endif
+       ierr=NF_GET_VAR_REAL(infid,tmpvarid,teta)
+       if (ierr.ne.NF_NOERR) then
+         stop "Error: Failed reading atmospheric temperature"
+       endif
+     
+  do itim=1,timelength
+    do ilev=1,altlength
+      do ilat=1,latlength
+        do ilon=1,lonlength
+          temp(ilon,ilat,ilev,itim)=teta(ilon,ilat,ilev,itim)*(press(ilon,ilat,ilev,itim)/ps(ilon,ilat,itim))**(.256793)
+        enddo
+      enddo
+    enddo
+  enddo
+
   endif
 else