zooming in on WRF plots

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
Received on Mon Jul 14 2008 - 11:14:13 MDT