Module Geometry - Revision history

Davidbrito: /* Modification Phase */

2018-03-12T12:41:57Z

‎Modification Phase

Davidbrito: /* Modification Phase */

2018-03-12T12:37:35Z

‎Modification Phase

Davidbrito: /* Modification Phase */

2018-03-12T12:34:02Z

‎Modification Phase

Davidbrito: /* Modification Phase */

2018-03-12T12:28:03Z

‎Modification Phase

Davidbrito: /* Modification Phase */

2018-03-12T12:23:55Z

‎Modification Phase

Davidbrito: /* Lagrangian Process */

2018-03-12T12:17:25Z

‎Lagrangian Process

Davidbrito: /* Domains Definition with Lagrangian Method */

2013-08-28T17:23:06Z

‎Domains Definition with Lagrangian Method

Davidbrito: /* Input data file */

2013-07-03T18:24:12Z

‎Input data file

Davidbrito: /* Sigma Domain on Top of Cartesian both with Lagrangian method */

2013-07-03T18:22:54Z

‎Sigma Domain on Top of Cartesian both with Lagrangian method

Davidbrito: /* Examples */

2013-07-03T18:21:26Z

‎Examples

@@ Line 114: / Line 114: @@
-At the end, is tested if the bottom layer (the last face to be processed is the bottom layer top face) is estimated with a thickness lower than the minimum. Two options are evaluated, testing if the water column height is lower than the sum of the layers minimum thichnesses:
+At the end, is tested if the bottom layer (the last face to be processed is the bottom layer top face) was estimated with a thickness lower than the minimum.
-# if yes, all layers are at its maximum compression and can not go thinner, so the process puts all the layers with an average thickness (water column / number of layers) to avoid very thin layers at the bottom that can create model instabilities.
+To avoid very thin layers at the bottom that can create model instabilities, two options are evaluated, testing if the water column height is lower than the sum of the layers minimum thichnesses:
 # if not, then minimum thickness is set to bottom layer and the process goes from bottom to top setting the minimum thickness where this condition is not met. Until it reaches a cell that has enough thickness (greater than its minimum).
-The above processes are not dependent on vertical velocity (trying to avoid thickness lower than user assigned) and can lead to stability problems. If that happens try to use more loose compression and expansion parameters (increase MINEVOLVELAYERTHICKNESS) and the process itself will be able to evolve to lower thicknesses.
+The above processes are not dependent on vertical velocity (trying to avoid thickness lower than user assigned) and ultimately can lead to stability problems. If that happens try to use more loose compression and expansion parameters (increase MINEVOLVELAYERTHICKNESS) and the process itself will be able to evolve to lower thicknesses.

@@ Line 116: / Line 116: @@
 At the end, is tested if the bottom layer (the last face to be processed is the bottom layer top face) is estimated with a thickness lower than the minimum. Two options are evaluated, testing if the water column height is lower than the sum of the layers minimum thichnesses:
-# if yes, all layers are at its maximum compression and can not go thinner, so the process puts all the layers with an average thickness (water column / number of layers)
+# if yes, all layers are at its maximum compression and can not go thinner, so the process puts all the layers with an average thickness (water column / number of layers) to avoid very thin layers at the bottom that can create model instabilities.
 # if not, then minimum thickness is set to bottom layer and the process goes from bottom to top setting the minimum thickness where this condition is not met. Until it reaches a cell that has enough thickness (greater than its minimum).

@@ Line 116: / Line 116: @@
 At the end, is tested if the bottom layer (the last face to be processed is the bottom layer top face) is estimated with a thickness lower than the minimum. Two options are evaluated, testing if the water column height is lower than the sum of the layers minimum thichnesses:
-# if yes, then the process puts all the layers with an average thichness (water column / number of layers)
+# if yes, all layers are at its maximum compression and can not go thinner, so the process puts all the layers with an average thickness (water column / number of layers)
-# if not, then minimum thicknesses are set from bottom to top until a cell has tickhness greater than its minimum
+# if not, then minimum thickness is set to bottom layer and the process goes from bottom to top setting the minimum thickness where this condition is not met. Until it reaches a cell that has enough thickness (greater than its minimum).
 The above processes are not dependent on vertical velocity (trying to avoid thickness lower than user assigned) and can lead to stability problems. If that happens try to use more loose compression and expansion parameters (increase MINEVOLVELAYERTHICKNESS) and the process itself will be able to evolve to lower thicknesses.

