Problems in make_clim.m

[lassad]

Hi all,

I'm using the ROMS TOOLS Package to produce my own initial and boundary file. I've already produced my grid file using seagrid and after that I succsessfully run make_forcing. When I tried to run make_clim to produce my boundary and initial files I've got an error message when the script tried to vertically interpolate the temperature field (transcripted below).
Does anyone knows what could be my source of error or what I'm doing wrong here?

Thanks for the help and attention,

Luiz Paulo

************************************************************************

Vertical interpolations

Temperature...
Warning: Divide by zero.
> In zlevs at 60
In vinterp_clm at 57
In make_clim at 135
Time index: 1 of total: 12
??? Subscript indices must either be real positive integers or logicals.

Error in ==> ztosigma at 62
z1=depth(thezlevs);

Error in ==> vinterp_clm at 114
nc{vname}(l,:,:,:)=ztosigma(flipdim(var,1),zroms,flipud(z));

Error in ==> make_clim at 135
vinterp_clm(clmname,grdname,oaname,'temp','tclm_time','Z',0,'r');

***********************************************************************

[miguel.solano]

Hi,
Since you have an error in:

Code: Select all

% Create S-coordinate system: based on model topography h(i,j),
% fast-time-averaged free-surface field and vertical coordinate
% transformation metrics compute evolving depths of of the three-
% dimensional model grid.
%    
hinv=1./h;

and:

Code: Select all

% Get the sigma depths
%
zroms=zlevs(h,0.*h,theta_s,theta_b,hc,N,'r');

My first guess is that make_clim cannot interpolate vertically in the way you set up your grid with Seagrid. I suggest you check your vertical grid parameters, specially (hc and hmin), and check the vertical grid with test_vgrid. You can also create your grid using make_grid and make sure hc < hmin.