Using hydrodynamic results recorded

MOHID Water forum. Questions and discussion around MOHID Water
guifranz@gmail.com
MOHID Beginner
Posts: 12
Joined: Thu Jan 01, 1970 00:00

Using hydrodynamic results recorded

Postby guifranz@gmail.com » Mon Jan 24, 2011 16:42

Hi,

I need to run the Lagrangian and Water quality modules using the results of a baroclinic hydrodynamic model that I already ran in a reservoir. I ran this model with temperature in the water property module, and the both results of hydrodynamic and temperature seem good.

After that, I managed to run the Lagrangian module with the recorded results of the hydrodynamic module. But I had problems with the temperature results (NaN). So, I need to fix it to be able to run the water quality module. I have not changed anything associated with temperature, only the Hydrodynamic, HydrodynamicFile, Model and Lagrangian modules. I want the same results of temperature in the new simulation.

Do I need to record the temperature results like I did to the hydrodynamic module? Is there something different that I have to do to run the water quality module with the previous results of the hydrodynamic module?

Here are some of the modules that I set in both simulations.

Model_1
START : 2008 1 1 0 0 0
END : 2008 1 30 0 0 0
DT : 6.
VARIABLEDT : 0
GMTREFERENCE : 0
SPLITTING : Double_Splitting

Model_2
START : 2008 1 2 0 0 0
END : 2008 1 30 0 0 0
DT : 86400.
VARIABLEDT : 0
GMTREFERENCE : 0
SPLITTING : No_Splitting
LAGRANGEANE : 1
LAGRANGIAN : 1

Hydrodynamic_1
ADV_METHOD_H : 4
TVD_LIMIT_H : 4
BAROCLINIC : 1
CONTINUOUS : 0
SUBMODEL : 0
TIDE : 0
WATER_DISCHARGES : 1
ENERGY : 0
TIME_SERIE : 1
TIDEPOTENTIAL : 0
WIND : 1
ATM_PRESSURE : 0
SURFACEWATERFLUX : 0
DATA_ASSIMILATION : 0
OUTPUT_TIME : 0 43200.
RECORDING : 1
TIME_SERIE_LOCATION : ..\\..\\GeneralData\\monitoring_stations3.dat

Hydrodynamic_2
EVOLUTION : Read_File
WATER_DISCHARGES : 1
OUTPUT_TIME : 0 3600.
TIME_SERIE_LOCATION : ..\\..\\GeneralData\\monitoring_stations3.dat


HydrodynamicFile_1
OUTPUT: 1
OUT_FIELD : C:\\Guilherme\\Mohid\\SaltoGrande-teste\\Simulation_2\\res\\Run1\\hidro.bin
TIME_INTEGRATION: 0
OUT_FILE_VERSION: 2
SPACE_INTEGRATION:0
DT_HYDROFILE: 86400

HydrodynamicFile_2
IN_FIELD : C:\\Guilherme\\Mohid\\SaltoGrande-teste\\Simulation_2\\res\\Run1\\hidro.bin
IN_FILE_TYPE : BeginEnd_type
IN_FILE_VERSION : 2

WaterProperties_1 = WaterProperties_2
ADV_METHOD_H : 4
TVD_LIMIT_H : 4
OUTPUT_TIME : 0. 43200.
TIME_SERIE_LOCATION : ..\\..\\GeneralData\\monitoring_stations3.dat
<beginproperty>
NAME : temperature
UNITS : ºC
DESCRIPTION : No description was given.
IS_COEF : 1
PARTICULATE : 0
INITIALIZATION_METHOD : CONSTANT
DEFAULTVALUE : 22.
DEFAULTBOUNDARY : 0.
BOUNDARY_INITIALIZATION : EXTERIOR
BOUNDARY_CONDITION : 2
OLD : 0
ADVECTION_DIFFUSION : 1
DISCHARGES : 1
DISCHARGES_TRACKING : 0
SURFACE_FLUXES : 0
DATA_ASSIMILATION : 0
ADVECTION_H_IMP_EXP : 0
TIME_SERIE : 1
BOX_TIME_SERIE : 0
OUTPUT_HDF : 1
<endproperty>


Thank you for your help.

frank.maretec@ist.utl.pt
MOHID Guru
Posts: 151
Joined: Thu Jan 01, 1970 00:00
Location: Technical University of Lisbon
Contact:

Postby frank.maretec@ist.utl.pt » Wed Jan 26, 2011 15:59

Hello,

from what I understand of your problem you should record also the temperature.

To me it seems that you problem is related with numerical stability. When you record you run with a DT of 6 secs (quiet small). Everything works fine.

When you run again, you calculate the transport process of temperature with a DT of 1 day, which most probably is to high.

So you should consider recording water temperature as well. How? I don't know from scratch but digging in the forum might bring up a solution.

Hope this helps,
Frank
In order to understand recursivity you must first understand recursivity.

guifranz@gmail.com
MOHID Beginner
Posts: 12
Joined: Thu Jan 01, 1970 00:00

Postby guifranz@gmail.com » Thu Jan 27, 2011 14:05

Hi Frank,

Thank you for your attention.

I solved the problem with temperature choosing DATA_ASSIMILATION : 1 in the WaterProperties module and defining the Assimilation module:
<beginproperty>
NAME: temperature
UNITS: °C
DIMENSION: 3D

<<begin_field>>
DEFAULTVALUE: 22
TYPE_ZUV: z
INITIALIZATION_METHOD: HDF
FILE_IN_TIME: HDF
FILENAME: C:\\Guilherme\\Mohid\\SaltoGrandeteste\\Simulation_2\\res\\WaterProperties_1.hdf5
<<end_field>>

<<begin_coef>>
DEFAULTVALUE: 0.001
TYPE_ZUV: z
FILE_IN_TIME: NONE
REMAIN_CONSTANT:1
INITIALIZATION_METHOD:constant
<<end_coef>>

<endproperty>

It works well and I had the same results for temperature of the first run. But when I try to run the WaterQuality module I had problems in the results of all properties. So now I am trying to record the hydrodynamic results more often, with smaller DT_HYDROFILE. I am using constant wind and discharges for tests.

I need to calibrate the water quality model for a long data series, and I don`t want to run the hydrodynamic model every time, because it spends a lot of time and will be impractical. I will work with some big reservoirs using 3D hydrodynamic, and variable wind, discharges, radiation, air temperature...

If someone has any experience in calibrating the water quality model with previous results of the hydrodynamic model and can help me, I`ll be thankful.

renanbribeiro@yahoo.com.br
MOHID User
Posts: 26
Joined: Thu Jan 01, 1970 00:00
Location: Universidade Santa Cecília - Santos, Brazil

Re: Using hydrodynamic results recorded

Postby renanbribeiro@yahoo.com.br » Thu Dec 22, 2011 02:41

Hello,

I am also trying to run Water quality modules using the results of a hydrodynamic model and also the results from temperature,salinity and coeshive sediments that I already ran.
When I run with the waterquality properties (oxygen, phyto, zoo, etc), the results from this properties are very bad (with spurius values)

I saved the hydrodymamic and waterproperties with a time step of 3600s, and I am using a DT = 3600s on the waterquality simulations.

Did anyone already use the previous results of the hydrodynamic model to simulate water quality? :D

Regards,
Renan


Return to “MOHID Water”

Who is online

Users browsing this forum: No registered users and 1 guest