[ncl-talk] gsn_csm_contour/ reversing x and y-axes/ overlay sampled points

Sam McClatchie smcclatchie at fishocean.info
Sun Sep 16 00:06:45 MDT 2018


Colleagues

I have two questions with regard to plotting a slice of data using 1-d 
array variables (temperature, depth and station):

1) I do not seem to be using the trYReverse and trXReverse resources 
correctly because my plot fails to show the filled contours when I 
uncomment these lines (see below). It plots fine before trying to 
reverse the axes.

2) I have failed to overlay the sample points using gsn_coordinates with 
a variable that has lat/ lon attributes. I used this method successfully 
when plotting surfaces rather than slices. My code does not modify the 
DrawOrder anywhere, so I assume the default PostDraw order.

Here is my code:
==== snip
  fname                      = 
"../data/CalCOFI_merged_line90_section_temperature.txt"
       lines                    = asciiread(fname,-1,"string")
       lines at _FillValue         = -999
       lines                    = 
str_sub_str(lines,",-999,",",lines at _FillValue,")
       Depthm               = tofloat(str_get_field(lines(1:),1," "))
       T_degC                = tofloat(str_get_field(lines(1:),2," "))
       T_degC at lat1d     = tofloat(str_get_field(lines(1:),3," "))  ; 
Attach as lat1d, lon1d
       T_degC at lon1d    = tofloat(str_get_field(lines(1:),4," "))  ; for 
plotting later
       Station                 = tofloat(str_get_field(lines(1:),6," "))

       ; plot options
       res at cnLevelSelectionMode = "ManualLevels"
       res at cnMinLevelValF            = 6.
       res at cnMaxLevelValF           = 14.
       res at cnLevelSpacingF          = 1.0 ; contour spacing
       res at cnFillOn                        = True
       res at tiMainString                 = "CalCOFI line 90 section"
       res at tiYAxisString                 = "Depth"       ; y-axis title
       res at tiXAxisString                = "Station"       ; x-axis title
       ; use a cmocean colormap
       res at cnFillPalette        = "cmocean_thermal"
       res at sfXArray             = Station
       res at sfYArray             = Depthm

       ; res at trYReverse         = True          ; reverses y-axis
       ; res at trXReverse         = True          ; reverses x-axis
       plot_TdegC_line90        = gsn_csm_contour(wks,T_degC,res)

       ; add stations
       gsn_coordinates(wks,plot_TdegC_line90,T_degC,pmres)

===== snip

And here is a snippet of my data

