| [464] | 1 | ### Wrapper for all gcm utilities. |
|---|
| 2 | ### |
|---|
| 3 | ### A. Colaitis -- LMD -- 08/11/2011 <-- Added ZRECAST for altitude mode, and MCS and TES modes for Pressure |
|---|
| 4 | ### A. Colaitis -- LMD -- 08/11/2011 <-- Added HRECAST for MCS and TES predefined grids |
|---|
| [466] | 5 | ### Added wrapper to streamfunction.e |
|---|
| 6 | ### Added wrapper to localtime.e |
|---|
| [464] | 7 | |
|---|
| 8 | def call_zrecast ( interp_mode = '4', \ |
|---|
| 9 | input_name = None, \ |
|---|
| 10 | fields = 'all', \ |
|---|
| 11 | predefined = None): |
|---|
| 12 | |
|---|
| 13 | import numpy as np |
|---|
| 14 | from myplot import separatenames |
|---|
| 15 | from os import system |
|---|
| 16 | pressure_axis_tes=[1658.152,1291.37,1005.72,783.2555,610.,475.0685,369.9837,288.1436,224.4065,174.7679,136.1094,106.0021,82.55452,64.29353,50.07185,38.99599,30.37011,23.65227,18.4204,14.34583,11.17254] |
|---|
| 17 | pressure_axis_mcs=[1878.9, 1658.2, 1463.3, 1291.4, 1139.6, 1005.7, 887.54, 783.26,691.22, 610., 538.32, 475.07, 419.25, 369.98, 326.51, 288.14, 254.29,224.41, 198.04, 174.77, 154.23, 136.11, 120.12, 106., 93.547, 82.555,72.854, 64.294, 56.739, 50.072, 44.188, 38.996, 34.414, 30.37, 26.802,23.652, 20.873, 18.42, 16.256, 14.346, 12.66, 11.173, 9.8597, 8.7012,7.6788, 6.7765, 5.9802, 5.2775, 4.6574, 4.1101, 3.6272, 3.201, 2.8249,2.4929, 2.2, 1.9415, 1.7134, 1.512, 1.3344, 1.1776, 1.0392, 0.9171,0.80934, 0.71424, 0.63031, 0.55625, 0.49089, 0.43321, 0.3823, 0.33738,0.29774, 0.26275, 0.23188, 0.20463, 0.18059, 0.15937, 0.14064, 0.12412,0.10953, 0.096661, 0.085303, 0.07528, 0.066434, 0.058628, 0.051739,0.04566, 0.040294, 0.03556, 0.031381, 0.027694, 0.02444, 0.021568,0.019034, 0.016797, 0.014824, 0.013082, 0.011545, 0.010188, 0.0089909,0.0079345, 0.0070021, 0.0061794, 0.0054533, 0.0048125, 0.004247] |
|---|
| 18 | system("rm -f zrecast.auto.def") |
|---|
| 19 | system("touch zrecast.auto.def") |
|---|
| 20 | indicefile=0 |
|---|
| 21 | outputfilename="" |
|---|
| 22 | f = open('zrecast.auto.def', 'w') |
|---|
| 23 | for zfile in input_name: |
|---|
| 24 | f.write(zfile+"\n") |
|---|
| 25 | for zvar in fields: |
|---|
| 26 | f.write(zvar+"\n") |
|---|
| 27 | f.write("\n") |
|---|
| 28 | if interp_mode == 4: |
|---|
| 29 | append="_S" |
|---|
| 30 | f.write("3"+"\n") |
|---|
| 31 | f.write("yes"+"\n") |
|---|
| 32 | f.write("0 150000"+"\n") |
|---|
| 33 | elif interp_mode == 2: |
|---|
| 34 | append="_P" |
|---|
| 35 | if predefined in ["TES","tes"]: |
|---|
| 36 | print "Using TES pressure levels" |
|---|
| 37 | f.write("1"+"\n") |
|---|
| 38 | f.write("no"+"\n") |
|---|
| 39 | f.write(str(len(pressure_axis_tes))+"\n") |
|---|
| 40 | for zp in pressure_axis_tes: |
|---|
| 41 | f.write(str(zp)+"\n") |
|---|
| 42 | f.write("\n") |
|---|
| 43 | elif predefined in ["MCS","mcs"]: |
|---|
| 44 | print "Using MCS pressure levels" |
|---|
| 45 | f.write("1"+"\n") |
|---|
| 46 | f.write("no"+"\n") |
|---|
| 47 | f.write(str(len(pressure_axis_mcs))+"\n") |
|---|
| 48 | for zp in pressure_axis_mcs: |
|---|
| 49 | f.write(str(zp)+"\n") |
|---|
| 50 | f.write("\n") |
|---|
| 51 | else: |
|---|
| 52 | f.write("1"+"\n") |
|---|
| 53 | f.write("yes"+"\n") |
|---|
| 54 | f.write("370 0.1"+"\n") |
|---|
| 55 | f.write("20"+"\n") |
|---|
| 56 | else: |
|---|
| 57 | print "zrecast interp option unsupported for now. Exiting." |
|---|
| 58 | exit() |
|---|
| 59 | f.close() |
|---|
| 60 | system("zrecast.e < zrecast.auto.def") |
|---|
| 61 | |
|---|
| 62 | if indicefile is 0: |
|---|
| 63 | outputfilename=input_name[indicefile][0:len(input_name[0])-3]+append+".nc" |
|---|
| 64 | else: |
|---|
| 65 | outputfilename=outputfilename+","+input_name[indicefile][0:len(input_name[0])-3]+append+".nc" |
|---|
| 66 | indicefile = indicefile +1 |
|---|
| 67 | |
|---|
| 68 | for i in range(len(outputfilename)): |
|---|
| 69 | output_name = separatenames(outputfilename) |
|---|
| 70 | |
|---|
| 71 | return output_name |
|---|
| 72 | |
|---|
| 73 | def call_hrecast ( input_name = None, \ |
|---|
| 74 | fields = 'all', \ |
|---|
| 75 | predefined = None): |
|---|
| 76 | |
|---|
| 77 | import numpy as np |
|---|
| 78 | from myplot import separatenames |
|---|
| 79 | from os import system |
|---|
| 80 | # Beware that latitude axis in mcs and tes files is reversed |
|---|
| 81 | # When comparing to MCS or TES data, one has to reverse the latitude axis and data along latitude axis. |
|---|
| 82 | # This is for exemple automatically done in mcs.py. |
|---|
| 83 | # This could also be done in a future tes.py. |
|---|
| [465] | 84 | latitudes_reversed_tes_mcs=[88.5, 85.5, 82.5, 79.5, 76.5, 73.5, 70.5, 67.5, 64.5, 61.5, 58.5, 55.5, 52.5, 49.5, 46.5, 43.5, 40.5, 37.5, 34.5, 31.5, 28.5, 25.5, 22.5, 19.5, 16.5, 13.5, 10.5, 7.5, 4.5, 1.5, -1.5, -4.5, -7.5, -10.5, -13.5, -16.5, -19.5, -22.5, -25.5, -28.5, -31.5, -34.5, -37.5, -40.5, -43.5, -46.5, -49.5, -52.5, -55.5, -58.5, -61.5, -64.5, -67.5, -70.5, -73.5, -76.5, -79.5, -82.5, -85.5, -88.5] |
|---|
| [464] | 85 | longitudes_tes_mcs=[-176.25, -168.75, -161.25, -153.75, -146.25, -138.75, -131.25,-123.75, -116.25, -108.75, -101.25, -93.75, -86.25, -78.75, -71.25, -63.75, -56.25, -48.75, -41.25, -33.75, -26.25, -18.75, -11.25, -3.75, 3.75, 11.25, 18.75, 26.25, 33.75, 41.25, 48.75, 56.25, 63.75, 71.25, 78.75, 86.25, 93.75, 101.25, 108.75, 116.25, 123.75, 131.25, 138.75, 146.25, 153.75, 161.25, 168.75, 176.25] |
|---|
| 86 | system("rm -f hrecast.auto.def") |
|---|
| 87 | system("touch hrecast.auto.def") |
|---|
| 88 | indicefile=0 |
|---|
| 89 | outputfilename="" |
|---|
| 90 | append="_h" |
|---|
| 91 | f = open('hrecast.auto.def', 'w') |
|---|
| 92 | for zfile in input_name: |
|---|
| 93 | f.write(zfile+"\n") |
|---|
| 94 | for zvar in fields: |
|---|
| 95 | f.write(zvar+"\n") |
|---|
| 96 | f.write("\n") |
|---|
| 97 | if predefined in ["MCS","mcs","TES","tes"]: |
|---|
| 98 | f.write(str(len(longitudes_tes_mcs))+"\n") |
|---|
| 99 | for lon in longitudes_tes_mcs: |
|---|
| 100 | f.write(str(lon)+"\n") |
|---|
| 101 | f.write(str(len(latitudes_reversed_tes_mcs))+"\n") |
|---|
| 102 | for lat in latitudes_reversed_tes_mcs: |
|---|
| 103 | f.write(str(lat)+"\n") |
|---|
| 104 | f.close() |
|---|
| 105 | system("hrecast.e < hrecast.auto.def") |
|---|
| 106 | |
|---|
| 107 | if indicefile is 0: |
|---|
| 108 | outputfilename=input_name[indicefile][0:len(input_name[0])-3]+append+".nc" |
|---|
| 109 | else: |
|---|
| 110 | outputfilename=outputfilename+","+input_name[indicefile][0:len(input_name[0])-3]+append+".nc" |
|---|
| 111 | indicefile = indicefile +1 |
|---|
| 112 | |
|---|
| 113 | for i in range(len(outputfilename)): |
|---|
| 114 | output_name = separatenames(outputfilename) |
|---|
| 115 | |
|---|
| 116 | return output_name |
|---|
| 117 | |
|---|
| [466] | 118 | def call_localtime ( input_name = None, \ |
|---|
| 119 | fields = 'all', \ |
|---|
| 120 | times = None): |
|---|
| 121 | |
|---|
| 122 | import numpy as np |
|---|
| 123 | from myplot import separatenames |
|---|
| 124 | from os import system |
|---|
| 125 | system("rm -f localtime.auto.def") |
|---|
| 126 | system("touch localtime.auto.def") |
|---|
| 127 | indicefile=0 |
|---|
| 128 | outputfilename="" |
|---|
| 129 | append="_LT" |
|---|
| 130 | f = open('localtime.auto.def', 'w') |
|---|
| 131 | for zfile in input_name: |
|---|
| 132 | f.write(zfile+"\n") |
|---|
| 133 | for zvar in fields: |
|---|
| 134 | f.write(zvar+"\n") |
|---|
| 135 | f.write("\n") |
|---|
| 136 | f.write(str(len(times))+"\n") |
|---|
| 137 | for t in times: |
|---|
| 138 | f.write(str(t)+"\n") |
|---|
| 139 | f.close() |
|---|
| 140 | system("localtime.e < localtime.auto.def") |
|---|
| 141 | |
|---|
| 142 | if indicefile is 0: |
|---|
| 143 | outputfilename=input_name[indicefile][0:len(input_name[0])-3]+append+".nc" |
|---|
| 144 | else: |
|---|
| 145 | outputfilename=outputfilename+","+input_name[indicefile][0:len(input_name[0])-3]+append+".nc" |
|---|
| 146 | indicefile = indicefile +1 |
|---|
| 147 | |
|---|
| 148 | for i in range(len(outputfilename)): |
|---|
| 149 | output_name = separatenames(outputfilename) |
|---|
| 150 | |
|---|
| 151 | return output_name |
|---|
| 152 | |
|---|
| 153 | def call_streamfunction ( input_name = None): |
|---|
| 154 | import numpy as np |
|---|
| 155 | from myplot import separatenames |
|---|
| 156 | from os import system |
|---|
| 157 | system("rm -f streamfunction.auto.def") |
|---|
| 158 | system("touch streamfunction.auto.def") |
|---|
| 159 | indicefile=0 |
|---|
| 160 | outputfilename="" |
|---|
| 161 | append="_stream" |
|---|
| 162 | f = open('streamfunction.auto.def', 'w') |
|---|
| 163 | for zfile in input_name: |
|---|
| 164 | f.write(zfile+"\n") |
|---|
| 165 | f.close() |
|---|
| 166 | system("streamfunction.e < streamfunction.auto.def") |
|---|
| 167 | |
|---|
| 168 | if indicefile is 0: |
|---|
| 169 | outputfilename=input_name[indicefile][0:len(input_name[0])-3]+append+".nc" |
|---|
| 170 | else: |
|---|
| 171 | outputfilename=outputfilename+","+input_name[indicefile][0:len(input_name[0])-3]+append+".nc" |
|---|
| 172 | indicefile = indicefile +1 |
|---|
| 173 | |
|---|
| 174 | for i in range(len(outputfilename)): |
|---|
| 175 | output_name = separatenames(outputfilename) |
|---|
| 176 | |
|---|
| 177 | return output_name |
|---|
| 178 | |
|---|
