Ocean Modeling Discussion

Thank you all for your replies and help. AVERAGES_DETIDE or using python seems to solve the problem. With regard to reading all data into memory, I think even in the MATLAB tidal analysis we can avoid that. I mean, we can extract the time series at a single node and store it in a dummy variable and ...

Thanks for the reply. As I mentioned, I was wondering about something more convenient than extracting the data from a netcdf file and then doing the tidal analysis both in FORTRAN.

Hello all, I have been using ROMS for (barotropic) tidal simulations. For high resolution tidal models, big data files should be first transferred from the computer cluster (which does not have MATLAB due to licence issue) to PC, and then locally analysed. This task is usually time consuming and ...

thanks John,
that was very helpful.

Hello All Referring to this page, https://www.myroms.org/wiki/index.php/SSW_BBL I think, there might be a mistake in the equation for computing the wave orbital velocities ub=Hs/(2sinh(kh)), at least the dimensions are not consistent. The relationship in linear wave theory [1] is ub=pi*H/(Tsinh(kh ...

Hi Kate
thanks for your reply. Is there any test case with spatially variable friction coefficient which implements this?

In COWAST model, the contribution of waves on the bottom stress can be incorporated by various methods such as SSW_BBL (Madsen, 1994). Additionally, the dissipation of the wave energy in the BBL results in a bottom stress which can be included by various switches for bottom streaming such as BOTTOM ...

Many thanks John That was very helpful. Therefore, wind induced stresses are actually wave induced stresses in varinfo.dat. I suppose if somebody wanted to compute wind induced stresses, he could run ROMS (decoupled ignoring waves) with and without wind forcing and subtract the resulting current ...

Hello all I have run a coupled ROMS-SWAN model and I used MB_BBL for bottom boundary layer formulation. I would like to extract the following bed stresses: > Current induced > Combined Current induced & wave induced ( tauCW=tauC*[1+1.2*(tauW/(tauC+tauW))^3.2] ) > Wind Induced > Total resulting ...

That was really helpful John. I compiled the COAWST version of the headland test case and it works fine without any issues.
Thank you very much

Reza

Hello All Has anybody run the headland test case successfully using the latest version of ROMS? I have compiled ROMS based on this test case and it blows up after about 80 time steps. I guess, the boundaries are not forced properly in my case. In contrast the inlet test case works fine. I was also ...

Further to my previous post, I realized that the header file for this test case in the Include directory of ROMS source code is different from the header file which comes in the ROMs example directory (whic includes all the necessary files for running this case). In terms of the ROMS version, I am ...

Many thanks John. That make sense. I have done this for other models as you say.
That solves the problem

Reza

Thanks. I added two lines in configure.wrf to define the MCT include and library files the compilation now finishes and makes coawstM file. Unless i do clean and reconfigure wrf, it is not overwritten in my case. In roms, I have defined these MCT paths in compiler fortran setting file. Do you mean ...

Hello all I have compiled both inlet test case and head land test case in ROMS to start studying wave current interactions. Inlet test case works fine. But, the head land test case blows up after 80 time steps due to nan values in the results. I have used the build.bash and test_head.h to build the ...

Thank you John. I am very sorry for the late reply, I was expecting to get an email if somebody replied to my post. I should have checked the forum instead. Anyway, I looked at the log file of the build and although I have put the right path for the MCT directories in ROMS compilation file, WRF ...

Hello all I am interested in turbulence properties of my model. I can see the following parameters in my history file by choosing Hout(idMtke) == T ! Write out turbulent kinetic energy. Hout(idMtls) == T ! Write out turbulent kinetic energy times length scale. However, when I tried to output other ...

Dear All Could you tell me about how can I prepare input data files for Caspian Sea? Caspian Sea is a great lakes but most of datasets are related to the global models and this area is blank. If I were you, I would look at the previous publications regarding the Caspian sea. They provide a Data ...

Hello I have been able to compile the COAWST code for UPWELLING case (Just ROMS) and then inlet test case which involves SWAN and ROMS coupling without any problem problem. Using the same buildcoawst.sh file, I started to compile the code for another applications which involves ROMS-SWAN and WRF ...

Hello everybody I am working on a coupled wave-tide model for Northwest European Shelf Seas. In the latest version of the ROMS that I have downloaded via svn, the version of SWAN in the ROMS\trunk\Waves folder is 40.51 while the latest version of the original SWAN code is 40.91. I was wondering, is ...