"Depthm" "T_degC" "Lat_Dec" "Lon_Dec" "Rpt_Line" "Rpt_Sta" "Cruise" 
"Year" "Month"
0 13.95 33.483333 -117.77 90 28 201203 2012 3
2 13.95 33.483333 -117.77 90 28 201203 2012 3
5 13.94 33.483333 -117.77 90 28 201203 2012 3
10 13.68 33.483333 -117.77 90 28 201203 2012 3
11 13.39 33.483333 -117.77 90 28 201203 2012 3
20 11.89 33.483333 -117.77 90 28 201203 2012 3
30 11.34 33.483333 -117.77 90 28 201203 2012 3
41 10.91 33.483333 -117.77 90 28 201203 2012 3
50 10.6 33.483333 -117.77 90 28 201203 2012 3
51 10.55 33.483333 -117.77 90 28 201203 2012 3
60 10.45 33.483333 -117.77 90 28 201203 2012 3
0 13.99 33.418333 -117.898333 90 30 201203 2012 3
2 13.99 33.418333 -117.898333 90 30 201203 2012 3
10 13.93 33.418333 -117.898333 90 30 201203 2012 3
11 13.89 33.418333 -117.898333 90 30 201203 2012 3
11 13.89 33.418333 -117.898333 90 30 201203 2012 3
20 13.74 33.418333 -117.898333 90 30 201203 2012 3
30 13.22 33.418333 -117.898333 90 30 201203 2012 3
31 12.79 33.418333 -117.898333 90 30 201203 2012 3
40 12.17 33.418333 -117.898333 90 30 201203 2012 3
50 11.25 33.418333 -117.898333 90 30 201203 2012 3
60 10.75 33.418333 -117.898333 90 30 201203 2012 3
70 10.59 33.418333 -117.898333 90 30 201203 2012 3
75 10.4 33.418333 -117.898333 90 30 201203 2012 3
85 10.26 33.418333 -117.898333 90 30 201203 2012 3
100 10.1 33.418333 -117.898333 90 30 201203 2012 3
119 9.93 33.418333 -117.898333 90 30 201203 2012 3
125 9.76 33.418333 -117.898333 90 30 201203 2012 3
141 9.56 33.418333 -117.898333 90 30 201203 2012 3
150 9.48 33.418333 -117.898333 90 30 201203 2012 3
170 9.32 33.418333 -117.898333 90 30 201203 2012 3
200 9.26 33.418333 -117.898333 90 30 201203 2012 3
201 9.22 33.418333 -117.898333 90 30 201203 2012 3
231 8.84 33.418333 -117.898333 90 30 201203 2012 3
250 8.72 33.418333 -117.898333 90 30 201203 2012 3
270 8.61 33.418333 -117.898333 90 30 201203 2012 3
300 8.57 33.418333 -117.898333 90 30 201203 2012 3
320 8.21 33.418333 -117.898333 90 30 201203 2012 3
381 7.74 33.418333 -117.898333 90 30 201203 2012 3
400 7.54 33.418333 -117.898333 90 30 201203 2012 3
441 6.98 33.418333 -117.898333 90 30 201203 2012 3
500 6.57 33.418333 -117.898333 90 30 201203 2012 3
516 6.46 33.418333 -117.898333 90 30 201203 2012 3
0 13.9 33.25 -118.255 90 35 201203 2012 3
2 13.9 33.25 -118.255 90 35 201203 2012 3
10 13.88 33.25 -118.255 90 35 201203 2012 3
10 13.88 33.25 -118.255 90 35 201203 2012 3
20 13.01 33.25 -118.255 90 35 201203 2012 3
30 11.66 33.25 -118.255 90 35 201203 2012 3
40 11.2 33.25 -118.255 90 35 201203 2012 3
50 10.79 33.25 -118.255 90 35 201203 2012 3
51 10.73 33.25 -118.255 90 35 201203 2012 3
60 10.58 33.25 -118.255 90 35 201203 2012 3
70 10.46 33.25 -118.255 90 35 201203 2012 3
75 10.38 33.25 -118.255 90 35 201203 2012 3
86 10.23 33.25 -118.255 90 35 201203 2012 3
100 10.06 33.25 -118.255 90 35 201203 2012 3
120 9.82 33.25 -118.255 90 35 201203 2012 3
125 9.82 33.25 -118.255 90 35 201203 2012 3
140 9.8 33.25 -118.255 90 35 201203 2012 3
150 9.64 33.25 -118.255 90 35 201203 2012 3
169 9.25 33.25 -118.255 90 35 201203 2012 3
199 8.96 33.25 -118.255 90 35 201203 2012 3
200 8.95 33.25 -118.255 90 35 201203 2012 3
231 8.6 33.25 -118.255 90 35 201203 2012 3
250 8.3 33.25 -118.255 90 35 201203 2012 3
270 8.16 33.25 -118.255 90 35 201203 2012 3
300 7.89 33.25 -118.255 90 35 201203 2012 3
320 7.79 33.25 -118.255 90 35 201203 2012 3
370 7.35 33.25 -118.255 90 35 201203 2012 3
0 13.92 33.183333 -118.391667 90 37 201203 2012 3
2 13.92 33.183333 -118.391667 90 37 201203 2012 3
10 13.93 33.183333 -118.391667 90 37 201203 2012 3
20 13.64 33.183333 -118.391667 90 37 201203 2012 3
21 13.56 33.183333 -118.391667 90 37 201203 2012 3
30 13.38 33.183333 -118.391667 90 37 201203 2012 3
41 12.89 33.183333 -118.391667 90 37 201203 2012 3
50 11.61 33.183333 -118.391667 90 37 201203 2012 3
51 11.47 33.183333 -118.391667 90 37 201203 2012 3
60 11.13 33.183333 -118.391667 90 37 201203 2012 3
71 10.71 33.183333 -118.391667 90 37 201203 2012 3
75 10.64 33.183333 -118.391667 90 37 201203 2012 3
85 10.46 33.183333 -118.391667 90 37 201203 2012 3
100 10.07 33.183333 -118.391667 90 37 201203 2012 3
122 9.76 33.183333 -118.391667 90 37 201203 2012 3
125 9.7 33.183333 -118.391667 90 37 201203 2012 3
141 9.6 33.183333 -118.391667 90 37 201203 2012 3
150 9.47 33.183333 -118.391667 90 37 201203 2012 3
170 9.12 33.183333 -118.391667 90 37 201203 2012 3
200 8.84 33.183333 -118.391667 90 37 201203 2012 3
201 8.83 33.183333 -118.391667 90 37 201203 2012 3
231 8.7 33.183333 -118.391667 90 37 201203 2012 3
250 8.5 33.183333 -118.391667 90 37 201203 2012 3
270 8.27 33.183333 -118.391667 90 37 201203 2012 3
300 7.96 33.183333 -118.391667 90 37 201203 2012 3
321 7.75 33.183333 -118.391667 90 37 201203 2012 3
381 7.19 33.183333 -118.391667 90 37 201203 2012 3
400 7.09 33.183333 -118.391667 90 37 201203 2012 3
439 6.76 33.183333 -118.391667 90 37 201203 2012 3
500 6.34 33.183333 -118.391667 90 37 201203 2012 3
515 6.28 33.183333 -118.391667 90 37 201203 2012 3
0 13.97 32.918333 -118.938333 90 45 201203 2012 3
2 13.97 32.918333 -118.938333 90 45 201203 2012 3
10 13.95 32.918333 -118.938333 90 45 201203 2012 3
10 13.95 32.918333 -118.938333 90 45 201203 2012 3

I have searched the examples and the maillist, but clearly have missed 
something.

Best fishes
Sam
-- 
email signature Sam McClatchie (fisheries oceanographer)
& Elena Turin (accounting & auditing)
FishOcean Enterprises
38 Upland Rd, Huia, Auckland 0604, New Zealand
cell: 027 752 8495
Internet <http://www.fishocean.info>


"The time has come", the tui said,
"to talk of many things:
Of songs - and ferns - and flowering flax,
of Pukekos and dreams ..."

(not Lewis Carroll)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180916/bd96f91f/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SamMcClatchie2_small.jpg
Type: image/jpeg
Size: 45362 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180916/bd96f91f/attachment.jpg>


More information about the ncl-talk mailing list