[ncl-talk] Why do I get "warning:gsnLeftString is not a valid resource in Plots00-P850_28_00Z_contour at this time"
Zilore Mumba
zmumba at gmail.com
Tue Apr 5 12:20:51 MDT 2022
I will try that. Thanks.
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
Virus-free.
www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
On Tue, Apr 5, 2022 at 8:07 PM Adam Phillips <asphilli at ucar.edu> wrote:
> Hi Zilore,
> I believe many if not all of the wrf plotting functions do not use the
> gsn_csm_* suite of plotting functions, and thus one cannot set subtitles
> with the gsn*String resources. (Looking at the code for wrf_contour
> (available in $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl), you can
> see that gsn_contour is being called.) According to the wrf_contour
> documentation
> <https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_contour.shtml>,
> you can set NoHeaderFooter = True to turn off the subtitles, did you try
> that? If that does not work, I would suggest you look at the
> WRFUserARW.ncl code (specifically look at the _SetMainTitle function) to
> see if you can turn off the labels manually.
>
> Hope that helps!
> Adam
>
> On Tue, Apr 5, 2022 at 3:55 AM Zilore Mumba via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> In my wrf plot script below, I am trying to suppress the left, center and
>> right string as well as the label at the bottom right. None of the wrf
>> script examples I have seen has these fields suppressed. My gsn and cn
>> resources (marked in blue in the code below) do not seem to have any
>> effect, hence the warning message.
>> When I print these resources (in purple below) they are captured ok.
>> Is it that with this kind of plot one cannot suppress these fields?
>>
>>
>> ; Example script to produce plots for a WRF real-data run,
>> ; with the ARW coordinate dynamics option.
>> ; Interpolating to specified pressure levels
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>>
>> begin
>> ;
>> ; The WRF ARW input file.
>> DATADir = "/home/zmumba/DA/OUTPUT/2022022800/noda/"
>> FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d01* ")
>> numFILES = dimsizes(FILES)
>>
>> files_all = addfiles(FILES+".nc","r")
>>
>> ; We generate plots, but what kind do we prefer?
>> ; type = "x11"
>> ; type = "pdf"
>> ; type = "ps"
>> ; type = "ncgm"
>> type = "png"
>> ; wks = gsn_open_wks(type,"plt_PressureLevel1")
>>
>> ; Set some Basic Plot options
>> res = True
>> res at MainTitle = "WRF MODEL OUTPUT"
>> res at Footer = False
>>
>> res at cnInfoLabelOn = False ; Suppress text "CONTOUR FROM # TO # by #"
>> res at gsnCenterString = "" ;"gsnCenterString ~C~ ~Z80~ (Default: None)"
>> res at gsnLeftString = "" ;"gsnRightString ~C~ ~Z80~ (Default:
>> Long_Name)"
>> res at gsnRightString = "" ;"gsnRightString ~C~ ~Z80~ (Default: Units)"
>>
>> pltres = True
>> mpres = True
>> mpres at mpFillOn = False ; turn off gray fill
>> mpres at mpOutlineBoundarySets = "National" ; turn on country
>> boundaries
>> mpres at mpGeophysicalLineColor = "Navy" ; color of cont.
>> outlines
>> mpres at mpGeophysicalLineThicknessF = 1.5 ; thickness of
>> outlines
>> mpres at mpNationalLineColor = "Black"
>> mpres at mpNationalLineThicknessF = 2.5 ; turn on country
>> boundaries
>>
>> mpres at mpDataBaseVersion = "MediumRes" ; choose higher
>> resolution
>> mpres at mpDataSetName = "Earth..4" ; choose most recent
>> boundaries
>>
>> mpres at mpMaxLatF = -4 ; choose subregion
>> mpres at mpMinLatF = -21
>> mpres at mpMaxLonF = 37
>> mpres at mpMinLonF = 21
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>
>> ; What times and how many time steps are in the data set?
>> times = files_all[:]->Times
>> ntimes = dimsizes(times) ; number of times in the file
>>
>> ; The specific pressure levels that we want the data interpolated to.
>> pressure_levels = (/ 850., 700., 500., 300./) ; pressure levels to
>> plot
>> nlevels = dimsizes(pressure_levels) ; number of pressure
>> levels
>>
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>
>> ; VTlab = (/ "T+00", T+00", "T+06", "T+12", "T+18", "T+24", "T+30",
>> "T+36", "T+42", "T+48"/) ; Validity for file name
>>
>> ; get time information and strip out the day and hour
>> times_in_file = files_all[:]->Times
>> dims = dimsizes (times)
>> ntimes = dimsizes (times(:,0))
>>
>> times_to_use = new(dims(0),string)
>>
>> do i=0,dims(0)-1
>> times_to_use(i) = chartostring(times_in_file(i,8:12))
>> end do
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> do it = 0,ntimes(0)-1,2 ; TIME LOOP
>>
>> print("Working on time: " + times_to_use(it) )
>> res at TimeLabel = times_to_use(it) ; Set Valid time to use on plots
>>
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> ; First get the variables we will need
>>
>> tc = wrf_user_getvar(files_all,"tc",it) ; T in C
>> u = wrf_user_getvar(files_all,"ua",it) ; u averaged to mass
>> points
>> v = wrf_user_getvar(files_all,"va",it) ; v averaged to mass
>> points
>> p = wrf_user_getvar(files_all, "pressure",it) ; pressure is our
>> vertical coordinate
>> z = wrf_user_getvar(files_all, "z",it) ; grid point height
>> rh = wrf_user_getvar(files_all,"rh",it) ; relative humidity
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> do level = 0,nlevels-1 ; LOOP OVER LEVELS
>> pressure = pressure_levels(level)
>>
>> tc_plane = wrf_user_intrp3d(tc,p,"h",pressure,0.,False)
>> z_plane = wrf_user_intrp3d( z,p,"h",pressure,0.,False)
>> rh_plane = wrf_user_intrp3d(rh,p,"h",pressure,0.,False)
>> u_plane = wrf_user_intrp3d( u,p,"h",pressure,0.,False)
>> v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
>>
>> spd = (u_plane*u_plane + v_plane*v_plane)^(0.5) ; m/sec
>> spd at description = "Wind Speed"
>> spd at units = "m/s"
>> u_plane = u_plane*1.94386 ; kts
>> v_plane = v_plane*1.94386 ; kts
>> u_plane at units = "kts"
>> v_plane at units = "kts"
>>
>> wks = gsn_open_wks(type,"Plots00-" + "P" + pressure_levels(level) +
>> "_" + times_to_use(it) + "Z")
>>
>> ; Plotting options for T
>> opts = res
>> print(opts))
>>
>> opts at cnLineColor = "Red"
>> opts at ContourParameters = (/ 5.0 /)
>> opts at cnInfoLabelOrthogonalPosF = 0.07 ; offset second label
>> information
>> opts at gsnContourLineThicknessesScale = 2.0
>> contour_tc = wrf_contour(files_all[it],wks,tc_plane,opts)
>> delete(opts)
>>
>> ; MAKE PLOT
>> if ( pressure .eq. 850 ) then ; plot temp, rh, height, wind
>> barbs
>> opts_z at ContourParameters = (/ 20.0 /)
>> contour_height = wrf_contour(files_all[it],wks,z_plane,opts_z)
>> plot =
>> wrf_map_overlays(files_all[it],wks,(/contour_rh,contour_tc,contour_height, \
>> vector/),pltres,mpres)
>> end if
>>
>> end do ; END OF LEVEL LOOP
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> end do ; END OF TIME LOOP
>> end
>>
>>
>>
>> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Virus-free.
>> www.avast.com
>> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
>> <#m_-2063917176174049324_m_-1170957471822237576_m_6353160373876870093_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> --
> Adam Phillips
> Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220405/686aedb3/attachment.html>
More information about the ncl-talk
mailing list