WRF subregion

On Fri, 18 Feb 2011 12:00:03 -0700, ncl-talk-request wrote
> Send ncl-talk mailing list submissions to
> ncl-talk@ucar.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> or, via email, send a message with subject or body 'help' to
> ncl-talk-request@ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner@ucar.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of ncl-talk digest..."
>
> Today's Topics:
>
> 1. how to plot a subregion with wrf_contour (WuLongtao)
> 2. Re: how to plot a subregion with wrf_contour (Dennis Shea)
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 18 Feb 2011 18:40:22 +0000
> From: WuLongtao <aos_rain@hotmail.com>
> Subject: how to plot a subregion with wrf_contour
> To: <ncl-talk@ucar.edu>
> Message-ID: <SNT134-w23E35090439D5590F20D5AF9D40@phx.gbl>
> Content-Type: text/plain; charset="gb2312"
>
> Hello,
>
> I would like to plot a subregion of my wrf output. I tried to set
> res@mpMinLatF = -50
> res@mpMaxLatF = 0
> res@mpMinLonF =-100
> res@mpMaxLonF =-70
>
> But it didn't work. The followings are my setting. Do you have any
> suggestion? Thanks!
>
> -Longtao
>
#######################################################################################

Longtao-

I have been trying to do the same thing & my colleague has come up with this method:

;zoom in on map
  mpres@mpLimitMode = 1
  minlat = 60
  maxlat = 62
  minlon = -152
  maxlon = -147

  mpres@mpMinLatF = minlat
  mpres@mpMaxLatF = maxlat
  mpres@mpMinLonF = minlon
  mpres@mpMaxLonF = maxlon

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

I have confirmed that it works, and that the final product after wrf_map_overlays()
is what I expect.

Comparing to yours, it looks like the magic bullet might be:

          mpres@mpLimitMode = 1

to get things off on the right foot.

PQO

> ********************************************************
> a = addfile("wrfout_d01_2006-09-22_12:00:00.nc","r")
> pblh = wrf_user_getvar(a,"PBLH",0)
> wks=gsn_open_wks("ps","pbl_wrf")
> res=True
> res@cnFillOn = True
> pltres = True
> mpres = True
> mpres@mpGeophysicalLineColor = "Black"
> mpres@mpNationalLineColor = "Black"
> mpres@mpUSStateLineColor = "Black"
> mpres@mpGridLineColor = "Black"
>
> gsn_define_colormap(wks,"rainbow")
> res@FieldTitle = "PBLH_WRF"
> contour=wrf_contour(a,wks,pblh,res)
> plot=wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
> ***************************************************************

    ___________________________________________________________
   |
   | Dr. Peter Q. Olsson,
   | Alaska State Climatologist, Alaska State Climate Center
   | Chief Scientist, Alaska Experimental Forecast Facility
   | University of Alaska Anchorage
   | 2811 Merrill Field Drive
   | Anchorage, AK 99501
   | voice: (907) 786-7234, fax (907) 786-7237|

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Feb 18 18:01:07 2011