[Thread Prev][Thread Next][Index]

RE: [ferret_users] regridding



Hi,
so I've tried to use rect_to_curv function put a rectilinear data set on a curvilinear grid but with no success:

let a = rect_to_curv(hi[l=1],tlon,tlat,2)
yes? shade a, tlon, tlat

Bailing out of external function "rect_to_curv":
    The range of source grid longitude may be no larger than tpi
 **ERROR: error in external function

and so details of the dataset

yes? sh dat
     currently SET data sets:
    1> ./R1009RBRcevp01a.cice.h.hi.nc  (default)
 name     title                             I         J         K         L         M         N
 TLON     T grid center longitude          1:1280    1:720     ...       ...       ...       ...
 TLAT     T grid center latitude           1:1280    1:720     ...       ...       ...       ...
 ULON     U grid center longitude          1:1280    1:720     ...       ...       ...       ...
 ULAT     U grid center latitude           1:1280    1:720     ...       ...       ...       ...
 TMASK    ocean grid mask                  1:1280    1:720     ...       ...       ...       ...
 TAREA    area of T grid cells             1:1280    1:720     ...       ...       ...       ...
 ANGLE    angle grid makes with latitude   1:1280    1:720     ...       ...       ...       ...
 HI       grid cell mean ice thickness     1:1280    1:720     ...       1:120     ...       ...

yes? sh grid hi
    GRID GEE2
 name       axis              # pts   start                end
 NI        X                 1280 r   1                    1280
 NJ        Y                  720 r   1                    720
 normal    Z
 TIME      TIME               120 i   01-FEB-1980 00:00    01-JAN-1990 00:00
 normal    E
 normal    F
yes? sh grid tlon
    GRID GEE1
 name       axis              # pts   start                end
 NI        X                 1280 r   1                    1280
 NJ        Y                  720 r   1                    720
 normal    Z
 normal    T
 normal    E
 normal    F

thanks,
Robert



From:
owner-ferret_users@xxxxxxxx [owner-ferret_users@xxxxxxxx] on behalf of Ansley C. Manke [ansley.b.manke@xxxxxxxx]

Sent: Monday, March 21, 2016 5:03 PM
To: ferret_users@xxxxxxxx
Subject: Re: [ferret_users] regriding

Hi Robert,
You would need to use the regridding functions to regrid, say, one of the fields to a rectilinear grid and from there to the other curvilinear grid; or regrid both of them to the same rectilinear grid and do your calculations there.

http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/appendix-a-external-functions/CURV-TO-RECT
http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/appendix-a-external-functions/CURV-TO-RECT-MAP
http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/appendix-a-external-functions/RECT-TO-CURV

We do have in the pipeline some development which will make handling curvilinear grids more seamless, and some regridding functions which are faster; but the fact remains that you'd need to get the data onto the same grids before being able to take differences.

Ansley


On 3/21/2016 3:53 PM, Osinski, Robert FORNATL, PL wrote:
Dear ferreters,
I need your help in interpolating/regridding model outputs onto satellite measurements. I have model results of ice thickness distribution on:
sh dat
    1> ./R1009RBRcaak01a.cice.h.2004-01.nc  (default)
 name     title                             I         J         K         L
 TLON     T grid center longitude          1:1280    1:720     ...       ...
 TLAT     T grid center latitude           1:1280    1:720     ...       ...
 ULON     U grid center longitude          1:1280    1:720     ...       ...
 ULAT     U grid center latitude           1:1280    1:720     ...       ...
 TMASK    ocean grid mask                  1:1280    1:720     ...       ...
 HI       grid cell mean ice thickness     1:1280    1:720     ...       1:1

and satellite measurements:
yes? sh dat
     currently SET data sets:
    1> ./kwok_icesat.nc  (default)
 name     title                             I         J         K         L         M         N
 CELL_AREA
          area of grid cell on stereograp  1:140     1:140     ...       ...       ...       ...
 H_FM04   February-March 2004 sea ice thi  1:140     1:140     ...       ...       ...       ...
 LATITUDE latitude                         1:140     1:140     ...       ...       ...       ...
 LONGITUDE
          longitude                        1:140     1:140     ...       ...       ...       ...
 MASK     ocean-land mask                  1:140     1:140     ...       ...       ...       ...

I need to have interpolated model's output onto measurements grid to have possibility to do some calculations between them for instance to get distribution of differences between them. Attached are figures of modeled and measured ice.


Thanks in advance,
Robert Osinski
NPS,
Monterey, CA


[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement