Ocean Modeling Discussion

If you have made changes to variables you can use the inbuilt function of matlab (nwrite) itself to store them in the existing variable of the netcdf file.

The function is c_grid.m. The function takes the following as arguments on input: % % Lp Number of RHO-points in the XI-direction % Mp Number of RHO-points in the ETA-direction % Gname Grid NetCDF file name (string) % NewFile Switch to create a new file (logical, OPTIONAL) % spherical Spherical grid ...

ROMS tools (Matlab based) have a romsgui function, an interactive visual toolbox. ( The link http://forge.ipsl.jussieu.fr/roms_locea ... order=name). Pl Check folder visualisation tools

The forcing file which is required is climatological (roms_frc_woa365.nc), which could be created easily using the Romstools toolbox (specifically for your grid) using
climatological World Ocean Atlas data (WOA2013).

Yes, adding the same in the header file and recompiling it.

Have you tried adding salinity correction term to the CPPDEFS key ?

The reason for this is probably the location of the file directory you provided, ROMS could not locate the restart file location. Pl check if the location you provided in .in file is correct or not.

Checking your log file it has been found that there exist NaN values in your U wind and V wind component (1.0000E+35) and also in many other variables of your forcing file. So you have to mask those bad values before running the model by setting the NaN flag in your forcing file. The forcing file is ...

I think you could use the roms_2d slice.m function in the same repository or use secondary software like Ferret (sample_xyt) to extract the data from the post-processed ROMS output data.

You would have to set the switch for intermediate file (OA ) ON (ie to 1) in the parameter file(romstools_param.m). And try re-running the script again.

Dear Kobashi, You may please note that ROMS has the capability to perform the spatial and temporal interpolation of 2D forcing fields. So it is not necessary that you need to interpolate your forcing files to your grid. But make sure that forcing file is either larger in dimension or smaller ( I ...

Dear Kobashi, I also had the same problem when I used same forcing files for nested grids(exception tidal forcing), the same was solved by repeating the forcing names as follows The following is an excerpt from my .in file NFFILES == 9 9 ! number of unique forcing files FRCNAME == input/ncmrwf/Input ...

The following is the syllabus of Department of Physical Oceanography CUSAT, splitted into theory and practical sections. Please find the attached figures giving details about the theory syllabus MODELING1.png The syllabus for practicals section is given by the following attachment Practical_modeling ...

Have you created a new contact grid file after you changed the domain ? As domain changes the contact point also changes.

Dear osua72
I have also met up with such instable values in salinity during my two way nesting and for one way it is almost fine. Would it be possible for you to check the configuration with one way nesting defined instead of two way?

Dear osua72
Was you able to resolve the blowup? Could you please state out what was wrong with the setup.

I have tried using a simple nesting set up (both one way and two way), when I run the two way nesting setup negative values arise and further it blows up[ cfl criterion is still valid] but for one way no such negative values arose. I would like to know the possible reasons for the same. The both ...

I have 2 grids which are running independently one of a larger region with coarse resolution and other within the larger region and smaller in domain and of higher resolution. The grids are needed to be nested, to create contact point file I used contact.m but error of non existence of contact point ...

My standalone ROMS set up for coarse(1/12) and fine(1/36) resolution is running well individually but when I use two way nesting the Cw value becomes greater than one for finer resolution grid. Is that problem due to the higher time step? I have used 300,100 for coarse and finer resolution ...

Hi kate, I need your help to resolve one of my doubt regarding Nesting, I have a coarser grid (1/16) and a refined grid(1/48) the refined grid is created from the coarse using coarse2refine.m matlab script. Rather than having a higher resolution bathymetry for the refined grid, the same bathymetry ...

Thank you for the reply Kate. I tried with multiple combinations of Ntilei and NtileJ keeping their ratio almost equal to Lm/Mm, as I reduced the total number of processors from 40x30 to 10x7 (40/30=1.33~Lm/Mm)(10/7=1.423~Lm/Mm) I was able to decrease the time taken for running the model ...

I have started to run a two way nesting simulation with a grid size of 1079*755 755*663, for running 2 minutes the time taken is 3 minutes approximately(slower as compared to standalone model), I am attaching my .in file for your reference. Kindly suggest a solution for this one. I have seen ...

Yes you are right, it was due to an error in the input script, I corrected it. Thank you.

Thank you for your suggestions, I have figured out that the problem was due to some syntax error in my ocean.in file, and have solved it. I have tried to use export use_debug but the same results in creation of oceanG (shared executable) where as I need oceanM. Is it possible to get oceanM (Mpi ...

I am trying to run a refined nesting case for 2 grids, after compiling with define NESTINg I am getting the following error report: Model Input Parameters: ROMS/TOMS version 3.7 Monday - June 3, 2019 - 4:32:54 PM ...

I suppose that you are trying to compile the test case of upwelling, i too had similiar problems while compiling the test case. Please try to rectify it by changing line 52 of Linux-gfortran.mk as follows ( code from unidata): Before LIBS := -L$(NETCDF_LIBDIR) -lnetcdf After LIBS := -L$(NETCDF ...

Thank you Andres

No I have just the global COADS atlas files for each month in the folder. Could you please specify where to get the files mentioned by you.

I would like to get data for running the script make_forcing.m whenever I run make_forcing.m in matlab session I am getting this as my output: ASHSW Model Read in the grid... Create the forcing file... Warning: NARGCHK will be removed in a future release. Use NARGINCHK or NARGOUTCHK instead. > In ...

I tried using make_grid.m of ROMS_tools to create grid, but I did not get any desired result. Kindly resolve the problem.