1 | ; |
---|
2 | ; Macros definitions (leave blanks around each code) |
---|
3 | ; |
---|
4 | ; [ ]macro[ = ]sign_1[ ]field_1[ ].... [ ]<legend>[ ][unit] |
---|
5 | ; or |
---|
6 | ; [ ]macro[ =: ]<name of idl function to call>[ ]<legend>[ ][unit] |
---|
7 | |
---|
8 | ; WORK |
---|
9 | |
---|
10 | soellips =: make_work <Ellipse> [-] |
---|
11 | |
---|
12 | ; HADAM |
---|
13 | totprep = lsrain + lssnow + cvrain + cvsnow <E-P> [kg m-2 s-1] |
---|
14 | -e-mp = evap - lsrain - lssnow - cvrain - cvsnow <E-P> [kg m-2 s-1] |
---|
15 | nrads = solar + longwave - sh - 2.501e6 * evap <Net surface heat flux> [W m-2] |
---|
16 | toabud = field200 - field201 - olr <Net Radiation at TOA> [W m-2] |
---|
17 | totcld = field31 + field32 + field33 <Total Could Amount> [ ] |
---|
18 | |
---|
19 | ; ECHAM4 |
---|
20 | -e-mp = - evap - tprep <E-P> [m/s] |
---|
21 | e-mp = - qevap - aprl - aprc - aprs <E-P> [m/s] |
---|
22 | emp = - evap - tprc <E-P> [m/s] |
---|
23 | -n-rads = ahfs + ahfl + srads + trads <Net surface heat flux> [W/m**2] |
---|
24 | curltau =: make_curltau <Wind Stress Curl> [N/m] |
---|
25 | |
---|
26 | ; OPA |
---|
27 | vozonbsf =: make_msf <Meridional Stream Function> [Sv] |
---|
28 | votkeavt =: make_votkeavt <Vertical Eddy Diffusivity> [m2/s] |
---|
29 | vodenpot =: make_eos <Potential Density> [kg/m3] |
---|
30 | vocrtmod =: make_crtm <Current Module> [m/s] |
---|
31 | vovem300 =: make_vem300 <Mean Vertical Velocity 300 m> [m/s] |
---|
32 | vozom300 =: make_zom300 <Mean Zonal Velocity 300 m> [m/s] |
---|
33 | vomem300 =: make_mem300 <Mean Meridional Velocity 300 m> [m/s] |
---|
34 | |
---|
35 | vodeptht =: make_depth <Depth> [m] |
---|
36 | vosigthi =: make_depth <Sigma Thickness> [m] |
---|
37 | vosigvol =: make_depth <Sigma Volume> [m3] |
---|
38 | |
---|
39 | vozotrax =: make_xtran <Zonal Transport> [Sv] |
---|
40 | vometray =: make_ytran <Meridional Transport> [Sv] |
---|
41 | vovetraz =: make_ztran <Vertical Transport> [Sv] |
---|
42 | |
---|
43 | ; GRADIENTS |
---|
44 | |
---|
45 | sosstgrady =: make_grad <Meridional SST Gradients> [C.(100km)-1] \sosstsst, y/ |
---|
46 | sosstgradx =: make_grad <Zonal SST Gradients> [C.(100km)-1] \sosstsst, x/ |
---|
47 | tauxgradx =: make_grad <Zonal Taux Gradients> [Pa.(m)-1] \sozotaux, x/ |
---|
48 | tauxgrady =: make_grad <Meridional Taux Gradients> [Pa.(m)-1] \sozotaux, y/ |
---|
49 | tsolgrady =: make_grad <Meridional TSOL Gradients> [C.(100km)-1] \tsol, y/ |
---|
50 | tsolgradx =: make_grad <Zonal TSOL Gradients> [C.(100km)-1] \tsol, x/ |
---|
51 | |
---|
52 | ; RATIO |
---|
53 | |
---|
54 | ratiotaux =: make_ratio <Ratio> [no units] \taux, ztauxx_SCi, ztemp_taux, ztauxx_SCi/ |
---|
55 | ratiotauy =: make_ratio <Ratio> [no units] \tauy, ztauyy_SCi, ztemp_tauy, ztauyy_SCi/ |
---|
56 | ratiosflx =: make_ratio <Ratio> [no units] \soshfoce_m, soshfoce_SCi, temp_soshfoce, soshfoce_SCi/ |
---|
57 | rationsflx =: make_ratio <Ratio> [no units] \sonsfoce_m, sonsfoce_SCi, temp_sonsfoce, sonsfoce_SCi/ |
---|
58 | |
---|
59 | ; LINEAR INTERP |
---|
60 | |
---|
61 | tauxSCi =: make_interp <Interp> [Nm-2] \ztauxx_SCa, 0.533333, ztauxx_SCb, 0.466667/ |
---|
62 | tauySCi =: make_interp <Interp> [Nm-2] \ztauyy_SCa, 0.533333, ztauyy_SCb, 0.466667/ |
---|
63 | sohefoce-f-m =: make_interp <Interp> [Wm-2] \sohefldo_f, 1.0, sonsfoce_m, -1.0, soshfoce_m, -1.0/ |
---|
64 | soshfoce-f-m =: make_interp <Interp> [Wm-2] \soshfldo_f, 1.0, soshfoce_m, -1.0/ |
---|
65 | |
---|
66 | ; NCEP HEAT FLUX |
---|
67 | |
---|
68 | nshf# = nlwrs + lhtfl + shtfl <Net Down Non Solar Heat Flux at Surf> [W/m2] |
---|
69 | shf# = nswrs <Net Down Solar Heat Flux at Surf> [W/m2] |
---|
70 | thf# = nswrs + nlwrs + lhtfl + shtfl <Net Down Total Heat Flux at Surf> [W/m2] |
---|
71 | |
---|
72 | ; ERA40 HEAT FLUX |
---|
73 | nshfs = str + slhf + sshf <Net Down Non Solar Heat Flux at Surf> [Wm-2] |
---|
74 | shfs = ssr <Net Down Solar Heat Flux at Surf> [W/m2] |
---|
75 | thfs = ssr + str + slhf + sshf <Net Down Total Heat Flux at Surf> [Wm-2] |
---|
76 | |
---|
77 | ; for stats, add a \field/ = base field to work on |
---|
78 | |
---|
79 | ; Deviation Standard |
---|
80 | sstdev =: make_stddev <SST Std Dev> [C] \sst/ |
---|
81 | sosstdev =: make_stddev <SST Std Dev> [C] \sosstsst/ |
---|
82 | tsoldev =: make_stddev <TSOL Std Dev> [C] \tsol/ |
---|
83 | airdev =: make_stddev <Air Temp Std Dev> [C] \air/ |
---|
84 | sotoxdev =: make_stddev <TAUX Oce Std Dev> [Nm-2] \sozotaux/ |
---|
85 | toxdev =: make_stddev <TAUX Atm Std Dev> [Nm-2] \taux_oce/ |
---|
86 | ztxdev =: make_stddev <TAUX Oce Std Dev> [Nm-2] \ztx/ |
---|
87 | sotoydev =: make_stddev <TAUY Oce Std Dev> [Nm-2] \sometauy/ |
---|
88 | toydev =: make_stddev <TAUY Atm Std Dev> [Nm-2] \tauy_oce/ |
---|
89 | mtydev =: make_stddev <TAUY Oce Std Dev> [Nm-2] \mty/ |
---|
90 | tempdev =: make_stddev <SST Std Dev> [C] \temp/ |
---|
91 | tsdev =: make_stddev <SST Std Dev> [C] \ts/ |
---|
92 | tosdev =: make_stddev <SST Std Dev> [C] \tos/ |
---|
93 | tasdev =: make_stddev <TAS Std Dev> [C] \tas/ |
---|
94 | tauudev =: make_stddev <TAUX Std Dev> [Pa] \tauu/ |
---|
95 | tauvdev =: make_stddev <TAUY Std Dev> [Pa] \tauv/ |
---|
96 | srftdev =: make_stddev <Surf. Temp. Std Deviation> [K] \srft/ |
---|
97 | airdev =: make_stddev <Air Temp Std Dev> [C] \air/ |
---|
98 | uflxdev =: make_stddev <Taux Std Dev> [Nm-2] \uflx/ |
---|
99 | vflxdev =: make_stddev <Tauy Std Dev> [Nm-2] \vflx/ |
---|
100 | nshfsdev =: make_stddev <Net Down Non Solar Heat Flx Std Dev> [Wm-2] \nshfs/ |
---|
101 | shfsdev =: make_stddev <Net Down Solar Heat Flx Std Dev> [Wm-2] \shfs/ |
---|
102 | thfsdev =: make_stddev <Net Down Total Heat Flx Std Dev> [Wm-2] \thfs/ |
---|
103 | sohefldodev =: make_stddev <Net Downwd Heat Flx Std Dev> [Wm-2] \sohefldo/ |
---|
104 | sonsfldodev =: make_stddev <Non Solar Heat Flx Std Dev> [Wm-2] \sonsfldo/ |
---|
105 | soshfldodev =: make_stddev <Solar Heat Flx Std Dev> [Wm-2] \soshfldo/ |
---|
106 | precipdev =: make_stddev <Precip Std Dev> [mm(day-1)] \precip/ |
---|
107 | topldev =: make_stddev <OLR Std Dev> [Wm-2] \topl/ |
---|
108 | |
---|
109 | ; Kurtosis et Skewness |
---|
110 | sosstskew =: make_skewness <SST Skewness> [C] \sosstsst/ |
---|
111 | sosstkurt =: make_kurtosis <SST Kurtosis> [C] \sosstsst/ |
---|
112 | tsskew =: make_skewness <SST Skewness> [C] \ts/ |
---|
113 | tasskew =: make_skewness <SST Skewness> [C] \tas/ |
---|
114 | tauuskew =: make_skewness <TAUX Atm Skewness> [C] \tauu/ |
---|
115 | tauvskew =: make_skewness <TAUY Atm Skewness> [C] \tauv/ |
---|
116 | tauxskew =: make_skewness <TAUX Oce Skewness> [C] \sozotaux/ |
---|
117 | tauyskew =: make_skewness <TAUY Oce Skewness> [C] \sometauy/ |
---|
118 | |
---|
119 | ; Ocean |
---|
120 | |
---|
121 | sowincur =: make_wcurl <Wind Stress Curl> [N/m] |
---|
122 | sowindiv =: make_wdiv <Wind Stress Divergence> [N/m] |
---|
123 | sowinmod =: make_wm <Wind Stress Module> [N/m2] |
---|
124 | so20cdep =: make_20C <Depth of 20C Isotherm> [m] |
---|
125 | sobarstf =: make_bsf <Barotropic Stream Function> [Sv] |
---|
126 | |
---|
127 | sodenflx =: make_denflx <Density Flux> [kg/m2/s] |
---|
128 | sodenflw =: make_denflw <Fresh Water Density Flux> [kg/m2/s] |
---|
129 | |
---|
130 | sonsfldo = sohefldo - soshfldo <Non Solar Heat Flux> [W/m2] |
---|
131 | sohefatm = sohefldo - sohefldp <Downward Heat Flux Atmos> [W/m2] |
---|
132 | sowatatm = sowaflup - sowafldp <Upward Water Flux Atmos> [Kg/m2/s] |
---|
133 | sohefcmp = q - qrp <Net Downward Heat Flux> [W/M2] |
---|
134 | |
---|
135 | sozonfht = sozonfha + sozonfhd <Total Heat Transport> [PW] |
---|
136 | sozanfht = sozanfha + sozanfhd <Atlantic Total Heat Transport> [PW] |
---|
137 | sozonfst = sozonfsa + sozonfsd <Total Salt Transport> [Giga.g/s] |
---|
138 | sozanfst = sozanfsa + sozanfsd <Atlantic Total Salt Transport> [Giga.g/s] |
---|
139 | sozdnfht = sozonfha + sozonfhd - sozanfha - sozanfhd <Indopac Total Heat Transport> [PW] |
---|
140 | sozdnfhg = sozonfhg - sozanfhg <Indopac Gyre Heat Transport> [PW] |
---|
141 | sozdnfho = sozonfho - sozanfho <Indopac Overturning Heat Transport> [PW] |
---|
142 | sozdnfhd = sozonfhd - sozanfhd <Indopac Diffusive Heat Transport> [PW] |
---|
143 | sozdnfst = sozonfsa + sozonfsd - sozanfsa - sozanfsd <Indopac Total Salt Transport> [Giga.g/s] |
---|
144 | sozdnfsg = sozonfsg - sozanfsg <Indopac Gyre Salt Transport> [PW] |
---|
145 | sozdnfso = sozonfso - sozanfso <Indopac Overturning Salt Transport> [PW] |
---|
146 | sozdnfsd = sozonfsd - sozanfsd <Indopac Diffusive Salt Transport> [PW] |
---|
147 | |
---|
148 | somltten = somlbtdz + somlbtdt + somlbtdv <T Total Vertical Entrainment> [C/s] |
---|
149 | somlttis = somlttdh + somlbtdv <T Total Isopycnal Diffusion> [C/s] |
---|
150 | somlttah = somlttax + somlttay <T Horizontal Advection> [C/s] |
---|
151 | somlttad = somlttax + somlttay + somlttaz <T Total Advection> [C/s] |
---|
152 | somltsen = somlbsdz + somlbsdt + somlbsdv <S Total Vertical Entrainment> [PSU/s] |
---|
153 | somltsis = somltsdh + somlbsdv <T Total Isopycnal Diffusion> [PSU/s] |
---|
154 | somltsah = somltsax + somltsay <S Horizontal Advection> [PSU/s] |
---|
155 | somltsad = somltsax + somltsay + somltsaz <S Total Advection> [PSU/s] |
---|
156 | |
---|
157 | ; ECMWF |
---|
158 | soheflec = ecmwf_146 + ecmwf_147 + ecmwf_176 + ecmwf_177 <Net Downward Heat Flux> [Wm-2] |
---|
159 | soshflec = ecmwf_176 <Solar Heat Flux> [Wm-2] |
---|
160 | sonsflec = ecmwf_146 + ecmwf_147 + ecmwf_177 <Non Solar Heat Flux> [Wm-2] |
---|
161 | sowaflec = ecmwf_182 + ecmwf_142 + ecmwf_143 <E-P> [mmj-1] |
---|
162 | |
---|
163 | ; LMDZ |
---|
164 | uwf = evap - rain <E-P> [mm/day] |
---|
165 | latent = bils - sens <Latent Heat Flux Atm> [Wm-2] |
---|
166 | |
---|