<div dir="auto">Dear Michael,<div dir="auto"><br></div><div dir="auto">Kindly check the chapter on postprocessing tools in the UserGuide of RegCM.</div><div dir="auto"><br></div><div dir="auto">Sample code from manual is below:</div><div dir="auto"><br></div><div dir="auto">nccopy -k 4 -d 9 -s <a href="http://test_001_SRF.1990060100.nc">test_001_SRF.1990060100.nc</a>  <span style="font-family:sans-serif"><a href="http://test_001_SRF.1990060100_compressed.nc">test_001_SRF.1990060100_compressed.nc</a></span></div><div dir="auto"><span style="font-family:sans-serif"><br></span></div><div dir="auto"><span style="font-family:sans-serif"><br></span></div><div dir="auto"><span style="font-family:sans-serif">Hope this is helpful.</span></div><div dir="auto"><span style="font-family:sans-serif"><br></span></div><div dir="auto"><span style="font-family:sans-serif">Regards,</span></div><div dir="auto"><span style="font-family:sans-serif"><br></span></div><div dir="auto"><span style="font-family:sans-serif">Peter.</span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 18 Jan 2018 15:28,  <<a href="mailto:ncl-talk-request@ucar.edu">ncl-talk-request@ucar.edu</a>> wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send ncl-talk mailing list submissions to<br>
        <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:ncl-talk-request@ucar.edu">ncl-talk-request@ucar.edu</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:ncl-talk-owner@ucar.edu">ncl-talk-owner@ucar.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of ncl-talk digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. Reducing file size (Michael Notaro)<br>
   2. Re: Reducing file size (Guido Cioni)<br>
   3. Re: ttest values not plotting (Sri Nandini)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Thu, 18 Jan 2018 15:02:19 +0000<br>
From: Michael Notaro <<a href="mailto:mnotaro@wisc.edu">mnotaro@wisc.edu</a>><br>
To: "<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Subject: [ncl-talk] Reducing file size<br>
Message-ID:<br>
        <<a href="mailto:DM5PR06MB3195ABB7F99D53454D04A6BFCAE80@DM5PR06MB3195.namprd06.prod.outlook.com">DM5PR06MB3195ABB7F99D53454D04<wbr>A6BFCAE80@DM5PR06MB3195.<wbr>namprd06.prod.outlook.com</a>><br>
