HydroGeoSphere/Reading an Existing 3-D Grid
In some cases, the grid generation step can be very time-consuming. If this is so, then the following instruction can be used to read a grid which was generated in a previous run:
Read 3D gridEdit
The grid whose prefix matches that of the prefix.grok file will be read in.
These instructions show how to set up the grid generation section of the prefix.grok file to use the Read 3d grid instruction:
! Generate the grid for first run ! skip on read gb 2d grid laurel generate layers from gb 2d grid .true. ! zone by layer? .false. ! base elevation constant? laurel.nprop.Bedrock 1 ! layer Whole domain 20 ! sublayer .false. ! top elevation constant? laurel.nprop.Topography write faces and segments !skip off skip on ! Read previously defined grid for subsequent runs read 3D grid skip off end grid definition
In the first run, one can read the slice and generate the layered grid. It is important that the instruction write faces and segments be included if you are using any of the choose face or choose segment instructions later in the prefix.grok file.
On subsequent runs, one can skip over the grid generation commands and use the read 3D grid instruction instead.
- • • •
- • • •