Changeset 1366 for trunk/UTIL

Timestamp:

Nov 7, 2014, 11:11:45 AM (10 years ago)

Author:

slebonnois

Message:

SL: zrecast corrected and tested to be compatible with Venus and Titan. Adopted.

Location:

trunk/UTIL

Files:

• compile (modified) (1 diff)
• zrecast.F90 (modified) (5 diffs)

trunk/UTIL/compile

@@ -4 +4 @@
 # > compile zrecast
 
-# pgf90 -Bstatic   $1.F90 \
-#-I/distrib/local/netcdf/pgi_7.1-6_32/include \
-#-L/distrib/local/netcdf/pgi_7.1-6_32/lib -lnetcdf  -o $1.e
+# pgf sur les machines du LMD:
+#-----------------------------
+#pgf95 -Bstatic $1.F90 \
+#-I/distrib/local/netcdf/pgi_7.1-6_64/include \
+#-L/distrib/local/netcdf/pgi_7.1-6_64/lib -lnetcdf -o $1.e
 
-#ifort $1.F90 \
-#pgf90 $1.F90 \
+# ifort
+#------
 ifort $1.F90 \
 -I$NETCDF/include \
 -L$NETCDF/lib -lnetcdf -o $1.e
 
+#-----------------------------------------------------------------
 # Before running that on you computer you might want to change :
+#-----------------------------------------------------------------
 # 1) replace "pgf90" with the name of your favorite compiler 
 #    (you may also add some non-agressive optimization options e.g. -O2)

trunk/UTIL/zrecast.F90

@@ -5 +5 @@
 real :: g0 ! gravity at a0  (Mars: g0=3.7257964 ; Lemoine et al. 2001)
 real :: R0 ! reduced gaz constant (Mars: R=191)
+character(len=20) :: planet
 end module planet_const
 
@@ -1714 +1715 @@
     stop
   endif
-else ! hybrid coordinates
+elseif (have_hyb) then  ! hybrid coordinates
   ierr=NF_PUT_VAR_REAL(outfid,aps_varid,aps)
   if (ierr.ne.NF_NOERR) then
@@ -2016 +2017 @@
 ! initialize planetary constants using the "controle" array in the file
 ! if "controle" array not found in file, look for it in "diagfi.nc"
+! if still no "controle" array nowhere, ask for planet name
+
 use planet_const
 use netcdf
@@ -2028 +2031 @@
 integer :: infid2 ! for an alternate input file
 real :: controle(100) ! to store "controle" array
+logical :: have_ctl
+
+planet="dunno"
+have_ctl=.false.
 
 ! look for "controle" in input file
@@ -2033 +2040 @@
 varname="controle" 
 status=nf90_inq_varid(infid,trim(varname),varid)
-if (status.ne.nf90_noerr) then
+
+if (status.eq.nf90_noerr) then
+  write(*,*) "OK, found ",trim(varname)
+  have_ctl=.true.
+else
   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
+  if (status.eq.nf90_noerr) then
+    status=nf90_inq_varid(infid2,trim(varname),varid)
+    if (status.eq.nf90_noerr) then
+      write(*,*) "OK, found ",trim(varname)
+      have_ctl=.true.
+    else
+      write(*,*) " Failed to find ",trim(varname)," in file!!"
+    endif
+  else
     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
+    if (status.eq.nf90_noerr) then
+      status=nf90_inq_varid(infid2,trim(varname),varid)
+      if (status.eq.nf90_noerr) then
+        write(*,*) "OK, found ",trim(varname)
+        have_ctl=.true.
+      else
+        write(*,*) " Failed to find ",trim(varname)," in file!!"
+      endif
     endif
   endif
-  status=nf90_inq_varid(infid2,trim(varname),varid)
+endif
+
+if (have_ctl) then
+  status=nf90_get_var(infid2,varid,controle)
   if (status.ne.nf90_noerr) then
-    write(*,*) " Failed to find ",trim(varname)," in file!!"
+    write(*,*) "init_planet_const: Failed to load ",trim(varname)
     stop
   endif
-  write(*,*) "OK, found ",trim(varname)
+  ! now that we have "controle"; extract relevent informations
+  a0=controle(5) ! = radius of planet (m)
+  g0=controle(7) ! = gravity (m.s-2) at a0
+  R0=1000.*8.314511/controle(8) ! controle(8)=mugaz = molar mass (g.mol-1) of atm.
+else
+   write(*,*) "No controle found." 
+   write(*,*) "You need to tell me on which planet we are:"
+   read(*,'(a20)') planet
+   write(*,*) ""
+   if (trim(planet).eq."titan") then
+     a0=2575.E3 ! = radius of planet (m)
+     g0=1.35 ! = gravity (m.s-2) at a0
+     R0=296.9 ! = molecular gas constant
+   else if (trim(planet).eq."venus") then
+     a0=6051.E3 ! = radius of planet (m)
+     g0=8.87 ! = gravity (m.s-2) at a0
+     R0=191.4 ! = molecular gas constant
+   else
+     write(*,*) "init_planet_const: Failed to init ",trim(planet), &
+                 " constants..."
+     stop 
+   endif
 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
-R0=1000.*8.314511/controle(8) ! controle(8)=mugaz = molar mass (g.mol-1) of atm.
 
 end subroutine init_planet_const