Ocean Modeling Discussion

Thank, I've got another problem that ROMS didn't read boundary sediment concentration if a number of classes > 2. For example, I've put 2 classes of non-cohensive mud in eastern and western boundaries both 0.1 mg/l. There are reasonable value for MUD_01 but nothing (zero) for MUD_02 in his, sta, rst ...

I used to run ROMS with RIVER_SEDIMENT (of course with UV_PSOURSE and TS_PSOURCE). But in present version (revision 643), RIVER_SEDIMENT seems to be no longer valid. Consequently, I can run ROMS for S and T but not sediments (mud, sand) from rivers.

Any ideas, Thanks

I've successfully install both MCT 2.2 and 2.6 by gfortran(/pgi/ifort) on UBUNTU x64. But I've just set

FC = gfortran (/pgf90/ifort).

the makefile will automatically find mpif90

I've got similar problem: + On my laptop dual 2 core (ubuntu 9.10, x64, MPICH2 and every related things made from gfortran 4.3,...) my problem described in https://www.myroms.org/forum/viewtopic.php?f=14&t=1743 + Yesterday, the above problem was disappeared and I've got fantastic results when I ...

Sorry, thanks for any suggestion,

I'm trying to run ROMS with SWAN_COUPLING, WET_DRY, UV_PSOURCE, TS_PSOURCE,SEDIMENT and SSW_BBL. I believed they're all toughest modules. I've spent almost two months on playing with ROMS but not yet been successful: - Without SWAN_COUPLING, I can get reasonable results - With SWAN_COUPLING, there ...

Hi,
I'm running ROMS with SWAN, SEDIMENT in 3D mode but it's quite time consuming (350x500 nodes, dt=6s).
Since my domain is shallow, I'm wonder whether the above combination can be run in 2D mode?
I supposed it can except SEDIMENT module.

Thanks

Thank for your quick reply. However I've set Tcline=hc=-2 m everywhere: *.ini, *.bry, *.in but ROMS took the value 0 somewhere and said 0 and -2 are "inconsistent value of variable....".

Dear all, In wet_dry simulation, can I set Hmin to negative (cells above MSL, say -2m) to see how the beach inundated by the tide of 3m high. If I can can do so, then, how to set Hc and Tcline. I've tried with many alternatives of Hc and Tcline without success BUT: NETCDF_CHECK_VAR - inconsistent ...

Thank you, Kate, for your hint. It's solved. I've made a mistake when generate netCDF file without taking care of the difference of the offsets of River.time between ROMS and Matlab Codes. By the way, I'm having some problem in reproduce inundation maps due to the combination of storm surge and ...

Thank you, Kate I've used the same configuration (header file) for both. But as fas as I remember, I've got this error before but I did not post it here and I've modified the source code instead (just to pass this error). Now, it comes back after upgrade. I didn't see any special thing in my input ...

I've just update ROMS to SVN Revision : 396M and got SET_NGFLD - current model time exceeds ending value for variable: river_transport TDAYS = 223.1389 Data Tmin = 2448638.0000 Data Tmax = 2448968.0000 Data Tstr = 2448968.0000 Data Tend = 2448968.0000 TINTRP1 = 0.0000 TINTRP2 = -100393198.0000 FAC1 ...

Hi! I didn't make ROMS run longer 1day and 15 hours. Our tidal amplitude is <1m. That means (if you said correctly "the water level to raise above the land" ) the water can come only on boundary cells with h>-1m. But my figure indicated that the water come from all boundary cells. Even though, as ...

Well, I didn't make any specific condition along all boundaries. They all start from initial condition (zeta=0) then driven by tidal constants. If BC at wet_dry cells was coded correctly then, Probably, I've got to check tidal constants that were interpolated form TPXO7.1. I'm not sure how the ...

wow! It works!!!
But I believed that there're still some bugs related to BC (when the water level is physically below bottom of the boundary cells).
The attached are output of HISTORY after 1 hour (in fact ROMS now can lives much more than that).

Thanks

Hi! It's me again, Since my domain was a bit big (200x600) so I was afraid of some uncontrolled grid cells. It took me sometime to prepared more ideal bathymetry as attached with the same size (so in thi run only netcdf grid file was different). ROMS died at the 2nd timestep (a step better than the ...

Thank you for the response, RESTART file did not follow any rule with differences of dt and all blow-up in the first step: - with dt of several second: patchy pattern of NaN and real value (*.E-6) without specific relation with topography; - with smaller dt (0.001 seconds): no NaN observed but one ...

I think probably compare Hmin to Tcline and hc is not a good idea since Hmin is elevation and the later are thickness.
I've commented out both Tcline and hc in mod_netcdf.F:

! exit_flag=5
! EXIT

but no improvement in ROMS.

Thanks for the information, but not much help. Using Vtransform == 1 ! transformation equation Vstretching == 1 ! stretching function with the combination hmin=-5m (as explained above), and (hc=-5m, Tcline=0m) or (hc=-5m, Tcline=-5m) I can run but: STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL ...

Sorry, with hc=.1m, hmin=-5m, I've got (SVN Revision : 342M): NETCDF_CHECK_VAR - inconsistent value of variable: hc 1.000000E-01 -5.000000E+00 in file: ROMS_FILES/roms_ini.nc with hc=0, I've got the same results. If I use hmin=.1m in SEAGRID, I can run ROMS but that's not what I want since SEAGRID ...

I've tried to use WET_DRY for a small island that its top is at about 2 m above MSL. To reserve its original topography, I've got to use hmin = -2 m (in SEAGRID). Start computing with zeta =0 and I 've expected that most of cells will be "WET" in flood tide. But I couldn't compute with any value of ...

Dear Jonh, In fact, update ROMS version 3.1 from trunks did not take much time to modify so we prefer this way. Unfortunately, I did not find any problem in the tidal forcing file: double tide_period(tide_period) ; tide_period:long_name = "Tide angular period" ; tide_period:units = "Hours" ; double ...

Thanks for you reply I'm using ROMS UCLA (trunk-r186, trunk-r195, trunk-r200, etc and they are all the same problem). I haven't modified any thing in varinfo.dat yet. Here is the related variable: ! !------------------------------------------------------------------------------ ! Tidal forcing ...

Dear all, Using ROMS-3.1 for tidal simulation in South China Sea, the output tidal amplitudes in all stations seemed to much higher than expected (says, ~ 4-5m compared to ~1.5-2m). The following CPPs were used: #define UV_ADV #define UV_COR #define UV_VIS2 #define MIX_S_UV #define TS_U3HADVECTION ...

Hi, I'm trying to simulate the variation of the tidal constants (and consequently, local Sea Level Rise) when MSL increases e.g. 1.0 meter. Using ROMS (now, trunk-r186) for tides at South China Sea with below CPP for all OBCs, I've got surface elevation ussually 50% higher than observation for most ...

Hi, 1. I'm trying to couple SWAN & ROMS on Fedora Core 7 with Intel Fortran 10.1.012, NETCDF4, MPICH2. I've attempted with the example of coupling_inlet_test.in however: * Trunk-r168 was failed in compilation: cd Build; /opt/mpich2/bin/mpif90 -c -ip -O3 -xW -I/opt/MCT/include ocean_coupler.f90 ...

By the way, did "# define SED_MORPH" was well tested by someone?
I've tried to use it for the river but ROMS 3.0 blows up after long enough time.
It's mostly seemed that at some corner nodes, the bottom layer was over eroded or sedimented.

Daer all, Before running wit SEDIMENT, I've deactivated SEDIMENT: ROMS worked fine! It seemed to be fine also with T & S (at least, no ridiculous values). You might notice that variable Dsrc, direction of point source discharge, was set to 1 in my case (which was not quite appropriate for sediment ...

I've defined it as:
#define TS_MPDATA
(details in the beginingmessage)

Do you have any suggestion?

Sorry for the bad explanation, What I' tried to say are: 2) 1325 value shows the sediment concentration as kg/m3 of Qbar (~5m3/s) at the point source. 3) if I put point source at (i=7, j=5) then i've received the horizontal distribution of sediment concentration as the below matrix: i= .....4 5 6 7 ...

Thank for the response 1) Yes, I have "is" and "k" loop 2) I mean the negative value next to (Isrc,Jsrc) and the 1325 value shows the sediment concentration as kg/m3 of Qbar (~5m3/s) In fact, in this simulation I've tried to figure out how the sea bottom is elevated by sediment dumping. quick look ...

Dear all, Using ROMS 3.0 to simulate transportation of sediment dumped from the ship (point source). There's always minus value of sediment next to the point source. In ANA_PSOURCE.H, I have modified Qshape, Qbar, etc. and: ...................................................... DO is=1,Nsrc Dsrc(is ...

Dear all, I could not run ROMS-2.1 with SEDIMENT option (with BIO_FASHAM defined). I believed there is a problem with this command in ANALYTICAL.F CALL exchange_r3d_tile (ng, Istr, Iend, Jstr, Jend, & & LBi, UBi, LBj, UBj, 1, MBEDP, & & bed) because bed is 4d array. I' ve corrected it to: CALL ...