Changeset 1120 for trunk/LMDZ.TITAN/Tools

Timestamp:

Dec 3, 2013, 3:05:40 PM (11 years ago)

Author:

slebonnois

Message:

SL: Titan and Venus modifications following a modif in dyn3d[par].

File:

• trunk/LMDZ.TITAN/Tools/psi.F90 (modified) (2 diffs)

trunk/LMDZ.TITAN/Tools/psi.F90

@@ -170 +170 @@
 allocate(za(lonlength,latlength,altlength,timelength))
 ! present only in _P regrided files
-! For others, using geop/g0
+! For others, using g0*a0*a0/(g0*a0-geop)-a0
 
 text="zareoid"
@@ -184 +184 @@
   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
+        za(ilon,ilat,ilev,itim) = (g0*a0*a0/(g0*a0-za(ilon,ilat,ilev,itim))-a0)/1000. ! in km
     else
         za(ilon,ilat,ilev,itim) = miss_val