<br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
I have about 10 TB of regional climate model output SRF files<br>
<br>
that I need to reduce in size probably to 1 TB.  For example,<br>
<br>
one monthly file, as dumped below, <a href="http://IPSL_SRF.1998010100.nc" rel="noreferrer" target="_blank">IPSL_SRF.1998010100.nc</a>,<br>
<br>
is 2.3 GB.  Any recommendations in NCL on how<br>
<br>
to effectively accomplish this task? (e.g. command<br>
<br>
to convert the contents to short, or way to<br>
<br>
compress the content, etc)<br>
<br>
Thanks, Michael<br>
<br>
<br>
<br>
[notaro@petenwell ~/processing]# ls -l -h *nc<br>
<br>
-rw-r--r-- 1 notaro notaro 2.3G Jan 18 08:35 <a href="http://IPSL_SRF.1998010100.nc" rel="noreferrer" target="_blank">IPSL_SRF.1998010100.nc</a><br>
<br>
[notaro@petenwell ~/processing]# ncdump -h <a href="http://IPSL_SRF.1998010100.nc" rel="noreferrer" target="_blank">IPSL_SRF.1998010100.nc</a><br>
<br>
netcdf IPSL_SRF.1998010100 {<br>
<br>
dimensions:<br>
<br>
jx = 217 ;<br>
<br>
iy = 141 ;<br>
<br>
kz = 28 ;<br>
<br>
time = UNLIMITED ; // (744 currently)<br>
<br>
m10 = 1 ;<br>
<br>
m2 = 1 ;<br>
<br>
soil_layer = 2 ;<br>
<br>
time_bounds = 2 ;<br>
<br>
variables:<br>
<br>
float jx(jx) ;<br>
<br>
jx:long_name = "x-coordinate in Cartesian system" ;<br>
<br>
jx:standard_name = "projection_x_coordinate" ;<br>
<br>
jx:units = "m" ;<br>
<br>
jx:axis = "X" ;<br>
<br>
jx:_CoordinateAxisType = "GeoX" ;<br>
<br>
float iy(iy) ;<br>
<br>
iy:long_name = "y-coordinate in Cartesian system" ;<br>
<br>
iy:standard_name = "projection_y_coordinate" ;<br>
<br>
iy:units = "m" ;<br>
<br>
iy:axis = "Y" ;<br>
<br>
iy:_CoordinateAxisType = "GeoY" ;<br>
<br>
float sigma(kz) ;<br>
<br>
sigma:long_name = "Sigma at half model layers" ;<br>
<br>
sigma:standard_name = "atmosphere_sigma_coordinate" ;<br>
<br>
sigma:units = "1" ;<br>
<br>
sigma:axis = "Z" ;<br>
<br>
sigma:positive = "down" ;<br>
<br>
sigma:formula_terms = "sigma: sigma ps: ps ptop: ptop" ;<br>
<br>
sigma:_CoordinateAxisType = "GeoZ" ;<br>
<br>
float ptop ;<br>
<br>
ptop:long_name = "Pressure at model top" ;<br>
<br>
ptop:standard_name = "air_pressure" ;<br>
<br>
ptop:units = "hPa" ;<br>
<br>
float xlon(iy, jx) ;<br>
<br>
xlon:long_name = "Longitude on Cross Points" ;<br>
<br>
xlon:standard_name = "longitude" ;<br>
<br>
xlon:units = "degrees_east" ;<br>
<br>
xlon:grid_mapping = "rcm_map" ;<br>
<br>
float xlat(iy, jx) ;<br>
<br>
xlat:long_name = "Latitude on Cross Points" ;<br>
<br>
xlat:standard_name = "latitude" ;<br>
<br>
xlat:units = "degrees_north" ;<br>
<br>
xlat:grid_mapping = "rcm_map" ;<br>
<br>
float mask(iy, jx) ;<br>
<br>
mask:long_name = "Land Mask" ;<br>
<br>
mask:standard_name = "land_binary_mask" ;<br>
<br>
mask:units = "1" ;<br>
<br>
mask:coordinates = "xlat xlon" ;<br>
<br>
mask:grid_mapping = "rcm_map" ;<br>
<br>
float topo(iy, jx) ;<br>
<br>
topo:long_name = "Surface Model Elevation" ;<br>
<br>
topo:standard_name = "surface_altitude" ;<br>
<br>
topo:units = "m" ;<br>
<br>
topo:coordinates = "xlat xlon" ;<br>
<br>
topo:grid_mapping = "rcm_map" ;<br>
<br>
float ps(time, iy, jx) ;<br>
<br>
ps:long_name = "Surface Pressure" ;<br>
<br>
ps:standard_name = "surface_air_pressure" ;<br>
<br>
ps:units = "hPa" ;<br>
<br>
ps:coordinates = "xlat xlon" ;<br>
<br>
ps:grid_mapping = "rcm_map" ;<br>
<br>
ps:cell_methods = "time: point" ;<br>
<br>
float drag(time, iy, jx) ;<br>
<br>
drag:long_name = "Surface drag stress coefficient in air" ;<br>
<br>
drag:standard_name = "surface_drag_coefficient_in_<wbr>air" ;<br>
<br>
drag:units = "1" ;<br>
<br>
drag:coordinates = "xlat xlon" ;<br>
<br>
drag:grid_mapping = "rcm_map" ;<br>
<br>
drag:cell_methods = "time: point" ;<br>
<br>
float ts(time, iy, jx) ;<br>
<br>
ts:long_name = "Ground surface temperature" ;<br>
<br>
ts:standard_name = "surface_temperature" ;<br>
<br>
ts:units = "K" ;<br>
<br>
ts:coordinates = "xlat xlon" ;<br>
<br>
ts:grid_mapping = "rcm_map" ;<br>
<br>
ts:cell_methods = "time: point" ;<br>
<br>
float tf(time, iy, jx) ;<br>
<br>
tf:long_name = "Foliage canopy temperature" ;<br>
<br>
tf:standard_name = "canopy_temperature" ;<br>
<br>
tf:units = "K" ;<br>
<br>
tf:coordinates = "xlat xlon" ;<br>
<br>
tf:grid_mapping = "rcm_map" ;<br>
<br>
tf:cell_methods = "time: point" ;<br>
<br>
tf:_FillValue = 1.e+20f ;<br>
<br>
float pr(time, iy, jx) ;<br>
<br>
pr:long_name = "Total precipitation flux" ;<br>
<br>
pr:standard_name = "precipitation_flux" ;<br>
<br>
pr:units = "kg m-2 s-1" ;<br>
<br>
pr:coordinates = "xlat xlon" ;<br>
<br>
pr:grid_mapping = "rcm_map" ;<br>
<br>
pr:cell_methods = "time: mean" ;<br>
<br>
float evspsbl(time, iy, jx) ;<br>
<br>
evspsbl:long_name = "Total evapotranspiration flux" ;<br>
<br>
evspsbl:standard_name = "water_evaporation_flux" ;<br>
<br>
evspsbl:units = "kg m-2 s-1" ;<br>
<br>
evspsbl:coordinates = "xlat xlon" ;<br>
<br>
evspsbl:grid_mapping = "rcm_map" ;<br>
<br>
evspsbl:cell_methods = "time: mean" ;<br>
<br>
float snv(time, iy, jx) ;<br>
<br>
snv:long_name = "Liquid water equivalent of snow thickness" ;<br>
<br>
snv:standard_name = "lwe_thickness_of_surface_<wbr>snow_amount" ;<br>
<br>
snv:units = "kg m-2" ;<br>
<br>
snv:coordinates = "xlat xlon" ;<br>
<br>
snv:grid_mapping = "rcm_map" ;<br>
<br>
snv:cell_methods = "time: mean" ;<br>
<br>
snv:_FillValue = 1.e+20f ;<br>
<br>
float hfss(time, iy, jx) ;<br>
<br>
hfss:long_name = "Sensible heat flux" ;<br>
<br>
hfss:standard_name = "surface_upward_sensible_heat_<wbr>flux" ;<br>
<br>
hfss:units = "W m-2" ;<br>
<br>
hfss:coordinates = "xlat xlon" ;<br>
<br>
hfss:grid_mapping = "rcm_map" ;<br>
<br>
hfss:cell_methods = "time: mean" ;<br>
<br>
float rsnl(time, iy, jx) ;<br>
<br>
rsnl:long_name = "Net upward longwave energy flux" ;<br>
<br>
rsnl:standard_name = "net_upward_longwave_flux_in_<wbr>air" ;<br>
<br>
rsnl:units = "W m-2" ;<br>
<br>
rsnl:coordinates = "xlat xlon" ;<br>
<br>
rsnl:grid_mapping = "rcm_map" ;<br>
<br>
rsnl:cell_methods = "time: mean" ;<br>
<br>
float rsns(time, iy, jx) ;<br>
<br>
rsns:long_name = "Net downward shortwave energy flux" ;<br>
<br>
rsns:standard_name = "net_downward_shortwave_flux_<wbr>in_air" ;<br>
<br>
rsns:units = "W m-2" ;<br>
<br>
rsns:coordinates = "xlat xlon" ;<br>
<br>
rsns:grid_mapping = "rcm_map" ;<br>
<br>
rsns:cell_methods = "time: mean" ;<br>
<br>
float rsdl(time, iy, jx) ;<br>
<br>
rsdl:long_name = "Surface downward longwave flux in air" ;<br>
<br>
rsdl:standard_name = "surface_downwelling_longwave_<wbr>flux_in_air" ;<br>
<br>
rsdl:units = "W m-2" ;<br>
<br>
rsdl:coordinates = "xlat xlon" ;<br>
<br>
rsdl:grid_mapping = "rcm_map" ;<br>
<br>
rsdl:cell_methods = "time: mean" ;<br>
<br>
float rsds(time, iy, jx) ;<br>
<br>
rsds:long_name = "Surface downward shortwave flux in air" ;<br>
<br>
rsds:standard_name = "surface_downwelling_<wbr>shortwave_flux_in_air" ;<br>
<br>
rsds:units = "W m-2" ;<br>
<br>
rsds:coordinates = "xlat xlon" ;<br>
<br>
rsds:grid_mapping = "rcm_map" ;<br>
<br>
rsds:cell_methods = "time: mean" ;<br>
<br>
float prc(time, iy, jx) ;<br>
<br>
prc:long_name = "Convective precipitation flux" ;<br>
<br>
prc:standard_name = "convective_rainfall_flux" ;<br>
<br>
prc:units = "kg m-2 s-1" ;<br>
<br>
prc:coordinates = "xlat xlon" ;<br>
<br>
prc:grid_mapping = "rcm_map" ;<br>
<br>
prc:cell_methods = "time: mean" ;<br>
<br>
float zmla(time, iy, jx) ;<br>
<br>
zmla:long_name = "Atmospheric Boundary Layer thickness" ;<br>
<br>
zmla:standard_name = "atmosphere_boundary_layer_<wbr>thickness" ;<br>
<br>
zmla:units = "m" ;<br>
<br>
zmla:coordinates = "xlat xlon" ;<br>
<br>
zmla:grid_mapping = "rcm_map" ;<br>
<br>
zmla:cell_methods = "time: point" ;<br>
<br>
float aldirs(time, iy, jx) ;<br>
<br>
aldirs:long_name = "Surface albedo to direct shortwave radiation" ;<br>
<br>
aldirs:standard_name = "surface_albedo_short_wave_<wbr>direct" ;<br>
<br>
aldirs:units = "1" ;<br>
<br>
aldirs:coordinates = "xlat xlon" ;<br>
<br>
aldirs:grid_mapping = "rcm_map" ;<br>
<br>
aldirs:cell_methods = "time: point" ;<br>
<br>
float aldifs(time, iy, jx) ;<br>
<br>
aldifs:long_name = "Surface albedo to diffuse shortwave radiation" ;<br>
<br>
aldifs:standard_name = "surface_albedo_short_wave_<wbr>diffuse" ;<br>
<br>
aldifs:units = "1" ;<br>
<br>
aldifs:coordinates = "xlat xlon" ;<br>
<br>
aldifs:grid_mapping = "rcm_map" ;<br>
<br>
aldifs:cell_methods = "time: point" ;<br>
<br>
float sund(time, iy, jx) ;<br>
<br>
sund:long_name = "Duration of sunshine" ;<br>
<br>
sund:standard_name = "duration_of_sunshine" ;<br>
<br>
sund:units = "s" ;<br>
<br>
sund:coordinates = "xlat xlon" ;<br>
<br>
sund:grid_mapping = "rcm_map" ;<br>
<br>
sund:cell_methods = "time: sum" ;<br>
<br>
float sndp(time, iy, jx) ;<br>
<br>
sndp:long_name = "Actual snow depth" ;<br>
<br>
sndp:standard_name = "snow_depth" ;<br>
<br>
sndp:units = "mm" ;<br>
<br>
sndp:coordinates = "xlat xlon" ;<br>
<br>
sndp:grid_mapping = "rcm_map" ;<br>
<br>
sndp:cell_methods = "time: mean" ;<br>
<br>
sndp:_FillValue = 1.e+20f ;<br>
<br>
float snfl(time, iy, jx) ;<br>
<br>
snfl:long_name = "Snowfall" ;<br>
<br>
snfl:standard_name = "snow_fall" ;<br>
<br>
snfl:units = "kg m-2 s-1" ;<br>
<br>
snfl:coordinates = "xlat xlon" ;<br>
<br>
snfl:grid_mapping = "rcm_map" ;<br>
<br>
snfl:cell_methods = "time: mean" ;<br>
<br>
float uas(time, m10, iy, jx) ;<br>
<br>
uas:long_name = "Anemometric zonal (westerly) wind component" ;<br>
<br>
uas:standard_name = "eastward_wind" ;<br>
<br>
uas:units = "m s-1" ;<br>
<br>
uas:coordinates = "xlat xlon" ;<br>
<br>
uas:grid_mapping = "rcm_map" ;<br>
<br>
uas:cell_methods = "time: point" ;<br>
<br>
float vas(time, m10, iy, jx) ;<br>
<br>
vas:long_name = "Anenometric meridional (southerly) wind component" ;<br>
<br>
vas:standard_name = "northward_wind" ;<br>
<br>
vas:units = "m s-1" ;<br>
<br>
vas:coordinates = "xlat xlon" ;<br>
<br>
vas:grid_mapping = "rcm_map" ;<br>
<br>
vas:cell_methods = "time: point" ;<br>
<br>
float tas(time, m2, iy, jx) ;<br>
<br>
tas:long_name = "Near surface air temperature" ;<br>
<br>
tas:standard_name = "air_temperature" ;<br>
<br>
tas:units = "K" ;<br>
<br>
tas:coordinates = "xlat xlon" ;<br>
<br>
tas:grid_mapping = "rcm_map" ;<br>
<br>
tas:cell_methods = "time: point" ;<br>
<br>
float qas(time, m2, iy, jx) ;<br>
<br>
qas:long_name = "Near surface air specific humidity" ;<br>
<br>
qas:standard_name = "specific_humidity" ;<br>
<br>
qas:units = "1" ;<br>
<br>
qas:coordinates = "xlat xlon" ;<br>
<br>
qas:grid_mapping = "rcm_map" ;<br>
<br>
qas:cell_methods = "time: point" ;<br>
<br>
float mrso(time, soil_layer, iy, jx) ;<br>
<br>
mrso:long_name = "Moisture content of the soil layers" ;<br>
<br>
mrso:standard_name = "soil_moisture_content_in_<wbr>layers" ;<br>
<br>
mrso:units = "kg m-2" ;<br>
<br>
mrso:coordinates = "xlat xlon" ;<br>
<br>
mrso:grid_mapping = "rcm_map" ;<br>
<br>
mrso:cell_methods = "time: point" ;<br>
<br>
mrso:_FillValue = 1.e+20f ;<br>
<br>
float mrro(time, soil_layer, iy, jx) ;<br>
<br>
mrro:long_name = "Runoff flux" ;<br>
<br>
mrro:standard_name = "runoff_flux" ;<br>
<br>
mrro:units = "kg m-2 s-1" ;<br>
<br>
mrro:coordinates = "xlat xlon" ;<br>
<br>
mrro:grid_mapping = "rcm_map" ;<br>
<br>
mrro:cell_methods = "time: mean" ;<br>
<br>
mrro:_FillValue = 1.e+20f ;<br>
<br>
float time(time) ;<br>
<br>
time:long_name = "time" ;<br>
<br>
time:standard_name = "time" ;<br>
<br>
time:units = "hours since 1949-12-01 00:00:00 UTC" ;<br>
<br>
time:calendar = "noleap" ;<br>
<br>
time:bounds = "time_bnds" ;<br>
<br>
float time_bnds(time, time_bounds) ;<br>
<br>
time_bnds:units = "hours since 1949-12-01 00:00:00 UTC" ;<br>
<br>
time_bnds:calendar = "noleap" ;<br>
<br>
char rcm_map ;<br>
<br>
rcm_map:grid_mapping_name = "lambert_conformal_conic" ;<br>
<br>
rcm_map:standard_parallel = 36., 52. ;<br>
<br>
rcm_map:longitude_of_central_<wbr>meridian = -97. ;<br>
<br>
rcm_map:latitude_of_<wbr>projection_origin = 45. ;<br>
<br>
rcm_map:_<wbr>CoordinateTransformType = "Projection" ;<br>
<br>
rcm_map:_CoordinateAxisTypes = "GeoX GeoY" ;<br>
<br>
<br>
// global attributes:<br>
<br>
:title = "ICTP Regional Climatic model V4" ;<br>
<br>
:institution = "ICTP" ;<br>
<br>
:source = "RegCM Model output file" ;<br>
<br>
:Conventions = "CF-1.4" ;<br>
<br>
:references = "<a href="http://gforge.ictp.it/gf/project/regcm" rel="noreferrer" target="_blank">http://gforge.ictp.it/gf/<wbr>project/regcm</a>" ;<br>
<br>
:model_revision = "tag 4.3.5.6" ;<br>
<br>
:history = "2015-01-25 04:40:03 : Created by RegCM RegCM Model program" ;<br>
<br>
:experiment = "IPSL" ;<br>
<br>
:projection = "LAMCON" ;<br>
<br>
:grid_size_in_meters = 25000. ;<br>
<br>
:latitude_of_projection_origin = 45. ;<br>
<br>
:longitude_of_projection_<wbr>origin = -97. ;<br>
<br>
:standard_parallel = 36., 52. ;<br>
<br>
:grid_factor = 0.696943758331507 ;<br>
<br>
:boundary_nspgx = 15 ;<br>
<br>
:boundary_nspgd = 15 ;<br>
<br>
:boundary_high_nudge = 3. ;<br>
<br>
:boundary_medium_nudge = 2. ;<br>
<br>
:boundary_low_nudge = 1. ;<br>
<br>
:model_is_restarted = "Yes" ;<br>
<br>
:model_simulation_initial_<wbr>start = "1977-06-01 00:00:00 UTC" ;<br>
<br>
:model_simulation_start = "1998-01-01 00:00:00 UTC" ;<br>
<br>
:model_simulation_end = "1999-01-01 00:00:00 UTC" ;<br>
<br>
:atmosphere_time_step_in_<wbr>seconds = 120. ;<br>
<br>
:surface_interaction_time_<wbr>step_in_seconds = 120. ;<br>
<br>
:radiation_scheme_time_step_<wbr>in_minuts = 30. ;<br>
<br>
:absorption_emission_time_<wbr>step_in_hours = 18. ;<br>
<br>
:lateral_boundary_condition_<wbr>scheme = 1 ;<br>
<br>
:boundary_layer_scheme = 1 ;<br>
<br>
:cumulus_convection_scheme = 2 ;<br>
<br>
:grell_scheme_closure = 2 ;<br>
<br>
:moisture_scheme = 1 ;<br>
<br>
:ocean_flux_scheme = 2 ;<br>
<br>
:zeng_ocean_roughness_formula = 1 ;<br>
<br>
:pressure_gradient_scheme = 0 ;<br>
<br>
:surface_emissivity_factor_<wbr>computed = 0 ;<br>
<br>
:lake_model_activated = 1 ;<br>
<br>
:chemical_aerosol_scheme_<wbr>activated = 0 ;<br>
<br>
:ipcc_scenario_code = "RF" ;<br>
<br>
:diurnal_cycle_sst_scheme = 0 ;<br>
<br>
:simple_sea_ice_scheme = 0 ;<br>
<br>
:seasonal_desert_albedo = 1 ;<br>
<br>
:convective_lwp_as_large_scale = 1 ;<br>
<br>
:rrtm_radiation_scheme_<wbr>activated = 0 ;<br>
<br>
:climatic_ozone_input_dataset = 0 ;<br>
<br>
:static_solar_constant_used = 1 ;<br>
<br>
:subex_bottom_level_with_no_<wbr>clouds = 1 ;<br>
<br>
:subex_maximum_cloud_fraction_<wbr>cover = 0.8 ;<br>
<br>
:subex_auto_conversion_rate_<wbr>for_land = 0.00025 ;<br>
<br>
:subex_auto_conversion_rate_<wbr>for_ocean = 0.00025 ;<br>
<br>
:subex_gultepe_factor_when_<wbr>rain_for_land = 0.4 ;<br>
<br>
:subex_gultepe_factor_when_<wbr>rain_for_ocean = 0.4 ;<br>
<br>
:subex_rh_with_fcc_one = 1.01 ;<br>
<br>
:subex_rh_threshold_for_land = 0.8 ;<br>
<br>
:subex_rh_threshold_for_ocean = 0.9 ;<br>
<br>
:subex_limit_temperature = 238. ;<br>
<br>
:subex_raindrop_evaporation_<wbr>rate = 0.0008 ;<br>
<br>
:subex_raindrop_accretion_rate = 3. ;<br>
<br>
:subex_cloud_fraction_maximum = 0.75 ;<br>
<br>
:subex_cloud_fraction_max_for_<wbr>convection = 0.25 ;<br>
<br>
:subex_cloud_liqwat_max_for_<wbr>convection = 5.e-05 ;<br>
<br>
:grell_min_shear_on_precip = 0.25 ;<br>
<br>
:grell_max_shear_on_precip = 0.5 ;<br>
<br>
:grell_min_precip_efficiency = 0.25 ;<br>
<br>
:grell_max_precip_efficiency = 0.5 ;<br>
<br>
:grell_min_precip_efficiency_o = 0.25 ;<br>
<br>
:grell_max_precip_efficiency_o = 0.5 ;<br>
<br>
:grell_min_precip_efficiency_x = 0.25 ;<br>
<br>
:grell_max_precip_efficiency_x = 0.5 ;<br>
<br>
:grell_min_shear_on_precip_on_<wbr>ocean = 0.25 ;<br>
<br>
:grell_max_shear_on_precip_on_<wbr>ocean = 0.5 ;<br>
<br>
:grell_min_precip_efficiency_<wbr>on_ocean = 0.25 ;<br>
<br>
:grell_max_precip_efficiency_<wbr>on_ocean = 0.5 ;<br>
<br>
:grell_min_precip_efficiency_<wbr>o_on_ocean = 0.25 ;<br>
<br>
:grell_max_precip_efficiency_<wbr>o_on_ocean = 0.5 ;<br>
<br>
:grell_min_precip_efficiency_<wbr>x_on_ocean = 0.25 ;<br>
<br>
:grell_max_precip_efficiency_<wbr>x_on_ocean = 0.5 ;<br>
<br>
:grell_max_depth_of_stable_<wbr>layer = 150. ;<br>
<br>
:grell_min_depth_of_cloud = 150. ;<br>
<br>
:grell_min_convective_heating = -250. ;<br>
<br>
:grell_max_convective_heating = 500. ;<br>
<br>
:grell_max_cloud_base_height = 0.4 ;<br>
<br>
:grell_FC_ABE_removal_<wbr>timescale = 30. ;<br>
<br>
:holtslag_critical_ocean_<wbr>richardson = 0.25 ;<br>
<br>
:holtslag_critical_land_<wbr>richardson = 0.25 ;<br>
<br>
}<br>
<br>
<br>
<br>
Michael Notaro<br>
Associate Director<br>
Nelson Institute Center for Climatic Research<br>
University of Wisconsin-Madison<br>
Phone: (608) 261-1503<br>
Email: <a href="mailto:mnotaro@wisc.edu">mnotaro@wisc.edu</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180118/cd13d908/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180118/cd13d908/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 18 Jan 2018 16:09:52 +0100<br>
From: Guido Cioni <<a href="mailto:guidocioni@gmail.com">guidocioni@gmail.com</a>><br>
To: Michael Notaro <<a href="mailto:mnotaro@wisc.edu">mnotaro@wisc.edu</a>><br>
Cc: "<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Reducing file size<br>
Message-ID: <<a href="mailto:1414CD3F-4265-4335-8AA4-20734B116CC2@gmail.com">1414CD3F-4265-4335-8AA4-<wbr>20734B116CC2@gmail.com</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Michael,<br>
first of all I would strongly suggest you to do that outside of NCL.<br>
If you try to read a 10 TB inside NCL and define some variables you will likely end up with a segmentation fault or extremely slow execution.<br>
You likely want to use fast-processing techniques that make use of C++/Fortran architectures, like CDO or NCO.<br>
<br>
Regarding your task, however, it is hard to give you a suggestion. If the original file was created like that it means that all variables have the correct type/units/attributes (I guess). Thus, I don't know what you can achieve by changing type, if not a strong headache in trying to read and write again the file :)<br>
<br>
Any chance you have the original GRIB file? You can try to convert it to GRIB, which will save you at most 50% of the space, but it is tricky....<br>
<br>
Anyway it would be nice to know what you want to achieve in the end. Do you need to process the data? Do a plot? You can still extract a single level, variable or even time step beforehand (CDO or NCO), and that will surely reduce the size and the time that you need to read it inside NCL.<br>
<br>
Cheers<br>
<br>
> On 18. Jan 2018, at 16:02, Michael Notaro via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
> I have about 10 TB of regional climate model output SRF files<br>
> that I need to reduce in size probably to 1 TB.  For example,<br>
> one monthly file, as dumped below, <a href="http://IPSL_SRF.1998010100.nc" rel="noreferrer" target="_blank">IPSL_SRF.1998010100.nc</a>,<br>
> is 2.3 GB.  Any recommendations in NCL on how<br>
> to effectively accomplish this task? (e.g. command<br>
> to convert the contents to short, or way to<br>
> compress the content, etc)<br>
> Thanks, Michael<br>
><br>
><br>
> [notaro@petenwell ~/processing]# ls -l -h *nc<br>
> -rw-r--r-- 1 notaro notaro 2.3G Jan 18 08:35 <a href="http://IPSL_SRF.1998010100.nc" rel="noreferrer" target="_blank">IPSL_SRF.1998010100.nc</a><br>
> [notaro@petenwell ~/processing]# ncdump -h <a href="http://IPSL_SRF.1998010100.nc" rel="noreferrer" target="_blank">IPSL_SRF.1998010100.nc</a><br>
> netcdf IPSL_SRF.1998010100 {<br>
> dimensions:<br>
> jx = 217 ;<br>
> iy = 141 ;<br>
> kz = 28 ;<br>
> time = UNLIMITED ; // (744 currently)<br>
> m10 = 1 ;<br>
> m2 = 1 ;<br>
> soil_layer = 2 ;<br>
> time_bounds = 2 ;<br>
> variables:<br>
> float jx(jx) ;<br>
> jx:long_name = "x-coordinate in Cartesian system" ;<br>
> jx:standard_name = "projection_x_coordinate" ;<br>
> jx:units = "m" ;<br>
> jx:axis = "X" ;<br>
> jx:_CoordinateAxisType = "GeoX" ;<br>
> float iy(iy) ;<br>
> iy:long_name = "y-coordinate in Cartesian system" ;<br>
> iy:standard_name = "projection_y_coordinate" ;<br>
> iy:units = "m" ;<br>
> iy:axis = "Y" ;<br>
> iy:_CoordinateAxisType = "GeoY" ;<br>
> float sigma(kz) ;<br>
> sigma:long_name = "Sigma at half model layers" ;<br>
> sigma:standard_name = "atmosphere_sigma_coordinate" ;<br>
> sigma:units = "1" ;<br>
> sigma:axis = "Z" ;<br>
> sigma:positive = "down" ;<br>
> sigma:formula_terms = "sigma: sigma ps: ps ptop: ptop" ;<br>
> sigma:_CoordinateAxisType = "GeoZ" ;<br>
> float ptop ;<br>
> ptop:long_name = "Pressure at model top" ;<br>
> ptop:standard_name = "air_pressure" ;<br>
> ptop:units = "hPa" ;<br>
> float xlon(iy, jx) ;<br>
> xlon:long_name = "Longitude on Cross Points" ;<br>
> xlon:standard_name = "longitude" ;<br>
> xlon:units = "degrees_east" ;<br>
> xlon:grid_mapping = "rcm_map" ;<br>
> float xlat(iy, jx) ;<br>
> xlat:long_name = "Latitude on Cross Points" ;<br>
> xlat:standard_name = "latitude" ;<br>
> xlat:units = "degrees_north" ;<br>
> xlat:grid_mapping = "rcm_map" ;<br>
> float mask(iy, jx) ;<br>
> mask:long_name = "Land Mask" ;<br>
> mask:standard_name = "land_binary_mask" ;<br>
> mask:units = "1" ;<br>
> mask:coordinates = "xlat xlon" ;<br>
> mask:grid_mapping = "rcm_map" ;<br>
> float topo(iy, jx) ;<br>
> topo:long_name = "Surface Model Elevation" ;<br>
> topo:standard_name = "surface_altitude" ;<br>
> topo:units = "m" ;<br>
> topo:coordinates = "xlat xlon" ;<br>
> topo:grid_mapping = "rcm_map" ;<br>
> float ps(time, iy, jx) ;<br>
> ps:long_name = "Surface Pressure" ;<br>
> ps:standard_name = "surface_air_pressure" ;<br>
> ps:units = "hPa" ;<br>
> ps:coordinates = "xlat xlon" ;<br>
> ps:grid_mapping = "rcm_map" ;<br>
> ps:cell_methods = "time: point" ;<br>
> float drag(time, iy, jx) ;<br>
> drag:long_name = "Surface drag stress coefficient in air" ;<br>
> drag:standard_name = "surface_drag_coefficient_in_<wbr>air" ;<br>
> drag:units = "1" ;<br>
> drag:coordinates = "xlat xlon" ;<br>
> drag:grid_mapping = "rcm_map" ;<br>
> drag:cell_methods = "time: point" ;<br>
> float ts(time, iy, jx) ;<br>
> ts:long_name = "Ground surface temperature" ;<br>
> ts:standard_name = "surface_temperature" ;<br>
> ts:units = "K" ;<br>
> ts:coordinates = "xlat xlon" ;<br>
> ts:grid_mapping = "rcm_map" ;<br>
> ts:cell_methods = "time: point" ;<br>
> float tf(time, iy, jx) ;<br>
> tf:long_name = "Foliage canopy temperature" ;<br>
> tf:standard_name = "canopy_temperature" ;<br>
> tf:units = "K" ;<br>
> tf:coordinates = "xlat xlon" ;<br>
> tf:grid_mapping = "rcm_map" ;<br>
> tf:cell_methods = "time: point" ;<br>
> tf:_FillValue = 1.e+20f ;<br>
> float pr(time, iy, jx) ;<br>
> pr:long_name = "Total precipitation flux" ;<br>
> pr:standard_name = "precipitation_flux" ;<br>
> pr:units = "kg m-2 s-1" ;<br>
> pr:coordinates = "xlat xlon" ;<br>
> pr:grid_mapping = "rcm_map" ;<br>
> pr:cell_methods = "time: mean" ;<br>
> float evspsbl(time, iy, jx) ;<br>
> evspsbl:long_name = "Total evapotranspiration flux" ;<br>
> evspsbl:standard_name = "water_evaporation_flux" ;<br>
> evspsbl:units = "kg m-2 s-1" ;<br>
> evspsbl:coordinates = "xlat xlon" ;<br>
> evspsbl:grid_mapping = "rcm_map" ;<br>
> evspsbl:cell_methods = "time: mean" ;<br>
> float snv(time, iy, jx) ;<br>
> snv:long_name = "Liquid water equivalent of snow thickness" ;<br>
> snv:standard_name = "lwe_thickness_of_surface_<wbr>snow_amount" ;<br>
> snv:units = "kg m-2" ;<br>
> snv:coordinates = "xlat xlon" ;<br>
> snv:grid_mapping = "rcm_map" ;<br>
> snv:cell_methods = "time: mean" ;<br>
> snv:_FillValue = 1.e+20f ;<br>
> float hfss(time, iy, jx) ;<br>
> hfss:long_name = "Sensible heat flux" ;<br>
> hfss:standard_name = "surface_upward_sensible_heat_<wbr>flux" ;<br>
> hfss:units = "W m-2" ;<br>
> hfss:coordinates = "xlat xlon" ;<br>
> hfss:grid_mapping = "rcm_map" ;<br>
> hfss:cell_methods = "time: mean" ;<br>
> float rsnl(time, iy, jx) ;<br>
> rsnl:long_name = "Net upward longwave energy flux" ;<br>
> rsnl:standard_name = "net_upward_longwave_flux_in_<wbr>air" ;<br>
> rsnl:units = "W m-2" ;<br>
> rsnl:coordinates = "xlat xlon" ;<br>
> rsnl:grid_mapping = "rcm_map" ;<br>
> rsnl:cell_methods = "time: mean" ;<br>
> float rsns(time, iy, jx) ;<br>
> rsns:long_name = "Net downward shortwave energy flux" ;<br>
> rsns:standard_name = "net_downward_shortwave_flux_<wbr>in_air" ;<br>
> rsns:units = "W m-2" ;<br>
> rsns:coordinates = "xlat xlon" ;<br>
> rsns:grid_mapping = "rcm_map" ;<br>
> rsns:cell_methods = "time: mean" ;<br>
> float rsdl(time, iy, jx) ;<br>
> rsdl:long_name = "Surface downward longwave flux in air" ;<br>
> rsdl:standard_name = "surface_downwelling_longwave_<wbr>flux_in_air" ;<br>
> rsdl:units = "W m-2" ;<br>
> rsdl:coordinates = "xlat xlon" ;<br>
> rsdl:grid_mapping = "rcm_map" ;<br>
> rsdl:cell_methods = "time: mean" ;<br>
> float rsds(time, iy, jx) ;<br>
> rsds:long_name = "Surface downward shortwave flux in air" ;<br>
> rsds:standard_name = "surface_downwelling_<wbr>shortwave_flux_in_air" ;<br>
> rsds:units = "W m-2" ;<br>
> rsds:coordinates = "xlat xlon" ;<br>
> rsds:grid_mapping = "rcm_map" ;<br>
> rsds:cell_methods = "time: mean" ;<br>
> float prc(time, iy, jx) ;<br>
> prc:long_name = "Convective precipitation flux" ;<br>
> prc:standard_name = "convective_rainfall_flux" ;<br>
> prc:units = "kg m-2 s-1" ;<br>
> prc:coordinates = "xlat xlon" ;<br>
> prc:grid_mapping = "rcm_map" ;<br>
> prc:cell_methods = "time: mean" ;<br>
> float zmla(time, iy, jx) ;<br>
> zmla:long_name = "Atmospheric Boundary Layer thickness" ;<br>
> zmla:standard_name = "atmosphere_boundary_layer_<wbr>thickness" ;<br>
> zmla:units = "m" ;<br>
> zmla:coordinates = "xlat xlon" ;<br>
> zmla:grid_mapping = "rcm_map" ;<br>
> zmla:cell_methods = "time: point" ;<br>
> float aldirs(time, iy, jx) ;<br>
> aldirs:long_name = "Surface albedo to direct shortwave radiation" ;<br>
> aldirs:standard_name = "surface_albedo_short_wave_<wbr>direct" ;<br>
> aldirs:units = "1" ;<br>
> aldirs:coordinates = "xlat xlon" ;<br>
> aldirs:grid_mapping = "rcm_map" ;<br>
> aldirs:cell_methods = "time: point" ;<br>
> float aldifs(time, iy, jx) ;<br>
> aldifs:long_name = "Surface albedo to diffuse shortwave radiation" ;<br>
> aldifs:standard_name = "surface_albedo_short_wave_<wbr>diffuse" ;<br>
> aldifs:units = "1" ;<br>
> aldifs:coordinates = "xlat xlon" ;<br>
> aldifs:grid_mapping = "rcm_map" ;<br>
> aldifs:cell_methods = "time: point" ;<br>
> float sund(time, iy, jx) ;<br>
> sund:long_name = "Duration of sunshine" ;<br>
> sund:standard_name = "duration_of_sunshine" ;<br>
> sund:units = "s" ;<br>
> sund:coordinates = "xlat xlon" ;<br>
> sund:grid_mapping = "rcm_map" ;<br>
> sund:cell_methods = "time: sum" ;<br>
> float sndp(time, iy, jx) ;<br>
> sndp:long_name = "Actual snow depth" ;<br>
> sndp:standard_name = "snow_depth" ;<br>
> sndp:units = "mm" ;<br>
> sndp:coordinates = "xlat xlon" ;<br>
> sndp:grid_mapping = "rcm_map" ;<br>
> sndp:cell_methods = "time: mean" ;<br>
> sndp:_FillValue = 1.e+20f ;<br>
> float snfl(time, iy, jx) ;<br>
> snfl:long_name = "Snowfall" ;<br>
> snfl:standard_name = "snow_fall" ;<br>
> snfl:units = "kg m-2 s-1" ;<br>
> snfl:coordinates = "xlat xlon" ;<br>
> snfl:grid_mapping = "rcm_map" ;<br>
> snfl:cell_methods = "time: mean" ;<br>
> float uas(time, m10, iy, jx) ;<br>
> uas:long_name = "Anemometric zonal (westerly) wind component" ;<br>
> uas:standard_name = "eastward_wind" ;<br>
> uas:units = "m s-1" ;<br>
> uas:coordinates = "xlat xlon" ;<br>
> uas:grid_mapping = "rcm_map" ;<br>
> uas:cell_methods = "time: point" ;<br>
> float vas(time, m10, iy, jx) ;<br>
> vas:long_name = "Anenometric meridional (southerly) wind component" ;<br>
> vas:standard_name = "northward_wind" ;<br>
> vas:units = "m s-1" ;<br>
> vas:coordinates = "xlat xlon" ;<br>
> vas:grid_mapping = "rcm_map" ;<br>
> vas:cell_methods = "time: point" ;<br>
> float tas(time, m2, iy, jx) ;<br>
> tas:long_name = "Near surface air temperature" ;<br>
> tas:standard_name = "air_temperature" ;<br>
> tas:units = "K" ;<br>
> tas:coordinates = "xlat xlon" ;<br>
> tas:grid_mapping = "rcm_map" ;<br>
> tas:cell_methods = "time: point" ;<br>
> float qas(time, m2, iy, jx) ;<br>
> qas:long_name = "Near surface air specific humidity" ;<br>
> qas:standard_name = "specific_humidity" ;<br>
> qas:units = "1" ;<br>
> qas:coordinates = "xlat xlon" ;<br>
> qas:grid_mapping = "rcm_map" ;<br>
> qas:cell_methods = "time: point" ;<br>
> float mrso(time, soil_layer, iy, jx) ;<br>
> mrso:long_name = "Moisture content of the soil layers" ;<br>
> mrso:standard_name = "soil_moisture_content_in_<wbr>layers" ;<br>
> mrso:units = "kg m-2" ;<br>
> mrso:coordinates = "xlat xlon" ;<br>
> mrso:grid_mapping = "rcm_map" ;<br>
> mrso:cell_methods = "time: point" ;<br>
> mrso:_FillValue = 1.e+20f ;<br>
> float mrro(time, soil_layer, iy, jx) ;<br>
> mrro:long_name = "Runoff flux" ;<br>
> mrro:standard_name = "runoff_flux" ;<br>
> mrro:units = "kg m-2 s-1" ;<br>
> mrro:coordinates = "xlat xlon" ;<br>
> mrro:grid_mapping = "rcm_map" ;<br>
> mrro:cell_methods = "time: mean" ;<br>
> mrro:_FillValue = 1.e+20f ;<br>
> float time(time) ;<br>
> time:long_name = "time" ;<br>
> time:standard_name = "time" ;<br>
> time:units = "hours since 1949-12-01 00:00:00 UTC" ;<br>
> time:calendar = "noleap" ;<br>
> time:bounds = "time_bnds" ;<br>
> float time_bnds(time, time_bounds) ;<br>
> time_bnds:units = "hours since 1949-12-01 00:00:00 UTC" ;<br>
> time_bnds:calendar = "noleap" ;<br>
> char rcm_map ;<br>
> rcm_map:grid_mapping_name = "lambert_conformal_conic" ;<br>
> rcm_map:standard_parallel = 36., 52. ;<br>
> rcm_map:longitude_of_central_<wbr>meridian = -97. ;<br>
> rcm_map:latitude_of_<wbr>projection_origin = 45. ;<br>
> rcm_map:_<wbr>CoordinateTransformType = "Projection" ;<br>
> rcm_map:_CoordinateAxisTypes = "GeoX GeoY" ;<br>
><br>
> // global attributes:<br>
> :title = "ICTP Regional Climatic model V4" ;<br>
> :institution = "ICTP" ;<br>
> :source = "RegCM Model output file" ;<br>
> :Conventions = "CF-1.4" ;<br>
> :references = "<a href="http://gforge.ictp.it/gf/project/regcm" rel="noreferrer" target="_blank">http://gforge.ictp.it/gf/<wbr>project/regcm</a> <<a href="http://gforge.ictp.it/gf/project/regcm" rel="noreferrer" target="_blank">http://gforge.ictp.it/gf/<wbr>project/regcm</a>>" ;<br>
> :model_revision = "tag 4.3.5.6" ;<br>
> :history = "2015-01-25 04:40:03 : Created by RegCM RegCM Model program" ;<br>
> :experiment = "IPSL" ;<br>
> :projection = "LAMCON" ;<br>
> :grid_size_in_meters = 25000. ;<br>
> :latitude_of_projection_origin = 45. ;<br>
> :longitude_of_projection_<wbr>origin = -97. ;<br>
> :standard_parallel = 36., 52. ;<br>
> :grid_factor = 0.696943758331507 ;<br>
> :boundary_nspgx = 15 ;<br>
> :boundary_nspgd = 15 ;<br>
> :boundary_high_nudge = 3. ;<br>
> :boundary_medium_nudge = 2. ;<br>
> :boundary_low_nudge = 1. ;<br>
> :model_is_restarted = "Yes" ;<br>
> :model_simulation_initial_<wbr>start = "1977-06-01 00:00:00 UTC" ;<br>
> :model_simulation_start = "1998-01-01 00:00:00 UTC" ;<br>
> :model_simulation_end = "1999-01-01 00:00:00 UTC" ;<br>
> :atmosphere_time_step_in_<wbr>seconds = 120. ;<br>
> :surface_interaction_time_<wbr>step_in_seconds = 120. ;<br>
> :radiation_scheme_time_step_<wbr>in_minuts = 30. ;<br>
> :absorption_emission_time_<wbr>step_in_hours = 18. ;<br>
> :lateral_boundary_condition_<wbr>scheme = 1 ;<br>
> :boundary_layer_scheme = 1 ;<br>
> :cumulus_convection_scheme = 2 ;<br>
> :grell_scheme_closure = 2 ;<br>
> :moisture_scheme = 1 ;<br>
> :ocean_flux_scheme = 2 ;<br>
> :zeng_ocean_roughness_formula = 1 ;<br>
> :pressure_gradient_scheme = 0 ;<br>
> :surface_emissivity_factor_<wbr>computed = 0 ;<br>
> :lake_model_activated = 1 ;<br>
> :chemical_aerosol_scheme_<wbr>activated = 0 ;<br>
> :ipcc_scenario_code = "RF" ;<br>
> :diurnal_cycle_sst_scheme = 0 ;<br>
> :simple_sea_ice_scheme = 0 ;<br>
> :seasonal_desert_albedo = 1 ;<br>
> :convective_lwp_as_large_scale = 1 ;<br>
> :rrtm_radiation_scheme_<wbr>activated = 0 ;<br>
> :climatic_ozone_input_dataset = 0 ;<br>
> :static_solar_constant_used = 1 ;<br>
> :subex_bottom_level_with_no_<wbr>clouds = 1 ;<br>
> :subex_maximum_cloud_fraction_<wbr>cover = 0.8 ;<br>
> :subex_auto_conversion_rate_<wbr>for_land = 0.00025 ;<br>
> :subex_auto_conversion_rate_<wbr>for_ocean = 0.00025 ;<br>
> :subex_gultepe_factor_when_<wbr>rain_for_land = 0.4 ;<br>
> :subex_gultepe_factor_when_<wbr>rain_for_ocean = 0.4 ;<br>
> :subex_rh_with_fcc_one = 1.01 ;<br>
> :subex_rh_threshold_for_land = 0.8 ;<br>
> :subex_rh_threshold_for_ocean = 0.9 ;<br>
> :subex_limit_temperature = 238. ;<br>
> :subex_raindrop_evaporation_<wbr>rate = 0.0008 ;<br>
> :subex_raindrop_accretion_rate = 3. ;<br>
> :subex_cloud_fraction_maximum = 0.75 ;<br>
> :subex_cloud_fraction_max_for_<wbr>convection = 0.25 ;<br>
> :subex_cloud_liqwat_max_for_<wbr>convection = 5.e-05 ;<br>
> :grell_min_shear_on_precip = 0.25 ;<br>
> :grell_max_shear_on_precip = 0.5 ;<br>
> :grell_min_precip_efficiency = 0.25 ;<br>
> :grell_max_precip_efficiency = 0.5 ;<br>
> :grell_min_precip_efficiency_o = 0.25 ;<br>
> :grell_max_precip_efficiency_o = 0.5 ;<br>
> :grell_min_precip_efficiency_x = 0.25 ;<br>
> :grell_max_precip_efficiency_x = 0.5 ;<br>
> :grell_min_shear_on_precip_on_<wbr>ocean = 0.25 ;<br>
> :grell_max_shear_on_precip_on_<wbr>ocean = 0.5 ;<br>
> :grell_min_precip_efficiency_<wbr>on_ocean = 0.25 ;<br>
> :grell_max_precip_efficiency_<wbr>on_ocean = 0.5 ;<br>
> :grell_min_precip_efficiency_<wbr>o_on_ocean = 0.25 ;<br>
> :grell_max_precip_efficiency_<wbr>o_on_ocean = 0.5 ;<br>
> :grell_min_precip_efficiency_<wbr>x_on_ocean = 0.25 ;<br>
> :grell_max_precip_efficiency_<wbr>x_on_ocean = 0.5 ;<br>
> :grell_max_depth_of_stable_<wbr>layer = 150. ;<br>
> :grell_min_depth_of_cloud = 150. ;<br>
> :grell_min_convective_heating = -250. ;<br>
> :grell_max_convective_heating = 500. ;<br>
> :grell_max_cloud_base_height = 0.4 ;<br>
> :grell_FC_ABE_removal_<wbr>timescale = 30. ;<br>
> :holtslag_critical_ocean_<wbr>richardson = 0.25 ;<br>
> :holtslag_critical_land_<wbr>richardson = 0.25 ;<br>
> }<br>
><br>
><br>
> Michael Notaro<br>
> Associate Director<br>
> Nelson Institute Center for Climatic Research<br>
> University of Wisconsin-Madison<br>
> Phone: (608) 261-1503<br>
> Email: <a href="mailto:mnotaro@wisc.edu">mnotaro@wisc.edu</a> <mailto:<a href="mailto:mnotaro@wisc.edu">mnotaro@wisc.edu</a>><br>
> ______________________________<wbr>_________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a> <mailto:<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a> <<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a>><br>
<br>
Guido Cioni<br>
<a href="http://guidocioni.altervista" rel="noreferrer" target="_blank">http://guidocioni.altervista</a> <<a href="http://guidocioni.altervista/" rel="noreferrer" target="_blank">http://guidocioni.altervista/</a><wbr>>.org<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180118/18da18c7/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180118/18da18c7/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 18 Jan 2018 15:28:33 +0000<br>
From: "Sri Nandini" <<a href="mailto:snandini@marum.de">snandini@marum.de</a>><br>
To: <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
Subject: Re: [ncl-talk] ttest values not plotting<br>
Message-ID: <<a href="mailto:0e969092ce17ab422667e2ce38b5cdfec18d21d6@marum.de">0e969092ce17ab422667e2ce38b5c<wbr>dfec18d21d6@marum.de</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hello<br>
<br>
It appears now the stippling is being added to the map, see pdf attached, however the the latlon seem to be shifted, am i wrong about this?<br>
Here is my modified script:<br>
<br>
<br>
load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/shea_util.ncl"<br>
<br>
?yrStrt = 0000<br>
?yrLast = 0100<br>
; ==============================<wbr>==============================<wbr>==<br>
; Open the file: Read only the user specified period<br>
; ==============================<wbr>==============================<wbr>==<br>
f= addfile("T2M_C.nc", "r") ;Model Control<br>
TIME?? = f->time<br>
? YYYY?? = cd_calendar(TIME,-1)/100??????<wbr>?????????? ; entire file<br>
? iYYYY? = ind(YYYY.ge.yrStrt .and. YYYY.le.yrLast)<br>
T41??? = f->TREFHT(iYYYY,:,:)<br>
;printVarSummary(T41)?????????<wbr>?????????????????????? ; (time, lat,lon)<br>
T4??? = lonFlip(T41)<br>
;printVarSummary(T4)??????????<wbr>????????????????????? ; (time, lat,lon)<br>
T4@_FillValue = -9.96921e+36?<br>
aveX = dim_avg_n_Wrap(T4, 0)?????????????????? ;; average over the 0th dim<br>
;printVarSummary(aveX)????????<wbr>??????????????????????? ; (lat,lon)<br>
varX = dim_variance_n_Wrap(T4,0) ; compute variance<br>
;printVarSummary(varX)????????<wbr>??????????????????????? ; (lat,lon)<br>
;=============================<wbr>==============================<wbr>===================<br>
<br>
f2= addfile("T2M_CC.nc", "r") ;Current Caspian<br>
TIME2? = f2->time<br>
? YYYY2?? = cd_calendar(TIME2,-1)/100?????<wbr>??????????? ; entire file<br>
? iYYYY2? = ind(YYYY2.ge.yrStrt .and. YYYY2.le.yrLast)<br>
air22?? = f2->TREFHT(iYYYY2,:,:)<br>
air2??? = lonFlip(air22)<br>
;printVarSummary(air2)????????<wbr>??????????????????????? ; (time, lat,lon)<br>
air2@_FillValue = -9.96921e+36??????????????????<wbr>???????????<br>
aveY = dim_avg_n_Wrap(air2,0)? ; average over the 0th dim<br>
;printVarSummary(aveY)????????<wbr>??????????????????????? ; (lat,lon)<br>
varY = dim_variance_n_Wrap(air2,0) ; compute variance<br>
;printVarSummary(varY)????????<wbr>??????????????????????? ; (lat,lon)<br>
;=============================<wbr>==============================<wbr>=====================<br>
;Use "ttest" to compute the probabilities? plot: ZeroNegDashLineContour and ShadeLtContour<br>
;?????????? are in shea_util.ncl.<br>
;=============================<wbr>==============================<wbr>======================<br>
? sigr = 0.05??????????????????????? ; critical sig lvl for r<br>
<br>
?xEqv = new(dimsizes(varY),typeof(<wbr>varY),varY@_FillValue)<br>
;printVarSummary(xEqv)????????<wbr>?????????????????????? ; (lat,lon)<br>
? xEqv = equiv_sample_size (T4(lat|:,lon|:,time|:), sigr,0)??<br>
?copy_VarMeta(varY,xEqv)<br>
;printVarSummary(xEqv)????????<wbr>?????????????????????? ; (lat,lon)<br>
<br>
?yEqv = new(dimsizes(varY),typeof(<wbr>varY),varY@_FillValue)<br>
;printVarSummary(yEqv)????????<wbr>??????????????????????? ; (lat,lon)<br>
? yEqv = equiv_sample_size (air2(lat|:,lon|:,time|:), sigr,0)<br>
copy_VarMeta(varY,yEqv)<br>
;printVarSummary(yEqv)????????<wbr>??????????????????????? ; (lat,lon)<br>
??????????????????<br>
;=============================<wbr>==============================<wbr>=======================<br>
?probt = new(dimsizes(varY),typeof(<wbr>varY),varY@_FillValue)<br>
<br>
?;probt = ttest(aveY,varY,xEqv,aveX,<wbr>varX,yEqv,True,False) ; values are between 0 and 1<br>
<br>
probt = 100.*(1. -ttest(aveY,varY,xEqv,aveX,<wbr>varX,yEqv,True,False)) ; all areas whose value > 95 are significant<br>
printVarSummary(probt)????????<wbr>??????? ; (lat,lon)<br>
<br>
copy_VarMeta(varY,probt)<br>
printVarSummary(probt)????????<wbr>??????? ; (lat,lon)<br>
<br>
?printMinMax(probt,0)?????????<wbr>?????????? ;? min=0?? max=97.2842<br>
;=============================<wbr>==============================<wbr>=<br>
;calculate the T2M difference<br>
;=============================<wbr>==============================<wbr>=<br>
?T2diff = new(dimsizes(aveY),typeof(<wbr>aveY),aveY@_FillValue)<br>
;printVarSummary(T2diff)??????<wbr>??????????????????????????<br>
copy_VarMeta(aveY,T2diff)<br>
;printVarSummary(T2diff)??????<wbr>??????????????????????????<br>
T2diff=aveY-aveX<br>
? ; printVarSummary(T2diff)???????<wbr>???????????????????????? ; (lat,lon)<br>
;=============================<wbr>==============================<wbr>=========<br>
;Here is a section of my code that draws a color fill plot, and then<br>
;overlays contours of statistical significance.<br>
;=============================<wbr>==============================<wbr>=========<br>
?? wks = gsn_open_wks("pdf","T2M_tt")<br>
?? gsn_define_colormap(wks,"<wbr>BlWhRe")<br>
<br>
? res = True<br>
? res@cnFillOn???????????? = True??????????????? ; turn on color<br>
? res@gsnSpreadColors????? = True??????????????? ; use full colormap<br>
<br>
? res@cnLinesOn??????????? = False?????????????? ; turn off contour lines<br>
? res@cnLevelSelectionMode = "ManualLevels"? ; set manual contour levels<br>
? res@cnMinLevelValF?????? = -2.???????????? ; set min contour level<br>
? res@cnMaxLevelValF?????? =? 2.???????????? ; set max contour level<br>
? res@cnLevelSpacingF????? =? 0.2??????????? ; set contour spacing<br>
<br>
? res@gsnDraw????????????? = False?????????? ; Do not draw plot<br>
? res@gsnFrame???????????? = False?????????? ; Do not advance frome<br>
<br>
? res@tiMainString???????? = "T2M Difference: 50yrs overlay with ttest"<br>
? ;res@gsnCenterString????? = "5% stippled"<br>
? ;res@gsnLeftString??????? = "K"<br>
???? res@cnFillMode?????? = "RasterFill"?????? ; Raster Mode<br>
?? res@cnRasterModeOn????????????<wbr>? = True????????????? ; Raster mode shows grid cells<br>
?? plot = gsn_csm_contour_map(wks,<wbr>T2diff,res)?<br>
?? plot = ZeroNegDashLineContour (plot)??????????? ; a shea_util.ncl function<br>
<br>
?aa=stat_dispersion(probt,<wbr>True)?? ; Make sure your data has values less than? or equal to 0.95 that can be plotted<br>
???<br>
? res2 = True??????????????????????????<wbr>? ; res2 probability plots<br>
? res2@gsnDraw???????????? = False?????? ; Do not draw plot<br>
? res2@gsnFrame??????????? = False?????? ; Do not advance frome<br>
? res2@cnLevelSelectionMode = "ExplicitLevels" ; set explicit cnlev<br>
? res2@cnLevels?? = (/.95/)??? ; only have 1 contour level<br>
? res2@cnInfoLabelOn?????? = False<br>
? res2@cnLinesOn?????????? = False?????? ; do not draw contour lines<br>
? res2@cnLineLabelsOn????? = False?????? ; do not draw contour labels<br>
? res2@cnFillScaleF??????? = 0.6???????? ; add extra density<br>
??????????????????????????????<wbr>?????????? ; add cyclic point<br>
? ;res@gsnAddCyclic?????????????<wbr>?? = True????????????? ; add cyclic point<br>
<br>
? plot2?? = gsn_csm_contour(wks,gsn_add_<wbr>cyclic_point(probt(:,:)), res2)<br>
? opt = True??????????????????? ; set up parameters for pattern fill<br>
? opt@gsnShadeFillType = "pattern"? ; specify pattern fill<br>
? opt@gsnShadeLow????? = 17???????? ; stipple pattern<br>
? plot2?? = gsn_contour_shade(plot2, 0.95, 3., opt) ; stipple all areas < 95 contour???????????????????????<wbr>???????<br>
? overlay (plot, plot2)<br>
<br>
? draw (plot)<br>
? frame(wks)<br>
<br>
;=============================<wbr>==============================<wbr>==========<br>
;T2diff contains the temperature difference, probt contains the ttest<br>
;values.<br>
;=============================<wbr>==============================<wbr>==========<br>
<br>
<br>
<br>
<br>
<br>
<br>
On Jan 17, 2018 9:36:03 PM, Adam Phillips wrote:<br>
<br>
> Hi Sri,> I think the issue is that you are letting NCL create the contour levels (for plot2), as you are setting?> res2@cnLevelSelectionMode = "ExplicitLevels" but you are?> not setting res2@cnLevels.?<br>
> > gsn_contour_shade uses the existing contour levels to modify the plot, it does not add them itself. I would suggest doing this:<br>
><br>
><br>
> ? stat_dispersion(probt,True)? ?;?Make sure your data has values less than? or equal to .05 that can be plotted<br>
> > ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ;?> <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/<wbr>Contributed/stat_dispersion.<wbr>shtml</a><br>
><br>
> ? res2 = True??????????????????????????<wbr>> ? ; res2 probability plots<br>
> ? res2@gsnDraw???????????? = False?????? ; Do not draw plot<br>
> ? res2@gsnFrame??????????? = False?????? ; Do not advance frome<br>
> ? res2@cnLevelSelectionMode = "ExplicitLevels" ; set explicit cnlev<br>
> ? res2@cnLevels?? = (/.05/)??? ; only have 1 contour level<br>
> ? ;res2@cnFillPatterns = (/-1,17/) ; don't fill <0.95, stipple >=0.95<br>
> ? res@cnFillOn = True<br>
><br>
> ? res2@cnInfoLabelOn?????? = False<br>
> ? res2@cnLinesOn?????????? = False?????? ; do not draw contour lines<br>
> ? res2@cnLineLabelsOn????? = False?????? ; do not draw contour labels<br>
> ? res2@cnFillScaleF??????? = 0.6???????? ; add extra density<br>
><br>
> ? plot2 = gsn_csm_contour(wks,probt,> res2)?????????????????????????<wbr>> ??????????? ; add cyclic point<br>
><br>
> ? opt = True??????????????????? ; set up parameters for pattern fill<br>
> ? opt@gsnShadeFillType = "pattern"? ; specify pattern fill<br>
> ? opt@gsnShadeLow????? = 17???????? ; stipple pattern<br>
> ? plot2?? = gsn_contour_shade(plot2, 0.05, 3., opt) ; stipple all areas <= 0.05 contour<br>
> ? overlay (plot, plot2)<br>
<br>
<br>
> > or, more simply, do not even use gsn_contour_shade:<br>
> > stat_dispersion(probt,True)? ?;?Make sure your data has values less than? or equal to .05 that can be plotted<br>
> > ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ;?> <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/<wbr>Contributed/stat_dispersion.<wbr>shtml</a><br>
><br>
> > ? res2 = True??????????????????????????<wbr>> ? ; res2 probability plots<br>
> ? res2@gsnDraw???????????? = False?????? ; Do not draw plot<br>
> ? res2@gsnFrame??????????? = False?????? ; Do not advance frome<br>
> ? res2@cnLevelSelectionMode = "ExplicitLevels" ; set explicit cnlev<br>
> ? res2@cnLevels?? = (/.05/)??? ; only have 1 contour level<br>
> ? res2@cnFillPatterns = (/17,-1/) ; don't fill <0.95, stipple >=0.95<br>
> ? res@cnFillOn = True<br>
><br>
> ? res2@cnInfoLabelOn?????? = False<br>
> ? res2@cnLinesOn?????????? = False?????? ; do not draw contour lines<br>
> ? res2@cnLineLabelsOn????? = False?????? ; do not draw contour labels<br>
> ? res2@cnFillScaleF??????? = 0.6???????? ; add extra density<br>
><br>
> ? plot2 = gsn_csm_contour(wks,probt,> res2)?????????????????????????<wbr>> ??????????? ; add cyclic point<br>
> ??> overlay (plot, plot2)<br>
><br>
<br>
> > Hope that helps. If you have any further questions please respond to the ncl-talk email list.<br>
> > Adam<br>
><br>
><br>
> On Wed, Jan 17, 2018 at 1:18 PM, Sri Nandini via ncl-talk > <> <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> >>  wrote:<br>
> > Hello<br>
> ><br>
> > I am trying to plot ttest values overlay onto the T2m differences, but it is not being plotted nor shaded region according to 95%confidence limits with ttest.<br>
> > Below is the script, i would be grateful if someone could let me know of my omission or deletion in the resources plotting attributes.<br>
> ><br>
> > ;Use NCL's named dimensions to reorder in time.<br>
> > ;Calculate the temporal means and variances using the dim_avg and dim_variance functions.<br>
> > ;Specify a critical significance level to test the lag-one auto-correlation coefficient and determine the (temporal) number of equivalent sample sizes in each grid point using equiv_sample_size.<br>
> > ;Specify a critical significance level for the ttest and test if the means are different at each grid point.<br>
> > ; This file still has to be loaded manually<br>
> > load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/shea_util.ncl"<br>
> ><br>
> > ?yrStrt = 0000<br>
> > ?yrLast = 0100<br>
> > ; ==============================<wbr>==============================<wbr>==<br>
> > ; Open the file: Read only the user specified period<br>
> > ; ==============================<wbr>==============================<wbr>==<br>
> > f= addfile("T2M_C.nc", "r") ;Model Control<br>
> > TIME?? = f->time<br>
> > ? YYYY?? = cd_calendar(TIME,-1)/100??????<wbr>?????????? ; entire file<br>
> > ? iYYYY? = ind(YYYY.ge.yrStrt .and. YYYY.le.yrLast)<br>
> > T41??? = f->TREFHT(iYYYY,:,:)<br>
> > printVarSummary(T41)??????????<wbr>????????????????????? ; (time, lat,lon)<br>
> > T4??? = lonFlip(T41)<br>
> > printVarSummary(T4)???????????<wbr>???????????????????? ; (time, lat,lon)<br>
> > T4@_FillValue = -9.96921e+36?<br>
> > aveX = dim_avg_n_Wrap(T4, 0)?????????????????? ;; average over the 0th dim<br>
> > printVarSummary(aveX)?????????<wbr>?????????????????????? ; (lat,lon)<br>
> > varX = dim_variance_n_Wrap(T4,0)?????<wbr>????????? ; compute variance<br>
> > printVarSummary(varX)?????????<wbr>?????????????????????? ; (lat,lon)<br>
> > ;=============================<wbr>==============================<wbr>===================<br>
> ><br>
> > f2= addfile("T2M_CC.nc", "r") ;Current Caspian<br>
> > TIME2? = f2->time<br>
> > ? YYYY2?? = cd_calendar(TIME2,-1)/100?????<wbr>??????????? ; entire file<br>
> > ? iYYYY2? = ind(YYYY2.ge.yrStrt .and. YYYY2.le.yrLast)<br>
> > air22?? = f2->TREFHT(iYYYY2,:,:)<br>
> > air2??? = lonFlip(air22)<br>
> > printVarSummary(air2)?????????<wbr>?????????????????????? ; (time, lat,lon)<br>
> > air2@_FillValue = -9.96921e+36??????????????????<wbr>???????????<br>
> > aveY = dim_avg_n_Wrap(air2,0)? ; average over the 0th dim<br>
> > printVarSummary(aveY)?????????<wbr>?????????????????????? ; (lat,lon)<br>
> > varY = dim_variance_n_Wrap(air2,0) ; compute variance<br>
> > printVarSummary(varY)?????????<wbr>?????????????????????? ; (lat,lon)<br>
> > ;=============================<wbr>==============================<wbr>=====================<br>
> > ;Use "ttest" to compute the probabilities<br>
> > ;=============================<wbr>==============================<wbr>======================<br>
> ><br>
> > ? sigr = 0.05??????????????????????? ; critical sig lvl for r<br>
> ><br>
> > ?xEqv = new(dimsizes(varY),typeof(<wbr>varY),varY@_FillValue)<br>
> > printVarSummary(xEqv)?????????<wbr>????????????????????? ; (lat,lon)<br>
> > ? xEqv = equiv_sample_size (T4(lat|:,lon|:,time|:), sigr,0)??<br>
> > ?copy_VarMeta(varY,xEqv)<br>
> > printVarSummary(xEqv)?????????<wbr>????????????????????? ; (lat,lon)<br>
> ><br>
> > ?yEqv = new(dimsizes(varY),typeof(<wbr>varY),varY@_FillValue)<br>
> > printVarSummary(yEqv)?????????<wbr>?????????????????????? ; (lat,lon)<br>
> > ? yEqv = equiv_sample_size (air2(lat|:,lon|:,time|:), sigr,0)<br>
> > copy_VarMeta(varY,yEqv)<br>
> > printVarSummary(yEqv)?????????<wbr>?????????????????????? ; (lat,lon)<br>
> > ??????????????????<br>
> > ;=============================<wbr>==============================<wbr>=======================<br>
> > ?probt = new(dimsizes(varY),typeof(<wbr>varY),varY@_FillValue)<br>
> ><br>
> > ?probt = ttest(aveY,varY,xEqv,aveX,<wbr>varX,yEqv,True,False) ; values are between 0 and 1<br>
> ><br>
> > ;probt = 100.*(1. -ttest(aveY,varY,xEqv,aveX,<wbr>varX,yEqv,True,False)) ; all areas whose value > 95 are significant<br>
> > printVarSummary(probt)????????<wbr>??????? ; (lat,lon)<br>
> > copy_VarMeta(varY,probt)<br>
> > printVarSummary(probt)????????<wbr>??????? ; (lat,lon)<br>
> > print(probt)?????????????????? ;check to make sure values are between 0 and 1.<br>
> > ;=============================<wbr>==============================<wbr>=<br>
> > ;calculate the T2M difference<br>
> > ;=============================<wbr>==============================<wbr>=<br>
> > ?T2diff = new(dimsizes(aveY),typeof(<wbr>aveY),aveY@_FillValue)<br>
> > printVarSummary(T2diff)???????<wbr>?????????????????????????<br>
> > copy_VarMeta(aveY,T2diff)<br>
> > printVarSummary(T2diff)???????<wbr>?????????????????????????<br>
> > T2diff=aveY-aveX<br>
> > ?? printVarSummary(T2diff)???????<wbr>???????????????????????? ; (lat,lon)<br>
> > ;=============================<wbr>==============================<wbr>=========<br>
> > ;Here is a section of my code that draws a color fill plot, and then<br>
> > ;overlays contours of statistical significance.<br>
> > ;=============================<wbr>==============================<wbr>=========<br>
> > ?? wks = gsn_open_wks("pdf","T2M_ttest"<wbr>)<br>
> > ?? gsn_define_colormap(wks,"<wbr>BlWhRe")<br>
> ><br>
> > ? res = True<br>
> > ? res@cnFillOn???????????? = True??????????????? ; turn on color<br>
> > ? res@gsnSpreadColors????? = True??????????????? ; use full colormap<br>
> ><br>
> > ? res@cnLinesOn??????????? = False?????????????? ; turn off contour lines<br>
> > ? res@cnLevelSelectionMode = "ManualLevels"? ; set manual contour levels<br>
> > ? res@cnMinLevelValF?????? = -2.???????????? ; set min contour level<br>
> > ? res@cnMaxLevelValF?????? =? 2.???????????? ; set max contour level<br>
> > ? res@cnLevelSpacingF????? =? 0.2??????????? ; set contour spacing<br>
> ><br>
> > ? res@gsnDraw????????????? = False?????????? ; Do not draw plot<br>
> > ? res@gsnFrame???????????? = False?????????? ; Do not advance frome<br>
> ><br>
> > ? res@tiMainString???????? = "T2M Difference: 50yrs overlay with ttest"<br>
> > ? res@gsnCenterString????? = "5% stippled"<br>
> > ? res@gsnLeftString??????? = "K"<br>
> ><br>
> > ?? res@cnRasterModeOn????????????<wbr>? = True????????????? ; Raster mode shows grid cells<br>
> > ?? plot = gsn_csm_contour_map_ce(wks,<wbr>T2diff,res)?<br>
> > ?? plot = ZeroNegDashLineContour (plot)<br>
> ><br>
> > ? res2 = True??????????????????????????<wbr>? ; res2 probability plots<br>
> > ? res2@gsnDraw???????????? = False?????? ; Do not draw plot<br>
> > ? res2@gsnFrame??????????? = False?????? ; Do not advance frome<br>
> > ? res2@cnLevelSelectionMode = "ExplicitLevels" ; set explicit cnlev<br>
> > ? ;res2@cnLevels?? = (/.95/)??? ; only have 1 contour level<br>
> > ? ;res2@cnFillPatterns = (/-1,17/) ; don't fill <0.95, stipple >=0.95<br>
> > ? res@cnFillOn = False<br>
> ><br>
> > ? res2@cnInfoLabelOn?????? = False<br>
> > ? res2@cnLinesOn?????????? = False?????? ; do not draw contour lines<br>
> > ? res2@cnLineLabelsOn????? = False?????? ; do not draw contour labels<br>
> > ? res2@cnFillScaleF??????? = 0.6???????? ; add extra density<br>
> ><br>
> > ? plot2 = gsn_csm_contour(wks,probt,<wbr>res2)?????????????????????????<wbr>??????????? ; add cyclic point<br>
> ><br>
> > ? opt = True??????????????????? ; set up parameters for pattern fill<br>
> > ? opt@gsnShadeFillType = "pattern"? ; specify pattern fill<br>
> > ? opt@gsnShadeLow????? = 17???????? ; stipple pattern<br>
> > ? plot2?? = gsn_contour_shade(plot2, 0.091, 3., opt) ; stipple all areas < 0.09 contour<br>
> > ? overlay (plot, plot2)<br>
> ><br>
> > ? draw (plot)<br>
> > ? frame(wks)<br>
> ><br>
<br>
> > ______________________________<wbr>_________________<br>
> ><br>
ncl-talk mailing list<br>
> > <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
> ><br>
List instructions, subscriber options, unsubscribe:<br>
> > <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
> ><br>
<br>
><br>
<br>
> --<br>
> Adam Phillips<br>
<br>
> > Associate Scientist,? > Climate and Global Dynamics Laboratory, NCAR<br>
<br>
><br>
> > <a href="http://www.cgd.ucar.edu/staff/asphilli/" rel="noreferrer" target="_blank">www.cgd.ucar.edu/staff/<wbr>asphilli/</a>> ?? > 303-497-1726<br>
><br>
<br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
> ><br>
><br>
<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180118/26db6d69/attachment.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180118/26db6d69/<wbr>attachment.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: T2M_tt.pdf<br>
Type: application/pdf<br>
Size: 171332 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180118/26db6d69/attachment.pdf" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180118/26db6d69/<wbr>attachment.pdf</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br>
<br>
------------------------------<br>
<br>
End of ncl-talk Digest, Vol 170, Issue 32<br>
******************************<wbr>***********<br>
</blockquote></div><br></div>