Search found 28 matches

by staalstrom
Fri Apr 24, 2015 1:33 pm
Forum: ROMS Discussion
Topic: Passive dye
Replies: 1
Views: 735

Passive dye

Hei I try to run with the option T_PASSIVE, but get the error message: ERROR: Abnormal termination: NetCDF INPUT. REASON: No error INQUIRE - unable to find requested variable: dye_01_sflux in files: subarctic-20km_AN_2007_unlim.nc subarctic-20km_FC_2007_unlim.nc tide.nc subarctic-20km_rivers_dye.nc ...
by staalstrom
Wed Dec 07, 2011 8:09 am
Forum: ROMS Discussion
Topic: Which GLS scheme is the best to model wind on surface water?
Replies: 3
Views: 1626

Re: Which GLS scheme is the best to model wind on surface wa

Hi I'm not an experienced turbulence closeure modeller. I have always used k-kl, but I have plans to try the others. I work in a fjord in Norway, the Oslofjord, with max depth of 200m and sill depth of 20m. It is clear that what kind of topography and vertical coordinates I use is important for the ...
by staalstrom
Wed Jul 13, 2011 10:18 am
Forum: ROMS Discussion
Topic: Hydrostatic consistency and the rx0 rx1 factors
Replies: 4
Views: 3226

Re: Hydrostatic consistency and the rx0 rx1 factors

Hi I study internal tides in a sill fjord (Oslofjord, Norway). The topography is complicated, and I use a horizontal resolution of 75m. In this application I get currents of up to 9 cm/s when I run a test case with a strong pycnocline at 20 m depth (same as the sill depth) and no horizontal gradient...
by staalstrom
Wed May 04, 2011 7:01 am
Forum: ROMS Discussion
Topic: Convert sigma coordinate to Z depth in ferret
Replies: 4
Views: 2777

Re: Convert sigma coordinate to Z depth in ferret

Hi It is necessary to calculate z_r (3D) and z_w (also 3D) from the stretching function (C_s and C_w). C_s and C_w is calculated in Utility/set_scoord.F and is constant in the horizontal directions. 4 different versions are available now (after roms-3.4). But the stretching functions was already sto...
by staalstrom
Fri Feb 04, 2011 8:32 am
Forum: ROMS Sediment
Topic: About the hc ,Tcline and hmin
Replies: 5
Views: 4538

Re: About the hc ,Tcline and hmin

Hello again

And one more question:
Is it necessary to recompile the code to use Vtransform=2?
Or is it enough to set it in ocean.in?

I use roms-3.2 by the way.

brds
Staalstrom
by staalstrom
Thu Feb 03, 2011 1:42 pm
Forum: ROMS Sediment
Topic: About the hc ,Tcline and hmin
Replies: 5
Views: 4538

Re: About the hc ,Tcline and hmin

Hello I have more questions conserning this subject. I want to set hc>hmin. hmin is now 2m and I will now try hc=20. I use Vtransform=2 and Vstreching=1. I make initial condition by using hc=20 and set Tcline=20 and hc=20 in this file. In ocean.in I set Tcline=20. But when I run my application I get...
by staalstrom
Tue Dec 07, 2010 11:35 am
Forum: User Applications
Topic: LTRANS: Larval TRANSport Lagrangian model
Replies: 39
Views: 33092

Re: LTRANS: Larval TRANSport Lagrangian model

Hi I have just downloaded LTRANS v.1 and try to compile it with gcc. I dont think I have any netcdf problems. This are some of the errors I got. Hydrodynamic hydrodynamic_module.f90:12.126: m,rho_nodes,u_nodes,v_nodes,max_rho_elements,max_u_elements,max_v_elements,rho_ 1 Error: Symbol 'rho_el' refer...
by staalstrom
Thu Sep 30, 2010 7:25 am
Forum: ROMS Ecosystem
Topic: how to add a continuous dye?
Replies: 4
Views: 3452

Re: how to add a continuous dye?

Hi One misprint in my last reply: I meant that you need to define T_PASSIVE to get a passive tracer. If you for instance spesify north_dye_01 in the boundary input file, I guess this will be the concentration at the boundary if the boundary condition works as it should. But you can have more timeste...
by staalstrom
Wed Sep 29, 2010 9:19 am
Forum: ROMS Ecosystem
Topic: how to add a continuous dye?
Replies: 4
Views: 3452

Re: how to add a continuous dye?

