Up to Transformations and interpolations
Dear Sophie, I would like to have your comment on the error occurred during a remap file generation (remap_conserv.F) for coupling of T63-ECHAM5 (1.8 degrees) to ORCA05-NEMO, on SX9. The error messages is below that may be related to North Pole: Using latitude bins to restrict search. Error sum wts map1:grid2_add,grid2_centroid_lat,norm_factor= 1 0.000000000000000E+000 1.000000000000000 F I use conserve remapping as below: CONSERV LR SCALAR LATITUDE 40 FRACAREA FIRST Do you have some idea how I can solve this problem? Best wishes, Wonsun
HI Wonsun, Regarding your error message, I know we get this for many cells using the SCRIP library. This is more a warning than a real error messages. However, I would advise you to look in details to the weights that the SCRIP provides. Sometimes the sum of the weights is equal to 1 but the weights themselves are completely wrong. For example, with 3 weights, one can obtain something like 369.6, -369., 0.4. The sum is indeed equal to 1 but the weights look terribly wrong. This is one of the reason that we are currently looking at other interpolation library, e.g. ESMF. I hope this helps!? With best regards, Sophie
Hi all, As Sophie already said, sometimes the sum of the weights is equal to 1 but the weights themselves are completely wrong with the SCRIP. The best way to verify the remapping quality is to use the toy oasis3-mct/examples/spoc/spoc_regridding with your grids, modifying the files data3_oasis/namcouple and run_regrid.sh. In this toy, an analytical field, calculated on the source grid, is interpolated on the target grid (following what is defined in the namcouple). The error of remapping is then calculated as the difference between the analytical field calculated on the target grid and the interpolated field. You can have an idea of the quality of the remapping by analyzing this error (file error.nc). Best regards, Laure
