% This file contains general diagnostics parameters. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Field information % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % When ever a new field is added to be analyzed using the diagnostics % package, it must be added here so that the packages knows what type of % field it is. fields2D = {'USTR' ,'VSTR' ,'TSR' ,'OLR' ,'SSR' ,... 'SLR' ,'SHF' ,'EVAP' ,'PRECNV','PRECLS',... 'CLOUDC','CLTOP' ,'CBMF' ,'DRAG' ,'aimV0' ,... 'aimT0' ,'aimQ0' ,'EnFxPr','albedo','dTsurf',... 'EmP' ,'TotP' ,'ETA' ,'ETAstd','Eta2' ,... 'Eta' ,... 'SLP','HF' ,'QSW','TX' ,'TY' ,'FW' ,... 'SFx','SIC','MXL','SST','SSS','vSq',... 'ice_fract' ,'ice_iceH' ,'ice_snowAge' ,... 'ice_snowH' ,'ice_snowPrc' ,'ice_Qice1' ,... 'ice_Qice2' ,'ice_Tice1' ,'ice_Tice2' ,... 'ice_Tsrf' ,... 'ICE_fract' ,'ICE_iceH' ,... 'ICE_fract-T' ,'ICE_iceH-T' ,... 'ice_fract_Ave' ,'ice_iceH_Ave' ,... 'ice_snowH_Ave' ,'ice_Tsrf_Ave' ,'ice_Tice1_Ave' ,... 'ice_Tice2_Ave' ,'ice_snowPr_Ave','ice_flx2oc_Ave',... 'ice_frw2oc_Ave','ice_salFx_Ave' ,'ice_flxAtm_Ave',... 'ice_frwAtm_Ave','ice_albedo_Ave','ICE_tMxL_Ave' ,... 'ICE_sMxL_Ave'}; fields3D = {'S','T','Temp','TT','U','UU','uVel','V','VV','vVel','W',... 'wVel','aim_RH','Tstd','KEpri','phiHyd','Psi','Bol','Res',... 'Conv','Test','ActT','ThetaEc'}; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % AIM physics parameters % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % This information is specific to the aim atmospheric physics package. The % data (unfortunately) comes in a single file with an array that must be % properly indexed to pull out the desired information. Here is the list % of fields and the indecies. Note that some fields do not have indecies % such as 'EmP' and 'TotP'. This is because these are compositions of % other aim physics fields. aimparameters = {'USTR' ,'VSTR' ,'TSR' ,'OLR' ,'SSR' ,... 'SLR' ,'SHF' ,'EVAP' ,'PRECNV','PRECLS',... 'CLOUDC','CLTOP' ,'CBMF' ,'DRAG' ,'aimV0' ,... 'aimT0' ,'aimQ0' ,'EnFxPr','albedo','dTsurf',... 'EmP' ,'TotP'}; USTR = 1; % U-stress (Pa) VSTR = 2; % V-stress (Pa) TSR = 3; % Top-of-atm. shortwave radiation (W/m2) OLR = 4; % Outgoing longwave radiation (W/m2) SSR = 5; % Surface shortwave radiation (W/m2) SLR = 6; % Surface longwave radiation (W/m2) SHF = 7; % Sensible heat flux (W/m2) EVAP = 8; % Evaporation [g/ (m2 s)] PRECNV = 9; % Convective precipitation [g/ (m2 s)] PRECLS = 10; % Large-scale precipitation [g/ (m2 s)] CLOUDC = 11; % Total cloud cover (fraction) CLTOP = 12; % Normalized pressure at cloud top CBMF = 13; % Cloud-base mass flux DRAG = 14; % Surface Drag term (= Cd*Rho*|V|) (land+sea combined) aimV0 = 15; % Surface wind speed [m/s] aimT0 = 16; % Surface air absolute temp. [K] aimQ0 = 17; % Surface air spec. humidity [g/kg] EnFxPr = 18; % Precip energy flux (snow, rain temp) [W/m^2] albedo = 19; % Surface albedo [0-1] dTsurf = 20; % Surf. Temp change from 1 iter to the next one (>0) [K] %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Coupled fields parameters % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% cplparameters = {'SLP','HF' ,'QSW','TX' ,'TY' ,'FW' ,... 'SFx','SIC','MXL','SST','SSS','vSq'}; SLP = 1; % Atmospheric Sea-Level pressure [N/m2] HF = 2; % Net surface heat-flux [W/m2] QSW = 3; % ??? [W/m2] TX = 4; % Surface stress in X [N/m2] TY = 5; % Surface stress in Y [N/m2] FW = 6; % Net fresh water flux (E-P-R) [m/s] SFx = 7; % Salt flux from sea-ice [psu.kg/m^2/s] SIC = 8; % Sea-ice mass [kg/m^2] MXL = 9; % Ocean mixed-layer depth [m] SST = 10; % Ocean surface temperature [degC] SSS = 11; % Ocean surface salinity [psu] vSq = 12; % Ocean surface velocity square [m^2/s^2] %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Ice time-average fields % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% iceparameters = {'ice_fract' ,'ice_iceH' ,'ice_snowAge' ,... 'ice_snowH' ,'ice_snowPrc' ,'ice_Qice1' ,... 'ice_Qice2' ,'ice_Tice1' ,'ice_Tice2' ,... 'ice_Tsrf' ,... 'ICE_fract' ,'ICE_iceH' ,... 'ICE_fract-T' ,'ICE_iceH-T' ,... 'ice_fract_Ave' ,'ice_iceH_Ave' ,... 'ice_snowH_Ave' ,'ice_Tsrf_Ave' ,'ice_Tice1_Ave' ,... 'ice_Tice2_Ave' ,'ice_snowPr_Ave','ice_flx2oc_Ave',... 'ice_frw2oc_Ave','ice_salFx_Ave' ,'ice_flxAtm_Ave',... 'ice_frwAtm_Ave','ice_albedo_Ave','ICE_tMxL_Ave' ,... 'ICE_sMxL_Ave'}; ice_fract_Ave = 1; % Ice fraction [0-1] ice_iceH_Ave = 2; % Ice thickness [m] ice_snowH_Ave = 3; % Snow thickness [m] ice_Tsrf_Ave = 4; % Surface Temperature [C] ice_Tice1_Ave = 5; % 1st level Temp. [C] ice_Tice2_Ave = 6; % 2nd level Temp. [C] ice_snowPr_Ave = 7; % Snow precipitation (+=down) [kg/m2/s] ice_flx2oc_Ave = 8; % Heat flux out of the ocean (+=up) [W/m2] ice_frw2oc_Ave = 9; % Fresh-water flux out off the ocean (E-P) [m/s] ice_salFx_Ave = 10; % Salt flux out of the ocean (+=up) [psu.kg/m2] ice_flxAtm_Ave = 11; % Net heat flux from Atmosphere (+=down) [W/m2] ice_frwAtm_Ave = 12; % Fresh-water flux from Atmos. (+=up) [kg/m2/s] ice_albedo_Ave = 13; % Sea-ice albedo [0-1] ICE_tMxL_Ave = 14; % Ocean mixed-layer Temp. [C] ICE_sMxL_Ave = 15; % Ocean mixed-layer salinity [psu] %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Experiment and configuration indecies % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Current experiment order with indecies defined: % { FieldName (fln), Experiment (trl), DataType (dat),... % DataDir (dad), GridDir (grd), Iterations (itr),... % TimeStep (tst), Fluid (flu), DataFormat (drm),... % Averaging (avg), SliceType (slc), PlotStyle (pst),... % Optional parameters} % % Field descriptions/options: % FieldName (fln): Variable name. % Experiment (trl): Experiment name. % DataType (dat): Data type, supports 'Tav', 'Ins'. % DataDir (dad): Data directory. % GridDir (grd): Grid data directory. % Iterations (itr): Iterations for analysis. % TimeStep (tst): Time step. % Fluid (flu): Fluid, supports 'A', 'O'. % DataFormat (dfm): Data format, supports 'MDS', 'MNC'. % Averaging (avg): 'Jan', 'Feb', 'Mar', 'Apr', 'May', 'Jun', % 'Jul', 'Aug', 'Sep', 'Oct', 'Nov', 'Dec', % 'DJF', 'MAM', 'JJA', 'SON', 'Ann', 'Ins', % SliceType (slc): 'Sur', 'Zon', 'i=#', 'j=#', 'k=#' % PlotStyle (pst): 'Grd', 'Int', 'Con', 'Cnf', 'Lin' ifln = 1; itrl = 2; idat = 3; idad = 4; igrd = 5; iitr = 6; itst = 7; iflu = 8; iddf = 9; igdf = 10; iavg = 11; islc = 12; ipst = 13; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % General grid parameters % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Often time cube sphere data must be interpolated onto a lon-lat grid. % These variables define the longitude (XL) and latitude (YL) values onto % which the grid sphere data is interpolated. XL = -179:2:180; YL = -89:2:90 ; % Zonal averaging parameters. Currently the zonal average function is set % to compute the grid information everytime. Though it only needs to do % this once (and then access the zonal average grid information file), it % takes little time, so it is currently set to do this every time. To set % it to do this only once, specify: kpr=3; nBas=-2; kwr=0; kpr=1; nBas=-2; kwr=1; % Resolution parameters. KEEP WITH OTHER GRID INFORMATION?!?!? hres = 32; faces = 6; % Constants g = 9.81; % Acceleration due to gravity [m/s2]. presrefA = 100000; % Atmospheric reference pressure [Pa]. RdA = 287.05; % Gas constant for dry air [J/kg/K]. RvA = 461; % Gas constant for ?moist? air [J/kg/K]. cpA = 1005; % Atmos. spec. heat at constant pres. [J/kg/K]. LHvapA = 2501000; % Atmos. latent heat of vaporization. A_CC = 611; % Clausius-Clapeyron A coefficient [Pa]. Beta_CC = 0.067; % Clausius-Clapeyron Beta coefficient [1/C]. % Conversions K2C = 273.15; % Kelvin = K2C + Celsius [K or C].