Hello This is something I also wonder about I would define continous dye as a continous source of any tracer. It is easy to amake a variable dye_01 in the INI fil. If you include UV_PSOURCE and TS_PSOURCE you transport the initial amount of dye around in your model area. To my knowledge a continous ...
by staalstrom
Tue Sep 28, 2010 7:54 am
Forum: ROMS Problems
Topic: M2FLATHER and M2CLM_NUDGIN
Replies: 0
Views: 791

M2FLATHER and M2CLM_NUDGIN

Hei I'm working on an application on the western coast of Norway and I have problems with the open boundary conditions. ubar and vbar have too high values close to the border. So I want to try nudging close to the boundary. I have made a nudging coefisient in ana_nudgcoef.h. I use these option on th...
by staalstrom
Tue May 25, 2010 9:44 am
Forum: ROMS Discussion
Topic: How to calculate vertical velocity of barotropic tide
Replies: 4
Views: 3757

Re: How to calculate vertical velocity of barotropic tide

Hi I think this is very interesting. I look at internal waves in a sill fjord. ROMS seems to create internal waves without any parametrization. The problem is to describe the sill without introducing to much pressure gradient errors because of the steep topography. What does the usual conversion for...
by staalstrom
Tue Dec 15, 2009 12:45 pm
Forum: ROMS Problems
Topic: Problems with stations.in
Replies: 4
Views: 1050

Re: Problems with stations.in

Hi again I managed to make it work in a horrible way. I went into ROMS/Utility/def_stations.F and changed DEFINE : IF (ldef) THEN CALL netcdf_create (ng, iNLM, TRIM(ncname), ncSTAid(ng)) to DEFINE : IF (.not.ldef) THEN CALL netcdf_create (ng, iNLM, TRIM(ncname), ncSTAid(ng)) and now the ocean_sta.nc...
by staalstrom
Tue Dec 15, 2009 8:22 am
Forum: ROMS Problems
Topic: Problems with stations.in
Replies: 4
Views: 1050

Re: Problems with stations.in

