Hi
When I run 'make_clim' I face following error:
Compute geostrophic currents
time index: 1 of total: 12
??? qhull precision warning:
The initial hull is narrow (cosine of min. angle is 1.0000000000000002).
A coplanar point may lead to a wide facet. Options 'QbB' (scale to unit box)
or 'Qbb' (scale last coordinate) may remove this warning. Use 'Pp' to skip
this warning. See 'Limitations' in http://www.qhull.org/html/qh-impre.htm
qhull precision error: initial simplex is not convex. Distance=-1.6e-14
While executing: | qhull d Qt Qbb Qc
Options selected for Qhull 2003.1 2003/12/30:
delaunay Qtriangulate Qbbound-last Qcoplanar-keep _pre-merge
_zero-centrum Pgood Qinterior-keep _max-width 28 Error-roundoff 1.5e-13
_one-merge 1.1e-12 Visible-distance 3e-13 U-coplanar-distance 3e-13
Width-outside 6e-13 _wide-facet 1.8e-12 _narrow-hull -2.2e-16
precision problems (corrected unless 'Q0' or an error)
1
flipped facets
The input to qhull appears to be less than 3 dimensional, or a
computation has overflowed.
Qhull could not construct a clearly convex simplex from points:
The center point is coplanar with a facet, or a vertex is coplanar
with a neighboring facet. The maximum round off error for
computing distances is 1.5e-13. The center point, facets and distances
to the center point are as follows:
facet
p176
p27612
p0
distance= -5.3e-15
facet
p27788
p27612
p0
distance= -1.2e-14
facet
p27788
p176
p0
distance= -5.3e-15
facet
p27788
p176
p27612
distance= -1.8e-15
These points either have a maximum or minimum x-coordinate, or
they maximize the determinant for k coordinates. Trial points
are first selected from points that maximize a coordinate.
The min and max coordinates for each dimension are:
0: 99 125 difference= 26
1: 2.166 30.03 difference= 27.87
2: 0 27.87 difference= 27.87
If the input should be full dimensional, you have several options that
may determine an initial simplex:
- use 'QJ' to joggle the input and make it full dimensional
- use 'QbB' to scale the points to the unit cube
- use 'QR0' to randomly rotate the input for different maximum points
- use 'Qs' to search all points for the initial simplex
- use 'En' to specify a maximum roundoff error less than 1.5e-13.
- trace execution with 'T3' to see the determinant for each point.
If the input is lower dimensional:
- use 'QJ' to joggle the input and make it full dimensional
- use 'Qbk:0Bk:0' to delete coordinate k from the input. You should
pick the coordinate with the least range. The hull will have the
correct topology.
- determine the flat containing the points, rotate the points
into a coordinate plane, and delete the other coordinates.
- add one or more points to make the input full dimensional.
This is a Delaunay triangulation and the input is co-circular or co-spherical:
- use 'Qz' to add a point "at infinity" (i.e., above the paraboloid)
- or use 'QJ' to joggle the input and avoid co-circular data
Error in ==> delaunayn at 81
t = qhullmx(x', 'd ', opt);
Error in ==> griddata>linear at 140
tri = delaunayn([x y]);
Error in ==> griddata at 109
zi = linear(x,y,z,xi,yi,opt);
Error in ==> geost_currents at 181
u_r(k,:,:)=griddata(lon(D),lat(D),u_r(k,D),lon,lat);
Error in ==> make_clim at 148
geost_currents(clmname,grdname,oaname,frcname,zref,ob
Everybody can help me tank so much
Best Regards
F.Daryabor