Impose water properties using a time serie at open boundary

[paulo.chambel.maretec@taguspark.pt]

Dear Paulo,

I would like to do is like, similar to the water discharges, to input a monthly value for some of the water properties, that could be read from a file. As an example one of my water properties would look like this:

<beginproperty>
NAME : temperature
UNITS : ºC
DESCRIPTION : temperature
IS_COEF : 1
PARTICULATE : 0
INITIALIZATION_METHOD : BOXES
FILENAME : C:\\MOHIDv1_3\\NewScheldt\\GeneralData\\Boxes4ordered.dat
BOXES_VALUES : 5.92 5.92 5.92 5.76 5.76 5.805 5.83 5.8 5.86 5.925 6.15 6.35
DEFAULTVALUE : 11.4
DEFAULTBOUNDARY : 11.4
BOUNDARY_INITIALIZATION : INTERIOR
BOUNDARY_CONDITION : 2
OLD : 0
ADVECTION_DIFFUSION : 1
DISCHARGES : 1
DISCHARGES_TRACKING : 1
VERTICAL_MOVEMENT : 0
SURFACE_FLUXES : 1
BOTTOM_FLUXES : 0
DATA_ASSIMILATION : 0
PARTITION : 0
WATER_QUALITY : 0
TIME_SERIE : 1
BOX_TIME_SERIE : 1
OUTPUT_HDF : 1
<endproperty>


What should I modify in the water properties to get this monthly input at the boundary? Thank you very much for the help, best regards,

Francisco

[paulo.chambel.maretec@taguspark.pt]

Dear Francisco,
it is not possible to define directly a value variable in time in the border. The only way is to active the assimilation module. The assimilation module allow the user to relax any property (flow or water) in any cell using the follow formula

dP/dt = (P(assimilation) - P) / Decay Time

If the decay time is very low the model impose in a cell the value P(assimilation) if the decay is very large the final result is equal to the model result. This type of methodology is use usually to impose boundary solution in a progressive way in space.

WaterProperties.dat

<beginproperty>
DATA_ASSIMILATION : 1

NAME : temperature
UNITS : ºC
DESCRIPTION : temperature
IS_COEF : 1
PARTICULATE : 0
INITIALIZATION_METHOD : BOXES
FILENAME : C:\\MOHIDv1_3\\NewScheldt\\GeneralData\\Boxes4ordered.dat
BOXES_VALUES : 5.92 5.92 5.92 5.76 5.76 5.805 5.83 5.8 5.86 5.925 6.15 6.35
DEFAULTVALUE : 11.4
DEFAULTBOUNDARY : 11.4
BOUNDARY_INITIALIZATION : INTERIOR
BOUNDARY_CONDITION : 2
OLD : 0
ADVECTION_DIFFUSION : 1
DISCHARGES : 1
DISCHARGES_TRACKING : 1
VERTICAL_MOVEMENT : 0
SURFACE_FLUXES : 1
BOTTOM_FLUXES : 0
DATA_ASSIMILATION : 0
PARTITION : 0
WATER_QUALITY : 0
TIME_SERIE : 1
BOX_TIME_SERIE : 1
OUTPUT_HDF : 1
<endproperty>

Assimilation.dat

<beginproperty>
NAME : temperature
UNITS : ºC
DIMENSION : 3D

<<begin_field>>
DEFAULTVALUE : 16.
TYPE_ZUV : z
FILE_IN_TIME : ASCII_FILE
FILENAME : TimeSerieTemperature.dat
<<end_field>>

<<begin_coef>>
DEFAULTVALUE : 1e32
TYPE_ZUV : z
FILE_IN_TIME : NONE
REMAIN_CONSTANT : 1
INITIALIZATION_METHOD : ASCII_FILE
FILENAME : DecayTimeZ.dat
<<end_coef>>

<endproperty>

TimeSerieTemperature.dat

Format time serie


DecayTimeZ.dat

Format Grid Data 3D

You need all the information use to define a grid that is in the Bahymetry plus
a block were you impose very low decay times in the boundary points. The temperature is always a 3D property heaven if the model has only one layer.

i j k Decay Time

<BeginGridData3D>
1 1 1 1e-6
2 1 1 1e-6
3 1 1 1e-6
4 1 1 1e-6
<EndGridData3D>