Re: WRF subregion issue

Hi Chris,

The error message "check_for_ylat_coord and check_for_lon_coord" is
usually a result of the units attribute not being set correctly or at
all for the latitude and longitude arrays. Change this:

    lat2d = f->XLAT
    lon2d = f->XLONG
    x_at_lat2d = f->XLAT(0,:,:)
    x_at_lon2d = f->XLONG(0,:,:)

to this:

    lat2d = f->XLAT(0,:,:)
    lon2d = f->XLONG(0,:,:)
    lat2d_at_units = "degrees_north"
    lon2d_at_units = "degrees_east"
    x_at_lat2d = lat2d
    x_at_lon2d = lon2d

That should set the units up correctly, and get rid of the error
message. I believe you should also comment out the WRF_map_c line, as
that sets resources to draw the entire WRF-file defined domain (which
you don't want)... Let me know if those 2 changes don't fix things...
Adam

Chris Misenis wrote:
> Adam and NCL,
>
> I tried attaching the info as attributes in the following way, if I
> understand you correctly:
>
> So my code looks like this:
>
> ;************************************************
> ; Read perturbation geopotential at all times and levels
> ;************************************************
> ; FIRST x for CHEM; SECOND x for METEOR
> ; x = f->PH ; (Time, bottom_top, south_north, west_east)
> ; x = f->hcho ; (Time, bottom_top, south_north, west_east)
> x = f->T2 ; (Time, south_north, west_east)
> u = f->U10 ;
> v = f->V10 ;
> lat2d = f->XLAT
> lon2d = f->XLONG
> x_at_lat2d = f->XLAT(0,:,:)
> x_at_lon2d = f->XLONG(0,:,:)
>
> ; printVarSummary(times)
> print(x_at_description)
> ; printVarSummary(znu)
> printVarSummary(u)
>
> I turned tfDoNDCOverlay = False.
>
> But now, the error I'm getting is that no data is being plotted, but the
> map is still zoomed. Also, I get the errors during generation that say
> check_for_ylat_coord and check_for_lon_coord, where it says the data
> does not contain a validy lat/lon array. Have I misplaced the info in
> the beginning? Frustrating.
>
> Thanks,
>
> Chris
>
> On 6/14/07, *Adam Phillips* <asphilli_at_cgd.ucar.edu
> <mailto:asphilli_at_cgd.ucar.edu>> wrote:
>
> Hi Chris,
>
> By setting tfDoNDCOverlay = True, you are telling NCL to not transform
> the data to a specific projection. This results in slightly faster
> plotting, but I believe it doesn't allow you to zoom in on a subregion
> of your grid.
>
> (This should be better explained on the NCL webpage.)
>
> To get around this, set tfDoNDCOverlay = False, and attach the 2D
> lats/lons to the data array as attributes:
>
> arr_at_lat2d = x->XLAT(0,:,:)
> arr_at_lon2d = x->XLON(0,:,:)
>
> Zooming in should be no problem then.
>
> Good luck,
> Adam
>
> Chris Misenis wrote:
>
>
>
> --
> Chris Misenis
> 919.513.4438
> M.S. Student, Atmospheric Science
> Air Quality Forecasting Lab
> North Carolina State University
> http://www.meas.ncsu.edu/aqforecasting
> http://www4.ncsu.edu/~camiseni <http://www4.ncsu.edu/~camiseni>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
--------------------------------------------------------------
Adam Phillips			             asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
ESSL/CGD/CAS                               fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000	  http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk