Personal tools

Difference between revisions of "Module LagrangianGlobal"

From MohidWiki

Jump to: navigation, search
(→‎Overview)
 
 
(50 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 
== Overview ==
 
== Overview ==
 +
 +
The LangrangianGlobal Module is a deep upgrade of the [[Module Lagrangian|Lagrangian Module]]. These upgrade has three main goals:
 +
 +
1 - Run only one lagrangian model in a multi-nesting implementation. This way particles can go from grid to grid without any problems ([[Media:SideBySideCloud.gif|see example of particles moving between grids]]). In the Lagrangian module particles are destroyed when they leave the grid where they were emitted. This is the reason way a lagrangian model is runned for each grid in the Lagrangian Module;
 +
 +
2 - Centralise all the interpolation and particle location (in a [[Grid|grid]]) methods in the [[Module HorizontalGrid|HorizontalGrid Module]] and (see source code [http://mohid.codeplex.com/SourceControl/changeset/view/64163#1208106]);
 +
 +
3 - Avoid redundances in the particle origin definitions. This way is possible to compress the input file (Lagrangian_x.dat).
 +
 +
 +
'''If the mohid user wants to activate this module it needs in the mohid compilation process to predefined a preprocessor symbol called _LAGRANGIAN_GLOBAL_.'''
 +
 +
== Related projects ==
 +
 +
This development was made by [http://www.Hidromod.com Hidromod] in the framework of the follow projects:
 +
 +
- [http://www.maretec.mohid.com/projects/easy/webgis/easy/easymap.aspx EASY] - Interreg project - Hidromod subcontract by [http://www.utl.ist.pt IST] - Debry tracjetory forecast in ocean/coastal domains
 +
 +
- [http://www.aquaplansistemas.es/es/soluciones/hidrometeorologia/cowama-gestion-calidad-aguas-de-ba%C3%B1o CoWaMa] Hidromod subcontract by [http://www.Clabsa.es Clabsa] and [http://www.lyonnaise-des-eaux.fr/ Lyonnaise-des-Eaux] - Forecast of fecal contamination in beaches
 +
 +
- [http://www.hidromod.com/Lenvis/LenvisMohidOnline/LENVIS/MohidOnline/ Lenvis] - FP7 - Hidromod [http://www.lenvis.eu/Partners.aspx consortium member] - Web model service focus in simulating the impact of emergency fecal contamination discharges;
 +
 +
- [http://www.argomarine.eu/public Argomarine] - FP7 project - Hidromod subcontract by [http://www.ualg.pt Algarve University] - FP7 - Implementation of a oil spill forecast model in the Toscanny region (Italy)
 +
 +
- [http://www.project-easy.info EASYCO] - Interreg project - Hidromod subcontract by [http://www.utl.ist.pt IST] - Data analysis and oil spill modelling via web
 +
 +
==Main upgrades==
 +
 +
=== Goal 1 - Run only one Lagrangian Model ===
 +
 
The LagrangianGlobal module is a module very similar to the [[Module Lagrangian|Lagrangian module]]. The input data file keywords are exactly the same. The modules give the same lagrangian result if the user only runs one model (no nesting).  
 
The LagrangianGlobal module is a module very similar to the [[Module Lagrangian|Lagrangian module]]. The input data file keywords are exactly the same. The modules give the same lagrangian result if the user only runs one model (no nesting).  
 
In the case of the first module when are run several nesting levels. The user only define one lagrangian input data file in the data file of the first nesting level (first model in the tree.dat). Each lagrangian tracer will use the hydrodynamic field with the higher priority depending on the tracer position. By default the priority is define  inverting the tree.dat order.  
 
In the case of the first module when are run several nesting levels. The user only define one lagrangian input data file in the data file of the first nesting level (first model in the tree.dat). Each lagrangian tracer will use the hydrodynamic field with the higher priority depending on the tracer position. By default the priority is define  inverting the tree.dat order.  
Line 9: Line 39:
 
  <EndModelPriority>
 
  <EndModelPriority>
  
−
The first model are the one with the higher priority.
+
The model are order by decreasing priority.
 +
 
 +
=== Goal 2 - Centralize all interpolation and particle location methods ===
 +
 
 +
In the [[Module Lagrangian|Lagrangian module]] the interpolation and particle location methods were developed in a independent way of similar methods implemented in the [[Module HorizontalGrid|HorizontalGrid Module]] for the [[Module Hydrodynamic|hydrodynamic]] and [[Module WaterProperties|water properties]] nesting. Now the lagrangian mdoule uses the same interpolation and location methods as the Eulerian modules. The subroutines from the HorizontalGrid Module used in the LagrangianGlobal Module are:
 +
 
 +
- GetXYInsideDomain
 +
 
 +
Check if a location X,Y (eg. particle location) is inside a specific model domain;
 +
 
 +
- GetXYCellZ
 +
 
 +
For a location X,Y (eg. particle location) and a specific model domain returns the correspondent cell (i,j) and even the relative position of the location inside the cell;
 +
 
 +
- GetCellZ_XY
 +
 
 +
Do the inverse process ofsuborutine  GetXYCellZ
 +
 
 +
- InterpolXYPoint
 +
 
 +
Interpolate a generic property for a location X,Y.
 +
 
 +
All this subroutines work in all [[Grid|grid]] types support by Mohid (uniform variable and constant spatial step, uniform rotate and curvilinear).
 +
 
 +
For the curvilinear case was necessary to develop a quite complex method to [[Relative position of a particle in a cell|calculate the relative position of a particle in a cell]].
 +
 
 +
 
 +
=== Goal 3 - Decrease redundances in Lagrangian_x.dat ===
 +
 
 +
To decrease this redundance the "clone origin" concept was created. There are several [[Module Lagrangian#Concepts|Lagrangian Origins]] options.
 +
 
 +
== Ilustrative examples of new options ==
 +
 
 +
=== Clone Origin concept ===
 +
 
 +
=== Fecal contamination variable discharges ===
 +
 
 +
The keywords used to controle the discretization (number of particles) in a operational system focus in the forecast of water quality based in Lagrangian discharges are the follow:
 +
 
 +
The maximum volume allowed for each particle (in this case 1 m3)
 +
<span style="color:#800000">POINT_VOLUME                  : 1</span>
 +
 
 +
Time between particle emission in a origin (in this case 100 s)
 +
<span style="color:#800000">DT_EMIT                        : 100</span>
 +
 
 +
Time the particle takes to double its volume due to turbulence(in this case 7200 s or 2 h)
 +
<span style="color:#800000">TVOL200                        : 7200.</span>
 +
 
 +
Factor that multiply by the particle initial volume gives the volume limit above which the particle is killed.
 +
<span style="color:#800000">VOLFAC                        : 100.</span>
 +
 
 +
True (or 1) if the particles volume is compute automaticly function of a flow time serie
 +
<span style="color:#800000">ESTIMATE_MIN_VOL              : 1</span>
 +
 
 +
In this case at least for each run for each origin at least 1000 particles will be emitted along the run.
 +
<span style="color:#800000">MAX_PART                      : 1000</span>
 +
 
 +
====Particle volume computation procedure per Origin ====
 +
 
 +
'''Step 1'''
 +
 
 +
If the option ESTIMATE_MIN_VOL is ON (ESTIMATE_MIN_VOL : 1) a first estimation of the particles volumes for each run  (eg. 1 day simulation) and for each origin is computed:
 +
 
 +
Particle Volume estimated = max (Total Discharge Volume in origin X along the run / MAX_PART, POINT_VOLUME)
 +
 
 +
'''Step 2'''
  
−
If the mohid user wants to activate this module it needs to predefined a preprocessor symbol called _LAGRANGIAN_GLOBAL_ in the compilation phase of the mohid.
+
After for each emission iteration is computed the volume to be discharge. The period of each emission iteration is from Now to Now + DT_EMIT
 +
 
 +
Along this period is compute the total volume emitted in a specific emission iteration (Total Volume by Iteration). This volume is computed based in the flow time serie associated to the origin.
 +
 
 +
'''Step 3'''
 +
 
 +
For each emission iteration is compute the actual particle volume emitted
 +
 
 +
'''Step 3.1'''
 +
 
 +
Number of particles emitted per iteration and per origin = int (Total Volume per Iteration and per Origin / Particle Volume Estimated) + 1
 +
 
 +
'''Step 3.2'''
 +
 
 +
Actual Particle Volume  Emitted = Total Volume per Iteration and per Origin / Number of particles emitted per iteration and per origin
  
−
== Relative position of tracer in a generic cell ==
 
  
−
A fundamental in the "MOHID system" lagrangian approach is to know the relative position of the tracer in the hydrodynamic cell to speed up the interpolation procedure. This is quite easy for square cells however this not the case for cells of curvilinear grids . In this case was develop a subroutine that is able to compute the generic position of a tracer in a generic cell (see subroutine RelativePosition4VertPolygon in module functions). The algorithm used is described below.
 
  
−
[[Image:FromGenericToSquare.jpg|850px|thumb|center|'''Relative position of tracer in a generic cell''']]
 
 
   
 
   
 +
The follow example used as default in the [http://www.aquaplansistemas.es/es/soluciones/hidrometeorologia/cowama-gestion-calidad-aguas-de-ba%C3%B1o CoWaMa] system.
 +
 +
 +
<span style="color:#006400"><BeginGeneralKeywords></span>
 +
  <span style="color:#006400">[[OUTPUT TIME|OUTPUT_TIME]]                    : 0 3600</span>
 +
 +
  <span style="color:#006400">[[RESTART FILE OUTPUT TIME|RESTART_FILE_OUTPUT_TIME]]      : 0 3600</span>
 +
  <span style="color:#006400">[[RESTART FILE OUTPUT TIME|RESTART_FILE_OVERWRITE]]        : 0</span>
 +
 +
  <span style="color:#006400">DT_PARTIC                      : 100</span>
 +
 +
  <span style="color:#006400">PARTIC_BOX                    : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\BOXRESTARTFILE_Lagrangian.dat</span>
 +
 +
  <span style="color:#006400">[[How to create MOHID timeseries outputs#Step 3 - Configure the model so as to output time series results|TIME_SERIE_LOCATION]]            : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\TimeSeriesLocationLag.dat</span>
 +
 +
  <span style="color:#006400">OUTPUT_CONC                    : 2</span>
 +
 +
  <span style="color:#006400">LW_PERCENTAGE                  : .4</span>
 +
  <span style="color:#006400">LW_EXTINCTION_COEF            : .3333333</span>
 +
  <span style="color:#006400">SW_PERCENTAGE                  : .6</span>
 +
  <span style="color:#006400">SW_EXTINCTION_TYPE            : 3</span>
 +
 +
  <span style="color:#006400">IGNORE_ON                      : 1</span>
 +
<span style="color:#006400"><EndGeneralKeywords></span>
  
−
<math> y_1  = a_1  + b_1 x_1  \wedge a_1  = {{x_d y_c  - x_c y_d } \over {x_d  - x_c }} \wedge b_1  = {{y_d  - y{}_c} \over {x_d  - x_c }} </math>
 
  
−
<math> y_2  = a_2  + b_2 x_2  \wedge a_2  = {{x_b y_a  - x_a y_b } \over {x_b  - x_a }} \wedge b_2  = {{y_b  - y{}_a} \over {x_b  - x_a }} </math>
 
  
−
<math> y_3  = a_3  + b_3 x_3  \wedge a_3  = {{x_a y_c  - x_c y_a } \over {x_a  - x_c }} \wedge b_3 = {{x_a - x_c } \over {y_a - y{}_c}} </math>  
+
<BeginModelPriority>
 +
  Level3ValidationEcoli
 +
  Level2Validation
 +
  <EndModelPriority>  
  
−
<math> y_4 = a_4 + b_4 x_4 \wedge a_4  = {{x_b y_d  - x_d y_b } \over {x_b  - x_d }} \wedge b_4  = {{x_b  - x_d } \over {y_b  - y{}_d}} </math>
+
<BeginOrigin>
 +
  GROUP_ID                      : 1
 +
  ORIGIN_NAME                    : Default
 +
  DEFAULT                        : 1
  
 +
EMISSION_SPATIAL              : Point
 +
EMISSION_TEMPORAL              : Continuous
 +
EMISSION_ON                    : 1
 +
OLD                            : 0
 +
POSITION_COORDINATES          : 0. 0.
 +
DEPTH_METERS                  : 0.
 +
MOVEMENT                      : SullivanAllen
 +
VARVELHX                      : 0.2
 +
VARVELH                        : 0.1
 +
TURB_V                        : Constant
 +
VARVELVX                      : 0.02
 +
VARVELV                        : 0.01
  
−
<math> y_g = {{a_3  - a_4 } \over {b_4  - b_3 }} \wedge x_g  = a_3  + b_3 y_g </math>
+
  NBR_PARTIC                    : 1
  
−
if (b1==b2) then
+
The maximum volume allowed for each particle (in this case 1 m3)
−
  <math> tg\alpha  = b_2</math> 
+
  <span style="color:#800000">POINT_VOLUME                  : 1</span>
−
else
 
−
<math> x_f  = {{a_1  - a_2 } \over {b_2  - b_1 }} \wedge y_f  = a_1  + b_1 x_f </math>  
 
−
<math> tg\alpha  = {{y_e  - y_f } \over {x_e  - x_f }} </math>  
 
−
endif
 
  
−
<math> y_h  = {{y_e  + \left( {a_3  - x_e } \right)Tg\alpha } \over {1 - b_3 Tg\alpha }} \wedge x_h  = a_3  + b_3 y_h </math>
 
  
−
if (b3==b4) then
+
Time between particle emission in a origin (in this case 100 s)
−
  <math> tg\beta  = b_3 </math>
+
  <span style="color:#800000">DT_EMIT                        : 100</span>
−
else
 
−
<math> x_g  = {{a_3  - a_4 } \over {b_4  - b_3 }} \wedge y_g  = a_3  + b_3 x_g </math>
 
−
<math> tg\beta  = {{y_e  - y_g } \over {x_e  - x_g }} </math>
 
−
endif
 
  
−
<math> y_i  = {{y_e  + \left( {a_1  - x_e } \right)Tg\beta } \over {1 - b_1 Tg\beta }} \wedge x_i = a_1  + b_1 y_i </math>  
+
Time the particle takes to double its volume due to turbulence(in this case 7200 s or 2 h)
 +
  <span style="color:#800000">TVOL200                        : 7200.</span>
  
 +
Factor that multiply by the particle initial volume gives the volume limit above which the particle is killed.
 +
<span style="color:#800000">VOLFAC                        : 100.</span>
  
−
<math> seg\_ac = \sqrt {\left( {x_a  - x_c } \right)^2 + \left( {y_a  - y_c } \right)^2 } </math>  
+
True (or 1) if the particles volume is compute automaticly function of a flow time serie
 +
  <span style="color:#800000">ESTIMATE_MIN_VOL              : 1</span>
  
−
<math> seg\_dc = \sqrt {\left( {x_d  - x_c } \right)^2  + \left( {y_d  - y_c } \right)^2 } </math>  
+
In this case at least for each run for each origin at least 1000 particles will be emitted along the run.
 +
<span style="color:#800000">MAX_PART                      : 1000</span>
  
−
<math> seg\_ic = \sqrt {\left( {x_i - x_c } \right)^2  + \left( {y_i - y_c } \right)^2 } </math>
+
  FLOW                          : -9e16
 +
FLOW_VARIABLE                  : 1
 +
FLOW_COLUMN                    : 2
 +
DISCHARGE_FILE                : *****.***
 +
!    T90_VAR_METHOD = Canteras            = 1
 +
  !    T90_VAR_METHOD =  Chapra              = 2
 +
!    T90_VAR_METHOD =  TimeSerie            = 3
  
−
<math> seg\_hc = \sqrt {\left( {x_h - x_c } \right)^2  + \left( {y_h - y_c } \right)^2 } </math>  
+
<<BeginProperty>>
 +
NAME                          : fecal coliforms
 +
UNITS                          : MPN/100ml
 +
CONCENTRATION                  : 0
 +
CONC_VARIABLE                  : 1
 +
CONC_COLUMN                    : 2
 +
DISCHARGE_FILE                : *****.***
 +
T90_VARIABLE                  : 0
 +
T90_VAR_METHOD                : 0
 +
T90                            : 10800
 +
T90_FILE                      : T902.dat
 +
T90_COLUMN                    : 2
 +
OUTPUT_HDF                    : 1
 +
TIME_SERIE                    : 0
 +
<<EndProperty>>
 +
<<BeginProperty>>
 +
NAME                          : escherichia coli
 +
UNITS                          : MPN/100ml
 +
CONCENTRATION                  : 0
 +
CONC_VARIABLE                  : 1
 +
CONC_COLUMN                    : 2
 +
DISCHARGE_FILE                : *****.***
 +
T90_VARIABLE                  : 0
 +
T90_VAR_METHOD                : 3
 +
T90                            : 10800
 +
T90_FILE                      : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\DailyT90.srw
 +
  T90_COLUMN                    : 2
 +
OUTPUT_HDF                    : 1
 +
TIME_SERIE                    : 0
 +
<<EndProperty>>
 +
<<BeginProperty>>
 +
NAME                          : temperature
 +
UNITS                          : ?C Summer temperature
 +
CONCENTRATION                  : 0
 +
OUTPUT_HDF                    : 1
 +
TIME_SERIE                    : 0
 +
<<EndProperty>>
 +
<<BeginProperty>>
 +
NAME                          : salinity
 +
UNITS                          : ppm
 +
CONCENTRATION                  : 0
 +
OUTPUT_HDF                    : 1
 +
TIME_SERIE                    : 0
 +
  <<EndProperty>>
 +
  <<BeginProperty>>
 +
NAME                          : cohesive sediment
 +
UNITS                          : mg/l
 +
CONCENTRATION                  : 0
 +
OUTPUT_HDF                    : 1
 +
TIME_SERIE                    : 0
 +
<<EndProperty>>
 +
<EndOrigin>
  
−
<math> x_e^' = {{Seg\_ic} \over {Seg\_dc}} </math>  
+
<BeginOrigin_Clone>
 +
CLONE                          : Default
 +
DEFAULT                        : 0
 +
GROUP_ID                      : 1
 +
ORIGIN_NAME                    : Discharge 1
 +
EMISSION_SPATIAL              : Point
 +
EMISSION_TEMPORAL              : Continuous
 +
DISCHARGE_FILE                : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\DISCHARGE_Discharge 1.dat
 +
POSITION_COORDINATES          : 1.3 41.11
 +
OLD                            : 0
 +
DEPTH_METERS                  : 1000000
 +
<<BeginProperty>>
 +
NAME                          : salinity
 +
UNITS                          : PSU
 +
CONCENTRATION                  : 0.2
 +
<<EndProperty>>
 +
<<BeginProperty>>
 +
NAME                          : temperature
 +
UNITS                          : ºC
 +
CONCENTRATION                  : 19.5
 +
<<EndProperty>>
 +
<<BeginProperty>>
 +
NAME                          : cohesive sediment
 +
UNITS                          : mg/l
 +
CONCENTRATION                  : 100
 +
<<EndProperty>>
 +
  <<BeginProperty>>
 +
NAME                          : fecal coliforms
 +
UNITS                          : MPN/100ml
 +
CONCENTRATION                  : 100000
 +
DISCHARGE_FILE                : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\DISCHARGE_FECAL_COLIFORMSDischarge 1.dat
 +
<<EndProperty>>
 +
<<BeginProperty>>
 +
NAME                          : escherichia coli
 +
UNITS                          : MPN/100ml
 +
CONCENTRATION                  : 100000
 +
DISCHARGE_FILE                : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\DISCHARGE_E_COLIDischarge 1.dat
 +
<<EndProperty>>
 +
<EndOrigin_Clone>
  
−
<math> y_e^'  = {{Seg\_hc} \over {Seg\_ac}} </math>
 
  
 
[[Category:Modules]]
 
[[Category:Modules]]
 
[[Category:MOHID Water]]
 
[[Category:MOHID Water]]

Latest revision as of 13:37, 22 March 2011

Overview

The LangrangianGlobal Module is a deep upgrade of the Lagrangian Module. These upgrade has three main goals:

1 - Run only one lagrangian model in a multi-nesting implementation. This way particles can go from grid to grid without any problems (see example of particles moving between grids). In the Lagrangian module particles are destroyed when they leave the grid where they were emitted. This is the reason way a lagrangian model is runned for each grid in the Lagrangian Module;

2 - Centralise all the interpolation and particle location (in a grid) methods in the HorizontalGrid Module and (see source code [1]);

3 - Avoid redundances in the particle origin definitions. This way is possible to compress the input file (Lagrangian_x.dat).


If the mohid user wants to activate this module it needs in the mohid compilation process to predefined a preprocessor symbol called _LAGRANGIAN_GLOBAL_.

Related projects

This development was made by Hidromod in the framework of the follow projects:

- EASY - Interreg project - Hidromod subcontract by IST - Debry tracjetory forecast in ocean/coastal domains

- CoWaMa Hidromod subcontract by Clabsa and Lyonnaise-des-Eaux - Forecast of fecal contamination in beaches

- Lenvis - FP7 - Hidromod consortium member - Web model service focus in simulating the impact of emergency fecal contamination discharges;

- Argomarine - FP7 project - Hidromod subcontract by Algarve University - FP7 - Implementation of a oil spill forecast model in the Toscanny region (Italy)

- EASYCO - Interreg project - Hidromod subcontract by IST - Data analysis and oil spill modelling via web

Main upgrades

Goal 1 - Run only one Lagrangian Model

The LagrangianGlobal module is a module very similar to the Lagrangian module. The input data file keywords are exactly the same. The modules give the same lagrangian result if the user only runs one model (no nesting). In the case of the first module when are run several nesting levels. The user only define one lagrangian input data file in the data file of the first nesting level (first model in the tree.dat). Each lagrangian tracer will use the hydrodynamic field with the higher priority depending on the tracer position. By default the priority is define inverting the tree.dat order. The user can specify the oder using a block define in the lagrangian input file where the user can define the model priority: 

<BeginModelPriority> 
Model name x
Model name y
<EndModelPriority>

The model are order by decreasing priority.

Goal 2 - Centralize all interpolation and particle location methods

In the Lagrangian module the interpolation and particle location methods were developed in a independent way of similar methods implemented in the HorizontalGrid Module for the hydrodynamic and water properties nesting. Now the lagrangian mdoule uses the same interpolation and location methods as the Eulerian modules. The subroutines from the HorizontalGrid Module used in the LagrangianGlobal Module are: 

- GetXYInsideDomain

Check if a location X,Y (eg. particle location) is inside a specific model domain; 

- GetXYCellZ

For a location X,Y (eg. particle location) and a specific model domain returns the correspondent cell (i,j) and even the relative position of the location inside the cell; 

- GetCellZ_XY

Do the inverse process ofsuborutine GetXYCellZ 

- InterpolXYPoint

Interpolate a generic property for a location X,Y.

All this subroutines work in all grid types support by Mohid (uniform variable and constant spatial step, uniform rotate and curvilinear).

For the curvilinear case was necessary to develop a quite complex method to calculate the relative position of a particle in a cell.


Goal 3 - Decrease redundances in Lagrangian_x.dat

To decrease this redundance the "clone origin" concept was created. There are several Lagrangian Origins options.

Ilustrative examples of new options

Clone Origin concept

Fecal contamination variable discharges

The keywords used to controle the discretization (number of particles) in a operational system focus in the forecast of water quality based in Lagrangian discharges are the follow:

The maximum volume allowed for each particle (in this case 1 m3) 

POINT_VOLUME                   : 1

Time between particle emission in a origin (in this case 100 s) 

DT_EMIT                        : 100

Time the particle takes to double its volume due to turbulence(in this case 7200 s or 2 h) 

TVOL200                        : 7200.

Factor that multiply by the particle initial volume gives the volume limit above which the particle is killed. 

VOLFAC                         : 100.

True (or 1) if the particles volume is compute automaticly function of a flow time serie 

ESTIMATE_MIN_VOL               : 1

In this case at least for each run for each origin at least 1000 particles will be emitted along the run. 

MAX_PART                       : 1000

Particle volume computation procedure per Origin

Step 1

If the option ESTIMATE_MIN_VOL is ON (ESTIMATE_MIN_VOL : 1) a first estimation of the particles volumes for each run (eg. 1 day simulation) and for each origin is computed:

Particle Volume estimated = max (Total Discharge Volume in origin X along the run / MAX_PART, POINT_VOLUME)

Step 2

After for each emission iteration is computed the volume to be discharge. The period of each emission iteration is from Now to Now + DT_EMIT

Along this period is compute the total volume emitted in a specific emission iteration (Total Volume by Iteration). This volume is computed based in the flow time serie associated to the origin.

Step 3

For each emission iteration is compute the actual particle volume emitted

Step 3.1

Number of particles emitted per iteration and per origin = int (Total Volume per Iteration and per Origin / Particle Volume Estimated) + 1

Step 3.2

Actual Particle Volume Emitted = Total Volume per Iteration and per Origin / Number of particles emitted per iteration and per origin



The follow example used as default in the CoWaMa system. 


<BeginGeneralKeywords>
 OUTPUT_TIME                    : 0 3600

 RESTART_FILE_OUTPUT_TIME       : 0 3600
 RESTART_FILE_OVERWRITE         : 0

 DT_PARTIC                      : 100

 PARTIC_BOX                     : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\BOXRESTARTFILE_Lagrangian.dat

 TIME_SERIE_LOCATION            : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\TimeSeriesLocationLag.dat

 OUTPUT_CONC                    : 2

 LW_PERCENTAGE                  : .4
 LW_EXTINCTION_COEF             : .3333333
 SW_PERCENTAGE                  : .6
 SW_EXTINCTION_TYPE             : 3

 IGNORE_ON                      : 1
<EndGeneralKeywords>



<BeginModelPriority>
Level3ValidationEcoli
Level2Validation
<EndModelPriority> 

<BeginOrigin>
GROUP_ID                       : 1
ORIGIN_NAME                    : Default
DEFAULT                        : 1

EMISSION_SPATIAL               : Point
EMISSION_TEMPORAL              : Continuous
EMISSION_ON                    : 1
OLD                            : 0
POSITION_COORDINATES           : 0. 0.
DEPTH_METERS                   : 0.
MOVEMENT                       : SullivanAllen
VARVELHX                       : 0.2
VARVELH                        : 0.1
TURB_V                         : Constant
VARVELVX                       : 0.02
VARVELV                        : 0.01

NBR_PARTIC                     : 1

The maximum volume allowed for each particle (in this case 1 m3) 

POINT_VOLUME                   : 1


Time between particle emission in a origin (in this case 100 s) 

DT_EMIT                        : 100

Time the particle takes to double its volume due to turbulence(in this case 7200 s or 2 h) 

TVOL200                        : 7200.

Factor that multiply by the particle initial volume gives the volume limit above which the particle is killed. 

VOLFAC                         : 100.

True (or 1) if the particles volume is compute automaticly function of a flow time serie 

ESTIMATE_MIN_VOL               : 1

In this case at least for each run for each origin at least 1000 particles will be emitted along the run. 

MAX_PART                       : 1000

FLOW                           : -9e16
FLOW_VARIABLE                  : 1
FLOW_COLUMN                    : 2
DISCHARGE_FILE                 : *****.***
!    T90_VAR_METHOD =  Canteras             = 1
!    T90_VAR_METHOD =  Chapra               = 2
!    T90_VAR_METHOD =  TimeSerie            = 3

<<BeginProperty>>
NAME                           : fecal coliforms
UNITS                          : MPN/100ml
CONCENTRATION                  : 0
CONC_VARIABLE                  : 1
CONC_COLUMN                    : 2
DISCHARGE_FILE                 : *****.***
T90_VARIABLE                   : 0
T90_VAR_METHOD                 : 0
T90                            : 10800
T90_FILE                       : T902.dat
T90_COLUMN                     : 2
OUTPUT_HDF                     : 1
TIME_SERIE                     : 0
<<EndProperty>>
<<BeginProperty>>
NAME                           : escherichia coli
UNITS                          : MPN/100ml
CONCENTRATION                  : 0
CONC_VARIABLE                  : 1
CONC_COLUMN                    : 2
DISCHARGE_FILE                 : *****.***
T90_VARIABLE                   : 0
T90_VAR_METHOD                 : 3
T90                            : 10800
T90_FILE                       : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\DailyT90.srw
T90_COLUMN                     : 2
OUTPUT_HDF                     : 1
TIME_SERIE                     : 0
<<EndProperty>>
<<BeginProperty>>
NAME                           : temperature
UNITS                          : ?C Summer temperature
CONCENTRATION                  : 0
OUTPUT_HDF                     : 1
TIME_SERIE                     : 0
<<EndProperty>>
<<BeginProperty>>
NAME                           : salinity
UNITS                          : ppm
CONCENTRATION                  : 0
OUTPUT_HDF                     : 1
TIME_SERIE                     : 0
<<EndProperty>>
<<BeginProperty>>
NAME                           : cohesive sediment
UNITS                          : mg/l
CONCENTRATION                  : 0
OUTPUT_HDF                     : 1
TIME_SERIE                     : 0
<<EndProperty>>
<EndOrigin>

<BeginOrigin_Clone>
CLONE                          : Default
DEFAULT                        : 0
GROUP_ID                       : 1
ORIGIN_NAME                    : Discharge 1
EMISSION_SPATIAL               : Point
EMISSION_TEMPORAL              : Continuous
DISCHARGE_FILE                 : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\DISCHARGE_Discharge 1.dat
POSITION_COORDINATES           : 1.3 41.11
OLD                            : 0
DEPTH_METERS                   : 1000000
<<BeginProperty>>
NAME                           : salinity
UNITS                          : PSU
CONCENTRATION                  : 0.2
<<EndProperty>>
<<BeginProperty>>
NAME                           : temperature
UNITS                          : ºC
CONCENTRATION                  : 19.5
<<EndProperty>>
<<BeginProperty>>
NAME                           : cohesive sediment
UNITS                          : mg/l
CONCENTRATION                  : 100
<<EndProperty>>
<<BeginProperty>>
NAME                           : fecal coliforms
UNITS                          : MPN/100ml
CONCENTRATION                  : 100000
DISCHARGE_FILE                 : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\DISCHARGE_FECAL_COLIFORMSDischarge 1.dat
<<EndProperty>>
<<BeginProperty>>
NAME                           : escherichia coli
UNITS                          : MPN/100ml
CONCENTRATION                  : 100000
DISCHARGE_FILE                 : ..\..\LEVEL1_VALIDATIONECOLI\..\Lagrangian\Level2Validation\data\DISCHARGE_E_COLIDischarge 1.dat
<<EndProperty>>
<EndOrigin_Clone>
Retrieved from "http://wiki.mohid.com/index.php?title=Module_LagrangianGlobal&oldid=4539"