@@ Line 108: / Line 108: @@
-First a new estimation of SZZ is computed. Starting at the top layer, the displacement of each face is estimated using the last vertical velocity computed (displacement [m] is vertical velocity[m/s] * dt [s]).
+First a new estimation of SZZ is computed. Surface and bottom faces SZZ's are fixed.
-At each layer, the estimation of SZZ based on vertical velocity is tested to check if they are within the maximum deformation allowed for the above and below cells of that face. In case the test fails, SZZ is corrected to the minimum of the deformations allowed.
+At each layer, the estimation of SZZ based on vertical velocity is tested to check if they are within the maximum deformation allowed for the above and below cells of that face. In case the test fails, SZZ displacement is corrected to the minimum of the deformations allowed.
 At the end, is tested if the bottom layer (the last face to be processed is the bottom layer top face) is estimated with a thickness lower than the minimum. Two options are evaluated, testing if the water column height is lower than the sum of the layers minimum thichnesses:
-# if yes, then the option is to put all the layers with an average thichness (water column / number of layers)
+# if yes, then the process puts all the layers with an average thichness (water column / number of layers)
 # if not, then minimum thicknesses are set from bottom to top until a cell has tickhness greater than its minimum

@@ Line 108: / Line 108: @@
-First a new estimation of SZZ is computed. Starting at the top layer, the displacement of each face is estimated using the last vertical velocity computed. And the new value of SZZ is evaluated.
+First a new estimation of SZZ is computed. Starting at the top layer, the displacement of each face is estimated using the last vertical velocity computed (displacement [m] is vertical velocity[m/s] * dt [s]).
-The new values of SZZ are tested to check if they are within the maximum deformation allowed for the above and below cells of that face and corrected if not.
+At each layer, the estimation of SZZ based on vertical velocity is tested to check if they are within the maximum deformation allowed for the above and below cells of that face. In case the test fails, SZZ is corrected to the minimum of the deformations allowed.
-At the end, is tested if the bottom layer (the last one to be processed) is estimated with a thickness lower than the minimum. Two options are evaluated testing if the water column height is lower than the sum of the layers minimum thichnesses:
+# if yes, then the option is to put all the layers with an average thichness (water column / number of layers)
-#1 if yes, then the option is to put all the layers with an average thichness (water column / number of layers)
+The above processes are not dependent on vertical velocity (trying to avoid thickness lower than user assigned) and can lead to stability problems. If that happens try to use more loose compression and expansion parameters (increase MINEVOLVELAYERTHICKNESS) and the process itself will be able to evolve to lower thicknesses.

@@ Line 104: / Line 104: @@
 ====Modification Phase====
-In case of Lagrangian evolution the grid is deformed at each time step. First a new estimation of SZZ is computed. Starting at the top layer the displacement of each layer is estimated using the last vertical velocity computed. And the new value of SZZ is evaluated.
+In case of Lagrangian evolution the grid is deformed at each time step.
 The new values of SZZ are tested to check if they are within the maximum deformation allowed and corrected if not. Using the final values of SZZ all the geometric properties are computed (distances, areas and volumes) and the vertical fluxes between cells are corrected to account for grid deformation (in ModuleHydrodynamic).
 == Bathymetry consistency diagnostic ==

@@ Line 201: / Line 201: @@
   LAGRANGIAN                : 1
   MINEVOLVELAYERTHICKNESS   : 0.5
   <enddomain>
@@ Line 218: / Line 217: @@
   LAGRANGIAN                : 1
   MINEVOLVELAYERTHICKNESS   : 0.5
   <enddomain>
@@ Line 235: / Line 233: @@
   LAGRANGIAN                : 1
   MINEVOLVELAYERTHICKNESS   : 0.5
   <enddomain>
@@ Line 247: / Line 244: @@
   LAGRANGIAN                : 1
   MINEVOLVELAYERTHICKNESS   : 0.5
   <enddomain>

@@ Line 148: / Line 148: @@
     DISPLACEMENT_LIMIT          : real           [1000]         !Maximum displacement in meters (if LAGRANGIAN : 1)
   <enddomain>
 ===Examples===
@@ Line 181: / Line 180: @@
   MININITIALLAYERTHICKNESS  : 0.2
   <enddomain>
 ====Domains Definition with Lagrangian Method====

@@ Line 224: / Line 224: @@
 =====Sigma Domain on Top of Cartesian both with Lagrangian method=====
-!General options
+ !General options
   MINIMUMDEPTH              : 0.1
   FACESOPTION               : 2