Search found 33 matches

by dhson
Mon Apr 15, 2013 2:39 am
Forum: ROMS Sediment
Topic: Sediment from river!
Replies: 19
Views: 28523

Re: Sediment from river!

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...
by dhson
Tue Jan 15, 2013 9:21 am
Forum: ROMS Sediment
Topic: Sediment from river!
Replies: 19
Views: 28523

Sediment from river!

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
by dhson
Fri Apr 02, 2010 9:16 am
Forum: ROMS Installation
Topic: MCT 2.6 issue
Replies: 3
Views: 3036

Re: MCT 2.6 issue

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
by dhson
Sat Feb 27, 2010 4:52 am
Forum: ROMS Installation
Topic: ROMS whit SWAN problem
Replies: 12
Views: 11445

Re: ROMS whit SWAN problem

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 ...
by dhson
Tue Feb 23, 2010 4:59 pm
Forum: ROMS Discussion
Topic: SWAN & ROMS coupling problem
Replies: 1
Views: 2301

Re: SWAN & ROMS coupling problem

Sorry, thanks for any suggestion,
by dhson
Tue Feb 23, 2010 4:58 pm
Forum: ROMS Discussion
Topic: SWAN & ROMS coupling problem
Replies: 1
Views: 2301

SWAN & ROMS coupling problem

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 w...
by dhson
Mon Feb 15, 2010 3:12 pm
Forum: ROMS Sediment
Topic: SEDIMENT in 2D mode
Replies: 1
Views: 2606

SEDIMENT in 2D mode

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
by dhson
Mon Feb 01, 2010 5:07 am
Forum: ROMS Problems
Topic: How to set Hmin (<0m), Tcline, Hc in WET_DRY
Replies: 3
Views: 3216

Re: How to set Hmin (<0m), Tcline, Hc in WET_DRY

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....".
by dhson
Sat Jan 30, 2010 4:13 pm
Forum: ROMS Problems
Topic: How to set Hmin (<0m), Tcline, Hc in WET_DRY
Replies: 3
Views: 3216

How to set Hmin (<0m), Tcline, Hc in WET_DRY

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 va...
by dhson
Wed Sep 16, 2009 8:36 am
Forum: ROMS Discussion
Topic: SET_NGFLD - current model time exceeds ending value for var
Replies: 5
Views: 3136

Re: SET_NGFLD - current model time exceeds ending value for var

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 tide...
by dhson
Tue Sep 15, 2009 1:51 am
Forum: ROMS Discussion
Topic: SET_NGFLD - current model time exceeds ending value for var
Replies: 5
Views: 3136

Re: SET_NGFLD - current model time exceeds ending value for var

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 f...
by dhson
Sun Sep 13, 2009 3:04 pm
Forum: ROMS Discussion
Topic: SET_NGFLD - current model time exceeds ending value for var
Replies: 5
Views: 3136

SET_NGFLD - current model time exceeds ending value for var

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 ...
by dhson
Tue Apr 14, 2009 2:22 pm
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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 tho...
by dhson
Sun Apr 12, 2009 3:05 am
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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 inlan...
by dhson
Sat Apr 11, 2009 6:52 am
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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
by dhson
Fri Apr 10, 2009 4:00 pm
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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 ...
by dhson
Wed Apr 08, 2009 2:30 am
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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 v...
by dhson
Tue Apr 07, 2009 11:48 am
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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.
by dhson
Mon Apr 06, 2009 3:29 pm
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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_ENR...
by dhson
Mon Apr 06, 2009 1:43 am
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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 w...
by dhson
Sun Apr 05, 2009 2:42 pm
Forum: ROMS Problems
Topic: potential mass conservation problem on dry cells
Replies: 20
Views: 9354

Re: potential mass conservation problem on dry cells

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...
by dhson
Tue Aug 05, 2008 2:31 am
Forum: ROMS Discussion
Topic: High output amplitutes in tidal simulation
Replies: 4
Views: 2731

Re: High output amplitutes in tidal simulation

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 = "H...
by dhson
Mon Aug 04, 2008 2:52 am
Forum: ROMS Discussion
Topic: High output amplitutes in tidal simulation
Replies: 4
Views: 2731

Re: High output amplitutes in tidal simulation

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 varia...
by dhson
Thu Jul 31, 2008 9:34 am
Forum: ROMS Discussion
Topic: High output amplitutes in tidal simulation
Replies: 4
Views: 2731

High output amplitutes in tidal simulation

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 #...
by dhson
Wed Jul 09, 2008 2:33 am
Forum: ROMS FAQ
Topic: Tidal forcing at the open boundary
Replies: 23
Views: 62483

Re: Tidal forcing at the open boundary

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 ...
by dhson
Tue Mar 25, 2008 4:53 am
Forum: ROMS Bugs
Topic: SWAN_COUPLING: Trunk-r165 vs Trunk-r168
Replies: 2
Views: 5917

SWAN_COUPLING: Trunk-r165 vs Trunk-r168

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 f...
by dhson
Thu Jul 05, 2007 1:41 am
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 21632

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.
by dhson
Wed Jul 04, 2007 1:37 am
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 21632

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 sedime...
by dhson
Tue Jul 03, 2007 10:00 am
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 21632

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

Do you have any suggestion?
by dhson
Tue Jul 03, 2007 7:40 am
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 21632

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 ...
by dhson
Tue Jul 03, 2007 2:04 am
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 21632

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 d...
by dhson
Mon Jul 02, 2007 5:10 am
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 21632

Minus values of sediment next to sediment point sources

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)...
by dhson
Tue Jul 06, 2004 7:32 am
Forum: ROMS Problems
Topic: I could not run ROMS2.1 with SEDIMENT defined
Replies: 0
Views: 2450

I could not run ROMS2.1 with SEDIMENT defined

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