Ocean Modeling Discussion

Hello, I tried running the default Kelvin wave test case for 10 days and the results show a growth in free surface amplitudes, starting from the eastern boundary. I am using the default boundary options that came with kelvin.h. #define RADIATION_2D #define EAST_FSRADIATION #define EAST_M2RADIATION ...

Hello, I am running ROMS on a unix cluster, which has an 8 hour walltime limit. My ROMS simulation is 270000 time steps long and the 8 hour walltime limit completes about 26000 time steps. So, I wrote a script(given below) moving the output files to another directory; copied the ocean_rst.nc to ...

I am using gfortran on my 64 bit dell machine running on redhat linux using the latest-trunk build(4.4) from the gfortran website. When I tried running the executable, I get the following error: ./oceanS < ocean.in ./oceanS: error while loading shared libraries: libgfortran.so.3: cannot open shared ...

Kate,

Thanks Kate!!! That is it.

Yes, USE mod_stepping had crept into my step2d_tile. I removed it and it solved the problem.


Sankar

Hello, I get the following compilation error, while trying to compile ROMS on a 64 bit linux using gfortran. Any ideas to sort this out would be appreciated. Thanks, Sankar In file step2d.f90:125 integer, intent(in) :: krhs, kstp, knew 1 Error: Name 'krhs' at (1) is an ambiguous reference to 'krhs ...

Hello, RAMP_TIDES does not seem to be activated in zetabc.f for SOUTH_FSCLAMPED. elif defined SOUTH_FSCLAMPED ! ! Southern edge, clamped boundary condition. ! DO i=Istr,Iend zeta(i,Jstr-1,kout)=BOUNDARY(ng)%zeta_south(i) # ifdef MASKING zeta(i,Jstr-1,kout)=zeta(i,Jstr-1,kout)* & & GRID(ng)%rmask(i ...

Hello, I am using the following southern boundary conditions for observed tidal height forcing on the southern boundary. Eventhough I use the RAMP_TIDES in my CPP options, it does not show up in the output log file. The effect of ramping is also not seen in the output. #define SOUTH_FSCLAMPED ...

Hello, I undefined the following CPP options and the 10 day 3D barotropic tidal simulation looks fine. #undef SALINITY #undef SOUTH_M3RADIATION When turn on SALINITY and SOUTH_M3RADIATION, I see large eddies near the southern boundary and the model blows up after 10 hours. The problem is with either ...

Thanks for the suggestion John. #define Masking was all that it needed: Now, I seem to have a problem with the tidal forcing. I ran the model for 10 days. But the current vectors are in all directions especially near the open boundary. (1) The maximum courant number is 0.227 (2) grid stiffness ratio ...

Thanks, John for the response. It was a great hint. I fixed the problem in the grid. Otherwise, I would be going through all seagrid code. Earlier, I had smoothed the bathymetry and while writing the smoothed bathymetry, had not followed nc format for writing the arrays. I have put the results in ...

Hello,

As a follow up to my last post, grid and mat files are available at

http://sankaranarayanans.tripod.com/prov_river_grid.nc
http://sankaranarayanans.tripod.com/prov_river_grid.mat

I would go ahead with a rectangular grid (with no inclination) and see what happens.

Thanks,

Sankar

Hello, I generated a non-rectangular grid ( rectangular grid inclined at an angle) using seagrid and ran a simple estuarine application in ROMS. The 10-day simulation using the grid ran, but the results do not look right. I have put the figures in the form of pdf in the following website. http ...

John, Yes, I had the wrong dimensions in the input file. I should have carefully read the glossary at the end of the input file. Lm Number of INTERIOR grid RHO-points in the XI-direction for ! each nested grid, [1:Ngrids]. If using NetCDF files as ! input, Lm=xi_rho-2 where "xi_rho" is the NetCDF ...

Hello, I getting the following error, while running ROMS for an estuary application, with tides and river flow. Any clues about where the problem is?. I generated the grid using seagrid and edited the grid using editmask. Thanks, Sankar OPENCDF - inconsistent size of dimension: xi_psi 59 61 ROMS ...

Thanks to John again for the fast response. The compilation went fine. Now, I am into the next step in ROMS. If my assumption is correct, for real cases, one would give the following s-coordinate parameters in ini_hydro file. How do I calculate values for these parameters, for a given number of ...

Hello, I am trying to build a real estuarine setup,closely following the estuary.h. I get the following error during the compilation. In file analytical.f90:19 END MODULE analytical_mod 1 Error: CONTAINS block at (1) is empty make: *** [Build/analytical.o] Error 1 The analytical.f90 in my build ...

Kate,

Thanks for your mail.

