<html>
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Colleagues<br>
<br>
Thanks for pointing out the gc_inout function, Rick. My regions are
currently simple shape files. <br>
<br>
The function gc_inout is working for me (please see test code and
test data at
<a class="moz-txt-link-rfc2396E" href="https://my.pcloud.com/publink/show?code=kZsanC7ZlKPCnfV5wnLFWm4C2yCww4mtUIYV"><https://my.pcloud.com/publink/show?code=kZsanC7ZlKPCnfV5wnLFWm4C2yCww4mtUIYV></a>).
But now I have what I think is a simple problem that I have not
solved. How to subset the data variables based on the logical
variable produced by gs_inout? I have tried to do so using a looped
conditional statement, but there is a problem with the indexed
assignment as I have written it. I'm also expecting that a loop like
this will be way to slow for the large data size I have. There must
be a smarter way to do this, but in reading the manuals and help
files, I'm still missing it. <br>
<br>
Any hints would be useful.<br>
<br>
Thanks, Sam<br>
<br>
<div class="moz-cite-prefix">On 28/12/18 2:44 AM, Rick Brownrigg
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:CAGKRhbGWcJ=AzGD=FO+YGp3rBDBw7wBEUbXVR0XGtP7gezB6ow@mail.gmail.com">
<meta http-equiv="content-type" content="text/html;
charset=windows-1252">
Hi Sam,
<div><br>
</div>
<div>You might take a look at the gc_inout function. Not sure how
you are defining your regions - shapefiles? Simple bounding
boxes?</div>
<div><br>
</div>
<div>Rick<br>
<br>
On Tuesday, December 25, 2018, Sam McClatchie <<a
href="mailto:smcclatchie@fishocean.info"
moz-do-not-send="true">smcclatchie@fishocean.info</a>>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF"> Colleagues<br>
<br>
I have some voluminous station data in ascii files that I
have written into a netcdf file with minimal pre-definition
of dimensions and attributes (please see the link below for
the data and code). My problem is that I want to subset the
global fishing_hours variable conveniently by region. <br>
<br>
The obvious solution would be to dimension the fishing_hours
variable properly into a 3-dimensional variable (time, lat,
lon) and use named coordinate subscripting. I realised when
I went back to rewrite the data with variable pre-definition
that I have a problem, because these data are not on any
kind of grid, so the sizes of the coordinate arrays are the
same size as the vector lengths for lat/ lon.<br>
<br>
My goal is to plot fishing_hours on maps where each point is
colour-coded by the magnitude of the fishing_hours variable
(as for example, here: <see <a
href="https://www.ncl.ucar.edu/Applications/Scripts/station_2.ncl"
target="_blank" moz-do-not-send="true">https://www.ncl.ucar.edu/<wbr>Applications/Scripts/station_<wbr>2.ncl</a>>).
Note that the data were not subset, I simply set the map
regions. Also, the plot only shows lat/lon variables, not
the fishing_hours variable because it has no coordinate
variables.<br>
<br>
The data, NCL code, and figure referred to in this message
can be downloaded from here:
<a
href="https://my.pcloud.com/publink/show?code=kZsanC7ZlKPCnfV5wnLFWm4C2yCww4mtUIYV"
target="_blank" moz-do-not-send="true"><https://my.pcloud.com/<wbr>publink/show?code=<wbr>kZsanC7ZlKPCnfV5wnLFWm4C2yCww4<wbr>mtUIYV></a><br>
I suggest downloading the directory as an archive and
unpacking it locally to a scratch directory.<br>
<br>
Any help on how to conveniently subset the global
fishing_hours data by lat/ lon would be appreciated.<br>
<br>
Best fishes<br>
Sam<br>
</div>
</blockquote>
</div>
</blockquote>
<br>
<div class="moz-signature">-- <br>
<meta content="text/html; charset=windows-1252"
http-equiv="content-type">
<title>email signature</title>
Sam McClatchie (fisheries oceanographer)<br>
& Elena Turin (accounting & auditing)<br>
FishOcean Enterprises<br>
38 Upland Rd, Huia, Auckland 0604, New Zealand<br>
cell: 027 752 8495<br>
<span style="text-decoration: underline;"><a
href="http://www.fishocean.info">Internet</a></span><br>
<img style="width: 150px; height: 149px;" alt=""
src="cid:part5.A070011D.DAFCFBA2@fishocean.info"><br>
<br>
"The time has come", the tui said,<br>
"to talk of many things:<br>
Of songs - and ferns - and flowering flax,<br>
of Pukekos and dreams ..."<br>
<br>
(not Lewis Carroll)
</div>
</body>
</html>