Direkt zum Inhalt | Direkt zur Navigation

Personal tools

Sektionen

SCRIP "Conservative" remapping

Up to Transformations and interpolations with the SCRIP library in OASIS3-MCT

SCRIP "Conservative" remapping

Posted by Anonymous at June 22. 2015

Hi,

If you look to the interpolated field plotted on figures O_Wind10_conserv.gif and O_Wind10_conserv_zoom.gif, you can see that there are negative values which is a problem as I interpolate the wind module which should be always positive. The most evident point with a problme is located at j=283, i=95.

What can I do ?

Thanks, Aurore

Re: SCRIP "Conservative" remapping

Posted by Laure Coquart at June 22. 2015

Figure O_Wind10_conserv.gif

Attachments

Re: SCRIP "Conservative" remapping

Posted by Laure Coquart at June 22. 2015

Figure O_Wind10_conserv_zoom.gif

Attachments

Re: SCRIP "Conservative" remapping

Posted by Anonymous at June 22. 2015

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

Reply to this
URLs will be automatically hyperlinked. Basic HTML tags are OK.
Powered by Ploneboard
© Copyright ENES Portal 2011