Changeset 829

Timestamp:

Nov 6, 2012, 12:02:45 PM (12 years ago)

Author:

emillour

Message:

Mars GCM utility:

update of zrecast utility: now the planetary constants (radius, gravity,

reduced gas constant) are read from the "controle" array in the file and
not hard coded. That way zrecast can also be used for the generic model.
Note that the "distance to center of Mars" and "MCD" options still contain
some hard coded Mars constants.

EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• util/zrecast.F90 (modified) (23 diffs)

trunk/LMDZ.MARS/README

@@ -1815 +1815 @@
    If ndynstep is not specified or is negative, nday is used. Otherwise nday value is discarded.
    The problem with nday alone is that one can only run the GCM for a decimal fraction of one sol.
+
+== 06/11/2012 == EM
+>> update of zrecast utility: now the planetary constants (radius, gravity,
+   reduced gas constant) are read from the "controle" array in the file and
+   not hard coded. That way zrecast can also be used for the generic model.
+   Note that the "distance to center of Mars" and "MCD" options still contain
+   some hard coded Mars constants.

trunk/LMDZ.MARS/util/zrecast.F90

@@ -1 @@
-
+module planet_const
+! planetary constants (set via init_planet_const routine)
+implicit none
+real :: a0 ! Mean planetary radius (a0=3396.E3 for Mars)
+real :: g0 ! gravity at a0  (Mars: g0=3.7257964 ; Lemoine et al. 2001)
+real :: Rmean ! reduced gaz constant (Mars: R=191)
+end module planet_const
 
 program zrecast
@@ -51 +57 @@
 ! EM 03/2011 : Add possibility to have outputs as distance to center
 !              of planet
+! EM 11/2012 : Adapted so it can be used on "generic" model outputs; planet
+!              constants (radius, R, etc.) are now read from file
 implicit none
 
@@ -154 +162 @@
 endif
 
+! load planet constants (radius, gravity, ...)
+
+call init_planet_const(infid)
+ 
 !===============================================================================
 ! 1.2 Get # and names of variables in input file
@@ -1743 +1755 @@
 end program
 
+!===============================================================================
+
+subroutine init_planet_const(infid)
+! initialize planetary constants using the "controle" array in the file
+! if "cointrole" array not found in file, look for it in "diagfi.nc"
+use planet_const
+use netcdf
+implicit none
+
+integer,intent(in) :: infid ! input file ID
+
+! local variables
+character(len=100) :: varname
+integer :: varid ! variable ID
+integer :: status
+integer :: infid2 ! for an alternate input file
+real :: controle(100) ! to store "controle" array
+
+! look for "controle" in input file
+infid2=infid ! initialization
+varname="controle" 
+status=nf90_inq_varid(infid,trim(varname),varid)
+if (status.ne.nf90_noerr) then
+  write(*,*) "init_planet_const: Failed to find ",trim(varname)
+  write(*,*) " looking for it in file diagfi.nc"
+  status=nf90_open("diagfi.nc",NF90_NOWRITE,infid2)
+  if (status.ne.nf90_noerr) then
+    write(*,*) " no diafi.nc file ... looking for diagfi1.nc"
+    status=nf90_open("diagfi1.nc",NF90_NOWRITE,infid2)
+    if (status.ne.nf90_noerr) then
+      write(*,*) "might as well stop here..."
+      stop
+    endif
+  endif
+  status=nf90_inq_varid(infid2,trim(varname),varid)
+  if (status.ne.nf90_noerr) then
+    write(*,*) " Failed to find ",trim(varname)," in file!!"
+    stop
+  endif
+  write(*,*) "OK, found ",trim(varname)
+endif
+status=nf90_get_var(infid2,varid,controle)
+if (status.ne.nf90_noerr) then
+  write(*,*) "init_planet_const: Failed to load ",trim(varname)
+  stop
+endif
+
+! now that we have "controle"; extract relevent informations
+a0=controle(5) ! = radius of planet (m)
+g0=controle(7) ! = gravity (m.s-2) at a0
+Rmean=1000.*8.314511/controle(8) ! controle(8)=mugaz = molar mass (g.mol-1) of atm.
+
+end subroutine init_planet_const
 
 !===============================================================================
@@ -1752 +1817 @@
 !          surface altitudes" corresponding to GCM atmospheric levels
 !==============================================================================
+use planet_const, only : a0,g0
 implicit none
 !==============================================================================
@@ -1785 +1851 @@
 
 ! Parameters needed to integrate hydrostatic equation:
-real,parameter :: g0=3.7257964
+!real,parameter :: g0=3.7257964
 !g0: exact mean gravity at radius=3396.km (Lemoine et al. 2001)
-real,parameter :: a0=3396.E3
+!real,parameter :: a0=3396.E3
 !a0: 'mean' Mars radius=3396.km
 real :: gz 
@@ -1871 +1937 @@
 !===============================================================================
 
-!#include"build_gcm_alt.F90"
 subroutine build_gcm_za(lonlength,latlength,altlength,timelength, &
                          phis,ps,press,temp,rho,za_gcm)
@@ -1878 +1943 @@
 !          altitudes" corresponding to GCM atmospheric levels
 !==============================================================================
+use planet_const, only : a0,g0
 implicit none
 !==============================================================================
@@ -1911 +1977 @@
 
 ! Parameters needed to integrate hydrostatic equation:
-real,parameter :: g0=3.7257964
+!real,parameter :: g0=3.7257964
 !g0: exact mean gravity at radius=3396.km (Lemoine et al. 2001)
