Ocean Modeling Discussion

Line 824-825 in distribute.f90:

CALL mpi_bcast (A, Npts, MPI_INTEGER, MyMaster, OCN_COMM_WORLD, &
& MyError)

Line 219 in initial.f90:

CALL mp_bcasti (ng, iNLM, exit_flag)

Hei I'm trying to start a ROMS model run, but get problem after reading all the data from the INI file. I think the code is about to read the BRY file. This is what I get in the log file: NLM: GET_STATE - Reading state initial conditions, 2017-07-24 12:00:00.00 (Grid 01, t = 17371.5000, File ...

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 ...

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 ...

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 ...

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 ...

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

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

Hi again

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

Andre Staalstrom

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.

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 ...

Hello again

Do anyone know the answer to this one?

Staalstrom

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

staalstrom

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 ...

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 ...

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 ...

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

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 ...

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 ...