Re: zooming in on WRF plots

Joe,

Mary is correct. This needs to be integers.
Please let me know if this does not solve your problem.

wrfhelp

On Jul 14, 2008, at 1:57 PM, Mary Haley wrote:

> Joe,
>
> I looked at the WRFUserARW.ncl file, and it appears it is expecting
> index integers for these start and end values, and not actual lat/lon
> values. I think you need to determine at what indices your desired
> lat/lon values fall at in your 2D lat/lon arrays, and use these values
> instead.
>
> Just FYI: if you try to put a floating point value in an NCL variable
> that's already an integer, you will get the error message about "can't
> be coerced to type of...".
>
> --Mary
>
>
> On Mon, 14 Jul 2008, Joe Grim wrote:
>
>> Hi,
>>
>> I am new to NCL and attended the NCL workshop last week. I have
>> been trying to find a way to zoom in on WRF plots, and according
>> to the information on the WRF graphics site (http://
>> www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/
>> NCL_functions.htm#map_overlay), it's just a simple matter of
>> adding five lines to the code:
>> mpres_at_ZoomIn = True
>> mpres_at_Xstart = -74.75 (here I tried using longitude and latitude,
>> but I also tried using the 0.0 - 1.0 plotting coordinates as well)
>> mpres_at_Xend = -72.25
>> mpres_at_Ystart = 39.50
>> mpres_at_Yend = 41.75
>>
>> However, when I try to run it, I get the following errors:
>> fatal:Assignment type mismatch, right hand side can't be coerced
>> to type of left hand side
>> fatal:Execute: Error occurred at or near line 1518 in file ./
>> WRFUserARW.ncl
>>
>> fatal:Execute: Error occurred at or near line 1727 in file ./
>> WRFUserARW.ncl
>>
>> fatal:Execute: Error occurred at or near line 133 in file SST_new.ncl
>>
>> Do you know what I'm doing wrong? If I comment out these five new
>> lines, the plot is made beautifully, except that it's not zoomed
>> in, of course.
>>
>> Below is my entire code, in case that would be helpful.
>>
>> I would appreciate very much any help you could give!
>>
>> Joe Grim
>>
>>
>>
>> << SST_new.ncl >>
>> ; Example script to produce plots for a WRF real-data run,
>> ; with the ARW coordinate dynamics option.
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> ;load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>> load "./WRFUserARW.ncl"
>>
>> begin
>>
>> start_month="04"
>> start_day="10"
>> start_hour="06"
>> end_day="10"
>> end_hour="21"
>> ;
>> ; The WRF ARW input file.
>> ; This needs to have a ".nc" appended, so just do it.
>>
>> cntrl_file="/d1/grim/NASA_NYC/WRF_output/cntrl/2006"+start_month
>> +start_day+start_hour+"/wrfout_d03_2006-"+start_month+"-"+end_day
>> +"_"+end_hour+":00:00.GRM_P+FCST.nc"
>> new_sst_file="/d1/grim/NASA_NYC/WRF_output/new_sst/
>> 2006"+start_month+start_day+start_hour+"/
>> wrfout_d03_2006-"+start_month+"-"+end_day+"_"+end_hour+":
>> 00:00.GRM_P+FCST.nc" a = addfile(cntrl_file,"r")
>> b = addfile(new_sst_file,"r")
>>
>> ; We generate plots, but what kind do we prefer?
>> ; type = "x11"
>> ; type = "pdf"
>> type = "ps"
>> ; type = "ncgm"
>> wks = gsn_open_wks(type,"SST_diff/SST_new")
>> gsn_define_colormap(wks,"rainbow")
>>
>> ; Set some basic resources
>> res = True
>> res_at_MainTitle = "REAL-TIME WRF"
>>
>> pltres = True
>> mpres = True
>> mpres_at_mpUSStateLineColor = "black"
>> mpres_at_mpUSStateLineThicknessF = 2
>> mpres_at_mpGeophysicalLineColor = "black"
>> mpres_at_mpGeophysicalLineThicknessF = 2
>> mpres_at_mpOutlineBoundarySets = "GeophysicalAndUSStates"
>> mpres_at_ZoomIn = True
>> mpres_at_Xstart = -74.75
>> mpres_at_Xend = -72.25
>> mpres_at_Ystart = 39.50
>> mpres_at_Yend = 41.75
>> ; mpres_at_Xstart = 0.25
>> ; mpres_at_Xend = 0.75
>> ; mpres_at_Ystart = 0.25
>> ; mpres_at_Yend = 0.75
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>
>> it=0
>> times = wrf_user_list_times(a) ; get times in the file
>> res_at_TimeLabel = times(it) ; Set Valid time to use on plots
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> ; First get the variables we will need slp = wrf_user_getvar
>> (a,"slp",it) ; slp
>> wrf_smooth_2d( slp, 3 ) ; smooth slp
>> u_old = wrf_user_getvar(a,"ua",it) ; 3D U at mass points
>> v_old = wrf_user_getvar(a,"va",it) ; 3D V at mass points
>> u_new = wrf_user_getvar(b,"ua",it) ; 3D U at mass points
>> v_new = wrf_user_getvar(b,"va",it) ; 3D V at mass points
>> u10_old = u_old(0,:,:) ; Use lowest level at time 0
>> v10_old = v_old(0,:,:)
>> u10_new = u_new(0,:,:) ; Use lowest level at time 0
>> v10_new = v_new(0,:,:)
>>
>> ; Extra variables I want to look at
>> Land_mask = wrf_user_getvar(a,"XLAND",it) ; Land Mask (1 FOR
>> LAND, 2 FOR WATER
>> SST_old = wrf_user_getvar(a,"SST",it) ; SST
>> SST_new = wrf_user_getvar(b,"SST",it) ; SST
>> SST_diff = SST_new - SST_old ; SST difference btwn the
>> sims
>> SST_old = SST_old-273.15
>> SST_old_at_units = "C"
>> SST_old_at_description = "Old Sea Surface Temperatures"
>> SST_new = SST_new-273.15 SST_new_at_units
>> = "C"
>> SST_new_at_description = "New Sea Surface Temperatures"
>> SST_old = mask(SST_old,Land_mask,2)
>> SST_new = mask(SST_new,Land_mask,2)
>> SST_diff = mask(SST_diff,Land_mask,2)
>> SST_diff_at_description = "SST Difference (new-old)"
>> SST_diff_at_units = "C"
>> u10_diff = u10_new - u10_old
>> v10_diff = v10_new - v10_old
>> u10_old_at_units = "m/s"
>> v10_old_at_units = "m/s"
>> u10_new_at_units = "m/s"
>> v10_new_at_units = "m/s"
>> u10_diff_at_units = "m/s"
>> v10_diff_at_units = "m/s"
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>
>> ; Plotting options for New SST Wind Vectors
>> opts = res
>> opts_at_FieldTitle = "Wind" ; overwrite Field Title
>> opts_at_NumVectors = 24 ; density of wind barbs
>> opts_at_vcGlyphStyle = "LineArrow"; type of wind arrow/barb
>> opts_at_vcRefAnnoOn = True
>> opts_at_vcRefLengthF = 0.03
>> opts_at_vcRefMagnitudeF = 5.0
>> vector_new = wrf_vector(a,wks,u10_new,v10_new,opts)
>> delete(opts)
>>
>> ; Plotting options for SST_new
>> opts = res
>> opts_at_cnFillOn = True
>> opts_at_cnLinesOn = False
>> opts_at_ContourParameters = (/ 0., 12., 2./)
>> opts_at_gsnSpreadColors = True
>> opts_at_gsnSpreadColorEnd = -3 ; End third from the last color in
>> color map
>> contour_SST_new = wrf_contour(a,wks,SST_new,opts)
>> delete(opts)
>>
>>
>> ; MAKE PLOTS
>>
>> ; cmap = (/(/1.,1.,1./),(/0.,0.,0./),(/1.,.5,1./),(/1.,0.,1./),\
>> ; (/.5,0.,1./),(/0.,0.,1./),(/0.,.5,1./),(/0.,1.,1./),\
>> ; (/1.,1.,1./),(/1.,1.,1./),(/.7,1.,0./),\
>> ; (/1.,1.,0./),(/1.,.5,0./),(/1.,0.,0./),(/.5,0.,0./),\
>> ; (/.3,.3,0./),(/.5,.5,0./),(/.7,.7,0./)/)
>>
>> ; gsn_define_colormap(wks, cmap)
>> ; gsn_define_colormap(wks,"ViBlGrWhYeOrRe")
>>
>> plot = wrf_map_overlays(a,wks,(/
>> contour_SST_new,vector_new/),pltres,mpres)
>>
>> cmd="convert SST_diff/SST_new.ps SST_diff/SST_new.gif"
>> cmd2="rm -f SST_diff/SST_new.ps"
>> system (cmd)
>> system (cmd2)
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>
>> end
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk_at_ucar.edu
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jul 15 2008 - 09:59:37 MDT