Re: station data / gridded data

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon, 17 Jul 2006 12:02:59 -0600 (MDT)

>I have 2180 tree ring sites with lat(0:2179) and lon(0:2179)
>and ringwidth(0:2179,0:99) for (site,time).
>I also have a gridded dataset with sfct(0:99,0:359,0:719)
>for (time,lat,lon) and lato(0:359) for the latitudes and lono(0:719)
>for the longitudes.
>For each tree ring site, I want to find within which lat/lon box in the
>gridded dataset that it lies and use that sfct. What is the
>fastest way to do this, avoiding too many do loops?

Mike ...

Adam's solution is close.

Would the following work?

It might not find the "nearest" but rather the lower
left corner of the box. Given a 1x1 grid, I don't think
it makes much difference.

   nTim = 100
   nTree = 2178
   SFCT = new ( (/nTree,ntim/) , typeof(ringwidth), getFillValue(sfct))

   do n=0,nTree-1
      SFCT(n,:) = (/ sfct(:,{lat(n)},{lon(n)}) /) ; don't transfer meta data
   end do
                           ; optional meta data
   SFCT!0 = "tree"
   SFCT!1 = "time"
   SFCT&time = sfct&time
   SFCT_at_long_name = sfct_at_long_name +: nearby grid point"


ncl-talk mailing list
Received on Mon Jul 17 2006 - 12:02:59 MDT

This archive was generated by hypermail 2.2.0 : Mon Jul 17 2006 - 14:09:20 MDT