wrf 2.1 subregion

Hi all,

I'm in a situation where I need to compared 12-km output to 4-km output and
I'd like to do so using the same subregion. To do so requires zooming in on
the 12-km data to the 4-km lat/lons. I've tried using the mpLimitmode =
"Corners" and such, but am only able to zoom the background map, not
transform the data. I've looked at other examples on NCL talk and the
website, but none seem to really answer the question. Below is a sample of
the code we're trying to use:

  wks = gsn_open_wks("ps" ,fnames(nf)+"_HFX") ;
ps,pdf,x11,ncgm,eps
  gsn_define_colormap(wks,"BlAqGrYeOrReVi200") ; select color map

  res = True ; plot mods desired
; res_at_gsnMaximize = True ; uncomment to maximize size
  res_at_gsnSpreadColors = True ; use full range of colormap
  res_at_gsnScalarContour = True ;the third field will be
drawn as a separate contour field
  res_at_cnFillOn = True ; color plot desired
  res_at_cnLinesOn = False ; turn off contour lines
  res_at_cnLineLabelsOn = False ; turn off contour labels
  res_at_vcRefMagnitudeF = 4.0 ; define vector ref mag
  res_at_vcRefLengthF = 0.04 ; define length of vec ref
  res_at_vcRefAnnoOrthogonalPosF = -0.09 ; move ref vector
  res_at_vcRefAnnoParallelPosF = 1.0
; res_at_vcRefAnnoSide = "Right" ; where to place the vector
reference annotation
  res_at_vcGlyphStyle = "CurlyVector" ; seems doesn't work, check
later
; res_at_tiXAxisString = "Longitude"
; res_at_tiYAxisString = "Latitude"
  res_at_lbOrientation = "Horizontal" ; Move labelbar
; res_at_cnLevelSelectionMode = "ManualLevels"
; res_at_cnMinLevelValF = 12.000
; res_at_cnMaxLevelValF = 48.000
; res_at_cnLevelSpacingF = 1.500
; vcres = True ; plot mods desired
  res_at_lbLabelStride = 4 ; plot every other colar bar label
; res_at_lbLabelStrings = sprintf("%4.2f",levels)
; res_at_lbLabelStrings = ("0.040",\
; ,"0.080","0.120"/)

;************************************************
; Use WRF_contributed procedure to set map resources
;************************************************
  WRF_map_c(f, res, 0) ; reads info from file
;************************************************
; set True for native projection (faster)
;************************************************
   res_at_tfDoNDCOverlay = True
;************************************************
; Turn on lat / lon labeling
;************************************************
  res_at_pmTickMarkDisplayMode = "Always" ; turn on tickmarks
  res_at_tmXTOn = False ; turn off top labels
  res_at_tmYROn = False ; turn off right labels
;************************************************
; Loop over all times and levels ( uncomment )
; Demo: one arbitrarily closen time and level
;************************************************
  dimx = dimsizes(x) ; dimensions of x
  ntim = dimx(0) ; number of time steps

  klev = dimx(1) ; number of "bottom_top"
levels

  nt = ntim/2 ; arbitrary time
  kl = 0 ; " level
  nt = 0
      res_at_tiMainFuncCode = "~"
       res_at_tiMainString = times
       res_at_gsnLeftString = x_at_description +" "
       plot = gsn_csm_contour_map(wks,x(nt,:,:),res)
  dummy = systemfunc ("convert "+fnames(nf)+"_HFX.ps
"+fnames(nf)+"_HFX.gif")
  end do
end

As I mentioned, I can get the zoomed in background map, but no further. I
did so by both going into WRF_contributed and editing, as well as putting it
in the code. Any ideas?

Thanks,

Chris

-- 
Chris Misenis
919.513.4438
M.S. Student, Atmospheric Science
Air Quality Forecasting Lab
North Carolina State University
http://www.meas.ncsu.edu/aqforecasting
http://www4.ncsu.edu/~camiseni

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jun 13 2007 - 09:58:19 MDT