Defining Cyclic water properties at open boundaries

[paulo.chambel.maretec@taguspark.pt]

Francisco,
in the hydrodynamic and water properties is possible to impose cyclic boundaries but is assumed that all open boundaries are cyclic.

Hydrodynamic.dat :

CYCLIC_BOUNDARY : 1

WaterProperties.dat

<beginproperty>
BOUNDARY_CONDITION : 8
NAME : salinity
UNITS : ºC
DESCRIPTION : salinity profile measure off the Nazare DEFAULTVALUE : 30
TYPE_ZUV : z
FILENAME : SalProfile.dat
INITIALIZATION_METHOD : PROFILE
FILE_IN_TIME : NONE
BOUNDARY_INITIALIZATION : INTERIOR
OLD : 0

<endproperty>

[paulo.chambel.maretec@taguspark.pt]

Francisco,
in the hydrodynamic and water properties is possible to impose cyclic boundaries but is assumed that all open boundaries are cyclic.

Hydrodynamic.dat :

CYCLIC_BOUNDARY : 1

WaterProperties.dat

<beginproperty>
!This is the keyword use to define the boundary condition in the water properties
BOUNDARY_CONDITION : 8
NAME : salinity
UNITS : ºC
DESCRIPTION : salinity profile measure off the Nazare DEFAULTVALUE : 30
TYPE_ZUV : z
FILENAME : SalProfile.dat
INITIALIZATION_METHOD : PROFILE
FILE_IN_TIME : NONE
BOUNDARY_INITIALIZATION : INTERIOR
OLD : 0

<endproperty>

[paulo.chambel.maretec@taguspark.pt]

Francisco,
in the hydrodynamic and water properties is possible to impose cyclic boundaries but is assumed that all open boundaries are cyclic.

Hydrodynamic.dat :

CYCLIC_BOUNDARY : 1

WaterProperties.dat

<beginproperty>
!This is the keyword use to define the boundary condition in the water properties
BOUNDARY_CONDITION : 8

NAME : salinity
UNITS : ºC
DESCRIPTION : salinity profile measure off the Nazare DEFAULTVALUE : 30
TYPE_ZUV : z
FILENAME : SalProfile.dat
INITIALIZATION_METHOD : PROFILE
FILE_IN_TIME : NONE

<endproperty>

[paulo.chambel.maretec@taguspark.pt]

Francisco,
in the hydrodynamic and water properties is possible to impose cyclic boundaries but is assumed that all open boundaries are cyclic.

Hydrodynamic.dat :

CYCLIC_BOUNDARY : 1

WaterProperties.dat

<beginproperty>
!This is the keyword use to define the boundary condition in the water properties
BOUNDARY_CONDITION : 8

NAME : salinity
UNITS : ºC
DESCRIPTION : salinity profile measure off the Nazare DEFAULTVALUE : 30
TYPE_ZUV : z
FILENAME : SalProfile.dat
INITIALIZATION_METHOD : PROFILE
FILE_IN_TIME : NONE

<endproperty>

[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>