Project

General

Profile

Forums » Support »

Generation weights between two SCRIP grids

Added by Vladimir Vovo about 11 years ago

Hello!
I am trying to build weights between two SCRIP grids.
Files contain only grids, no data.
The quiestion is - how can i explain it to CDO, like this:

cdo genbil,grid1.nc grid2.nc w.nc

Instead of:

cdp genbil,grid1.nc ifile w.nc.

Because at initialization i have no data 'ifile', i have only SCRIP-grids.

Thanks.


Replies (2)

RE: Generation weights between two SCRIP grids - Added by Uwe Schulzweida about 11 years ago

You have to create a temporary data field, e.g. with the operator random or const:

cdo genbil,grid1.nc -random,grid2.nc w.nc
cdo genbil,grid1.nc -const,1,grid2.nc w.nc

RE: Generation weights between two SCRIP grids - Added by Vladimir Vovo about 11 years ago

Thanks, it helps. But than very strange thing happened..

I decided to rewrite my SCRIP weights generation(using original Fortran90 package), because CDO is much convenient.
I decided to compare results. And there are absolutely different!

First of all, I noticed strange thing - CDO 'gencon' change fields of source grid - for example imask field is different for source grid and for result weights file.
And as a result - I have different remap matrix.

I made test grids in CDO-SCRIP and original SCRIP formats(the same, but arrays are unrolled to vectors).

Build weights with command:
cdo gencon,grid2_cdo.nc -const,1,grid1_cdo.nc w_cdo.nc

You can see, that, for example, src_grid_center_lat is different for w_cdo.nc and grid1.nc.

Help me to understand it, please.

Download all files
grid1_cdo.nc (66.5 KB) grid1_cdo.nc
grid2_cdo.nc (9.05 KB) grid2_cdo.nc
grid1_scrip.nc (76.7 KB) grid1_scrip.nc
grid2_scrip.nc (19.3 KB) grid2_scrip.nc
    (1-2/2)