Hi Dennis
I had done already and it did not work (I get a white band on the right
side of the plot). A colleague figured out how to fix the problem. The
projection parameters are wrongly specified in the WRF file. For a
polar stereographic projection the center latitude should be 90., that is:
res_at_mpCenterLatF = 90.
and not what is specified with
res_at_mpCenterLatF = ff_at_TRUELAT1
Perhaps the next version of "WRFUserARW.ncl" will contain the fix.
Thanks,
Andrea
Dennis Shea wrote:
> I think the following will work:
>
> landmask_at_lat2d = xlat
> landmask_at_lon2d = xlon
>
> map = gsn_csm_contour_map(wks,landmask,res)
>
> Andrea Hahmann wrote:
>> Hi!
>>
>> I have been having problems trying to plot a land mask (or land use)
>> field on a stereographic projection. The data comes from a WRF "geo"
>> file. When I use res_at_tfDoNDCOverlay = True, then the coastline of the
>> data does not match the drawn geographic coastline. When I turn it off,
>> the coastlines line up but I get a strange white band on the right side
>> of the plot (also a small one on the bottom). I have experimented with
>> various parameters to set the projection (CEN_LAT, CEN_LON, and
>> STAND_LON), but nothing seems to help.
>> Any other ideas of what might be going wrong?
>>
>> The data is on an ftp site
>> (ftp://ftp.rap.ucar.edu/pub/hahmann/geo_em.d01.nc) if anybody is
>> interested in trying the script. The script and plot are attached.
>>
>> Thank you so much!
>> Andrea
>>
>>
>>
>> ------------------------------------------------------------------------
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk_at_ucar.edu
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-- ---------------------------------------------------------------- Andrea N. Hahmann, Ph.D. Research Applications Laboratory Natl. Center for Atmospheric Research Phone: 1-303-497-8383 PO BOX 3000 Fax: 1-303-497-8401 Boulder, CO 80301 hahmann_at_ucar.edu ---------------------------------------------------------------- _______________________________________________ ncl-talk mailing list ncl-talk_at_ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Mon Oct 15 2007 - 08:20:25 MDT
This archive was generated by hypermail 2.2.0 : Tue Oct 16 2007 - 09:18:51 MDT