Ocean Modeling Discussion

Hi all, Is there anyone here ever adjust the ROMS CODE to change the TNUDG / M3NUDG for different side (North, South, East, West). Is it will affect the volume conservation in the model domain? The case is, in my model domain the North Equatorial Counter Current cant freely radiate out during the ...

Thanks for your suggestion, Kate. I tried to increase the frequency of atmospheric fields and it helps. Now, I am setting up to increase the frequency of boundary data.

Dear all, I wonder if anyone here ever experienced the unmatch "lateral boundary tracer (T and S) data at the surface layer" to the "calculated data which is affected by atmospheric forcing", and how to solve this problem? I found this in my model, higher surface layer temperature than available ...

If you turn on climatology nudging, it will change the OBC time scales. If you don't want that, change the open boundary code to do what you want it to do - that's what I did in my branch. Dear Kate, I read this from https://www.myroms.org/projects/src/ticket/627. I want to ask: 1. How to set the ...

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

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

Thank you so much for your suggestion. I'll try reduce bottom drag and change the tracer/momentum advection scheme before nudging to climatology. Thanks a lot. I faced a similar issue in my application. The problem occurred only in the slope region where the bathymetry was steeper. My suspicious at ...

Dear all, I have a problem on my vertical profile of velocity in almost all area in my model domain. The vertical profile of meridional velocity has a counter current below the thermocline layer. I can sure the initial condition of T and S is quietly good. In the initial simulation, the velocity ...

Thank you Kate.

Dear all, I want to know, if is it a must activate the VolCons if I am using RadNud boundary condition for 3D U,V,T,and S? I have 4 open boundary conditions, with boundary value from reanalysis ocean model. And I am forcing the model with surface heat flux, freshwater flux (calculated using bulk ...

Where are you getting your net longwave heat flux? I've had issues with unrealistic temperatures in shallow regions when applying heat fluxes. The problem is that the heat flux depends on the temperature difference between the air and the ocean, but when we take heat fluxes from some other model ...

Ok Kate, I will try to turn off the BODYFORCE. I used surface heat flux and freshwater flux in the surface boundary layer.

I did bathymetry smoothing using B.C Saphiro constant available at GridBuilder. Sorry Kate, I think I made a wrong copy. I used this CPP #define UV_ADV #define UV_COR #define UV_QDRAG #define DJ_GRADPS #define SPLINES_VDIFF #define SPLINES_VVISC #define TS_U3HADVECTION #define TS_C4VADVECTION ...

Dear All I have a problem with my model. The bottom temperature is intruding to the shallow water area very rapidly (the bathymetry range is 60 - 5000m, with 30 vertical levels). For the longer simulation, the cooler temperature is spreading towards the surface . I activated the following CPP. 1 ...

Your values look like you should use 'second-1' as the units. 'day-1' is also a valid option. See get_nudgcoef.F to check that. You say you want 1 day at the edge and one year in the middle. I think if it rather as strong (1 dayish) at the edge to infinite time (zero coefficient) over some width ...

Yes, you have to have a boundary file as well as a climatology file for the options you have picked. The LBC choices demand boundary values from a boundary file or from ana_xxx.h files. For the same number of records, the boundary files should be much smaller than the climatology file, but in ...

Thank you very much for your useful comments. I found above problem was caused by "mask" itself. when I created "grid file" with new version of Matlab (2014a) that problem is gone. I still do not know why it happen because when I check mask value of both "creating by old version of Matlab, 2013a ...

Thank you for your replies. Isn't it has similar data format? Similar, yes, but the script is expecting certain names and conventions for the coordinates and dimensions. Code doesn't do similar . By the way, I interpolated the WOA data (1 x 1 degree) to my grid (5 x 5 km), which might be the ...

Looks to me like you're using a World Ocean Atlas climatology as input. Tfile = 'woa18_ann_temp.nc'; Sfile = 'woa18_ann_salt.nc'; This script is designed to process Mercator-Océan data from Copernicus.eu. Thank you for your replies. Isn't it has similar data format? By the way, I interpolated the ...

% % D_MERCATOR2ROMS: Driver script to create a ROMS initial conditions % % This a user modifiable script that can be used to prepare ROMS % initial conditions NetCDF file from Mercator dataset. It sets-up % all the necessary parameters and variables. USERS can use this % as a prototype for their ...

One of those is easy - change the timestep and see if the time of the trouble changes. It's been a while since I've had an unambiguous blowup from steep topography. Steep topography causes unphysical currents, but not a blowup unless you have done no smoothing at all (there are those in here who ...

Thank you Kosa, I have solved my problem.

Thank you Kate, I have solved my problem.

Build_romsG/propagator.o /public/home/corry/model/project/upwelling/Build_romsG/master.o /public/home/corry/model/project/upwelling/Build_romsG/roms_export.o /public/home/corry/model/project/upwelling/Build_romsG/coupler.o /public/home/corry/model/project/upwelling/Build_romsG/esmf_roms.o /public ...

ROMS/TOMS - Output NetCDF summary for Grid 01: Analytical header files used: ROMS/Functionals/ana_btflux.h ROMS/Functionals/ana_diag.h ROMS/Functionals/ana_grid.h ROMS/Functionals/ana_initial.h ROMS/Functionals/ana_smflux.h ROMS/Functionals/ana_stflux.h Found Error: 02 Line: 350 Source: ROMS/Utility ...