3d interpolation is not simple - like any interpolation. You have to make assumptions on the rules how to find values between source grid points. May be "trilinear" interpolation. But how do you justify this? The ocean is highly anisotropic between horizontal and vertical. Stratification is formed by mostly vertical processes, mixing is isoneutral in many cases. Hence ocean dynamics is not following linear "interpolation rules". So why to invest much work in this direction. Also, in some areas the WOA is a climatology that is supported by only a few measurements anyway. So seeking accuracy from best interpolation onto the model grid may be hopeless.
So I just use combined 2d-interpolation horizontally (bilinear with @ave) and 1d-interpolation vertically. The steps are:
- extend the original field into land by fill_xy before interpolation to the model grid. This is a 2d nearest-neighbour averaging filler. Please consult the manual. The aim is to improve interpolation results for coastal points. Nevertheless, WOA is not well supported by data points near coasts in some cases.
- map horizontally to the model grid by bilinear interpolation with @ave
- apply your model land mask
- interpolate vertically
- fill missing model (mostly bottom and surface) cells by the 1d-nearest neighbour filler of ferret @fnr. This is an extrapolation.
I am aware that your initial total nitrate content of the model may be incorrect this way. But the WOA itself has uncertaincy. Stratification from inclined isolines may be disturbed from this method. But this should be corrected during the models spin up by the model dynamics itself. At least the model should run into its own equilibrium during some time and should "forget" initial errors due to nitrate surfcace flux, denitrification, nitrification flux to sediments ...
If you are interested I can send you my scripts, just let me know.
I hope this helps a little bit,
Martin
Am 30.11.2015 09:48, schrieb Pearse Buchanan:
Hi ferreters,
I’m in the process of building some new routines within a global biogeochemical model. This means adding a climatology of Nitrate into the model’s restart file.
But I’m having some issues. Essentially,
· I want to take the World Ocean Atlas global annual mean of Nitrate, which has the following lat,lon,depth coordinates: [180,360,102]
· And I want to regrid the data into the model’s coordinate system: [112,128,21]
I have looked everywhere online for some information on 3D regridding!!!!!!!! I’ve tried downloading the Earth System Modelling Framework toolbox for regridding, but to no avail.
So, I turn to this amazing ferret community for help. Does anyone have any experience doing 3-dimensional regridding!? How can I do this seemingly simple thing?
Thank you,
Pearse J. Buchanan
PhD Candidate / CSIRO-UTAS Quantitative Marine Science
Institute for Marine and Antarctic Studies (IMAS), University of Tasmania
~~~ “Ocean deoxygenation: A Palaeo-Modelling Perspective” ~~~
Parts of the Ocean are predicted to lose oxygen in the coming century. But why?
Looking into the past, how does oxygen in the ocean change across climate transitions?
University of Tasmania Electronic Communications Policy (December, 2014).
This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.