Changeset 980 for trunk/LMDZ.TITAN/Tools

Timestamp:

Jun 6, 2013, 2:41:55 PM (12 years ago)

Author:

slebonnois

Message:

SL: encore une correction sur psi.F90 (posttraitement Venus, Titan)

File:

• trunk/LMDZ.TITAN/Tools/psi.F90 (modified) (1 diff)

trunk/LMDZ.TITAN/Tools/psi.F90

@@ -169 +169 @@
 
 allocate(za(lonlength,latlength,altlength,timelength))
+! present only in _P regrided files
+! For others, using geop/g0
 
 text="zareoid"
 call get_var4d(infid,lonlength,latlength,altlength,timelength,text,za,ierr1,ierr2)
-if (ierr1.ne.NF_NOERR) stop "Error: Failed to get za ID"
-if (ierr2.ne.NF_NOERR) stop "Error: Failed reading zareoid"
+if (ierr1.ne.NF_NOERR) then
+  write(*,*) "  looking for geop instead... "
+  text="geop"
+  call get_var4d(infid,lonlength,latlength,altlength,timelength,text,za,ierr1,ierr2)
+  if (ierr1.ne.NF_NOERR) stop "Error: Failed to get geop ID"
+do itim=1,timelength
+do ilon=1,lonlength
+ do ilat=1,latlength
+  do ilev=1,altlength
+    if (za(ilon,ilat,ilev,itim).ne.miss_val) then
+        za(ilon,ilat,ilev,itim) = za(ilon,ilat,ilev,itim)/(g0*1000.) ! in km
+    else
+        za(ilon,ilat,ilev,itim) = miss_val
+    endif
+  enddo
+ enddo
+enddo
+enddo
+endif
+if (ierr2.ne.NF_NOERR) stop "Error: Failed reading zareoid/geop"
 
 !===============================================================================