Changeset 1356 for trunk/LMDZ.TITAN/Tools

Timestamp:

Oct 8, 2014, 9:26:28 AM (10 years ago)

Author:

slebonnois

Message:

SL: update to newstart/start2archive tools in Venus+Titan / additional diagnostics in radiative fluxes for Titan

Location:

trunk/LMDZ.TITAN/Tools

Files:

• angmom.F90 (modified) (3 diffs)
• compile_pgf (modified) (1 diff)
• planet.h (modified) (1 diff)

trunk/LMDZ.TITAN/Tools/angmom.F90

@@ -242 +242 @@
 call get_var3d(infid,lonlength,latlength,timelength,text,ps,ierr1,ierr2)
 if (ierr1.ne.NF_NOERR) then
-  write(*,*) "  looking for psol instead... "
-  text="psol"
+  write(*,*) "  looking for ps instead... "
+  text="ps"
   call get_var3d(infid,lonlength,latlength,timelength,text,ps,ierr1,ierr2)
-  if (ierr1.ne.NF_NOERR) stop "Error: Failed to get psol ID"
+  if (ierr1.ne.NF_NOERR) then
+    write(*,*) "  looking for psol instead... "
+    text="psol"
+    call get_var3d(infid,lonlength,latlength,timelength,text,ps,ierr1,ierr2)
+    if (ierr1.ne.NF_NOERR) stop "Error: Failed to get psol ID"
+  endif
 endif
 if (ierr2.ne.NF_NOERR) stop "Error: Failed reading surface pressure"
@@ -467 +472 @@
   do ilev=1,altlength
 ! Need to be consistent with GCM computations
-!    if (za(ilon,ilat,ilev,itim).ne.miss_val) then
+    if (vitu(ilon,ilat,ilev,itim).ne.miss_val) then
      rayon(ilon,ilat,ilev,itim) = a0
 !     rayon(ilon,ilat,ilev,itim) = za(ilon,ilat,ilev,itim) + a0
       grav(ilon,ilat,ilev,itim) = g0*a0*a0 &
                  /(rayon(ilon,ilat,ilev,itim)*rayon(ilon,ilat,ilev,itim))
-!    else
-!     rayon(ilon,ilat,ilev,itim) = miss_val
-!      grav(ilon,ilat,ilev,itim) = miss_val
-!    endif
+    else
+     rayon(ilon,ilat,ilev,itim) = miss_val
+      grav(ilon,ilat,ilev,itim) = miss_val
+    endif
   enddo
  enddo
@@ -526 +531 @@
  do ilat=1,latlength
   do ilev=1,altlength
+    if (rayon(ilon,ilat,ilev,itim).ne.miss_val) then
       oaam(itim) = oaam(itim) &
        + osam(ilon,ilat,ilev,itim)/ hadday * dmass(ilon,ilat,ilev,itim)
       raam(itim) = raam(itim) &
        + rsam(ilon,ilat,ilev,itim)/ hadday * dmass(ilon,ilat,ilev,itim)
+    endif
   enddo
  enddo

trunk/LMDZ.TITAN/Tools/compile_pgf

@@ -1 +1 @@
 # path for netcdf should be adapted to your configuration !
 
-\rm planet.h
-ln -s $2.h planet.h
+\cp -f $2.h planet.h
  pgf95 -Bstatic cpdet.F90 moyzon.F moyzon2.F moytim.F dx_dp.F epflux.F90 \
     io.F90 dmass.F90 reverse.F90 $1.F90 \

trunk/LMDZ.TITAN/Tools/planet.h

@@ -1 +1 @@
 ! Parameters needed to integrate hydrostatic equation:
 
-real,parameter :: g0=8.87
-!g0: exact mean gravity at radius=6051.km
+real,parameter :: g0=1.35
+!g0: exact mean gravity at radius=2575.km
 
-real,parameter :: a0=6051.E3
-!a0: 'mean' radius=6051.km
+real,parameter :: a0=2575.E3
+!a0: 'mean' radius=2575.km
 
-real,parameter :: R0=191.4 ! molecular gas constant
+real,parameter :: R0=296.9 ! molecular gas constant
 
-real,parameter :: psref=9.2e6 ! reference pressure at surface (Pa)
+real,parameter :: psref=1.4e5 ! reference pressure at surface (Pa)
 
-real,parameter :: omega=2.992677e-7 ! angular rotation speed (s-1)
+real,parameter :: omega=4.5238899E-06 ! angular rotation speed (s-1)
 
-real,parameter :: localday=1.0087e7 ! local day (s)
+real,parameter :: localday=1.37889e6 ! local day (s)
 
-character (len=5),parameter :: planet="Venus"
+character (len=5),parameter :: planet="Titan"
 
-real,parameter :: cp0=1000. !doit etre egal a cpp (dyn) et RCPD (phy)
-real,parameter :: t0=460.
-real,parameter :: nu=0.35
-
+real,parameter :: cp0=1039. !doit etre egal a cpp (dyn) et RCPD (phy)
+real,parameter :: t0=0.
+real,parameter :: nu=0.