[ferret_users] Nearest coordinate regridding NRST

[Brockmann Patrick <Patrick.Brockmann@xxxxxx>]

Hi all,

I have a problem using the dynamic regidding NRST.
It works with a variable with a 2 year time axis but not
with a variable that has a 1 year time axis.

The goal is to extend yearly values to a variable defined monthly.

Here are the tests made.

$ ferret -script nrst.jnl var_2years.nc
use ($01)
!-> SET DAT/FORM=CDF var_2years.nc

list GPP_YEAR[gt=MONTHLY@NRST]
            VARIABLE : GPP[GT=YEARLY@AVE]
                       regrid: on T@NRST
            FILENAME : var_2years.nc
            SUBSET   : 24 points (TIME)
            LONGITUDE: 56.3E
            LATITUDE : 51.7N
                        56.3E 
                         1
16-JAN-2005 12 /  1:  0.09577
15-FEB-2005 00 /  2:  0.09577
16-MAR-2005 12 /  3:  0.09577
16-APR-2005 00 /  4:  0.09577
16-MAY-2005 12 /  5:  0.09577
16-JUN-2005 00 /  6:  0.09577
16-JUL-2005 12 /  7:  0.09577
16-AUG-2005 12 /  8:  0.09577
16-SEP-2005 00 /  9:  0.09577
16-OCT-2005 12 / 10:  0.09577
16-NOV-2005 00 / 11:  0.09577
16-DEC-2005 12 / 12:  0.09577
16-JAN-2006 12 / 13:  0.00999
15-FEB-2006 00 / 14:  0.00999
16-MAR-2006 12 / 15:  0.00999
16-APR-2006 00 / 16:  0.00999
16-MAY-2006 12 / 17:  0.00999
16-JUN-2006 00 / 18:  0.00999
16-JUL-2006 12 / 19:  0.00999
16-AUG-2006 12 / 20:  0.00999
16-SEP-2006 00 / 21:  0.00999
16-OCT-2006 12 / 22:  0.00999
16-NOV-2006 00 / 23:  0.00999
16-DEC-2006 12 / 24:  0.00999


$ ferret -script nrst.jnl var_1year.nc
use ($01)
!-> SET DAT/FORM=CDF var_1year.nc

list GPP_YEAR[gt=MONTHLY@NRST]
            VARIABLE : GPP[GT=YEARLY@AVE]
                       regrid: on T@NRST
            FILENAME : var_1year.nc
            SUBSET   : 12 points (TIME)
            LONGITUDE: 59.8E(59.8)
            LATITUDE : 53.4N
                      60E
                       1
16-JAN-2006 12 /  1:....
15-FEB-2006 00 /  2:....
16-MAR-2006 12 /  3:....
16-APR-2006 00 /  4:....
16-MAY-2006 12 /  5:....
16-JUN-2006 00 /  6:....
16-JUL-2006 12 /  7:....
16-AUG-2006 12 /  8:....
16-SEP-2006 00 /  9:....
16-OCT-2006 12 / 10:....
16-NOV-2006 00 / 11:....
16-DEC-2006 12 / 12:....

As you see, all values are set to missing
whereas I expect a variable filled by the yearly value.
Time axis seems to be correct and are defined the same ways.
Tested with ferret 6.67

Patrick

--
LSCE/IPSL, Laboratoire CEA-CNRS-UVSQ
Data Analysis and Visualization Engineer
ICMC - IPSL Climate Modelling Centre
--

use ($01)

list GPP_YEAR[gt=MONTHLY@NRST]

Attachment: var_1year.nc
Description: Binary data

Attachment: var_2years.nc
Description: Binary data