Re: areahi2lores and plotting

Adam, Thank you.
I didn't realize that my use of lonFlip was causing problems. I erased
all of the lonFlip commands that were causing discrepancies, but I
still got the same messy result.But when I tried your suggestion at
the end ...I got the plot that I needed.
Most helpful
Thanks again.
-Erik

On Mon, May 10, 2010 at 4:53 PM, Adam Phillips <asphilli@ucar.edu> wrote:
> Hi Erik,
> Whenever you see that type of plotting effect after regridding you
> should immediately take a good look at the input/output
> latitudes/longitudes and data arrays in your regridding call to make
> sure that they span the same longitude range and are in the same direction.
>
> In your case:
> WRF_prate_day = model->rain
> WRF_lat = model->south_north
> WRF_lon = model->west_east
> WRF_prate_day = lonFlip(WRF_prate_day)
> WRF_lat=WRF_lat(::-1)
>
> (snip)
>
> latGRID        = fspan(-90,  90, 73)
> latGRID!0      = "latGRID"
> latGRID@units  = "degrees_north"
> latGRID&latGRID = latGRID
> NLAT = dimsizes(latGRID)
>
> lonGRID        = fspan( 0, 357.5, 144)
> lonGRID!0      = "lonGRID"
> lonGRID@units  = "degrees_east"
> lonGRID&lonGRID = lonGRID
> MLON = dimsizes(lonGRID)
>
> (snip)
>
> newWRF_prate_day = area_hi2lores_Wrap (WRF_lon,WRF_lat, WRF_prate_day,
> False, 1, lonGRID, latGRID, opt)
>
> The WRF_lon array, which was never lonFlip'd, is input as the 1st
> argument. WRF_prate_day, input as the 3rd argument, was lonFlip'd. I
> think this discrepancy is the cause of your plot cutting off at the
> prime meridian.
>
> Also, you flip the latitudes of WRF_lat, but you don't flip the
> latitudes of the WRF_prate_day array. (Unless I am missing something.)
> You should either flip both of those or flip neither of them, as WRF_lat
> is input as the 2nd argument in your call to area_hi2lores_Wrap.
>
> My advice: Do not use lonFlip or flip the latitudes, alter latGRID to go
> from 90:-90, and lonGRID to go from -180 to 180, and regrid. After
> regridding, lonFlip and flip the latitudes to match your NCPR2 data and
> to perform the difference.
>
> Good luck,
> Adam
>
>
>
> On 05/10/2010 09:41 AM, Erik Noble wrote:
>>Dear NCL,
> I am using the areahi2lores function to plot hi-res & lo-res
> precipitation data onto the same grid. Then I create a panel plot of
> each data set (2 time vs. longitude plots along with a difference
> plot).
>
> My plots will not "plot" past the prime meridian (see picture).
> Both of my datasets "should" be plotting past this point since they
> contain data at all of these points.
> I am using the lonFlip command in my script so I don't think this is
> my problem. Could there be something else causing this?
> (I'm not trying to ask others to debug my program; helpful suggestions
> would be nice.)
> Thank you in advance.
> Sincerely,
> Erik Noble
>
> __________________________________________________
> Adam Phillips
> asphilli@ucar.edu
> National Center for Atmospheric Research   tel: (303) 497-1726
> Climate and Global Dynamics Division         fax: (303) 497-1333
> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon May 10 20:43:53 2010