Ocean Modeling Discussion

****SOLVED**** In the coawst.bash I was simply pointing to another directory and *.h file. ooooops! ************************** Dear John, I already hav all of those set to 1, but still I have the problem. By the way I'm using the Sandy *.in as "template" so maybe there is some other CPP op...

Dears, which CPP option (and/or *.in variables) requires the NGCNAME variable to actually point a file? I mean, I put in my *.h CPP option "#undef nesting" BUT the program still looks for a "contact" *.nc file. I'm trying to setup an application by modifying a nested domains appl...

Sorry for the late reply, fix all my errors took me a long time ! I did as you advised me, I turned on USE_DEBUG and I finally found where my errors came from, thank you ! If someone has the same problem I can advice first to update the varinfo.dat file and then, my main error was in my forcing fil...

Update: I was able to reduce such value of error in BICONJ by increasing number of iterations of Nbico in s4dvar.in but it seems not affect the result, so maybe it is not the responsible. I have to explore in some other direction. What is strange to me is that the nonlinear is stable, while the firs...

of Hi everybody, When trying to run I4DVAR for an area of Western Med (with an Island in the middle of the domain),assimilating just ADT-AVISO (corrected in respect of the bias with model by subtacting and adding respectively their averages and model averages for zeta), I have the following symptoms...

Dear Thanos,
I'm also interested on this point.
Do you have news on that?

I think I will try this excercise ASAP, when I will find some time to do it-
Cheers
Antonio

For salinity and temp it could be that you are imposing every day heavy "winter" conditions, i.e. cold winds cooling the sea surface and tends to concentrate salt? The effect could cumulate in time as it is "fixed", despite you relax to climatology. Zeta values however suggest th...

For example, if this is netCDF file including two days data, time array is 0 12 24 36 48 3 6 9 15 18 21 27 30 33 39 42 45. Thank you very much to hbzong. I was struggling from a while with internal times of ERA-interim netcdf's extracted from the data portal and now finally I have understood. For a...

The COARE algorithm is clear that the wind vector needs to be referenced to the surface current. This is straightforward to do in bulk_fluxes because the wind is already expressed in local roms coordinate directions. Schematically, we modify uwind,vwind: Uwind(i,j) = Uwind(i,j) - 0.5*(u(i,j,N)+u(i+...

Dear All.
Any news about that issue?

First of all, sorry for resurrecting such an old 3d :roll: I'm interested to this issue about combined use of RadNud and Bulk fluxes. I understood that there could be a (complex) problem about heat gain or loss at boundaries by using such a combination. I have some question about that, hoping that s...

I found another typo in the d_ecmwf2roms.m script, here: F( 6).Vname = {'Uwind', 'v10v'}; and F( 7).Vname = {'Vwind', 'v10u'}; U and V are inverted. Further v10u and v10v are no (longer?) the exact names of variables as downloaded from ERA website in netcdf files, while the names are u10 and v10. So...

First of all, sorry for catch again this very old post. In any case: this is exactly the error I have just now. I'm just trying to port my test-case on a new pc, just copied *.in, *.h and build.sh and installed new version of roms from svn. Upwelling compile and execute correctly, but when I try to ...

In the matlab/netcdf routines just downloaded from svn, in nc_read.m there are lines 429-430 (about) like that case {nc_int, nc_short, nc_byte} f = double(f).*scale+offset; wich give me an error about undefined variable nc_int when launching from d_ecmwf2roms.m I fixed just commenting them, but mayb...

OK, now it's fixed. There was a disagreement between "model time" (the ocean_time in the IC) and the time in the BC. Further, there was some problem with the old BC netcdf file (I don't know where extactly): I rewrote BC on a simpler manner (less attributes, variables ...etc...less stuff) ...

Hi kate, records are not reversed. I had days but now I turn all the stuff in seconds. I think I found a first error: I just realized that my "ocean_time" in initial conditions was wrong of 588 days (!!) and didn't correspond to the first record of *_time variables in bry.nc. I fixed this ...

Hi everybody, I'm newby to roms. I'm trying to provide IC and BC from a coarser model to my own implementation, with ana surface forcings. While initial are passed correctly (I think), when bounds are checked the following error occurrs: SET_NGFLD - current model time exceeds ending value for variab...

Hi everybody, I have a doubt on the editmask tool. I generated my grid by using seagrid but initially setting all the points as seapoints. Then I generated a land mask for rho_points just using bathymetry as filter (not in seagrid). After that, to refine the mask, I used the editmask tool, that work...

Thank you very much kate,
I thought for mistake the heat fluxes were defined by default.
Cheers,
Antonio

Hello everybody, As you will argue I'm new with ROMS.... I wonder why in the windbasin test case, which uses analytical forcings and closed boundaries, the temperature does not change, even if ana_smflux is defined and some wind-stress is applied. Temp is initialized constant allover the basin to 20...

I am afraid I do not follow you; If you for example copy your original roms grid to say new_grid.nc (in order to have original to compare with) and than smooth variable h and put it back inside that new file. Than in ocean.in you have to change GRDNAME to point to the new_grid.nc; run roms and chec...

Dear Ivica and all. I have some trouble with my rx0. I smoothed the grid as initially I've got 0.9 of rx0. It now fall to 0.3. But, without changing grid or vertical levels, but only the configuration (ana forcing, boundaries etc) I saw my rx0 coming back to somthing like 0.9. Is this possible? What...

Sorry,
I've found the answer in this viewtopic.php?f=14&t=2033 post.
Cheers

Hello everybody,
I have I think a dummy question but...here it is:
I read through the post that the h (and maybe even hraw) should not have 0 values. But what about the LAND points? Even the Land should not have zeros or what?

Thanks in advance,
Cheers,
Antonio

There are several idealized and realistic cases to download and play with it. They are in the test repository! This information is in :arrow: WikiROMS . This will keep you busy for very long time. There is a lot to learn. We will add the LATTE test case to that repository when we get some free time...

Hello, I'm a new user too and I'd like too to run the latte_c realistic testcase. I ran successfully upwelling test case but latte it seems to have also compiling problems due to some issue on old *.h files and other stuff I don't know. Here is the print of the last lines of compilation error (compi...