Re: Errors: wrf_map_overlay: and fatal:ContourPlotPreDraw

Cindy, Mary,

Thanks. I'll try these now. I'm having an interesting set of problems
getting pnetcdf-enabled WRF to work, so I can work on post-processing
while I ponder next steps!

gerry

wrfhelp wrote:
> Gerry,
>
> In addition to the comments below - have you tried this scripts below?
> http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/BASIC_SFC/wrf_Surface1.htm
>
> It makes use of the newer functions, which remove the need to create a
> map for your plots.
>
> For very large datastes, I have found using RasterFill which Mary
> mentions very helpful.
>
> Cindy
>
> On Jan 29, 2010, at 4:22 PM, Mary Haley wrote:
>
>> Hi Gerry,
>>
>> The wrf_map_overlay function should still work, but the
>> wrf_map_overlays function will be the one that continues to get
>> updated in the future.
>>
>> One difference is that with wrf_map_overlay, you had to call wrf_map
>> to create the map. If you use wrf_map_overlays, it effectively calls
>> wrf_map for you, and does the overlay.
>>
>> I would try bumping up wsMaximumSize. You can do this in your
>> .hluresfile in your home directory if you have one. The line would
>> look like this:
>>
>> *wsMaximumSize : 100000000
>>
>> I see you also have it in your NCL script, which is also fine.
>>
>> I usually start by bumping this up by a factor of 10. If that fixes
>> the problem, then I bump it back down to half the size.
>>
>> The other thing you might try is turning on raster fill for your
>> wrf_contour plot:
>>
>> opts@cnFillMode = "RasterFill"
>>
>> This can be indistinguishable from the default contour fill, if you
>> have lots of points. You can also add this:
>>
>> opts@cnRasterSmoothingOn = True
>>
>> --Mary
>>
>>
>> On Jan 28, 2010, at 3:18 PM, Gerald Creager wrote:
>>
>>> Hi, Mary,
>>>
>>> I'm calling wrf_map_overlay, but it tells me that's obsolete. I'll
>>> see what I can find
>>> Here's the error... was working on this now, in fact.
>>>
>>> (0) /data/gerry/wind2/output/
>>> (0) Working on FILE:
>>> /data/gerry/wind2/output/wrfout_d03_2010-01-12_12:00:00
>>> (0) Working on time: 2010011212
>>> (0) wrf_map_overlay: Warning: This procedure is obsolete. Consider
>>> (0) using wrf_map_overlays instead.
>>> (0) Working on FILE:
>>> /data/gerry/wind2/output/wrfout_d03_2010-01-12_15:00:00
>>> (0) Working on time: 2010011215
>>> (0) wrf_map_overlay: Warning: This procedure is obsolete. Consider
>>> (0) using wrf_map_overlays instead.
>>>
>>> Let me see what I can find regarding the differences and I'll try to
>>> fix it on my own. I'm also attaching the code in question that's
>>> resulting in both errors, including:
>>>
>>> fatal:ContourPlotPreDraw: Workspace reallocation would exceed maximum
>>> size 32556688
>>> fatal:ContourPlotDraw: draw error
>>> warning:WorkstationDeactivate: workstation not active or not opened
>>>
>>> I can put several of the wrfout somewhere if needed for test.
>>>
>>> Thanks, gerry
>>>
>>>
>>> Mary Haley wrote:
>>>> On Jan 26, 2010, at 3:46 PM, Gerald Creager wrote:
>>>>> I'm looking at documentation now to see what I need to do, but it
>>>>> appears it's to simply change the name to wrf_map_overlays...
>>>>> except doing that yielded an error.
>>>>>
>>>> Hi Gerry,
>>>> [I see you have another question posted, so maybe answering this one
>>>> is a moot point...]
>>>> I didn't see in your code where you called wrf_map_overlays, or did
>>>> you rename it back to wrf_map_overlay?
>>>> The two functions are actually quite different. wrf_map_overlay
>>>> wants the map object to be passed, while wrf_map_overlays wants the
>>>> file pointer passed. wrf_map_overlays creates the map for you, based
>>>> on parameters set in the file.
>>>> What was the full error message?
>>>> Also, I've seen problems with just passing "True" to the wrf_xxxx
>>>> functions. To be safe, create a variable, like "mres = True", and
>>>> pass this to wrf_map_overlays.
>>>> --Mary
>>>>> On the second error, I tried up-sizing the wsMaximumSize to
>>>>> 40000000 but that didn't help much. Am I missing an update or
>>>>> similar command?
>>>>>
>>>>> Note that the WRF filesize is on the order of 1.5GB.
>>>>>
>>>>> The NCL script in question is attached.
>>>>>
>>>>> Thanks, Gerry
>>>>> --
>>>>> Gerry Creager -- gerry.creager@tamu.edu
>>>>> Texas Mesonet -- AATLT, Texas A&M University
>>>>> Cell: 979.229.5301 Office: 979.458.4020 FAX: 979.862.3983
>>>>> Office: 1700 Research Parkway Ste 160, TAMU, College Station, TX 77843
>>>>> ;### 20091001 gc
>>>>> ;### TAMU_SFC_sst.ncl
>>>>> ;### plot SSTs, wind, pressure...
>>>>>
>>>>>
>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
>>>>>
>>>>> begin
>>>>>
>>>>> ; Requires dataTime = 00z, 06z, 12z or 18z
>>>>> ; Requires initDate for plt name
>>>>> ; pltRoot = "Refract_"+ initDate + "T" + dataTime
>>>>> pltRoot = "Refract-test"
>>>>>
>>>>> it=0
>>>>> ; pltName = pltRoot + "."+ sprinti("%0.2i", it) ; create name
>>>>>
>>>>> ;
>>>>> ; The WRF ARW input file.
>>>>> ; This needs to have a ".nc" appended, so just do it.
>>>>> ; DATADir = "/jobs/wrfuser/wrf/runs/SCOOP12km/wrfprd/" + dataTime
>>>>> +"/" + initDate + "/"
>>>>> ;print (DATADir)
>>>>>
>>>>> ; FILES = systemfunc("ls " + DATADir + "wrfout_d01_*00")
>>>>> ; numFILES = dimsizes(FILES)
>>>>> ; a = addfiles(FILES+".nc","r")
>>>>>
>>>>> DATADir = "/data/gerry/wind2/output/"
>>>>> FILES = systemfunc("ls " + DATADir + "wrfout_d03_*00")
>>>>> numFILES = dimsizes(FILES)
>>>>> a = addfiles(FILES+".nc","r")
>>>>>
>>>>> ;
>>>>> ; Moved the next lines up here to try to get it to not overwrite
>>>>> the ps and pdf stuff.
>>>>> ; We generate plots, but what kind do we prefer?
>>>>> ; type = "x11"
>>>>> type = "pdf"
>>>>> ; type = "ps"
>>>>> ; type = "ncgm"
>>>>> ; type = "png"
>>>>> wks = gsn_open_wks(type,pltRoot)
>>>>>
>>>>> setvalues NhlGetWorkspaceObjectId()
>>>>> "wsMaximumSize" : 40000000
>>>>> end setvalues
>>>>> ; What times and how many time steps are in this data set?
>>>>> ; times = wrf_user_list_times(a) ; get times in the file
>>>>> Times = a[:]->Times
>>>>> times = wrf_times_c(Times,3)
>>>>> ntimes = dimsizes(times) ; number of times in the file
>>>>> ; Set some Basic Plot options
>>>>> ARWres = True
>>>>> ARWres@MainTitle = "WRF 12km Atlantic Basin: Surface level
>>>>> refractivity index"
>>>>>
>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>> do ifil = 0, numFILES-2,3 ; BIG FILES LOOP
>>>>> print("Working on FILE: " +FILES(ifil) )
>>>>>
>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>>
>>>>> ; do it = 0,ntimes-1,1 ; TIME LOOP
>>>>> it = 0
>>>>> ; Bodily insert the MMM WRF_SURFACE1 script here
>>>>>
>>>>> ; Example script to produce plots for a WRF real-data run,
>>>>> ; with the ARW coordinate dynamics option.
>>>>>
>>>>> print("Working on time: " + times(ifil) )
>>>>> ARWres@TimeLabel = times(ifil) ; Set Valid time to use on plots
>>>>>
>>>>> mpres = True ; Create map background
>>>>> mpres@mpFillOn = False
>>>>> mpres@mpFillColors = (/-1, -1, 8, -1/)
>>>>> mpres@mpFillDrawOrder = "postdraw"
>>>>> ;###
>>>>> gsnPaperOrientation = "landscape"
>>>>> lbOrientation = "vertical"
>>>>> mpres@gsnMaximize = True
>>>>> ;###
>>>>> map = wrf_map(wks,a[ifil],mpres)
>>>>>
>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>> ; First get the variables we will need
>>>>>
>>>>> slp = wrf_user_getvar(a[ifil],"slp",0) ; slp
>>>>> wrf_smooth_2d( slp, 3 ) ; smooth slp
>>>>> if ( it .eq. 0 ) then
>>>>> ; tc = wrf_user_getvar(a[ifil],"tc",0) ; 3D tc
>>>>> ; sst = wrf_user_getvar(a[ifil],"SST",0) ; SST
>>>>> ; sst=sst-273.16 ; SST2 in C
>>>>> ; td = wrf_user_getvar(a[ifil],"td",0) ; 3D td
>>>>> u = wrf_user_getvar(a[ifil],"ua",0) ; 3D U at mass points
>>>>> v = wrf_user_getvar(a[ifil],"va",0) ; 3D V at mass points
>>>>> ; tc2 = tc(0,:,:) ; Use lowest T at time zero
>>>>> ; sst2 = sst(0,:,:)
>>>>> ; td2 = td(0,:,:) ; Use lowest Td at time zero
>>>>> u10 = u(0,:,:) ; Use lowest level at time 0
>>>>> v10 = v(0,:,:)
>>>>> tc2 = wrf_user_getvar(a[ifil],"T2",0) ; T2 in Kelvin
>>>>> else
>>>>> ; sst = wrf_user_getvar(a[ifil],"SST",0) ; SST
>>>>> ; sst2 = wrf_user_getvar(a[ifil],"SST",0) ;sst2 in K
>>>>> ; tc2=a->T2
>>>>> tc2 = wrf_user_getvar(a[ifil],"T2",0) ; T2 in Kelvin
>>>>> ; sst=sst-273.16 ; SST2 in C
>>>>> ;/print(sst)
>>>>> ; tc2 = tc2-273.16 ; T2 in C
>>>>> ; td2 = wrf_user_getvar(a[ifil],"td2",0) ; Td2 in C
>>>>> u10 = wrf_user_getvar(a[ifil],"U10",0) ; u at 10 m, mass point
>>>>> v10 = wrf_user_getvar(a[ifil],"V10",0) ; v at 10 m, mass point
>>>>> end if
>>>>>
>>>>> ; tf2 = 1.8*tc2+32. ; Turn temperature into
>>>>> Fahrenheit
>>>>> ; tf2@description = "Surface Temperature"
>>>>> ; tf2@units = "F"
>>>>> ; sst@description = "Sea Surface Temperature"
>>>>> ; sst@units = "C"
>>>>>
>>>>> ; td_f = 1.8*td2+32. ; Turn temperature into
>>>>> Fahrenheit
>>>>> ; td_f@description = "Surface Dew Point Temp"
>>>>> ; td_f@units = "F"
>>>>> u10 = u10*1.94386 ; Turn wind into knots
>>>>> v10 = v10*1.94386
>>>>> u10@units = "kts"
>>>>> v10@units = "kts"
>>>>> tc2@units = "K"
>>>>>
>>>>> e = 10^((10.79574 * (1-273.16/tc2) - (5.02800 * log10(tc2/273.16))
>>>>> + (1.50475 * 10^(-4)* (1 - 10^(-8.2969*(tc2/273.16 - 1)))) + (
>>>>> 0.42873 * 10^(-3) ) * (10^(4.76955 * (1-273.16/tc2))-1)+.78614))
>>>>>
>>>>> RI = (77.6 * (slp/tc2))+(3.73*(10^5)^(e/tc2^2))
>>>>> RI@description = "Refractivity index"
>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>>
>>>>> ; Plotting options for T
>>>>> opts = ARWres
>>>>> opts@cnFillOn = True
>>>>> opts@ContourParameters = (/ 260., 300., 1./)
>>>>> ; opts@gsnSpreadColorEnd = -3 ; End third from the last color
>>>>> in color map
>>>>> opts@gsnSpreadColorEnd = -2 ; End third from the last color in
>>>>> color map
>>>>> contour_RI = wrf_contour(a[ifil],wks,RI,opts)
>>>>> delete(opts)
>>>>> delete(RI) ;because we're running outta memory!
>>>>>
>>>>> ; Plotting options for Td
>>>>> ; opts = ARWres
>>>>> ; opts@cnFillOn = True
>>>>> ; opts@cnLinesOn = True
>>>>> ; opts@cnLineLabelsOn = True
>>>>> ; opts@ContourParameters = (/ -20., 90., 5./)
>>>>> ; opts@cnLineLabelBackgroundColor = -1
>>>>> ; opts@gsnSpreadColorEnd = -3 ; End third from the last color
>>>>> in color map
>>>>> ; contour_td = wrf_contour(a[ifil],wks,td_f,opts)
>>>>> ; delete(opts)
>>>>> ; delete(td_f)
>>>>>
>>>>> ; Plotting options for SLP
>>>>> ; opts = ARWres
>>>>> ; opts@cnLineColor = "NavyBlue"
>>>>> ; opts@cnHighLabelsOn = True
>>>>> ; opts@cnLowLabelsOn = True
>>>>> ; opts@ContourParameters = (/ 960., 1040., 4. /)
>>>>> ; opts@cnLineLabelBackgroundColor = -1
>>>>> ; opts@gsnContourLineThicknessesScale = 2.0
>>>>> ; contour_psl = wrf_contour(a[ifil],wks,slp,opts)
>>>>> ; delete(opts)
>>>>> ; delete(slp)
>>>>>
>>>>> ; Plotting options for Wind Vectors
>>>>> opts = ARWres
>>>>> opts@FieldTitle = "Wind" ; overwrite Field Title
>>>>> opts@NumVectors = 47 ; density of wind barbs
>>>>> opts@vcGlyphStyle = "CurlyVector"
>>>>> vector = wrf_vector(a[ifil],wks,u10,v10,opts)
>>>>> delete(opts)
>>>>> delete(u10)
>>>>> delete(v10)
>>>>>
>>>>> ; Maximize plot
>>>>> ; psres = True
>>>>> ; maximize_output(wks,psres)
>>>>> ; MAKE PLOTS
>>>>> wrf_map_overlay(wks,map,(/contour_RI,vector/),True)
>>>>> ; wrf_map_overlay(wks,map,(/contour_td,vector/),True)
>>>>>
>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>> ; delete(wks)
>>>>> ; end do ; end of the time loop
>>>>>
>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>> end do ; END IF BIG FILE LOOP
>>>>>
>>>>> end
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>> --
>>> Gerry Creager -- gerry.creager@tamu.edu
>>> Texas Mesonet -- AATLT, Texas A&M University
>>> Cell: 979.229.5301 Office: 979.458.4020 FAX: 979.862.3983
>>> Office: 1700 Research Parkway Ste 160, TAMU, College Station, TX 77843
>>> ;### 20091001 gc
>>> ;### TAMU_SFC_sst.ncl
>>> ;### plot SSTs, wind, pressure...
>>>
>>>
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW_add.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
>>>
>>> begin
>>>
>>> ; Requires dataTime = 00z, 06z, 12z or 18z
>>> ; Requires initDate for plt name
>>> ; pltRoot = "Refract_"+ initDate + "T" + dataTime
>>> pltRoot = "Refract-test"
>>>
>>> it=0
>>> ; pltName = pltRoot + "."+ sprinti("%0.2i", it) ; create name
>>>
>>> ;
>>> ; The WRF ARW input file.
>>> ; This needs to have a ".nc" appended, so just do it.
>>> ; DATADir = "/jobs/wrfuser/wrf/runs/SCOOP12km/wrfprd/" + dataTime
>>> +"/" + initDate + "/"
>>> ;print (DATADir)
>>>
>>> ; FILES = systemfunc("ls " + DATADir + "wrfout_d01_*00")
>>> ; numFILES = dimsizes(FILES)
>>> ; a = addfiles(FILES+".nc","r")
>>>
>>> DATADir = "/data/gerry/wind2/output/"
>>> FILES = systemfunc("ls " + DATADir + "wrfout_d03_*00")
>>> numFILES = dimsizes(FILES)
>>> a = addfiles(FILES+".nc","r")
>>>
>>> ;
>>> ; Moved the next lines up here to try to get it to not overwrite the
>>> ps and pdf stuff.
>>> ; We generate plots, but what kind do we prefer?
>>> type = "x11"
>>> ; type = "pdf"
>>> ; type = "ps"
>>> ; type = "ncgm"
>>> ; type = "png"
>>> wks = gsn_open_wks(type,pltRoot)
>>>
>>> setvalues NhlGetWorkspaceObjectId()
>>> "wsMaximumSize" : 40000000
>>> end setvalues
>>> ; What times and how many time steps are in this data set?
>>> ; times = wrf_user_list_times(a) ; get times in the file
>>> Times = a[:]->Times
>>> times = wrf_times_c(Times,3)
>>> ntimes = dimsizes(times) ; number of times in the file
>>> ; Set some Basic Plot options
>>> ARWres = True
>>> ARWres@MainTitle = "WRF 12km Atlantic Basin: Surface level
>>> refractivity index"
>>>
>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>> do ifil = 0, numFILES-2,3 ; BIG FILES LOOP
>>> print("Working on FILE: " +FILES(ifil) )
>>>
>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>
>>> ; do it = 0,ntimes-1,1 ; TIME LOOP
>>> it = 0
>>> ; Bodily insert the MMM WRF_SURFACE1 script here
>>>
>>> ; Example script to produce plots for a WRF real-data run,
>>> ; with the ARW coordinate dynamics option.
>>>
>>> print("Working on time: " + times(ifil) )
>>> ARWres@TimeLabel = times(ifil) ; Set Valid time to use on plots
>>>
>>> mpres = True ; Create map background
>>> mpres@mpFillOn = False
>>> mpres@mpFillColors = (/-1, -1, 8, -1/)
>>> mpres@mpFillDrawOrder = "postdraw"
>>> ;###
>>> gsnPaperOrientation = "landscape"
>>> lbOrientation = "vertical"
>>> mpres@gsnMaximize = True
>>> ;###
>>> map = wrf_map(wks,a[ifil],mpres)
>>>
>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>> ; First get the variables we will need
>>>
>>> slp = wrf_user_getvar(a[ifil],"slp",0) ; slp
>>> wrf_smooth_2d( slp, 3 ) ; smooth slp
>>> if ( it .eq. 0 ) then
>>> ; tc = wrf_user_getvar(a[ifil],"tc",0) ; 3D tc
>>> ; sst = wrf_user_getvar(a[ifil],"SST",0) ; SST
>>> ; sst=sst-273.16 ; SST2 in C
>>> ; td = wrf_user_getvar(a[ifil],"td",0) ; 3D td
>>> u = wrf_user_getvar(a[ifil],"ua",0) ; 3D U at mass points
>>> v = wrf_user_getvar(a[ifil],"va",0) ; 3D V at mass points
>>> ; tc2 = tc(0,:,:) ; Use lowest T at time zero
>>> ; sst2 = sst(0,:,:)
>>> ; td2 = td(0,:,:) ; Use lowest Td at time zero
>>> u10 = u(0,:,:) ; Use lowest level at time 0
>>> v10 = v(0,:,:)
>>> tc2 = wrf_user_getvar(a[ifil],"T2",0) ; T2 in Kelvin
>>> else
>>> ; sst = wrf_user_getvar(a[ifil],"SST",0) ; SST
>>> ; sst2 = wrf_user_getvar(a[ifil],"SST",0) ;sst2 in K
>>> ; tc2=a->T2
>>> tc2 = wrf_user_getvar(a[ifil],"T2",0) ; T2 in Kelvin
>>> ; sst=sst-273.16 ; SST2 in C
>>> ;/print(sst)
>>> ; tc2 = tc2-273.16 ; T2 in C
>>> ; td2 = wrf_user_getvar(a[ifil],"td2",0) ; Td2 in C
>>> u10 = wrf_user_getvar(a[ifil],"U10",0) ; u at 10 m, mass point
>>> v10 = wrf_user_getvar(a[ifil],"V10",0) ; v at 10 m, mass point
>>> end if
>>>
>>> ; tf2 = 1.8*tc2+32. ; Turn temperature into
>>> Fahrenheit
>>> ; tf2@description = "Surface Temperature"
>>> ; tf2@units = "F"
>>> ; sst@description = "Sea Surface Temperature"
>>> ; sst@units = "C"
>>>
>>> ; td_f = 1.8*td2+32. ; Turn temperature into
>>> Fahrenheit
>>> ; td_f@description = "Surface Dew Point Temp"
>>> ; td_f@units = "F"
>>> u10 = u10*1.94386 ; Turn wind into knots
>>> v10 = v10*1.94386
>>> u10@units = "kts"
>>> v10@units = "kts"
>>> tc2@units = "K"
>>>
>>> e = 10^((10.79574 * (1-273.16/tc2) - (5.02800 * log10(tc2/273.16)) +
>>> (1.50475 * 10^(-4)* (1 - 10^(-8.2969*(tc2/273.16 - 1)))) + ( 0.42873
>>> * 10^(-3) ) * (10^(4.76955 * (1-273.16/tc2))-1)+.78614))
>>>
>>> RI = (77.6 * (slp/tc2))+(3.73*(10^5)^(e/tc2^2))
>>> RI@description = "Refractivity index"
>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>
>>> ; Plotting options for T
>>> opts = ARWres
>>> opts@cnFillOn = True
>>> opts@ContourParameters = (/ 260., 300., 1./)
>>> ; opts@gsnSpreadColorEnd = -3 ; End third from the last color in
>>> color map
>>> opts@gsnSpreadColorEnd = -2 ; End third from the last color in
>>> color map
>>> contour_RI = wrf_contour(a[ifil],wks,RI,opts)
>>> delete(opts)
>>> delete(RI) ;because we're running outta memory!
>>>
>>> ; Plotting options for Td
>>> ; opts = ARWres
>>> ; opts@cnFillOn = True
>>> ; opts@cnLinesOn = True
>>> ; opts@cnLineLabelsOn = True
>>> ; opts@ContourParameters = (/ -20., 90., 5./)
>>> ; opts@cnLineLabelBackgroundColor = -1
>>> ; opts@gsnSpreadColorEnd = -3 ; End third from the last color
>>> in color map
>>> ; contour_td = wrf_contour(a[ifil],wks,td_f,opts)
>>> ; delete(opts)
>>> ; delete(td_f)
>>>
>>> ; Plotting options for SLP
>>> ; opts = ARWres
>>> ; opts@cnLineColor = "NavyBlue"
>>> ; opts@cnHighLabelsOn = True
>>> ; opts@cnLowLabelsOn = True
>>> ; opts@ContourParameters = (/ 960., 1040., 4. /)
>>> ; opts@cnLineLabelBackgroundColor = -1
>>> ; opts@gsnContourLineThicknessesScale = 2.0
>>> ; contour_psl = wrf_contour(a[ifil],wks,slp,opts)
>>> ; delete(opts)
>>> ; delete(slp)
>>>
>>> ; Plotting options for Wind Vectors
>>> opts = ARWres
>>> opts@FieldTitle = "Wind" ; overwrite Field Title
>>> opts@NumVectors = 47 ; density of wind barbs
>>> opts@vcGlyphStyle = "CurlyVector"
>>> vector = wrf_vector(a[ifil],wks,u10,v10,opts)
>>> delete(opts)
>>> delete(u10)
>>> delete(v10)
>>>
>>> ; Maximize plot
>>> ; psres = True
>>> ; maximize_output(wks,psres)
>>> ; MAKE PLOTS
>>> wrf_map_overlay(wks,map,(/contour_RI,vector/),True)
>>> ; wrf_map_overlay(wks,map,(/contour_td,vector/),True)
>>>
>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>> ; delete(wks)
>>> ; end do ; end of the time loop
>>>
>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>> end do ; END IF BIG FILE LOOP
>>>
>>> end
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Feb 3 05:40:10 2010