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