Version 4 (modified by omamce, 7 years ago) (diff) 

Bug on linear interpolation
Normally, the wind stress is interpolated with the bilinear interpolation scheme of the Script library included in OASISMCT.
In IPSLCM6 at resolution ORCA1 x LMD 144x144, the bilinear interpolation of the wind stress is bugged. Two NEMO points around 0°N/180°W takes there wind stress values from LMDZ points around 0°N/0°E, and the weights, base on the distance between the source and destination points, are totally wrong..
The bilinear interpolation of OASISMCT compute weights and adresses that are written in a file : rmp_tlmd_to_uorc_BILINEAR.nc, rmp_tlmd_to_torc_BILINEAR.nc or rmp_tlmd_to_torc_BILINEAR.nc.
A python program has been developed. It makes a few corrections of adresses and weights in a copy of the file.
The program will be available soon.
To run the correction, the steps are:
 Run the model with bilinear interpolation of the wind stress ;
 Check the wind stress fields on the NEMO grid ;
 If they are wrong, get the weigths files rmp_tlmd_to_[uv]orc_BILINEAR.nc ;
 Correct them with the program. You may need to adapt the program to the bug you observed :(
 Put them in the IGCM repository where the weights file are ;
 Use MOZAIC interpolation for the wind stress, using these new weights files.
7 Check the result ;)
Attachments (1)

BugOasisBilinear.png
(109.2 KB) 
added by omamce 7 years ago.
Wind stress on ocean grid, U point
Download all attachments as: .zip