Changeset 466

Timestamp:

Dec 9, 2011, 5:10:15 PM (13 years ago)

Author:

acolaitis

Message:

Added localtime.e and streamfunction.e to the gcm wrapper.

Location:

trunk/UTIL/PYTHON

Files:

• README.MCS (modified) (5 diffs)
• gcm_transformations.py (modified) (2 diffs)

trunk/UTIL/PYTHON/README.MCS

@@ -27 +27 @@
 -------
 -------
-Ok, but after all, I want to interpolate these fields on the vertical grid used by MCS... (make sure to have rho in
-your file for a better interpolation. You also need aps, bps, and ps)
+Ok, but after all, I want to interpolate these fields on the vertical grid used in the MCS data... (make sure to have rho in
+your file for a better interpolation. You also need aps, bps, ps and phisinit, at least in diagfi.nc or phisinit.nc.)
 
 command -> mcs.py -f $YourSimulation/diagfi.nc -m /d5/lmlmd/MCS_processeddata/MCSdata_binned_MY29.nc --var temp,co2ice -i
@@ -34 +34 @@
 -------
 -------
-Ok, but I dont care for the intermediary diagfi_P.nc file, it takes to much space. I want to remove it automatically
+Ok, but I dont care for the intermediary diagfi_P.nc file, it takes too much space. I want to remove it automatically
 after the computation
 
@@ -44 +44 @@
 
 command -> mcs.py -f $YourSimulation/diagfi.nc -m /d5/lmlmd/MCS_processeddata/MCSdata_binned_MY29.nc --var temp,co2ice -i --ditch -H
-output -> diagfi_P_h_MCS.nc
+output -> diagfi_h_P_MCS.nc
 -------
 -------
@@ -51 +51 @@
 
 command -> mcs.py -f $YourSimulation/diagfi.nc -m /d5/lmlmd/MCS_processeddata/MCSdata_binned_MY29.nc --var temp,co2ice -i -H
-output -> diagfi_P_h_MCS.nc
+output -> diagfi_h_P_MCS.nc
 -------
 -------
@@ -57 +57 @@
 
 command -> mcs.py -f $YourSimulation/diagfi.nc -m /d5/lmlmd/MCS_processeddata/MCSdata_binned_MY29.nc --var temp,co2ice -i -H --override
-output -> diagfi_P_h_MCS.nc
+output -> diagfi_P.nc , diagfi_h_P.nc , diagfi_h_P_MCS.nc
 -------

trunk/UTIL/PYTHON/gcm_transformations.py

@@ -3 +3 @@
 ### A. Colaitis -- LMD -- 08/11/2011 <-- Added ZRECAST for altitude mode, and MCS and TES modes for Pressure
 ### A. Colaitis -- LMD -- 08/11/2011 <-- Added HRECAST for MCS and TES predefined grids
-###
-
-### TODO
-### add wrapper to streamfunction.e
-### add wrapper to localtime.e
+###                                      Added wrapper to streamfunction.e
+###                                      Added wrapper to localtime.e
 
 def call_zrecast (  interp_mode   = '4', \
@@ -119 +116 @@
     return output_name
 
+def call_localtime (  input_name = None, \
+                    fields  = 'all', \
+                    times = None):
+
+    import numpy as np
+    from myplot import separatenames
+    from os import system
+    system("rm -f localtime.auto.def")
+    system("touch localtime.auto.def")
+    indicefile=0
+    outputfilename=""
+    append="_LT"
+    f = open('localtime.auto.def', 'w')
+    for zfile in input_name:
+        f.write(zfile+"\n")
+        for zvar in fields:
+             f.write(zvar+"\n")
+        f.write("\n")
+    f.write(str(len(times))+"\n")
+    for t in times:
+        f.write(str(t)+"\n")
+    f.close()
+    system("localtime.e < localtime.auto.def")
+
+    if indicefile is 0:
+          outputfilename=input_name[indicefile][0:len(input_name[0])-3]+append+".nc"
+    else:
+          outputfilename=outputfilename+","+input_name[indicefile][0:len(input_name[0])-3]+append+".nc"
+    indicefile = indicefile +1
+
+    for i in range(len(outputfilename)):
+        output_name = separatenames(outputfilename)
+
+    return output_name
+
+def call_streamfunction (  input_name = None):
+    import numpy as np
+    from myplot import separatenames
+    from os import system
+    system("rm -f streamfunction.auto.def")
+    system("touch streamfunction.auto.def")
+    indicefile=0
+    outputfilename=""
+    append="_stream"
+    f = open('streamfunction.auto.def', 'w')
+    for zfile in input_name:
+        f.write(zfile+"\n")
+    f.close()
+    system("streamfunction.e < streamfunction.auto.def")
+
+    if indicefile is 0:
+          outputfilename=input_name[indicefile][0:len(input_name[0])-3]+append+".nc"
+    else:
+          outputfilename=outputfilename+","+input_name[indicefile][0:len(input_name[0])-3]+append+".nc"
+    indicefile = indicefile +1
+
+    for i in range(len(outputfilename)):
+        output_name = separatenames(outputfilename)
+
+    return output_name
+