[ncl-talk] X axis tick mark with Explicit Mode
Adam Phillips
asphilli at ucar.edu
Fri Sep 13 17:07:19 MDT 2019
Hi Prasad,
Please do not write to individuals who respond to ncl-talk and not to
ncl-talk itself, as that way others can assist and see the resolution of
the issue. My original email still stands: The range of values specified in
tmXBValues needs to match the range of x-axis coordinate values of the
plotted array. In your case, from your printVarSummary output, you have
latitudes that are around 30N:
lat: [29.84782..33.84531]
so you will want to set tmXBValues to have values between 29.84 and 33.84.
For example:
resx at tmXBMode = "Explicit"
resx at tmXBValues = fspan(29.5,34,10)
resx at tmXBLabels = resx at tmXBValues
Adam
On Thu, Sep 12, 2019 at 1:12 AM Prasad Kunjeer <pskunjeer at gmail.com> wrote:
> Sir,
> As suggested by you, I have drawn in default mode without specifying the
> explicit x labels. The output is enclosed as an attachment. The figure
> shows only one lat value of 30N. The output of the printVarSummary(slice1)
> is
>
> Variable: slice1
> Type: float
> Total Size: 19328 bytes
> 4832 values
> Number of Dimensions: 2
> Dimensions and sizes: [lev | 32] x [lat | 151]
> Coordinates:
> lev: [996.9073..5.045747]
> lat: [29.84782..33.84531]
> Number Of Attributes: 8
> lon : 76.08612
> Time : 2018-09-20_15:00:00
> FieldType : 104
> MemoryOrder : XYZ
> description : x-wind component
> units : m s-1
> stagger :
> coordinates : XLONG XLAT
>
> The Variable does not have longitude value, is this giving the problem in
> an explicit mode where I am specifying the longitude value.
>
> I have tried explicit mode without longitudinal value but this also gives
> no x value labels.
>
> Regards
> प्रसाद कुंजीर/ Prasad Kunjeer
> वैज्ञानिक - सी / Scientist - C
> केन्द्रीय जल और विद्युत अनुसंधान शाला / Central Water and Power Research
> Station
> पुणे 411 024/ Pune 411 024
>
>
> On Wed, Sep 11, 2019 at 9:44 PM Adam Phillips <asphilli at ucar.edu> wrote:
>
>> Hi Prasad,
>> Whenever you are trying to set explicit tickmarks and they disappear,
>> then it is likely that the domain of your set tm??Values is outside the
>> corresponding range of the data coordinate variable (in this case the
>> latitudes of your slice1 array). Try commenting out the resources that
>> start with res at tmXB*, and labels should appear (although they likely
>> will not be what you want). Once you see the range of the default tick mark
>> labels, that indicates what the range should be for your res at tmXBValues
>> resource.
>>
>> If the above does not help please respond to ncl-talk with the output
>> from printVarSummary(slice1)
>> Adam
>>
>>
>>
>>
>>
>> On Wed, Sep 11, 2019 at 12:30 AM Prasad Kunjeer via ncl-talk <
>> ncl-talk at ucar.edu> wrote:
>>
>>> Sir,
>>> I am using the following script to plot U wind from WRF output file. I
>>> am able to plot the figure, but the X-axis tick marks are not appearing.
>>> The script is follow
>>>
>>> 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/WRF_contributed.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
>>>
>>> begin
>>>
>>> fName = "/home/emsuser/Build_WRF/wrfprd/results/do3/t/do3.nc"
>>> f = addfile(fName, "r")
>>>
>>> U1 = f->U
>>> U = wrf_user_unstagger(U1,U1 at stagger)
>>> lat = f->XLAT(0,:,0)
>>> lon = f->XLONG(0,0,:)
>>> znu = f->ZNU(0,:) ; eta levels
>>> lev = znu*1000 ; [996.9073..5.045747]
>>> Time = chartostring( f->Times(:,:) )
>>>
>>> U!1 = "lev"
>>> lev at units= "hPa"
>>> U!0 = "Time"
>>> U!2 = "lat"
>>> lat at units = "degree_north"
>>> U!3 = "lon"
>>> lon at units = "degree_east"
>>>
>>> U&lev = lev
>>> U&Time = Time
>>> U&lat = lat
>>> U&lon = lon
>>>
>>> wgt = (/ 1., 3., 4., 3., 1./) ; wgts for temporal
>>> smooth
>>> wgt = wgt/sum(wgt) ; normalize
>>>
>>> lat2d = f->XLAT(0,:,:)
>>> lon2d = f->XLONG(0,:,:)
>>>
>>> slice = U(lev|:,Time|:,lat|:,lon|:) ; trick to copy
>>> meta data
>>> slice1 = U(lev|:,Time|1,lat|84:234,lon|175)
>>> printVarSummary(slice1)
>>>
>>>
>>> wks = gsn_open_wks("png" ,"WRF_U_Wind_New_1") ; open
>>> workstation
>>>
>>> ;+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
>>> ; Create resources for cross section(s)
>>> ;+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
>>> resx = True
>>> resx at gsnMaximize = True
>>> resx at cnFillOn = True ; turn on color fill
>>> resx at cnLinesOn = False ; turn lines on/off
>>> ; True is default
>>> resx at cnLineLabelsOn = False ; turn line labels
>>> on/off ; True is default
>>> resx at cnFillPalette = "ncl_default"
>>>
>>> nLabels = 8 ; arbitrary
>>> resx at tmXBMode = "Explicit"
>>> resx at tmXBLabels = new(nLabels,"string")
>>> resx at tmXBValues := toint(fspan(84,234,nLabels))
>>> do i=0,nLabels-1
>>> x = lon2d(resx at tmXBValues(i),176)
>>> y = lat2d(resx at tmXBValues(i),176)
>>> resx at tmXBLabels(i) = sprintf("%5.1f", y)+"~C~"+sprintf("%5.1f",
>>> x)
>>> print(resx at tmXBLabels(i))
>>> end do
>>>
>>> plot = gsn_csm_pres_hgt(wks, slice1, resx)
>>>
>>> end
>>>
>>> I have attached the plot for ready reference. Any help is highly
>>> appreciated.
>>>
>>> Regards
>>> प्रसाद कुंजीर/ Prasad Kunjeer
>>> वैज्ञानिक - सी / Scientist - C
>>> केन्द्रीय जल और विद्युत अनुसंधान शाला / Central Water and Power
>>> Research Station
>>> पुणे 411 024/ Pune 411 024
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://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>
>>
>
--
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: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190913/2f6d0b56/attachment.html>
More information about the ncl-talk
mailing list