
<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>&nbsp; while ( y &lt;= 29 )<br>&nbsp;&nbsp;&nbsp; x=81<br>&nbsp;&nbsp;&nbsp; while ( x &lt;= 102 )<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'set y 'y<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'set x 'x<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 'set defval pmask 'x' 'y' 1'<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; x=x+1<br>&nbsp;&nbsp;&nbsp; say ' Domain =&nbsp; lat = 'y',&nbsp;&nbsp; lon = 'x<br>&nbsp;&nbsp;&nbsp; endwhile<br>&nbsp;&nbsp;&nbsp; y=y+1 <br>&nbsp; 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-&gt; 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>&nbsp;-57.52 -53.81 -50.1 -46.39 -42.68 -38.97 -35.26 -31.54 -27.83 -24.12<br>&nbsp;-20.41 -16.7 -12.99 -9.28 -5.57 -1.86 1.86 5.57 9.28 12.99<br>&nbsp;16.7 20.41 24.12 27.83 31.54 35.26 38.97 42.68 46.39 50.1<br>&nbsp;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&nbsp; 1&nbsp; 99&nbsp; surface temperature&nbsp; [degk]<br>endvars<br>****************************************************************<br>ga-&gt; q config&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>Config:

 v2.0.a9 little-endian readline printim grib2 netcdf hdf4-sds hdf5 

opendap-grids,stn geotiff 

shapefile&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;

 <br>Grid Analysis and Display System (GrADS) Version 2.0.a9&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>Copyright (c) 1988-2010 by Brian Doty and the&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>Institute for Global Environment and Society (IGES)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>This program is distributed WITHOUT ANY WARRANTY&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>See file COPYRIGHT for more information.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br><br>Built Fri Sep&nbsp; 3 11:29:45 MDT 2010 for i686-pc-linux-gnu<br><br>This version of GrADS has been configured with th

e following options:<br>&nbsp; o Built on a LITTLE ENDIAN machine<br>&nbsp; o Athena Widget GUI DISABLED<br>&nbsp; o Command line editing ENABLED<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://tiswww.case.edu/php/chet/readline/rltop.html<br>&nbsp; o printim command for image output ENABLED<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://www.zlib.net<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://www.libpng.org/pub/png/libpng.html<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://www.libgd.org/Main_Page<br>&nbsp; o GRIB2 interface ENABLED<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://www.ijg.org<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://www.ece.uvic.ca/~mdadams/jasper<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; g2clib-1.1.8<br>&nbsp; o NetCDF interface ENABLED<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://www.unidata.ucar.edu/software/netcdf<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; netcdf 4.1.1-rc2 of Mar 26 2010 20:33:15 $<br>&nbsp; o OPeNDAP gridded data interface ENABLED<br>&nbsp; o OPeNDAP station data interface ENABLED<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://iges.org/grads/gadoc/supplibs.html<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; libgadap 2.0<br>&nbsp; o HDF4 and HDF5 interfaces ENABLED<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://hdfgroup.org<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; HDF 4.2r3<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; HDF5 1.8.4<br>&nbsp; o GeoTIFF and KML/TIFF output ENABLED<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://www.libtiff.org<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://geotiff.osgeo.org<br>&nbsp; o KML contour output ENABLED<br>&nbsp; o Shapefile interface ENABLED<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; http://shapelib.maptools.org<br><br>For additional information please consult http://iges.org/grads                                               </div></body>

</html>