-real,parameter :: a0=3396.E3
+!real,parameter :: a0=3396.E3
 !a0: 'mean' Mars radius=3396.km
 real :: gz 
@@ -2062 +2128 @@
 !          surface altitudes" corresponding to GCM atmospheric levels
 !==============================================================================
+use planet_const, only : a0,g0,Rmean
 implicit none
 !==============================================================================
@@ -2086 +2153 @@
 !===============================================================================
 real,dimension(:),allocatable :: sigma ! GCM sigma levels
-real,parameter :: R=191 ! molecular gas constant
+!real,parameter :: R=191 ! molecular gas constant
 !real,dimension(:,:,:,:),allocatable :: zlocal ! local above surface altitude
 real :: Tmean ! "mean" value of temperature for a given level
@@ -2092 +2159 @@
 
 ! Parameters needed to integrate hydrostatic equation:
-real,parameter :: g0=3.7257964
+!real,parameter :: g0=3.7257964
 !g0: exact mean gravity at radius=3396.km (Lemoine et al. 2001)
-real,parameter :: a0=3396.E3
+!real,parameter :: a0=3396.E3
 !a0: 'mean' Mars radius=3396.km
 real :: gz 
@@ -2114 +2181 @@
       sigma(1)=press(iloop,jloop,1,tloop)/ps(iloop,jloop,tloop)
       zs_gcm(iloop,jloop,1,tloop)=-log(sigma(1))* &
-                                               R*temp(iloop,jloop,1,tloop)/g0
+                                           Rmean*temp(iloop,jloop,1,tloop)/g0
 !      za_gcm(iloop,jloop,1,tloop)=zlocal(iloop,jloop,1,tloop)+ &
 !                                  phis(iloop,jloop)/g0
@@ -2139 +2206 @@
         ! compute zs_gcm(iloop,jloop,kloop,tloop)
         zs_gcm(iloop,jloop,kloop,tloop)=zs_gcm(iloop,jloop,kloop-1,tloop)- &
-                log(sigma(kloop)/sigma(kloop-1))*R*Tmean/gz
+                log(sigma(kloop)/sigma(kloop-1))*Rmean*Tmean/gz
         
         ! compute za_gcm(kloop)
@@ -2157 +2224 @@
 !===============================================================================
 
-!#include"crude_gcm_alt.F90"
 subroutine crude_gcm_za(lonlength,latlength,altlength,timelength, &
                          phis,ps,press,temp,za_gcm)
@@ -2164 +2230 @@
 !          altitudes" corresponding to GCM atmospheric levels
 !==============================================================================
+use planet_const, only : a0,g0,Rmean
 implicit none
 !==============================================================================
@@ -2188 +2255 @@
 !===============================================================================
 real,dimension(:),allocatable :: sigma ! GCM sigma levels
-real,parameter :: R=191 ! molecular gas constant
+!real,parameter :: R=191 ! molecular gas constant
 real,dimension(:,:,:,:),allocatable :: zlocal ! local above surface altitude
 real :: Tmean ! "mean" value of temperature for a given level
@@ -2194 +2261 @@
 
 ! Parameters needed to integrate hydrostatic equation:
-real,parameter :: g0=3.7257964
+!real,parameter :: g0=3.7257964
 !g0: exact mean gravity at radius=3396.km (Lemoine et al. 2001)
-real,parameter :: a0=3396.E3
+!real,parameter :: a0=3396.E3
 !a0: 'mean' Mars radius=3396.km
 real :: gz 
@@ -2216 +2283 @@
       sigma(1)=press(iloop,jloop,1,tloop)/ps(iloop,jloop,tloop)
       zlocal(iloop,jloop,1,tloop)=-log(sigma(1))* &
-                                               R*temp(iloop,jloop,1,tloop)/g0
+                                            Rmean*temp(iloop,jloop,1,tloop)/g0
       za_gcm(iloop,jloop,1,tloop)=zlocal(iloop,jloop,1,tloop)+ &
                                   phis(iloop,jloop)/g0
@@ -2239 +2306 @@
         ! compute zlocal(kloop)
         zlocal(iloop,jloop,kloop,tloop)=zlocal(iloop,jloop,kloop-1,tloop)- &
-                log(sigma(kloop)/sigma(kloop-1))*R*Tmean/gz
+                log(sigma(kloop)/sigma(kloop-1))*Rmean*Tmean/gz
         
         ! compute za_gcm(kloop)
@@ -2332 +2399 @@
 !===============================================================================
 
-!#include"za_coord_interp.F90"
 subroutine z_coord_interp(lonlen,latlen,altlen,tlen,newaltlen, &
                           missing,z_gcm,gcmdata,flag,z_new,newdata)
@@ -2462 +2528 @@
 !       if flag=0, and extrapolated (exponentially) if flag=1.
 !==============================================================================
+use planet_const, only : a0,g0,Rmean
 implicit none
 !==============================================================================
@@ -2503 +2570 @@
 real :: x,y ! temporary variables
 integer :: iloop,jloop,kloop,tloop
-real,parameter :: g0=3.7257964
+!real,parameter :: g0=3.7257964
 !g0: exact mean gravity at radius=3396.km (Lemoine et al. 2001)
-real,parameter :: a0=3396.E3
+!real,parameter :: a0=3396.E3
 !a0: 'mean' Mars radius=3396.km
-real,parameter :: Rmean=191 ! molecular gas constant
+!real,parameter :: Rmean=191 ! molecular gas constant
 real :: gz 
 ! gz: gravitational acceleration at a given (above areoid) altitude