Re: reading data with different increment

Hi Ahmed,
The plot looks good. As you found out, you can either use lonFlip, or
modify your input coordinate values (j1,h2,i1,i2) so that they are in the
same range as your data's coordinate variables (lat,lon)... Adam

> thanks
> when i used the function longFlib i got the following message
> lonflip: longitude dimension size must be even: mlon=541
> i donot what i have to do with this function have i to cut one point of
> the
> data.
>
> but when i used the following command it works fine
> elv = f->z({j1-90:j2:20},{i1-180:i2:20})
> i make it add the missed part of the coordinate.
> i have attached the image produced .
>
> On Thu, Oct 6, 2011 at 10:38 PM, Adam Phillips <asphilli@ucar.edu> wrote:
>
>> Hi Ahmed,
>> I should have added this to my initial reply:
>> Make sure to check your data before and after you subscript, to make
>> sure you selected the data you wanted. Note also that the orientation of
>> your data comes into play here. If, say, your longitudes go from
>> -180:180, then coordinate subscripting that you want {0:240} will only
>> grab data from 0:180... (I am guessing that this is happening in your
>> case.)
>>
>> So, look at the latitude/longitude span of your data, and make sure your
>> specified coordinate subscripts are in the same range. If your
>> longitudes go from -180:180, and you want them to go from 0:360, then
>> you can use the function lonFlip:
>> http://www.ncl.ucar.edu/Document/Functions/Contributed/lonFlip.shtml
>> Adam
>>
>>
>> On 10/06/2011 02:30 PM, ahmed lasheen wrote:
>> > thanks for your help
>> > when i use the bracket , i found that the image is produced for
>> > only quarter of the earth (the quarter has positive lat and lon only).
>> > i have attached the image produced.
>> >
>> >
>> > On Thu, Oct 6, 2011 at 10:15 PM, Adam Phillips <asphilli@ucar.edu
>> > <mailto:asphilli@ucar.edu>> wrote:
>> >
>> > Hi Ahmed,
>> > There are two ways to subscript: by using coordinate values, or by
>> using
>> > index values. I believe you are confusing the two methods here.
>> >
>> > Coordinate subscripting uses { } to denote that you are specifying
>> > coordinate values (latitude/longitude values in your case)
>> >
>> > elv = f->z({j1:j2:20},{i1:i2:20})
>> > elv = f->z({j1:j2},{i1:i2})
>> > elv = f->z({j1:},{:i2})
>> >
>> > Index subscripting doe not:
>> > elv = f->z(0:100:10,5:125:5})
>> > elv = f->z(::10,::5})
>> > elv = f->z(5::10,5::5})
>> >
>> > Note that in either case, the specified stride value must be an
>> > integer..
>> >
>> > In your case, I think you want this:
>> > elv = f->z({j1:j2:20},{i1:i2:20})
>> >
>> > Hope that helps..
>> > Adam
>> >
>> > On 10/06/2011 02:04 PM, ahmed lasheen wrote:
>> > > Hello all
>> > > I am trying to read and plot elevation data in which it's
>> > resolution is
>> > > 5 minute , and when i plot it keeping in mind maximizing the
>> > workstation
>> > > size and increasing the cnLevelSpacingF .
>> > > when i use the com elv = f->z(::20,::20) , then it works well
>> > but when
>> > > i use the command elv = f->z(j1:j2:20,i1:i2:20) where j1 is
>> the
>> > first
>> > > lat ,j2 is the second lat , i1 is the first lon and i2 is the
>> > second lon
>> > > , then the image produces is wrong .anybody know why , i
>> > > have attached the scripts.
>> > > thanks
>> > >
>> > > --
>> > > ===============
>> > > Ahmed Lasheen
>> > > Junior researcher at Cairo Numerical Weather Prediction Center
>> > (CNWPC)
>> > > Egyptian Meteorological Authority(EMA)
>> > > Cairo,Egypt
>> > > ===============
>> > >
>> > >
>> > >
>> > > _______________________________________________
>> > > ncl-talk mailing list
>> > > List instructions, subscriber options, unsubscribe:
>> > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>> >
>> > --
>> > ______________________________________________________________
>> > Adam Phillips asphilli@ucar.edu <mailto:asphilli@ucar.edu>
>> > NCAR/Climate and Global Dynamics Division (303) 497-1726
>> > P.O. Box 3000
>> > Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>> > _______________________________________________
>> > ncl-talk mailing list
>> > List instructions, subscriber options, unsubscribe:
>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>> >
>> >
>> >
>> >
>> > --
>> > ===============
>> > Ahmed Lasheen
>> > Junior researcher at Cairo Numerical Weather Prediction Center (CNWPC)
>> > Egyptian Meteorological Authority(EMA)
>> > Cairo,Egypt
>> > ===============
>> >
>> >
>> >
>> > _______________________________________________
>> > ncl-talk mailing list
>> > List instructions, subscriber options, unsubscribe:
>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>> --
>> ______________________________________________________________
>> Adam Phillips asphilli@ucar.edu
>> NCAR/Climate and Global Dynamics Division (303) 497-1726
>> P.O. Box 3000
>> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
>
> --
> ===============
> Ahmed Lasheen
> Junior researcher at Cairo Numerical Weather Prediction Center (CNWPC)
> Egyptian Meteorological Authority(EMA)
> Cairo,Egypt
> ===============
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

__________________________________________________
Adam Phillips asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division
P.O. Box 3000 (303) 497-1726
Boulder, CO 80307-3000
http://www.cgd.ucar.edu/cas/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Oct 7 09:16:22 2011