Dear all,I put:oro_d1b.ncl (script)WRF_oro_d1b.ncgm (ncgm)wrfout_d01_2005-12-12_06:00:00 (WRF output used in the script,it really will be available in 30 minutes because of the slow connession...) in the direcory you suggested.Anyway I reduced the shifted value (no more 50 km), but coastlinesare always not consistent with land/sea mask (see for instance the Sardinia island).Thanks,Andrea
|
|
To[IMAGE] |
[IMAGE]
andrea.storto@uniroma1.it |
|
cc[IMAGE] |
|
|
bcc[IMAGE] |
|
|
Subject[IMAGE] |
[IMAGE]
Re: [ncl-talk] shifted coastlines |
|
|
|
Buongiorno AndreaI think we have to see the full script and the data.Could you please ftp your script and data file toftp ftp.cgd.ucar.eduanonymousemailcd incomingput ... data file ...put ... ncl script ...quitAfter doing this please sent us an email telling usthe name of the script and file. [For securityreasons we can not look ar the 'incoming' directory.]or make them availabale via the WWW.THXDennis SheaOn Mon, 16 Jan 2006 andrea.storto@uniroma1.it wrote:>> Hi everyone,>> I run a two-nested WRF simulation,>> when I try to visualize horography>> (and consequently all the variables),>> I got the right domain and coastlines for the nested>> (finest resolution domain) while for>> the coarser one, coastlines are shifted>> of about 50 km westernward (and not>> in North-South direction). I'll try>> to change @mpLimitMode, the database version>> and the longitude of the center of projection,>> but it gets worst.>> The code I use (for both the domains) is:>> ....>> x = f->HGT(0,:,:)>> lat2d = f->XLAT(0,:,:)>> lon2d = f->XLONG(0,:,:)>> ....>> dimll = dimsizes(lat2d)>> nlat = dimll(0) ; to get dimension of height matrix>> mlon = dimll(1)>> ....>> res@mpProjection = "Sterographic" ; as set in WRF namelist>> res@mpDataBaseVersion = "HighRes">> res@mpLimitMode = "Corners">> res@mpLeftCornerLatF = lat2d(0,0)>> res@mpLeftCornerLonF = lon2d(0,0)>> res@mpRightCornerLatF = lat2d(0,0)>> res@mpRightCornerLonF = lon2d(0,0)>> res@mpCenterLonF = 12.50 ; as set in WRF namelist>> res@mpOutlineDrawOrder = "PostDraw">> res@mpOutlineBoundarySets = "GeophysicalAndUSStates">> res@tfDoNDCOverlay = True>> res@gsnAddCyclic = False>> ....>> plot = gsn_csm_contour_map(wks, x, res)>>>> Has anyone an idea why the d01 (coarser) domain>> has shifted coastlines and not the d02?>> Thanks in advance,>> Andrea>>>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jan 16 2006 - 09:40:30 MST