Re: MSLP: how to mask Tibet

Hi,

You might take a look at the shapefile_5.ncl example at:

http://ncl.ucar.edu/Applications/shapefiles.shtml

It reads a shapefile of a national boundary, and uses gc_inout() to
mask out areas outside the border.

Hope that helps...
Rick

On Tue, 24 Aug 2010 15:46:41 +0530
  Sahidul <sahidul@gmail.com> wrote:
> Hello NCL users,
>
> I am using following scripts to cal sea level pressure from WRF
>outputs. WRF
> outputs are saved in 1 our interval.
> I am able to plot it correctly.
> Now I want MASK TIBET area. Could you please help me how to do it.?
> (I am having Tibet boundary as latS = 28; latN = 40 ; lonW = 75
>; lonE
> = 104)
>
> ;;;;;;;;;;;;;; following is the script =============
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>
> begin
>
> datadir = "/WRF/output/2010062200/"
>FILES = systemfunc ("ls -1 " + datadir + "wrfout_d02*")
> numFILES = dimsizes(FILES)
>
> res = True
> res@InitTime = False
> res@Footer = False
> res@Footer = False
> pltres = True
> mpres = True
>
> sopts = True
> sopts@cnLineColor = "Magenta4"
> sopts@ContourParameters =
>(/950.,1050.,2./)
> sopts@cnInfoLabelOn = False
> sopts@cnLineLabelBackgroundColor = -1
> sopts@cnInfoLabelOrthogonalPosF = 0.07 ; offset
>second
> label information
> sopts@gsnContourLineThicknessesScale = 1.5
> mpres@mpGeophysicalLineColor = "Black"
> mpres@mpNationalLineColor = "black"
> mpres@mpOutlineOn = True
> mpres@mpNationalLineThicknessF = 1.8
> mpres@mpGeophysicalLineThicknessF = 1.8
>
> do ifil = 0,numFILES-1
> a = addfile(FILES(ifil)+".nc","r")
> times = wrf_user_list_times(a)
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> day = systemfunc("echo "+times +"|cut -c9-10")
> yr = systemfunc("echo "+times +"|cut -c1-4")
> mon = systemfunc("echo "+times +"|cut -c6-7")
> hr = systemfunc("echo "+times +"|cut -c12-13")
>
> if(mon .eq. "08")then
> mon="AUG"
> end if
>
> wks = gsn_open_wks("ps","slp_surf"+"_"+day+mon+yr+"_"+hr+"Z") ;;
>to save
> hourly plots in different names
> gsn_define_colormap(wks,"gui_default")
>
> ; PLOT for SLP
> slp = wrf_user_getvar(a,"slp",0)
> wrf_smooth_2d( slp,3 )
> contour_slp = wrf_contour(a,wks,slp,sopts)
> plot = wrf_map_overlays(a,wks,contour_slp,pltres,mpres)
> end do
> end
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> Please help me

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Aug 24 09:02:46 2010