Up to Transformations and interpolations
Hi, I interpolate the wind module which should be always positive but I have some negative values in the interpolated field which is a problem. The most evident point with a problme is located at j=283, i=95. What can I do ? Thanks, Aurore
Hi Aurore, The problem seems to be located along a constant latitude near to the North pole. In oasis3-mct/lib/scrip/src, in remap_conserv.F, could you replace north_thresh=1.45_dbl_kind by north_thresh=2.00_dbl_kind ? ie replaced : real (kind=dbl_kind), parameter :: & north_thresh = 1.45_dbl_kind, ! threshold for coord transf. & south_thresh =-2.00_dbl_kind ! threshold for coord transf. by : real (kind=dbl_kind), parameter :: & north_thresh = 2.00_dbl_kind, ! threshold for coord transf. & south_thresh =-2.00_dbl_kind ! threshold for coord transf. In the SCRIP library, if north_thresh = 1.45_dbl_kind, that means that there will be a projection for latitudes superior to lat= 1.45 rad. If you put north_thresh=2.00_dbl_kind, there will be no more projection. Let me know if it solves your problem. Best regards, Laure
