Ocean Modeling Discussion

Hi,

I'm attaching a python script here so it may be of use to others.

Good luck.

Thiago.

I have recently build a ROMS grid using RGFGRID and QUICKPLOT, from Delft3D's Graphical User Interface (GUI), and pyroms ( https://github.com/ESMG/pyroms ). The main idea is to convert a delft3d grid (.grd) created with RGFGRID into .mat using QUICKPLOT. This will give you the grid vertices in a .ma...

I did try to use TS_A4HADVECTION, TS_A4VADVECTION with UV_ADV and UV_C4ADVECTION but it leads on blowing up within several months of simulations due to a discontinuity on Temperature and Salinity. Is it need a smaller timestep to use Akima horizontal and vertical advection? Or any other things shou...

I faced a similar issue in my application. The problem occurred only in the slope region where the bathymetry was steeper. My suspicious at the time was that a strong Ekman layer was being formed, inducing a change in velocity direction. So I gradually reduced the bottom drag (if using UV_QDRAG, red...

Hi, the BLK_ZW parameter refers to height of your surface wind. Mine is set to 10.0 m as this is the height of the wind momentum components I used to build my forcing files. I used this dataset: "NCEP Climate Forecast System Reanalysis (CFSR) Selected Hourly Time-Series Products, January 1979 t...

I don't remember the exact steps I took, but the flag WIND_MINUS_CURRENT (see: https://www.myroms.org/projects/src/ticket/758 ) along with a fine tweak in viscosity (nl_visc2) led me to better results. Besides comparing to global models, I also made some comparisons with in situ data (with appropria...

Thank you very much for the reply. You are doing something wrong or something is weird with the compiler. You need to read carefully the following :arrow: trac ticket. Notice that even if SINGLE_PRECISION is activated, there are several variables that are still in double precision (dp 64-bit kind pa...

Hello all, I'm not sure this is whether a bug or a particular issue with my enviroment or setup. Anyway, I think it's worth to report. I have been running my ROMS IS4DVAR application for a while using revision 881. Now I decided to test the new flag SINGLE_PRECISION of revision 928, which is availab...

Hi, I'm running a realistic application spanning 6 years (2008 - 2014) and I was able to execute the entire simulation using ROMS version svn783. However, after upgrading to svn881, I am getting the following error: 44160 2009-01-01 00:00:00.00 2.761819E-03 1.905423E+04 1.905423E+04 6.367074E+15 (09...

Hi, I'm running ROMS (svn 836) in a realistic application in southeast of Brazil and I've noticed what appears to be an overestimation of surface currents, possibly due to wind stress. I'am using bulk fluxes and hourly fields from CFSR reanalyses as forcing and 1/12 degree resolution grid nested to ...

Update 3: I feel crazy answer my own post several times but this might be useful for others. It is possible to use the 3rd-order upstream split tracer advection scheme (TS_U3ADV_SPLIT) for the nonlinear model with IS4DVAR, or any other option, as explained in here https://www.myroms.org/forum/viewto...

Update 2:
After compiling and trying to run, I found that TS_U3ADV_SPLIT is indeed not supported in adjoint as it is not coded yet.

Update: I managed to compile by making a minor modification in "ad_t3dmix4_geo.h" from ./ROMS/Adjoint/. I included "USE mod_ncparam" statement in subroutines "ad_t3dmix4" and "ad_t3dmix4_tile" definition so "isTvar" becomes accessible. Cheers. Thiago.

I'm trying to compile an IS4DVAR application using TS_U3ADV_SPLIT as tracer advection scheme in NLM but it throws me the following error: ad_t3dmix.f90(268): error #6404: This name does not have a type, and must have an explicit type. [ISTVAR] IF (ad_LBC(iwest,isTvar(itrc),ng)%closed) THEN It seems ...