[gradsusr] Sub: Mask out error
Adv
advita6 at gmail.com
Fri Feb 6 17:34:31 EST 2015
Hi,
>
Can someone help me to mask file in grads....
Thank you ....
> Out put of ncdump -c of my netcdf file
> ncdump -c pr_Amon_EC-EARTH_decadal1960_r3i1p1_196011-196012.nc
> netcdf pr_Amon_EC-EARTH_decadal1960_r3i1p1_196011-196012 {
> dimensions:
> time = UNLIMITED ; // (2 currently)
> lat = 160 ;
> lon = 320 ;
> bnds = 2 ;
> variables:
> double time(time) ;
> time:bounds = "time_bnds" ;
> time:units = "days since 1960-1-1" ;
> time:calendar = "gregorian" ;
> time:axis = "T" ;
> time:long_name = "time" ;
> time:standard_name = "time" ;
> double time_bnds(time, bnds) ;
> double lat(lat) ;
> lat:bounds = "lat_bnds" ;
> lat:units = "degrees_north" ;
> lat:axis = "Y" ;
> lat:long_name = "latitude" ;
> lat:standard_name = "latitude" ;
> double lat_bnds(lat, bnds) ;
> double lon(lon) ;
> lon:bounds = "lon_bnds" ;
> lon:units = "degrees_east" ;
> lon:axis = "X" ;
> lon:long_name = "longitude" ;
> lon:standard_name = "longitude" ;
> double lon_bnds(lon, bnds) ;
> float pr(time, lat, lon) ;
> pr:standard_name = "precipitation_flux" ;
> pr:long_name = "Precipitation" ;
> pr:comment = "at surface; includes both liquid and solid phases
> from all types of clouds (both large-scale and convective)" ;
> pr:units = "kg m-2 s-1" ;
> pr:original_name = "TP=LSP+CP" ;
> pr:cell_methods = "time: mean (interval: 3 hours)" ;
> pr:cell_measures = "area: areacella" ;
> pr:history = "2012-04-02T08:18:24Z altered by CMOR: replaced
> missing value flag (1e+28) with standard missing value (1e+20).
> 2012-04-02T08:18:24Z altered by CMOR: Converted type from \'d\' to \'f\'.
> 2012-04-02T08:18:24Z altered by CMOR: Inverted axis: lat." ;
> pr:missing_value = 1.e+20f ;
> pr:_FillValue = 1.e+20f ;
> pr:associated_files = "baseURL:
> http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile:
> gridspec_atmos_fx_EC-Earth_decadal1960_r0i0p0.nc areacella:
> areacella_fx_EC-Earth_decadal1960_r0i0p0.nc" ;
>
> // global attributes:
> :institution = "EC-Earth (European Earth System Model)" ;
> :institute_id = "ICHEC" ;
> :experiment_id = "decadal1960" ;
> :source = "EC-Earth 2.3 (2011); atmosphere: IFS
> (cy31R1+modifications, T159L62); ocean: NEMO (version2+modifications,
> ORCA1-42lev); sea ice: LIM2; land: HTessel" ;
> :model_id = "EC-EARTH" ;
> :forcing = "GHG, Oz, SA, Sl, Vl, (GHG = CO2, N2O, CH4, CFCs)" ;
> :parent_experiment_id = "N/A" ;
> :parent_experiment_rip = "N/A" ;
> :branch_time = 0. ;
> :contact = "Alastair McKinstry <alastair.mckinstry at ichec.ie>" ;
> :history = "IC3 Decadal Experiment 2012-04-02T08:18:19Z CMOR
> rewrote data to comply with CF standards and CMIP5 requirements." ;
> :comment = "IC3 contact f.doblas-reyes at ic3.cat" ;
> :references = "Du, H., F.J. Doblas-Reyes, J. Garcia-Serrano, V.
> Guemas, Y. Soufflet, B. Wouters, 2012. Sensitivity of decadal predictions
> to the initial atmospheric and oceanic perturbations. Climate Dynamics, in
> press" ;
> :initialization_method = 1 ;
> :physics_version = 1 ;
> :tracking_id = "f6c0da05-03d9-4129-a131-d72903099305" ;
> :product = "output" ;
> :experiment = "10- or 30-year run initialized in year 1960" ;
> :frequency = "mon" ;
> :creation_date = "2012-04-02T08:18:24Z" ;
> :Conventions = "CF-1.4" ;
> :project_id = "CMIP5" ;
> :table_id = "Table Amon (26 July 2011)
> b26379e76858ab98b927917878a63d01" ;
> :title = "EC-Earth model output prepared for CMIP5 10- or 30-year
> run initialized in year 1960" ;
> :parent_experiment = "N/A" ;
> :modeling_realm = "atmos" ;
> :realization = 3 ;
> :cmor_version = "2.8.0" ;
> data:
>
> time = 320, 350.5 ;
>
> lat = -89.14152, -88.02943, -86.91077, -85.79063, -84.66992, -83.54895,
> -82.42782, -81.30659, -80.18531, -79.06398, -77.94262, -76.82124,
> -75.69984, -74.57843, -73.45701, -72.33558, -71.21414, -70.09269,
> -68.97124, -67.84978, -66.72833, -65.60686, -64.4854, -63.36393,
> -62.24246, -61.12099, -59.99952, -58.87804, -57.75657, -56.63509,
> -55.51361, -54.39214, -53.27066, -52.14917, -51.02769, -49.90621,
> -48.78473, -47.66325, -46.54176, -45.42028, -44.29879, -43.17731,
> -42.05582, -40.93434, -39.81285, -38.69137, -37.56988, -36.44839,
> -35.32691, -34.20542, -33.08393, -31.96244, -30.84096, -29.71947,
> -28.59798, -27.47649, -26.355, -25.23351, -24.11203, -22.99054,
> -21.86905, -20.74756, -19.62607, -18.50458, -17.38309, -16.2616,
> -15.14011, -14.01862, -12.89713, -11.77564, -10.65415, -9.53266,
> -8.41117, -7.28968, -6.16819, -5.0467, -3.92521, -2.80372, -1.68223,
> -0.56074, 0.56074, 1.68223, 2.80372, 3.92521, 5.0467, 6.16819,
> 7.28968,
> 8.41117, 9.53266, 10.65415, 11.77564, 12.89713, 14.01862, 15.14011,
> 16.2616, 17.38309, 18.50458, 19.62607, 20.74756, 21.86905, 22.99054,
> 24.11203, 25.23351, 26.355, 27.47649, 28.59798, 29.71947, 30.84096,
> 31.96244, 33.08393, 34.20542, 35.32691, 36.44839, 37.56988, 38.69137,
> 39.81285, 40.93434, 42.05582, 43.17731, 44.29879, 45.42028, 46.54176,
> 47.66325, 48.78473, 49.90621, 51.02769, 52.14917, 53.27066, 54.39214,
> 55.51361, 56.63509, 57.75657, 58.87804, 59.99952, 61.12099, 62.24246,
> 63.36393, 64.4854, 65.60686, 66.72833, 67.84978, 68.97124, 70.09269,
> 71.21414, 72.33558, 73.45701, 74.57843, 75.69984, 76.82124, 77.94262,
> 79.06398, 80.18531, 81.30659, 82.42782, 83.54895, 84.66992, 85.79063,
> 86.91077, 88.02943, 89.14152 ;
>
> lon = 0, 1.125, 2.25, 3.375, 4.5, 5.625, 6.75, 7.875, 9, 10.125, 11.25,
> 12.375, 13.5, 14.625, 15.75, 16.875, 18, 19.125, 20.25, 21.375, 22.5,
> 23.625, 24.75, 25.875, 27, 28.125, 29.25, 30.375, 31.5, 32.625, 33.75,
> 34.875, 36, 37.125, 38.25, 39.375, 40.5, 41.625, 42.75, 43.875, 45,
> 46.125, 47.25, 48.375, 49.5, 50.625, 51.75, 52.875, 54, 55.125, 56.25,
> 57.375, 58.5, 59.625, 60.75, 61.875, 63, 64.125, 65.25, 66.375, 67.5,
> 68.625, 69.75, 70.875, 72, 73.125, 74.25, 75.375, 76.5, 77.625, 78.75,
> 79.875, 81, 82.125, 83.25, 84.375, 85.5, 86.625, 87.75, 88.875, 90,
> 91.125, 92.25, 93.375, 94.5, 95.625, 96.75, 97.875, 99, 100.125,
> 101.25,
> 102.375, 103.5, 104.625, 105.75, 106.875, 108, 109.125, 110.25,
> 111.375,
> 112.5, 113.625, 114.75, 115.875, 117, 118.125, 119.25, 120.375, 121.5,
> 122.625, 123.75, 124.875, 126, 127.125, 128.25, 129.375, 130.5,
> 131.625,
> 132.75, 133.875, 135, 136.125, 137.25, 138.375, 139.5, 140.625,
> 141.75,
> 142.875, 144, 145.125, 146.25, 147.375, 148.5, 149.625, 150.75,
> 151.875,
> 153, 154.125, 155.25, 156.375, 157.5, 158.625, 159.75, 160.875, 162,
> 163.125, 164.25, 165.375, 166.5, 167.625, 168.75, 169.875, 171,
> 172.125,
> 173.25, 174.375, 175.5, 176.625, 177.75, 178.875, 180, 181.125,
> 182.25,
> 183.375, 184.5, 185.625, 186.75, 187.875, 189, 190.125, 191.25,
> 192.375,
> 193.5, 194.625, 195.75, 196.875, 198, 199.125, 200.25, 201.375, 202.5,
> 203.625, 204.75, 205.875, 207, 208.125, 209.25, 210.375, 211.5,
> 212.625,
> 213.75, 214.875, 216, 217.125, 218.25, 219.375, 220.5, 221.625,
> 222.75,
> 223.875, 225, 226.125, 227.25, 228.375, 229.5, 230.625, 231.75,
> 232.875,
> 234, 235.125, 236.25, 237.375, 238.5, 239.625, 240.75, 241.875, 243,
> 244.125, 245.25, 246.375, 247.5, 248.625, 249.75, 250.875, 252,
> 253.125,
> 254.25, 255.375, 256.5, 257.625, 258.75, 259.875, 261, 262.125,
> 263.25,
> 264.375, 265.5, 266.625, 267.75, 268.875, 270, 271.125, 272.25,
> 273.375,
> 274.5, 275.625, 276.75, 277.875, 279, 280.125, 281.25, 282.375, 283.5,
> 284.625, 285.75, 286.875, 288, 289.125, 290.25, 291.375, 292.5,
> 293.625,
> 294.75, 295.875, 297, 298.125, 299.25, 300.375, 301.5, 302.625,
> 303.75,
> 304.875, 306, 307.125, 308.25, 309.375, 310.5, 311.625, 312.75,
> 313.875,
> 315, 316.125, 317.25, 318.375, 319.5, 320.625, 321.75, 322.875, 324,
> 325.125, 326.25, 327.375, 328.5, 329.625, 330.75, 331.875, 333,
> 334.125,
> 335.25, 336.375, 337.5, 338.625, 339.75, 340.875, 342, 343.125,
> 344.25,
> 345.375, 346.5, 347.625, 348.75, 349.875, 351, 352.125, 353.25,
> 354.375,
> 355.5, 356.625, 357.75, 358.875 ;
> }
>
> ncdump -c of my newly created mask file....
ncdump -c mask.nc
netcdf mask {
dimensions:
time = UNLIMITED ; // (1 currently)
lat = 160 ;
lon = 320 ;
variables:
double time(time) ;
time:axis = "T" ;
time:standard_name = "time" ;
time:delta_t = "0000-00-00 00:00:00" ;
time:actual_range = 66410., 66410. ;
time:long_name = "Time" ;
time:units = "days since 1800-1-1 00:00:0.0" ;
float lat(lat) ;
lat:units = "degrees_north" ;
float lon(lon) ;
lon:units = "degrees_east" ;
float msk(time, lat, lon) ;
msk:long_name = "Land Sea Mask" ;
msk:valid_range = 0s, 1s ;
msk:units = "none" ;
msk:dataset = "NOAA Optimum Interpolation (OI) SST V2" ;
msk:var_desc = "Land-sea mask" ;
msk:level_desc = "Surface" ;
msk:statistic = "Other" ;
msk:add_offset = 0.f ;
msk:scale_factor = 1.f ;
msk:missing_value = 32766s ;
msk:_FillValue = 32766.f ;
// global attributes:
:title = "msk" ;
data:
time = 66410 ;
lat = -89.14152, -88.02945, -86.91738, -85.80531, -84.69324, -83.58117,
-82.4691, -81.35703, -80.24496, -79.13289, -78.02082, -76.90874,
-75.79668, -74.68461, -73.57254, -72.46046, -71.3484, -70.23633,
-69.12426, -68.01219, -66.90012, -65.78805, -64.67598, -63.56391,
-62.45184, -61.33977, -60.2277, -59.11563, -58.00356, -56.89149,
-55.77942, -54.66735, -53.55528, -52.44321, -51.33114, -50.21907,
-49.107, -47.99493, -46.88286, -45.77079, -44.65872, -43.54665,
-42.43458, -41.32251, -40.21044, -39.09837, -37.9863, -36.87423,
-35.76216, -34.65009, -33.53802, -32.42595, -31.31388, -30.20181,
-29.08974, -27.97767, -26.8656, -25.75353, -24.64146, -23.52939,
-22.41732, -21.30525, -20.19318, -19.08111, -17.96904, -16.85697,
-15.7449, -14.63283, -13.52076, -12.40869, -11.29662, -10.18455,
-9.072479, -7.960411, -6.848343, -5.736267, -4.624199, -3.512131,
-2.400063, -1.287987, -0.1759186, 0.9361496, 2.048218, 3.160286,
4.272362, 5.38443, 6.496498, 7.608566, 8.720642, 9.83271, 10.94478,
12.05685, 13.16892, 14.28099, 15.39306, 16.50513, 17.6172, 18.72927,
19.84134, 20.95341, 22.06548, 23.17755, 24.28962, 25.40169, 26.51376,
27.62583, 28.7379, 29.84997, 30.96204, 32.07411, 33.18618, 34.29825,
35.41032, 36.52238, 37.63446, 38.74653, 39.8586, 40.97066, 42.08274,
43.1948, 44.30688, 45.41895, 46.53101, 47.64309, 48.75517, 49.86723,
50.9793, 52.09136, 53.20344, 54.31551, 55.42757, 56.53965, 57.65171,
58.76379, 59.87586, 60.98792, 62.1, 63.21207, 64.32413, 65.43621,
66.54827, 67.66035, 68.77242, 69.88448, 70.99656, 72.10862, 73.2207,
74.33277, 75.44483, 76.55691, 77.66898, 78.78104, 79.89312, 81.00518,
82.11726, 83.22933, 84.34139, 85.45347, 86.56554, 87.6776 ;
lon = 0, 1.125, 2.25, 3.375, 4.5, 5.625, 6.75, 7.875, 9, 10.125, 11.25,
12.375, 13.5, 14.625, 15.75, 16.875, 18, 19.125, 20.25, 21.375, 22.5,
23.625, 24.75, 25.875, 27, 28.125, 29.25, 30.375, 31.5, 32.625, 33.75,
34.875, 36, 37.125, 38.25, 39.375, 40.5, 41.625, 42.75, 43.875, 45,
46.125, 47.25, 48.375, 49.5, 50.625, 51.75, 52.875, 54, 55.125, 56.25,
57.375, 58.5, 59.625, 60.75, 61.875, 63, 64.125, 65.25, 66.375, 67.5,
68.625, 69.75, 70.875, 72, 73.125, 74.25, 75.375, 76.5, 77.625, 78.75,
79.875, 81, 82.125, 83.25, 84.375, 85.5, 86.625, 87.75, 88.875, 90,
91.125, 92.25, 93.375, 94.5, 95.625, 96.75, 97.875, 99, 100.125,
101.25,
102.375, 103.5, 104.625, 105.75, 106.875, 108, 109.125, 110.25,
111.375,
112.5, 113.625, 114.75, 115.875, 117, 118.125, 119.25, 120.375, 121.5,
122.625, 123.75, 124.875, 126, 127.125, 128.25, 129.375, 130.5,
131.625,
132.75, 133.875, 135, 136.125, 137.25, 138.375, 139.5, 140.625, 141.75,
142.875, 144, 145.125, 146.25, 147.375, 148.5, 149.625, 150.75,
151.875,
153, 154.125, 155.25, 156.375, 157.5, 158.625, 159.75, 160.875, 162,
163.125, 164.25, 165.375, 166.5, 167.625, 168.75, 169.875, 171,
172.125,
173.25, 174.375, 175.5, 176.625, 177.75, 178.875, 180, 181.125, 182.25,
183.375, 184.5, 185.625, 186.75, 187.875, 189, 190.125, 191.25,
192.375,
193.5, 194.625, 195.75, 196.875, 198, 199.125, 200.25, 201.375, 202.5,
203.625, 204.75, 205.875, 207, 208.125, 209.25, 210.375, 211.5,
212.625,
213.75, 214.875, 216, 217.125, 218.25, 219.375, 220.5, 221.625, 222.75,
223.875, 225, 226.125, 227.25, 228.375, 229.5, 230.625, 231.75,
232.875,
234, 235.125, 236.25, 237.375, 238.5, 239.625, 240.75, 241.875, 243,
244.125, 245.25, 246.375, 247.5, 248.625, 249.75, 250.875, 252,
253.125,
254.25, 255.375, 256.5, 257.625, 258.75, 259.875, 261, 262.125, 263.25,
264.375, 265.5, 266.625, 267.75, 268.875, 270, 271.125, 272.25,
273.375,
274.5, 275.625, 276.75, 277.875, 279, 280.125, 281.25, 282.375, 283.5,
284.625, 285.75, 286.875, 288, 289.125, 290.25, 291.375, 292.5,
293.625,
294.75, 295.875, 297, 298.125, 299.25, 300.375, 301.5, 302.625, 303.75,
304.875, 306, 307.125, 308.25, 309.375, 310.5, 311.625, 312.75,
313.875,
315, 316.125, 317.25, 318.375, 319.5, 320.625, 321.75, 322.875, 324,
325.125, 326.25, 327.375, 328.5, 329.625, 330.75, 331.875, 333,
334.125,
335.25, 336.375, 337.5, 338.625, 339.75, 340.875, 342, 343.125, 344.25,
345.375, 346.5, 347.625, 348.75, 349.875, 351, 352.125, 353.25,
354.375,
355.5, 356.625, 357.75, 358.875 ;
>
> 2)Grads version i use: (GrADS) Version 2.0.a9
> 3) Linux version :
> Linux M141647-Atmos 2.6.32-431.23.3.el6.x86_64
>
> My grads script:
>
reinit
'sdfopen pr_Amon_EC-EARTH_decadal1960_r3i1p1_196101-197012.nc'
'sdfopen mask.nc'
'define mask=lterp(msk.2,pr.1)'
'd aave(maskout(pr,mask(t=1)),lon=-125,lon=-70,lat=24,lat=50)'
error i get Data Request Warning: Request is completely outside file limits
> Thank you,
> Adv
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150206/822e9afd/attachment-0001.html
More information about the gradsusr
mailing list