<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>
Dear grads users,<br><br>My dataset in T30 Gaussian grids is written from longitude 0 (x=1) to longitude 360 (X=96). Now I want to "maskout" so that I'll have data only for: lon = -60 20 and lat = -45 10 i.e. mainly for South Atlantic Ocean. So the data origin splits my box into two.<br><br>y=12<br> while ( y <= 29 )<br> x=81<br> while ( x <= 102 )<br> 'set y 'y<br> 'set x 'x<br> 'set defval pmask 'x' 'y' 1'<br> x=x+1<br> say ' Domain = lat = 'y', lon = 'x<br> endwhile<br> y=y+1 <br> endwhile<br><br>I have tried the above, changing x again and again but it's either I get only the western or the eastern segment of the box, not a full continuous box. My ctl and GrADS info are included below. <br><br>Thanks in advance for any suggestions.<br><br>Hyacinth<br><br>********************************************************************************************<br><br>ga-> q ctlinfo<br>dset noaa_anom_1950_2010_mean1950_2010.t30.grd<br>title Sea surface temperature on T30 gg<br>undef -9.99e+08<br>xdef 96 linear 0 3.75<br>ydef 48 levels -87.16 -83.47 -79.78 -76.07 -72.36 -68.65 -64.94 -61.23<br> -57.52 -53.81 -50.1 -46.39 -42.68 -38.97 -35.26 -31.54 -27.83 -24.12<br> -20.41 -16.7 -12.99 -9.28 -5.57 -1.86 1.86 5.57 9.28 12.99<br> 16.7 20.41 24.12 27.83 31.54 35.26 38.97 42.68 46.39 50.1<br> 53.81 57.52 61.23 64.94 68.65 72.36 76.07 79.78 83.47 87.16<br>zdef 1 linear 1013 1<br>tdef 732 linear 00Z01JAN1950 1mo<br>vars 1<br>ssta 1 99 surface temperature [degk]<br>endvars<br>****************************************************************<br>ga-> q config &nb
sp; <br>Config: v2.0.a9 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile <br>Grid Analysis and Display System (GrADS) Version 2.0.a9 <br>Copyright (c) 1988-2010 by Brian Doty and the <br>Institute for Global Environment and Society (IGES)
<br>This program is distributed WITHOUT ANY WARRANTY <br>See file COPYRIGHT for more information. <br><br>Built Fri Sep 3 11:29:45 MDT 2010 for i686-pc-linux-gnu<br><br>This version of GrADS has been configured with the following options:<br> o Built on a LITTLE ENDIAN machine<br> o Athena Widget GUI DISABLED<br> o Command line editing ENABLED<br> http://tiswww.case.edu/php/chet/readline/rltop.html<br> o printim command for image output ENABLED<br> http://www.zlib.net<br> http://www.libpng.org/pub/png/libpng.html<br> http://www.libgd.org/Main_Page<br> o GRIB2 interface ENABLED<br> http://www.ijg.org<br> http://www.ece.uvic.ca/~mdadams/jasper<br> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2<br> g2clib-1.1.8<br> o NetCDF interface ENABLED<br> http://www.unidata.ucar.edu/software/netcdf<br> netcdf 4.1.1-rc2 of Mar 26 2010 20:33:15 $<br> o OPeNDAP gridded data interface ENABLED<br> o OPeNDAP station data interface ENABLED<br> http:/
/iges.org/grads/gadoc/supplibs.html<br> libgadap 2.0<br> o HDF4 and HDF5 interfaces ENABLED<br> http://hdfgroup.org<br> HDF 4.2r3<br> HDF5 1.8.4<br> o GeoTIFF and KML/TIFF output ENABLED<br> http://www.libtiff.org<br> http://geotiff.osgeo.org<br> o KML contour output ENABLED<br> o Shapefile interface ENABLED<br> http://shapelib.maptools.org<br><br>For additional information please consult http://iges.org/grads<br><br><br><br><br> </div></body>
</html>