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">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">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"><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>-- <br>
      
      
      Sam McClatchie (fisheries oceanographer)<br>
      & Elena Turin (accounting & auditing)<br>
      FishOcean Enterprises<br>
      <a href="https://maps.google.com/?q=38+Upland+Rd,+Huia,+Auckland+0604,+New+Zealand&entry=gmail&source=g">38 Upland Rd, Huia, Auckland 0604, New Zealand</a><br>
      cell: 027 752 8495<br>
      <span style="text-decoration:underline"><a href="http://www.fishocean.info" target="_blank">Internet</a></span><br>
      <img style="width:150px;height:149px" alt="" src="cid:part2.B64C7967.C7905AA9@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>
  </div>

</blockquote></div>