ini bnd forc clim...create

Discussion of how to use ROMS on different regional and basin scale applications.

Moderators: arango, robertson

Post Reply
Message
Author
mistral
Posts: 23
Joined: Wed Jan 16, 2008 11:17 pm
Location: Privat

ini bnd forc clim...create

#1 Unread post by mistral »

hi
What is the most simple and complete tools for MATLAB to create these files???
tutorials???

gracias

slvester
Posts: 32
Joined: Mon Jun 01, 2009 12:59 pm
Location: second institution of oceanography,state oceanic administration China

Re: ini bnd forc clim...create

#2 Unread post by slvester »


mistral
Posts: 23
Joined: Wed Jan 16, 2008 11:17 pm
Location: Privat

Re: ini bnd forc clim...create

#3 Unread post by mistral »

I also use rtn but I have several problems.
for example when I create analytical grid:








% EDIT: set a name to recognize the grid
nameit='test'; % 8 km grid

% EDIT: set path where you want to store the grid.
grdfile=[nameit,'-grid.nc'];

dx=500; % m in X direction
dy=500; % m in Y direction
Lp=96; % number of points in X direction
Mp=24; % number of points in Y direction


rgrd_AnaGrid(dx,dy,Lp,Mp,grdfile);

% EDIT add new grid in the rnt_gridinfo.m
configfile=which('rnt_gridinfo');

% Set bathymetry from analytical function
%----------------------------------------
val1=(Mp/2)*dx;
val2=dx*(Mp);
Creating grid file ...
Computing (X,Y) coordinates on C-Grid
Adding Mask == 1
Computing coordinate transformation metrics 1/dx 1/dy
Computing Coriolis
Seting Topograhpy to constant 1000 m
... modify topo later as you like.
>> grd=rnt_gridload(nameit);
??? Error: File: subsref.m Line: 128 Column: 9
"subs" previously appeared to be used as a function or command, conflicting with its use here
as the name of a variable.
A possible cause of this error is that you forgot to initialize the
variable, or you have initialized it implicitly using load or eval.


Error in ==> netcdf.subsref at 52
result = subsref(result, s);

Error in ==> rnt_gridload at 59
grd.lonr=nc{'lon_rho'}(:)';lonr=nc{'lon_rho'}(:)';

>>

I modified the file rnt_gridload:


...
case 'test'
gridindo.id = gridid;
gridindo.name = 'test';
gridindo.grdfile = '/.../.../test-grid.nc';
gridindo.N = 1;
gridindo.thetas = 3;
gridindo.thetab = 0.4;
gridindo.cstfile = which('CoastlineNONE.mat');
...



how you doing?
Having problems?

sorry but i am newbie
Last edited by mistral on Thu Dec 03, 2009 2:42 pm, edited 1 time in total.

mistral
Posts: 23
Joined: Wed Jan 16, 2008 11:17 pm
Location: Privat

Re: ini bnd forc clim...create

#4 Unread post by mistral »

ok...
i solved (for previous error , not for clima.nc file)
:lol:

rnt for president :wink:

Post Reply