seagrid colapse!

[perezro]

>Dear Mr. Arango
I do have a problem with seagrid, even if I am able to create ROM
> files and GETM files for small regions (gulf of Finland, union of
> north and Baltic seas, etc..) when I try to compute the whole Baltic
> (upper latitude westernmost longitude easternmost longitude lower
> latitude ) with 1x1 bathymetry, my computer collapse...
> I ask the IT responsible of the building and said that maybe if I get
> more ram for my laptop (it only have 1Gb) it could work, as you know
> the program maybe you could give me some advice?? I attach the error
> message in case it helps:
>

It is possible that you don't have enough memory in your laptop to execute larger grids. I have been using seagrid for year and it always worked for me. However, I always use Linux desktops with a lot of memory and never a laptop.

Recently, a colleague was able to get seagrid in a laptop but under Linux. There is to much limitations in the window environment. His
grid was 258x182 which is a medium size grid. His laptop has
2 Gb of memory.

The fact that you have been able to run smaller grids tells me that everything is installed correctly in your laptop. Also tells me that you know how to use seagrid.

> ??? Error using ==> qhullmx
>
> 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 0.48
> Error-roundoff 2.1e-015 _one-merge 1.4e-014 Visible-distance 4.1e-015
> U-coplanar-distance 4.1e-015 Width-outside 8.3e-015 _wide-facet
> 2.5e-014
>
> Last point added to hull was p184595.
> Last merge was #1465098.
>
> Qhull has finished constructing the hull.
>
>
>
> At error exit:
>
>
>
> precision problems (corrected unless 'Q0' or an error)
>
> 1465098
> coplanar horizon facets for new vertices
>
>
>
> Error in ==> delaunayn at 117
> t = qhullmx(x', 'd ', opt);
>
> Error in ==> griddata>linear at 150
> tri = delaunayn([x y]);
>
> Error in ==> griddata at 119
> zi = linear(x,y,z,xi,yi,opt);
>
> Error in ==> seagrid.dobathymetry at 79
> zi = griddata(x, y, z, xi, yi, theGriddingMethod); % Grid z.
>
> Error in ==> seagrid.doupdate at 110
> self = dobathymetry(self, needsBathymetry); % <== Step #5 <==
>
> Error in ==> seagrid.doevent at 93
> doupdate(self, needsUpdate)
>
> Error in ==> psevent at 25
> doevent(ps(gcbf), theEvent, theMessage);
>
> ??? Error while evaluating uimenu Callback
>

The first error are internal to matlab itself and not seagrid.
The delauna triangulation is failing when calling griddata. There is nothing that we can do about it. This is Matlab internal interpolation routines. This interpolation is very complicated.
One thing that you can check is if the arguments to these routines has a lot of NaNs. If so, it is possible that you
are trying to do weird stretching in the grid. Just do
rectangular grid and the use land/sea mask to represent the coastlines. I do not recommend to users to do extreme curvilinear grids.

> Please, I know you are busy but I would be so gratefull...
> Thanks in advance: R


Good luck, H

[perezro]

Thanks a lot for your reply!!!!
The fact is that my grid is small, 120x100 but the area is big...I mean I try grid size of 300x280 but for a smaller region and it works!! so... I don't know... maybe it is better to convince my supervisor to get matlab for the cluster
Cheers! R