Hi Yes, according to ocean.out LDEFOUT=F. But in ocean.in LDEFOUT is set to True. There is no problems creating a new history file. I see from some old runs without STATIONS that LDEFOUT=F even if it is set to T in ocean.in. Is it set to F somewhere in the code? Best regards Andre Staalstrom (ans@ni...
by staalstrom
Mon Dec 14, 2009 1:13 pm
Forum: ROMS Problems
Topic: Problems with stations.in
Replies: 4
Views: 1050

Problems with stations.in

Hello I want to write data from 2 positions to station output files. I have defined stations in my cpp options. The right staions.in file is spesified in ocean.in But ROMS thinks ocean_sta.nc should already exist! How do I make one? Here's the output: STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_...
by staalstrom
Fri Dec 04, 2009 4:36 pm
Forum: ROMS Installation
Topic: Problems compiling ROMS NetCDF problem
Replies: 2
Views: 1281

Problems compiling ROMS NetCDF problem

Hello I try to compile ROMS on a Dell Laptop with Ubuntu Linux. netcdf-3.6.2 is installed When I try to compile I get the following error: makefile:227: INCLUDING FILE Build/make_macros.mk WHICH CONTAINS APPLICATION-DEPENDENT MAKE DEFINITIONS cd Build; /bin/g95 -c -fno-second-underscore -O3 -ffast-m...
by staalstrom
Mon Aug 17, 2009 2:28 pm
Forum: ROMS Problems
Topic: Blow up due to stratification
Replies: 0
Views: 1090

Blow up due to stratification

Hei I am running a fjord application, with one open boundary in the west. The horisontal grid size is 100m, an the model has 20 layers. One TS profile is used as initial condition for salt and temp (the stratification is homogeneous in the horisontal). u, v, ubar, vbar and zeta is set to zero in the...
by staalstrom
Mon May 25, 2009 7:00 am
Forum: ROMS Problems
Topic: start time for river is greater than current model time??
Replies: 8
Views: 1832

Re: start time for river is greater than current model time??

Hei For some reason I prefer to edit the cdl-files manually. Probably because I never took time to look into all the Roms matlab tools. Anyway, I think I understood from what you wrote that there should be a something called cycle_length in the the river_time variable in the forcing nc file. So I ad...
by staalstrom
Fri May 22, 2009 1:07 pm
Forum: ROMS Problems
Topic: start time for river is greater than current model time??
Replies: 8
Views: 1832

Re: start time for river is greater than current model time??

Hi again

Downgrading to ROMS3.1 did not solve the problem!
In my varinfo.dat there is variable called "river".

Andre Staalstrom
by staalstrom
Fri May 22, 2009 12:27 pm
Forum: ROMS Problems
Topic: start time for river is greater than current model time??
Replies: 8
Views: 1832

Re: start time for river is greater than current model time??

Hi
I also have this problem.
I also have the variable "river" specified in the forcing file. I have 16 rivers so the river identification numbers goes form 1 to 16.
But I still get this problem.
by staalstrom
Mon May 04, 2009 9:21 am
Forum: ROMS Usage
Topic: Eigenvalues and eigenvectors
Replies: 0
Views: 840

Eigenvalues and eigenvectors

Hello I have noticed that in ocean.in: ! Number of eigenvalues (NEV) and eigenvectors (NCV) to compute for the ! Lanczos/Arnoldi problem in the Generalized Stability Theory (GST) ! analysis. NCV must be greater than NEV (see documentation below). NEV = 12 ! Number of eigenvalues NCV = 10 ! Number of...
by staalstrom
Tue Apr 14, 2009 8:16 am
Forum: ROMS Problems
Topic: Forcing in ROMS3.0
Replies: 4
Views: 1080

Re: Forcing in ROMS3.0

Hello again

Do anyone know the answer to this one?

Staalstrom
by staalstrom
Thu Apr 02, 2009 1:03 pm
Forum: ROMS Installation
Topic: Compiling ROMS-3.2 with modifications in ana_nudgcoef.h
Replies: 2
Views: 1185

Re: Compiling ROMS-3.2 with modifications in ana_nudgcoef.h

Thanks.
I used and outdated version of ana_nudgcoef.h
When I do my edits in an updated version, everything works fine.

staalstrom
by staalstrom
Wed Apr 01, 2009 11:40 am
Forum: ROMS Installation
Topic: Compiling ROMS-3.2 with modifications in ana_nudgcoef.h
Replies: 2
Views: 1185

Compiling ROMS-3.2 with modifications in ana_nudgcoef.h

I have compiled ROMS-3.2. The application is a fjord (The Oslofjord) and has one open boundary in the south and I want to use enhanced climatological nudging near the boundary and no nudging in the interior. So I want to make some modifications in the ana_nudgcoef.h. #elif defined IOSLO Lh=user(1) !...
by staalstrom
Thu Jan 15, 2009 9:30 am
Forum: ROMS Problems
Topic: Forcing in ROMS3.0
Replies: 4
Views: 1080

Forcing in ROMS3.0

Hi, yesterday i finally checked out roms3.0 I have used roms2.2 up till now. I managed to compile the code. In ocean.in I have specified 2 forcing files like this NFFILES == 2 ! number of forcing files FRCNAME == oslo150m_tide_zero.nc \ oslo150m_wind_zero.nc But I got this message when i run the cod...
by staalstrom
Thu Jan 15, 2009 8:10 am
Forum: ROMS Discussion
Topic: convert sigma coordinate to depth coordinate
Replies: 5
Views: 6487

Re: convert sigma coordinate to depth coordinate

Hallo Here's a matlab script to convert from sigma to z coordinates. % USAGE: z=zdepth_roms(h,hc,theta,b,N) % z z-coordinate center of layer (m) % h depth (m) % hc = 20; % theta = 5; % b = 0.6; % N number of layers function z=zdepth_roms(h,hc,theta,b,N) % Center of layers ds=1/N; s = -1+ds/2:ds:0-ds...
by staalstrom
Mon Aug 18, 2008 7:44 am
Forum: ROMS Discussion
Topic: Segmentation Fault
Replies: 8
Views: 1558

Re: Segmentation Fault

Hei
I dont know much about memory and computer architecture.
All I know is that I run my ROMS applications on a Linux machine with Red Hat 9 that is siting somewhere in the basement.
I always get a segmentation fault if I dont write

ulimit -s unlimited

Andre Staalstrom
by staalstrom
Fri May 30, 2008 10:48 am
Forum: ROMS Source
Topic: Point sources of passive tracer
Replies: 3
Views: 3449

Point sources of passive tracer

Hello I run roms 2.2. The model is working with passive tracer dye_01 initiated in the ini NetCDF file. But I want the dye to be put into the model with a river. TS_MPDATA, TS_DIF2, TS_PSOURCE and UV_PSOURCE is defined inn cppdefs.h. I added this to varinfo.dat 'river_dye_' ! Input 'river runoff dye...
by staalstrom
Mon Feb 04, 2008 4:57 pm
Forum: ROMS Source
Topic: Passive tracers
Replies: 2
Views: 4711

Passive tracers

Hi I'm a new user to this forum and also a new ROMS user. We have set up the model for a part of the Norwegian Western coast, and we want to find out the relation of the water between different fish farms. So I want to simulate a passive tracer. I guess i define T_PASSIVE, set NPT=1 (one passive tra...