It looks like my nc{'hraw'}(:), nc{'bath'}(:) are both empty and this causing the problem. I can not seem to get the hraw filled.

Sankar

Kate,

I am doing something wrong somewhere while generating the grid using seagrid.

I create the grid and then compute the depth and land mask and save it as a mat file. Then, when I exit seagrid and open seagrid and load the saved grid, I do not see the masks.

Any hints?

Thanks,

Sankar

Hi Kate,

Thanks for your hint.

I still get the same error. I can see the masks in the seagrid, when I do compute depth and land mask. Still, my header file looks the same as before. I also still get the error message, ?? call to ncvarget failed, when try to editmask the grid.

Sankar

Hello all, I first loaded the amazon bathymetry, generated a coastline, and generated a sample grid. Then I used seagrid2roms to generate the .nc file. The header looks like the following: netcdf amazont_test { dimensions: xi_psi = 9 ; xi_rho = 10 ; xi_u = 9 ; xi_v = 10 ; eta_psi = 9 ; eta_rho = 10 ...

I get the following segmentation fault error while running the soliton case in cygwin. It prints all the parameters DEF_HIS - creating history file WRT_HIS - wrote history files -----------. 30 (main) ROMS_soliton 2596 _cygtls:: handle_exceptions: Error while dumping state (probably corrupted stack ...

Many thanks to Mark and John for your help with my netcdf build in cygwin. I have finally figured out that my successful netcdf build was using g77 compiler. I used g95 for compiling ROMS, removing the no second underscore option. I was successfully able to reproduce this in my desktop also. Finally ...

Mark, Many thanks for the help. I was able to run the upwelling case. Now I need to visualize the output. Is there any description of the test case and how the results look like? Do I have to compile each time for a new test case?. Now, another netcdf question: I copied the succesful libnetcdf.a and ...

Mark, Many thanks for the help. I was able to run the upwelling case. Now I need to visualize the output. Is there any description of the test case and how the results look like? Do I have to compile each time for a new test case?. Now, another netcdf question: I copied the succesful libnetcdf.a and ...

m.hadfield, Thanks for your encouragement and support. As you had suggested, I removed the -fno-second-underscore from FFLAGS and used the old libnetcdf.a. It worked and created the ROMS.exe. As I ran the ROMS.exe, it is running the default upwelling case. I get the standard ROMS/TOMS version 3.0 ...

Yes, something is broken.

I did install g77, in Cygwin setup. I tried g77 and man g77 and both worked fine.

Is there another easier route for running ROMS in pc than cygwin-g95 combination? I do have a visual fortran compiler for windows.

Thanks,
Sankar

Thanks m.hadfield. I implemented what you had suggested. When I did the make test, I get the following error message given below. I also get similar error messages in ROMS. When I do a nm on libnetcdf I get U ___assert with three underscores instead of __assert in the error message. I wrote a simple ...

I deleted everything. Then I typed m.hadfield's command as given by CC=gcc FC=g95 CPPFLAGS=-DpgiFortran FFLAGS=-fno-second-underscore FCFLAGS=-fno-second-underscore ./configure --disable-f90 --disable-cxx Then the double underscore problem in libnetcdf.a went away. However the make test for the ...

John, Thanks for your patience. If I do nm on the libnetcdf.a, I get the following. 000004b0 T _nf_abort__ 000004d0 T _nf_close__ 00000050 T _nf_create__ 000004f0 T _nf_delete__ ******more Similarly, If I do a close_io.o, I get the following: 00000000 b .bss 00000000 d .data 00000000 r .rdata ...

thanks to m.hadfield and John for their response. I tried: CC=gcc FC=g95 CPPFLAGS=-DpgiFortran FFLAGS=-fno-second-underscore FCFLAGS=-fno-second-underscore ./configure --disable-f90 --disable--cxx Now, I also see another file libnetcdff.a in addition to libnetcdf.a. Still, I do not seem to be ...

John, many thanks for your respone. You are correct, it could not find netcdf.inc file. I fixed the problem the I get the following error, it looks like it has again something to do with netcdf. ar: creating Build/libADM.a ar: creating Build/libRPM.a ar: creating Build/libTLM.a ar: creating Build ...

Hello, Many thanks to m.hadfield, jcwarner, and arango for their responses to my earlier email on CYGWIN in ROMS Now, another question from this novice. (1) I am trying to build netcdf in cygwin using g95 compiler. When I tried to configure it using ./configure I got the following error message ...

Hello, I am a very new ROMS and CYGWIN user. I am trying to compile ROMS in CYGWIN by, Make -f CYGWIN-g95.mk and I get the following message: Make: *** No targets. Stop. I beginning to learn about cygwin, cpp, and make, this week. For eg in th e cppdef.h file, it says The application CPP flag to run ...