Changeset 282 for trunk/LMDZ.MARS/util

Timestamp:

Sep 6, 2011, 12:10:28 PM (13 years ago)

Author:

emillour

Message:

Ajout utilitaire "extract" pour sorties GCM Mars.
EM

Location:

trunk/LMDZ.MARS/util

Files:

• README.exec (modified) (1 diff)
• extract.F90 (added)
• extract.points.def (added)
• extract.profile.def (added)

trunk/LMDZ.MARS/util/README.exec

@@ -129 +129 @@
      (if infile is not specified, "startfi.nc" is used as default)
      (if outfile is not specified, it is built as "infile_ex.nc")
+
+--------------------------------------------------------------------
+6) extract.e
+--------------------------------------------------------------------
+
+This program extracts (ie: interpolates) pointwise values of an atmospheric 
+variable from a 'zrecast'ed diagfi file (works if altitude is geometrical 
+height or a pressure vertical coordinates)
+user has to specify:
+- name of input file
+- date (in sols) offset wrt the input file (e.g. if the input file "begins"
+  at Ls=0, then the offset is 0; if the input file begins at Ls=30, the
+  offset date corresponding to the first 3 months is 61+66+66=193 sols, etc.)
+- the "extraction mode": 
+     1: extract individual values; user will specify values of
+        lon lat alt Ls LT (all on a same line)
+        on as many lines as there are sought values
+     2: extract a profile: user will specify on a first line the values of
+        lon lat Ls LT (all on a same line)
+        and then only specify values of altitudes (m or Pa depending on the
+        coordinate in the input file), one per line, at which values are
+        sought
+- output values are given on standard output; with extraction mode 1 as
+  lines of "lon lat alt Ls LT value" an for a profile (extraction mode 2)
+  as lines of "alt value"
+  NB: If there is no data to do an appropriate interpolation to extract
+      the sought value, then a "missing_value" (taken from the variable's
+      attribute in the input file, most likely -9.99E33) is returned.
+
+See files "extract.points.def" and "extract.profile.def" for examples of
+input format that would simply need to be redirected to extract.e, e.g.:
+
+extract.e < extract.points.def
+
+