Levitus_Error

[QAYETBEY@GMAIL.COM]

Hello MOHIDers
I have got this message:
Reading constructions...
Type of Interpolation : Bilinear
Constructing father grid...
Constructing new grid...
Constructing aux grid...
Constructing communication between grids...
Number of instants in available data : 12
reading salinity fields
Interpolating : salinity_00001 2
Interpolating : salinity_00001 2
Interpolating : salinity_00001 2
Interpolating : salinity_00001 2
Interpolating : salinity_00001 2
.
.
.
Interpolating : salinity_00001 2

All Triangulation Nodes are collinear
VerifyCollinearPoints-ModuleTriangulation-ERR01

Any idea about that ?!

[davidbrito]

Hello,

Please specify what tool are you using and what are you trying to do.
If you send the input data file also may help.


Best regards,
David

[QAYETBEY@GMAIL.COM]

Hi David
Thank you for your response.
I wanna feed MOHID with temperature and salinity values from Levitus as HDF5 file using ( ConverToHDF5_release_double ) Tool.
The program reads the Levitus files well , but gives the above error message on beginning interpolation. I use high resolution (1/4 degree) Temperature and Salinity Analyses of the World's Oceans.,Version 2 ( Gridded Fields) and the ConvertToHDFAction file looks like this :

<begin_file>
ACTION : CONVERT LEVITUS FORMAT

OUTPUTFILENAME : Levitus.hdf5
OUTPUT_GRID_FILENAME : LevitusGrid.dat
OUTPUT_GEOMETRY_FILENAME : LevitusGeometry.dat

PERIODICITY : monthly
SPATIAL_RESOLUTION : 0.25
FILL_VALUE : -99.9999

LOWER_LEFT_CORNER : 32.1 26.76
UPPER_RIGHT_CORNER : 34.41 30.092

<<beginfield>>
NAME : salinity
ANNUAL_FILE : D:\Mohid_Wks\ConvertToHDF5\LEVITUS\s000AV

<<<begin_input_files>>>
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s001
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s002
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s003
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s004
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s005
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s006
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s007
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s008
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s009
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s010
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s011
D:\Mohid_wks\ConvertToHDF5\LEVITUS\s012
<<<end_input_files>>>
<<endfield>>

<<beginfield>>
NAME : temperature
ANNUAL_FILE : D:\Mohid_Wks\ConvertToHDF5\LEVITUS\t000AV


<<<begin_input_files>>>
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t001
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t002
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t003
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t004
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t005
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t006
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t007
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t008
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t009
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t010
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t011
D:\Mohid_wks\ConvertToHDF5\LEVITUS\t012
<<<end_input_files>>>
<<endfield>>
<end_file>


<begin_file>
ACTION : INTERPOLATE GRIDS

TYPE_OF_INTERPOLATION : 1

FATHER_FILENAME : Levitus.hdf5
FATHER_GRID_FILENAME : LevitusGrid.dat

OUTPUTFILENAME : LeviTusAllPointsWithData.hdf5
NEW_GRID_FILENAME : LevitusGrid.dat

START : -9999 1 1 0 0 0
END : -9999 12 1 0 0 0

INTERPOLATION3D : 1
FATHER_GEOMETRY : LevitusGeometry.dat
NEW_GEOMETRY : LevitusGeometry.dat
AUX_GRID_FILENAME : LevitusGrid.dat
AUX_OUTPUTFILENAME : AuxLeviTusAllPointsWithData.hdf5

POLI_DEGREE : 3
DO_NOT_BELIEVE_MAP : 1
EXTRAPOLATE_2D : 2
<end_file>


<begin_file>
ACTION : INTERPOLATE GRIDS

TYPE_OF_INTERPOLATION : 1
FATHER_FILENAME : LeviTusAllPointsWithData.hdf5
OUTPUTFILENAME : Levitus_out.hdf5
FATHER_GRID_FILENAME : LevitusGrid.dat
NEW_GRID_FILENAME : LevitusGrid_out.dat

START : -9999 1 1 0 0 0
END : -9999 12 1 0 0 0

INTERPOLATION3D : 1
FATHER_GEOMETRY : LevitusGeometry.dat
NEW_GEOMETRY : Geometry_1.dat
AUX_OUTPUTFILENAME : AuxGoSMonthlyLevitus.hdf5
AUX_GRID_FILENAME : Aux12km.dat

POLI_DEGREE : 3
DO_NOT_BELIEVE_MAP : 1
<end_file>

[davidbrito]

Hello,

So as a start this is not a Levitus Error, but interpolation. That can be related to the file or not.
the first thing to do is separate the processes so that you know in which it failed (it was the first interpolation that failed?)

If you look to the code error it means that was an error on routine VerifyCollinearPoints inside ModuleInterpolation. And the error number inside that routine is 1. You can search this way every error code you got.

And it seems that you have collinear points. So we have to see that after it is separated the process that gives the problem.


Best regards,
David

[QAYETBEY@GMAIL.COM]

Hello David
Thank You for Your Reply.
I already separated processes and it stopped upon interpolating (second process). ok, I am going to check it as you suggest but ConvertToHDF5_release_double is not present in the repository !
Do not you have a facility for debugging ?

[davidbrito]

Hello,

the ConvertToHDF5 code is is in Codeplex in software.

So in your first interpolation the model found collinear points to interpolate.
Have you made other interpolations successfull?
If yes try those settings first.

Also you may try other interpolation method (not triangulation) to see if you get it done (and the problem may be in triangulation) or if no methos works and maybe the problem is the files you are trying to interpolate.

Best regards,
David

[QAYETBEY@GMAIL.COM]

Hello David
Thank You for Your Reply.There is another problem appeared
ConstructFatherGridLocation-ModuleHorizontalGrid-ERR10
Any Idea?!!

[davidbrito]

Hello,

This error appeared in the same run or after you changed the interpolation method?


But just some remarks about your interpolation:
- you are performing interpolations from one grid FATHER_GRID_FILENAME to the NEW_GRID_FILENAME. And you are using the same file LevitusGrid.dat. Interpolation in this case (in ConvertToHDF5) is a process to get HDF5 values in one grid and interpolate them to other grid (coarser or finer). From one grid to the same you do not need to do nothing. Probably that is why it gave collinear points error.
- You define year of start and end with -9999 I do not know what the model does with that. The number -99 or -9999 only means something if it is a grid value. In grids -99 will not be considered but in dates I do not think the model will use as it. it will use year -9999.

So if you can go back to the time where the error was collinear points that proabbly is because you are not interpolating between grids.
If you want interpolation to do something you have to provide two grids, the original and the new.

Probably you do not need Interpolation (you do not have a new grid) and you are done!


Best regards,

David

[QAYETBEY@GMAIL.COM]

Dear David
Sorry, I forget to tell you that this error appeared while nesting a model.

[davidbrito]

ok but you are giving the program the father and son the same grid.
and they need to be different or interpolation does not make sense.

Best regards,
David