Re: Monotonically smaller displays with time

Hello

You are using wrf_* functions from the "WRFUserARW.ncl"
library. The version being used is some modification of this library.

*all* questions generated from using wrf_* should be
sent to wrfhelp_at_ucar.edu

They can best answer the wrf_* related questions.

Good luck

jwsmith_at_ucar.edu wrote:
> Hello,
>
> I am using the attached script for plotting surface ozone. There is a
> display for each time, but they are getting monotonically smaller with
> time. It only seems to be a problem with this surface variable script.
>
> Jonathan
>
> ----------------------------------------------------------------
> This message was sent using IMP, the Internet Messaging Program.
>
> ------------------------------------------------------------------------
>
> ; 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/wrf/WRFUserARW.ncl"
> load "/ptmp/jwsmith/NCLscripts/RON_BROWN/WRFUserARW_mary.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> begin
> ;
> ; The WRF ARW input file.
> ; This needs to have a ".nc" appended, so just do it.
> dirname = "/ptmp/jwsmith/wrfout_30min_printouts/" ;WRFchem20081111_fire_megan/run/" wrfhist_out_Jan07/"
> fname = "wrfout_d01_2006-06-07_00:00:00.nc"
> a = addfile(dirname + fname, "r")
>
>
> ; We generate plots, but what kind do we prefer?
> ; type = "x11"
> type = "pdf"
> ; type = "ps"
> ; type = "ncgm"
> wks = gsn_open_wks(type,"o3_surface_June7_30min_printout")
>
>
> ; Set some Basic Plot options
> res = True
> res_at_MainTitle = "WRF-Chem Surface "
>
> pltres = True
> mpres = True
> mpres_at_mpGeophysicalLineColor = "Black"
> mpres_at_mpNationalLineColor = "Black"
> mpres_at_mpUSStateLineColor = "Black"
> mpres_at_mpAllBoundariesLineColor = "Black"
> mpres_at_mpGridLineColor = "Black"
> mpres_at_mpLimbLineColor = "Black"
> mpres_at_mpPerimLineColor = "Black"
> mpres_at_mpOutlineBoundarySets = "National"
> mpres_at_mpOutlineOn = True
> mpres_at_mpOutlineBoundarySets = "AllBoundaries"
> mpres_at_mpGridLineThicknessF = 2.0
> mpres_at_mpLimbLineThicknessF = 2.0
> mpres_at_mpNationalLineThicknessF = 1.0
> mpres_at_mpUSStateLineThicknessF = 2.0
> mpres_at_mpAllBoundariesLineColor = 2.0
> map = wrf_map(wks,a,mpres)
>
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ; What times and how many time steps are in the data set?
> times = wrf_user_list_times(a) ; get times in the file
> ntimes = dimsizes(times) ; number of times in the file
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> gsn_define_colormap(wks,"gui_default") ;"wxpEnIR")
>
> do it = 0,ntimes-1 ; TIME LOOP
>
> print("Working on time: " + times(it) )
> 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
> o3 = wrf_user_getvar(a,"o3",it) ; 3D Carbon Monoxide
> o3 = o3*1000.
> o3_at_description = "Ozone Surface Concentration"
> o3_at_units = "ppbv"
> ; u = wrf_user_getvar(a,"ua",it) ; 3D U at mass points
> ; v = wrf_user_getvar(a,"va",it) ; 3D V at mass points
> ; u = u*1.94386 ; Turn wind into knots
> ; v = v*1.94386
> ; u_at_units = "kts"
> ; v_at_units = "kts"
>
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ; Plotting options for CO
> opts = res
> opts_at_cnFillOn = True
> opts_at_cnLinesOn = True
> ; opts_at_cnLevelSelectionMode = "AutomaticLevels"
> opts_at_cnMinLevelValF = min(o3)
> opts_at_cnMaxLevelValF = max(o3)
> opts_at_cnLevelSpacing = 10
> ;opts_at_ContourParameters = (/ -50., 90., 5./)
> ; opts_at_cnFillColors = (/"BlueDarkRed18"/)
> ; opts_at_gsnSpreadColorEnd = -3 ; End third from the last color in color map
> contour_o3 = wrf_contour(a,wks,o3(0,:,:),opts)
> delete(opts)
>
>
> ; Plotting options for SLP
> ; opts = res
> ; opts_at_cnLineColor = "Blue"
> ; opts_at_cnHighLabelsOn = True
> ; opts_at_cnLowLabelsOn = True
> ; opts_at_ContourParameters = (/ 900., 1100., 4. /)
> ; opts_at_cnLineLabelBackgroundColor = -1
> ; opts_at_gsnContourLineThicknessesScale = 2.0
> ; contour_psl = wrf_contour(a,wks,slp,opts)
> ; delete(opts)
>
> ; Plotting options for Wind Vectors
> ; opts = res
> ; opts_at_FieldTitle = "Wind" ; overwrite Field Title
> ; opts_at_NumVectors = 15 ; density of wind barbs
> ; vector = wrf_vector(a,wks,u(0,:,:),v(0,:,:),opts)
> ; delete(opts)
>
>
> ; MAKE PLOTS
> wrf_map_overlay(wks,map,(/contour_o3/),True)
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> end do ; END OF TIME LOOP
>
> end
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk