Stations and sediment variables

[jprinehimer]

How do I get sediment bottom or bed parameters written out to the stations files?

I turned on #define STATIONS, created a stations.in file, and set some of the Sout(idBott) parameters. I get water column values and I also get mud_01, mud_02, but I don't get any bed surface parameters. I also can't seem to figure out where in the code it turns on writing out bed parameters like fraction and porosity. I have these turned on in sediment.in, but do you also turn them on in station.in?

Is there a typo in inp_par.F? In the stations section, a line reads:

Code: Select all

 ELSE IF (TRIM(KeyWord).eq.'Hout(idBott)') THEN

Shouldn't that be 'Sout(idBott)', the same as in the stations file?

[arango]

Yes, it needs to be Sout(idBott). Good Catch. I will correct the svn repository.

[jprinehimer]

So I tried getting sediment bed properties (like bed age, porosity, thickness, mass) into the stations file again. No luck. Surface properties (active layer thickness at least) are being output to the stations file.

However, there's nowhere to set any Sout flags for bed properties in the stations.in file. I tried setting:

Code: Select all

Sout(iMmass)=T

to get bed sediment mass in the stations file, but that didn't work either.

Has this feature not been implemented yet or am I doing something wrong?

J.Paul