Ocean Modeling Discussion

CBian

When I set NAVG ==1 in the .in file. The u,v results in output file ocean_avg.nc are not right. The u, v will accumulate with time-step.

CBian

Dear ROMS developers, I set the sediment parameters in the sediment.in file: ! Particle settling velocity (mm/s). SAND_WSED == 0.06 0.60 6.00 0.06 0.60 6.00 12.0 24.4 However, in the output history file, it said: double Wsed(NST) ; Wsed:long_name = "sediment particle settling velocity" ; W...

CBian

Dear ROMS users, I put cohesive mud and non-cohesive sand into rivers, however the log file output show : GET_NGFLD - river runoff suspended sediment concentratit = 345 00:00:00 (Rec=0012, Index=2, File: byecs-river.nc) (Tmin= 15.0000 Tmax= 345.0000) (Min = 0.00000000E+00 Max = 0.00000000E+00) ........

CBian

Dear ROMS users, I used ROMS model to simulate the Yangtze River sediment transport in East China Sea. I confront a problem these days: All the sediment particles in the model are non-cohesive, and most of the river discharged fine-grained particles (sand_size = 0.004 mm, 0.008mm) suspended in the w...

CBian

I think you have a bad continutation of the line in your input file. It is probably repeating the last value on the first line and not seeing the remaining entries. Did you edit the file in Windows? If so try dos2unix on your sediment input file. In any case, 42 sediment classes is a lot. It actual...

CBian

Recently, I want to add 42 sediment classes into model. In the ocean.in file, I set: NNS = 42 ! Number of non-cohesive (sand) sediment tracers In the sediment.in file I set: SAND_SD50 == 0.004d0 0.008d0 0.016d0 0.032d0 0.064d0 0.128d0 0.256d0 0.004d0 0.008d0 0.016d0 0.032d0 0.064d0 0.128d0 0.256d0 0...

CBian

[quote="CBian"]Dear ROMS users, *** 194 Could not allocate in ALLOCATE PROG=mod_ocean.allocate_ocean ELN=94(400500c88) Called from mod_arrays ELN=63(400451ba4) Called from ocean_control_mod.roms_initialize ELN=101(400001ac8) Called from ocean ELN=97(400000e68) [quote] The NEC computer manu...

CBian

When I #undef SEDIMENT, the model can run one step, after that there is some message: STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME C => (i,j,k) Cu Cv Cw Max Speed 0 0 00:00:00 2.361273E-14 2.194795E+04 2.194795E+04 5.944509E+15 (360,250,19) 1.367113E-06 0.000000E+00 0.000000E+00 3...

CBian

Dear ROMS users, I can run my model case successfully in a little cluster using gfortran. However, when I move my model code and data to a NEC computer which used sxmpif90 fortran compiler, I have an error: *** 194 Could not allocate in ALLOCATE PROG=mod_ocean.allocate_ocean ELN=94(400500c88) Called...

CBian

Very strange!! ZQ told me to move river point seaward to see what will happen. I tried it, same abnormal salt result. However, I tried to move river point landward. Now the model runs nearly one month, the salt is good and the river flow into sea too. But this means I put the river cell in land. I ...

CBian

Thanks, kate, you said very clear. But my result is different with what you said. mathieu said the river_transport is negative when river flows westward or southward. After change the sign, I can add all river successfully. However, one of my river is southward flow. I set the river point at river f...

CBian

Thank you very much, mathieu, this is really helpful!

I never know the river_transport could be negative. When I change the sign of the river_transport, I can add all the rivers in model domain.

I cannot find map_rivers.m in the web, could you attach this code?

CBian

Actually this is the same problem I have post in ROMS Bugs: Minus sediment concentration in river sediment case . But I think this is not a sediment concentration problem any more. So I make a new post so that more people can have a look. These days I tried to figure out how the river point was adde...

CBian

I still cannot understand how the river location was set in the model. So I tried to put all the rivers in my domain to the model. There are 10 rivers in this area. After lots of tests, I added 7 of them to the model successfully. However, other 3 rivers cannot be accepted near their location. These...

CBian

coelho wrote:Hi CBian,

Did you try to put any sediment with the new river location?

I still didn't put sediment into river, but I will do it as soon as possible and see if I still get negetive concentration near river mouth.

I am very glad to know the river point can be located at "land".

CBian

Very strange!! ZQ told me to move river point seaward to see what will happen. I tried it, same abnormal salt result. However, I tried to move river point landward. Now the model runs nearly one month, the salt is good and the river flow into sea too. But this means I put the river cell in land. I r...

CBian

i am still not clear about the direction of the river flow. If you set direction = 0, this just means it is in the XI ('x') direction. The magnitude of the flow can be positive or negative, which makes the transport be to the left or to the right. If you are specifying flow at a u point, and the fl...

CBian

I change the river position to the cell which has abnormal large salt. The result shows that this cell(now the river cell) have normal salt value, however, the cell next to this cell have more than 100 salt value. So abnormal salt always appeared next to the river cell.

CBian

is this flow correct? it seems that the net flow would be out of the domain, into the river. A In my opinion, the flow is correct. Because I check the current and tide harmonic constant results in this domain, the result looks well. http://b67.photo.store.qq.com/http_imgload.cgi?/rurl4_b=107c9cde58...

CBian

