Ocean Modeling Discussion

ocecept

Hi Mahmud; You have water entering the gulf in one side and leaving it in the other side. This horizontal shear can induce the eddy. I have seen it when studding the Spencer Gulf, Australia but it is a much larger gulf. The inverse circulation in Spencer Gulf is in a geostrophic adjustment which ...

ocecept

Dear Colleagues and IWMO Friends, We are very pleased to present the second announcement for the 10th International Workshop on Modeling the Ocean (IWMO2018) to be held in Santos, Brazil, on June 25-28 2018. The IWMO 2018 website is now open for registrations and abstract submissions: https ...

ocecept

Hi John; Yes. I missed the placeholder in my ana_passive. Thanks I was not sure if the rivers should have river_dye_01 and river_dye_01_age or river_dye_01 _02 _03 _04. Now it makes sense and the models is running fine. Thanks again; Carlos Look at the placeholder code in ana_passive.h # ifdef AGE ...

ocecept

Hi all; John, first of all, thanks for setting the AGE_MEAN option. I have a question about river forcing. Let's say I have two dye in my river forcing, do I need to have dye_01_age and dye_02_age in the river forcing? If I set NPT = 4 I get the following error: INQUIRE - unable to find requested ...

ocecept

Hi Milena; Thanks for that. I will read the paper. No, I'm not Carlos Lentini. He is a frind of mine who works at Universidade da Bahia. I'm Carlos Teixeira. Thanks again Carlos Teixeira Hi Carlos, I used the T and S diagnostics (in ROMS Agrif, but there shouldn't be a difference with Rutgers ROMS ...

ocecept

Can anyone suggest me a paper where the diagnostics for temperature and salinity has been used?

I couldn't find any.

Manu thanks

Carlos

ocecept

Hi John, thanks for your replay. I looked into output.F when I got the zeros in my output and realize that the code didn´t check for NTSAVG, but I didn´t know if it was done in elsewhere. Thanks for explaining. Anyway, since I want the averages results after two years of simulation I think it´s be ...

ocecept

Dear All: I set NTSAVG in order to accumulate the output time-averaged and diagnostics data after two years. However, whatever number of NTSAVG is used, the output of average data is always recorded since the first time step every NAVG and filled with zeros. Then, after NTSAVG (two years in my case ...

ocecept

Regarding the first problem, I would recommend you to save the output every time step to see where the problem occurs. Also post the configurations you have used in the model (your *.h file), it will let us to help you. Regarding the second problem, the model is telling you the problem: /home/nils ...

ocecept

Hi Ban; I'm wondering how many grid points (i x j) do you have? Probably you already know that, but have a look on the Herman comments about partitions and number of cores at https://www.myroms.org/forum/viewtopic.php?t=1979 Your question came in a good time, I'm going to buy a computer with Intel ...

ocecept

Herman;

I haven't changed set_avg.F or set_diags.F. I don't know what is happening

. I got the results that I need using DIAGNOSTICS. Lets leave it and see if someone else have the same kind of problem in the future. Thanks a lot.

ocecept

Herman; I think I have found the problem, but not the solution: :roll: If I have DIAGNOSTICS_UV and DIAGNOSTICS_TS (and AVERAGES) defined it works and I get nonzero for the 3D field of the potential vorticity. If these options are NOT defined (only AVERAGES is defined) I don't get the 3D potential ...

ocecept

Thanks Herman;

I will do some tests with different applications, to see where is the problem.

Cheers;

Carlos

ocecept

Hi guys; Anyone have used the new logic to get the Potential Vorticity on the average files (see Herman post https://www.myroms.org/projects/src/ticket/421). I updated the ROMS code yesterday, my *.in and varinfo.dat files, but the my 3D potential vorticity results are all zeros. The 2D field looks ...

ocecept

You can find the Mexnc and SNCTOOLS for Matlab 2009a (and later) at: http://mexcdf.sourceforge.net/downloads/2008b.php And some instructions on how to install it at: https://www.myroms.org/wiki/index.php/MEXNC Install and test it. If everything is Ok, try to use Seagrid. Note that now, Matlab is ...

ocecept

Hi Raul;

Have you read the information at https://www.myroms.org/wiki/index.php/Tidal_Forcing?

I had followed it when creating my tidal forcing with no big problems.

Cheers;

Carlos Teixeira

ocecept

Hi Angelo; In your configuration the model can't find the varinfo.dat file. You need to tell ROMS where to find this file. Have you changed the path (VARNAME) to your varinfo.dat into your *.in file? In my case I had to change VARNAME to VARNAME = /mnt/msa20/homes/teixec01/Oceanography/CARLOS/myroms ...

ocecept

Thanks for the corrections, you fixed one of the problems (the typo problem on step_3d_t.F). We still have something wrong . Now when I try to run the model I get the following error: forrtl: severe (408): fort: (2): Subscript #3 of the array DIAU2RHS has value 11 which is greater than the upper ...

ocecept

Hermann, Thanks for the improvements on the diagnostic output. I had problems running the model after the upgrade. First, it seems there is a typo on step_3d_t.F line 1043 Instead of DO idiaf=1,NDT , you should have DO idiag=1,NDT I fixed that but it is still not working. Running the model in the ...

ocecept

Hi; If you haven't done yet, have a look at https://www.myroms.org/wiki/index.php/Tidal_Forcing . It will help you. First, what do you call a large grid? What are you using to force the model? TPOX? Which version? They have a new version TPOX7.2 that came out few days ago ( http://www.coas ...

ocecept

Hi Bardur, here is one example of my *.h for a 2D - tidal model. I'm using the OSU / TPX 7.2 global model (see https://www.myroms.org/wiki/index.php/Tidal_Forcing for detais) to force the model. That's different from what you intend to do, but maybe it will help: #define UV_COR #define UV_LDRAG ...

ocecept

Hi John; Thanks for your replay. That's the problem that I told you about in the ROMS workshop at Sydney. I used the Matlab routines found in http://marine.rutgers.edu/~hunter/roms/tides/otps/ to create my tidal forcing. That's the details of the forcing file. netcdf tides_extended_gulfs_grd ...

ocecept

I'm trying to use the AVERAGES_DETIDE option, but I got the following error in the first time-step. DEF_VAR - unable to define variable: CosWCosW in NetCDF file: tides_extended_gulfs_grd.nc ROMS/TOMS - Output NetCDF summary for Grid 01: ROMS/TOMS - Output error ............ exit_flag: 3 ERROR ...

ocecept

Hi Maria; Have a look in the following wiki page https://www.myroms.org/wiki/index.php/Tidal_Forcing . It will answer most of your questions. You will find the software to create the tidal forcing from global (with luck regional) models. I'm using the OSU Tidal Prediction and the results are quite ...

ocecept

Hi Leon; Your problem is probable related with your time step since your Maximum barotropic Courant Number = 2.43502746E+00. Try to reduce it to get Maximum barotropic Courant Number < 0 and satisfy the CFL criteria. Have a look in arangos comments at https://www.myroms.org/forum/viewtopic.php?f=1&t ...

ocecept

You can use passive/inert tracers, but I'm not sure about the decay rate.

Reading the following posts will help:

viewtopic.php?f=4&t=789&hilit=passive

viewtopic.php?f=14&t=1051&hilit=passive

Cheers

Carlos

ocecept

Hi Robertson/ Kate; thanks for your comments. Sorry for the late replay I was out for some days. Yes. The NECTCDF4 does take longer to run, but it reduces the size of my output files a lot. I did some tests... For the upwelling test case using netcdf3 the ocean_his file is 29 MB, if I use netcdf4 it ...

ocecept

Robertson, thanks for your replay. The problem was in my Linux-ifort.mk, it's not a ROMS default file and it was not updated. :oops: I had added in the netcdf4 and hdf5 paths , but I missed something else. I started it again following Darwin-ifort.mk and now everything works fine. Sorry I'm always a ...

ocecept

Hi guys; I'm trying to have the NETCFD4 library working in serial and I get some kind of problem. That's what I did: I installed the NETCFD4/HDF5 libraries for the ifort compiler with no problems, all tests worked and I used h5dump to double check if I have a netcdf4 file. I added netcdf4 and hdf5 ...

ocecept

Thanks for your replay; Some times I look part of the code and the solution looks simple, but I'm sure the Kernel is a lot more complex than it :oops: . That's ok.. I will try another solution 8) . Now let me explain my idea: I'm modeling an inverse estuary (Spencer Gulf, South Australia), where we ...

ocecept

Hi guys; I'm running a test using TS_FIXED to get the free-surface and momentum that balances a specific temperature and salinity field. That work's fine. I tried to used T_PASSIVE to get some dyes working. Unfortunately TS_FIXED fix all the tracers and not only salinity and temperature as I ...

ocecept

Hi Daryabor; >>[z,sc,Cs]=scoord(40,6,0.4,10,15,0,0,1,1) >>??? Error using ==> set >>Bad value for axes property: 'XLim' >>Values must be increasing and non-NaN. Instead of 40 you have to enter your depth matrix. [z,sc,Cs]=scoord( h ,theta_s,theta_b,Tcline,N,kgrid, column , index ,plt); if column = 1 ...

ocecept

Hi Daryabor; I started using ROMs few months ago, so I'm not an expert. I will tell you what I did and people are free to correct it. First read these posts https://www.myroms.org/forum/viewtopic.php?f=14&t=76&hilit=scoord ; https://www.myroms.org/forum/viewtopic.php?f=14&t=612&hilit=haney and https ...

ocecept

Hi Daryabor;

I have used scoord.m to test these coordinates parameters. It's located in http://www.myroms.org/software/Processi ... lab.tar.gz.

Cheers

ocecept

Hi all. First thanks for EASYGRID. I was able to create the grid to my area, smooth the bathymetry and mask the coast with no problems. My grid has 955 Km Width, 755 Km Length, 6 Km resolution and is rotated -33 degrees. Now, I'm using the EASYGRID default "spherical" equals "F". I'm forcing a ...

ocecept

Sorry Diego, In my previous message, I was supposed to answer the question not repeat it. I had the same problem reported by user "eastwind11" and the solution as Diego said is: fillseg(coast_new,[0.5 0.5 0.5],[0 0 0]) It will result gray in land areas and black coastline. First, make sure you have ...

ocecept

Hi all; The search engine is not working here as well. I did some tests, that maybe will help to find the problem. I used Firefox 3 and Explorer7 in the tests and the words sponge and arango: First; I tried the "Search for" button in the header: It doesn't work when I'm logged or not logged at https ...

ocecept

Hi Wilkin; Thanks for these very good tools. I'm able to use some of then, but with some options the program complain about specific functions or sub-functions. For example: Using: roms_zview('ocean_his.nc','umag',1,1,grd) I get: ??? Undefined function or method 'av2' for input arguments of type ...

ocecept

Hi Kate; Thanks for your help. I did: nm libnetcdf.a | grep NF_CLOSE And it couldn't find NF_CLOSE . than I did: nm libnetcdf.a And there's a mix of lower and upper case names (without underscores) (here) . For exemplo nf_close is in lower case but other functions are in upper case. I think this is ...

ocecept

Hi Guys; I'm running ROMS in a Opteron Machine, running SUSE Linux 10 and Absoft 64-bit Fortran 95 10.0.6. I installed the netcdf-3.6.2 in two different folders using gfortran and Absoft compilers and linked it in my Linux_gfortran.mk and Linux_f90.mk files. I was able to get ROMs working using ...

ocecept

Hi Diego; That's the m_map package. Sorry for the incomplete information. If you REALLY want your program working for a large area with my solution (and m_map), change: ((pi/180) * cos(lat) * normradius) ); by (1000*m_lldist([lon lon+1],[lat lat]))); at lines 145;163; 183 and 202 I hope the program ...

ocecept

Hi Diego; First I would like to thank you for what you have been doing here. I tried to use your roms_make_simple_grid.m to create a grid for South Australia, without success. My grid is larger (930 Km x 846 Km) than your test case and located in the south hemisphere. The problem appears to be ...

Ocean Modeling Discussion

Search found 42 matches

Re: Average Surface Flow Through The Strait of Hormuz

10th International Workshop on Modeling the Ocean (IWMO2018)

Re: AGE_MEAN

AGE_MEAN

Re: Paper using Diagnostics for tracers

Paper using Diagnostics for tracers

Re: NTSAVG and NTSDIAG

NTSAVG and NTSDIAG

Re: Problems with basic model

Re: Intel’s new i7 980x CPU gives disappointing speedup

Re: Potential Vorticity

Re: Potential Vorticity

Re: Potential Vorticity

Potential Vorticity

Re: Problem with the MATLAB R2009a version (7.8.0) and Netcdf

Re: Tidal forcing file

Re: Problem with varinfo.dat

Re: Problem with ticket 405 - Diagnostic terms review

Problem with ticket 405 - Diagnostic terms review

Re: Input tidal forcing in a large grid

Re: 2D tidal modelling and options in .h and .in files

Re: Problems with AVERAGES_DETIDE

Problems with AVERAGES_DETIDE

Re: how build tidal forcing

Re: how does roms caculate the value of zeta?

Re: Chemical Tracers

Re: NETCDF4

Re: NETCDF4

NETCDF4

Re: TS_FIXED with TS_PASSIVE defined

TS_FIXED with TS_PASSIVE defined

Re: test_vgrid.m

Re: test_vgrid.m

Re: test_vgrid.m

Re: EASYGRID: Grid generation for ROMS

Re: EASYGRID: Grid generation for ROMS

Re: Forum "search" problem...

Absoft Fortran x Linux Error