Doesn't use river value

[changjw214]

Hello friends

I have a problem with my river discharge. I have one big river in my estuary.
and I set my own river discharge file.
But when I run my model(river discharge + tide),it doesn't show any change in salinity though it read my river values according to my log file.
Is there any problem in my river input or header???
I really don't have any idea...

Thanks for reading my problem~~

[kate]

Using ncview, the rivers file looks fine, though the timeseries of transport looks oddly stripy, alternating between zero and non-zero in time. Actually, ROMS is reading it:

Code: Select all

    GET_NGFLD   - river runoff  XI-positions at RHO-points
                   (Min =  1.97000000E+02 Max =  1.97000000E+02)
    GET_NGFLD   - river runoff ETA-positions at RHO-points
                   (Min =  1.20000000E+02 Max =  1.26000000E+02)
    GET_NGFLD   - river runoff direction
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - river runoff mass transport vertical profile
                   (Min =  2.00000000E-01 Max =  2.00000000E-01)
    GET_NGFLD   - river runoff mass transport,               t =     0 00:00:00
                   (Rec=0000001, Index=2, File: geum_river.nc)
                   (Tmin=          0.0000 Tmax=   34214400.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - river runoff potential temperature,        t =     0 00:00:00
                   (Rec=0000001, Index=2, File: geum_river.nc)
                   (Tmin=          0.0000 Tmax=   34214400.0000)
                   (Min =  8.89000000E+00 Max =  8.89000000E+00)
    GET_NGFLD   - river runoff salinity,                     t =     0 00:00:00
                   (Rec=0000001, Index=2, File: geum_river.nc)
                   (Tmin=          0.0000 Tmax=   34214400.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)

then

Code: Select all

   GET_NGFLD   - river runoff mass transport,               t = ***** 00:00:00
                   (Rec=0000002, Index=1, File: geum_river.nc)
                   (Tmin=          0.0000 Tmax=   34214400.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - river runoff potential temperature,        t = ***** 00:00:00
                   (Rec=0000002, Index=1, File: geum_river.nc)
                   (Tmin=          0.0000 Tmax=   34214400.0000)
                   (Min =  8.42700000E+00 Max =  8.42700000E+00)
    GET_NGFLD   - river runoff salinity,                     t = ***** 00:00:00
                   (Rec=0000002, Index=1, File: geum_river.nc)
                   (Tmin=          0.0000 Tmax=   34214400.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)

I think the time reporting assumes your time is in days, not seconds.

[changjw214]

I think, the reason why my transport shows alternatively zero and non-zero is the presence of dam.
And your comment means that my model read river_time value as in days???

[kate]

Sometimes ROMS can convert units - but it has to specifically check and do the conversion. The reason you get the five stars for the time is because it thinks the time should fit into five characters in the formatting. Do you think it should have read more than the first two time records for the rivers? Then you can try converting the river time to days - or you can investigate the code to see if it can read the units and behave accordingly. Personally, I like to watch this stuff in the debugger to see what the model is doing, but debugging seems to be a lost art.

[ckharris]

HI- I suggest trying to use a constant (steady) river input first until you get the time-stamps right, and see that the freshwater is really entering the estuary.