is the river flow coming into the model or out of the model? my river_direction is 0. I add tide in the model. http://b70.photo.store.qq.com/http_imgload.cgi?/rurl4_b=e0a48a1f6296ad7961d7adf695f736285e88837409b61a82a7aa81a4bc2ec2d19308d073d05c6b2bfd735c81a5baab80fbf59bbc32cdfc8851f93c1034847f9e522d...

CBian

Hi, I refined my domain, and the Dx, Dy are changed from 18km to 5km. I didn't add sediment to the model this time, just want to see the salt result near river mouth. Both of two rivers can be found this time. However, one water cell next to the river cell always has a largest salt value(salt >100)....

CBian

Seems that everything is fine. Try to put all the 20 levels of river_Vshape = 0.05(attached) and see if you find any sand coming from the rivers. Your river_sand_03 is always zero. Is that correct? There are 3 kind of sand in my model. I just 2 kind of sand from rivers, the third kind of sand doesn...

CBian

GET_NGFLD - river runoff XI-positions at RHO-points (Min = 3.40000000E+01 Max = 3.40000000E+01) GET_NGFLD - river runoff ETA-positions at RHO-points (Min = 7.20000000E+01 Max = 7.20000000E+01) The Min and Max values for XI-positions and ETA-positions are the same. Are you sure you are using two riv...

CBian

I expended my model domain to include more area. However, in my new model case (still include above two rivers), the river have no influence on the model result. It seems no river added to the model. The salt, temperature, sediment concentration have no difference with the model results without rive...

CBian

I mean that the river should be near to land point. I had some problems with rivers located not exactly beside a land point. Try to move the location of the river that you are getting negative values and see what happens. Thanks, coelho. Yes, I always make my river point next to land point. I tried...

CBian

Did you check if your river is located beside a land point? Maybe the negative sediment values are because of that. The river is located beside two land points. Did you mean the river should not near to land point? But another river is located beside three land points, its concentration is positive.

CBian

Here is an example: IF the flow is out of the grid, and the flow is going to take more sediment out of a cell than actually exists in the cell, then the concentrations will go negative. Dear Warner, thanks for your advice. I am confused about this. I have 3 kind of sediment in my model(sand_01,sand...

CBian

can u check if the flows are correctly coming out of the river? so look at the history file and compute the flow out of the river. Is the flow direction tidal or does it always flow into the grid? I check the flows at river point, it is tidal, sometimes flow into the grid, sometimes flow out. At fi...

CBian

Dear linzhenhua,

I always used TS_MPDATA in my river model. But the sediment concentration is always minus near one river mouth. Thank you for your suggestion.

CBian

Hello! I add two rivers in my model. When I check the sediment concentration, I found, near the river mouth, one river have normal sediment concentration and another always have minus sediment concentration. And both river have positive sediment mass. Is it a bug? The sediment concentration cannot b...

CBian

Dear jcwarner, I tried use Nbed=1,and Nbed=5 in my model, other configurations are same, but the sediment mass transport results are very different. Even use same Nbed=5, I tried different bed layer thickness, after 4 months, the mass transport results are also have obvious different. I use the late...

CBian

Dear ZhaoQiang, Thanks for your suggestion, it's really helpful. I made a stupid mistake, the vertical viscosity near coast is about 0.01 N/m2 not 0.001 N/m2. This doesn't seem too small. In my opinion, when there is a very strong northward wind (eg. 10 m/s), the wind-induced current near the coast(...

CBian

Dear Qiang, Thanks for your suggestion. I check the frc_uvstress.cdl: short sustr(sms_time, eta_u, xi_u) ; sustr:long_name = "surface u-momentum stress" ; sustr:units = "Newton meter-2" ; sustr:time = "sms_time" ; sustr:scale_factor = 0.0005 ; short svstr(sms_time, eta_...

CBian

Dear Roms users, I cannot get a good wind-driven current result near the coast, so I output the vertical viscosity (AKv) from model. The AKv results are very strange. I check some papers, and the AKv should have a bigger value in mixed coastal area than in deep water. However, when I used MY25_MIXIN...

CBian

Dear ROMS users, I want to user ROSM simulate the sediment transportation from coast to deep ocean. My domain is the whole East China Sea(117E-133E， 24N-41N)， and the resolution is 10/60 degree. I know the importance of wave in suspending sediments at coastal area, so I need couple the SWAN model to...

CBian

Hi, ROMS users, I change my server to a NEC SX-9 computer recently and download the newest version of ROMS (revision of 513). When I compile the code, there is error: f90: distribute.f90, distribute_mod: There are 8 errors and 28 warnings. sxf90 fatal : /SX/usr/lib/f90com command error : 1 make: ***...

CBian

Stone, John, cdw, thanks for your kindly suggestions! My classmate told me that the quadratic drag coefficent is 0.0025 in the deep sea but is 0.001 near the coast. And I checked the realistic ocean data in my domain, that's true. So I should use a changable drag coefficent in my domain. I believe t...

CBian

Hello, ROMS users! When I simulated the tide in Bohai and Yellow Sea, I confront a problem: the amplitude of the tide is about 30% less than the real amplitude near the coast （fig. a）. First, I thought maybe the error was caused by the topography. Therefore, I try the Etopo5 topography data. However...

CBian

but if I define ana_mask and define masking together, the model output are all 1e37; if I just define ana_mask only, model results looked reasonable. Could we defined ana_mask and define masking together? If MASKING is defined, the model will read mask_rho,mask_u,mask_v,mask_psi from our grid.nc da...

Search found 39 matches