22 |
'ice_Qice2' ,'ice_Tice1' ,'ice_Tice2' ,... |
'ice_Qice2' ,'ice_Tice1' ,'ice_Tice2' ,... |
23 |
'ice_Tsrf' ,... |
'ice_Tsrf' ,... |
24 |
'ICE_fract' ,'ICE_iceH' ,... |
'ICE_fract' ,'ICE_iceH' ,... |
25 |
'ICE_fract-T' ,'ICE_iceH-T' ,... |
'ICE_fract-T' ,'ICE_iceH-T' ,... |
26 |
'ice_fract_Ave' ,'ice_iceH_Ave' ,... |
'ice_fract_Ave' ,'ice_iceH_Ave' ,... |
27 |
'ice_snowH_Ave' ,'ice_Tsrf_Ave' ,'ice_Tice1_Ave' ,... |
'ice_snowH_Ave' ,'ice_Tsrf_Ave' ,'ice_Tice1_Ave' ,... |
28 |
'ice_Tice2_Ave' ,'ice_snowPr_Ave','ice_flx2oc_Ave',... |
'ice_Tice2_Ave' ,'ice_snowPr_Ave','ice_flx2oc_Ave',... |
32 |
|
|
33 |
fields3D = {'S','T','Temp','TT','U','UU','uVel','V','VV','vVel','W',... |
fields3D = {'S','T','Temp','TT','U','UU','uVel','V','VV','vVel','W',... |
34 |
'wVel','aim_RH','Tstd','KEpri','phiHyd','Psi','Bol','Res',... |
'wVel','aim_RH','Tstd','KEpri','phiHyd','Psi','Bol','Res',... |
35 |
'Conv','Test'}; |
'Conv','Test','ActT','MoiPTc'}; |
36 |
|
|
37 |
|
|
38 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
133 |
% Experiment and configuration indecies % |
% Experiment and configuration indecies % |
134 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
135 |
|
|
|
% REWRITE THE FOLLOWING PARAGRAPH!!! |
|
|
% In the incoming data from the user (fields), each entry (field) will have |
|
|
% an experiment cell array and an configuration cell array that must be in |
|
|
% the proper order. The following defines the order in which the |
|
|
% parameters are passed for each of the cell arrays. The field array had |
|
|
% the following structure: |
|
|
% |
|
|
% {Experiment1,... {Configuration1,... |
|
|
% field = { Experiment2,... , Configuration2,... , TitleStr } |
|
|
% ComparisonType} Configuration3} |
|
|
% |
|
|
% Note that while often there can be any number of experiments and |
|
|
% configurations, there are certain types of comparisons that set some |
|
|
% limitations on this. The Experiment# and Configuration# arrays have the |
|
|
% following order of fields: |
|
|
% |
|
|
% SHOULD BE GOOD FROM HERE!!! |
|
136 |
% Current experiment order with indecies defined: |
% Current experiment order with indecies defined: |
137 |
% { FieldName (fln), Experiment (exp), DataType (dat),... |
% { FieldName (fln), Experiment (trl), DataType (dat),... |
138 |
% DataDir (dad), GridDir (grd), Iterations (itr),... |
% DataDir (dad), GridDir (grd), Iterations (itr),... |
139 |
% TimeStep (tst), Fluid (flu), DataFormat (drm),... |
% TimeStep (tst), Fluid (flu), DataFormat (drm),... |
140 |
% Averaging (avg), SliceType (slc), PlotStyle (pst)} |
% Averaging (avg), SliceType (slc), PlotStyle (pst),... |
141 |
|
% Optional parameters} |
142 |
% |
% |
143 |
% Field descriptions/options: |
% Field descriptions/options: |
144 |
% FieldName (fln): Variable name. |
% FieldName (fln): Variable name. |
145 |
% Experiment (exp): Experiment name. |
% Experiment (trl): Experiment name. |
146 |
% DataType (dat): Data type, supports 'Tav', 'Int'. |
% DataType (dat): Data type, supports 'Tav', 'Ins'. |
147 |
% DataDir (dad): Data directory. |
% DataDir (dad): Data directory. |
148 |
% GridDir (grd): Grid data directory. |
% GridDir (grd): Grid data directory. |
149 |
% Iterations (itr): Iterations for analysis. |
% Iterations (itr): Iterations for analysis. |
156 |
% SliceType (slc): 'Sur', 'Zon', 'i=#', 'j=#', 'k=#' |
% SliceType (slc): 'Sur', 'Zon', 'i=#', 'j=#', 'k=#' |
157 |
% PlotStyle (pst): 'Grd', 'Int', 'Con', 'Cnf', 'Lin' |
% PlotStyle (pst): 'Grd', 'Int', 'Con', 'Cnf', 'Lin' |
158 |
ifln = 1; |
ifln = 1; |
159 |
iexp = 2; |
itrl = 2; |
160 |
idat = 3; |
idat = 3; |
161 |
idad = 4; |
idad = 4; |
162 |
igrd = 5; |
igrd = 5; |
189 |
|
|
190 |
% Resolution parameters. KEEP WITH OTHER GRID INFORMATION?!?!? |
% Resolution parameters. KEEP WITH OTHER GRID INFORMATION?!?!? |
191 |
hres = 32; |
hres = 32; |
|
vres = 5; |
|
192 |
faces = 6; |
faces = 6; |
193 |
|
|
|
g = 9.81; |
|
194 |
|
% Constants |
195 |
|
g = 9.81; % Acceleration due to gravity [m/s2]. |
196 |
|
presrefA = 100000; % Atmospheric reference pressure [Pa]. |
197 |
|
RdA = 287.05; % Gas constant for dry air [J/kg/K]. |
198 |
|
RvA = 461; % Gas constant for ?moist? air [J/kg/K]. |
199 |
|
cpA = 1005; % Atmos. spec. heat at constant pres. [J/kg/K]. |
200 |
|
LHvapA = 2501000; % Atmos. latent heat of vaporization. |
201 |
|
A_CC = 611; % Clausius-Clapeyron A coefficient [Pa]. |
202 |
|
Beta_CC = 0.067; % Clausius-Clapeyron Beta coefficient [1/C]. |
203 |
|
|
204 |
|
% Conversions |
205 |
|
K2C = 273.15; % Kelvin = K2C + Celsius [K or C]. |