Syntax:
write_isurf group-ID Nx Ny Nz filename ablateID keyword args ...
precision arg = int or double
Examples:
write_isurf block 100 100 200 isurf.material.* ablation
Description:
Write a grid corner point file in binary format describing the current corner point values which define the current set of implicit surface elements. See the read_isurf command for a definition of implicit surface elements and how they are defined from grid conner point values. The surface file can be used for later input to a new simulation or for post-processing and visualization.
The specified group-ID is the name of a grid cell group, as defined by the group grid command, which contains a set of grid cells, all of which are the same size, and which comprise a contiguous 3d array, with specified extent Nx by Ny by Nz. These should be the same parameters that were used by the read_isurf command, when the original grid corner point values were read in and used to define a set of implicit surface elements. For 2d simulations, Nz must be specified as 1, and the group must comprise a 2d array of cells that is Nx by Ny. These are the grid cells that contain implicit surfaces.
Similar to dump files, the filename can contain a "*" wildcard character. The "*" character is replaced with the current timestep value. For example isurf.material.0 or isurf.material.100000.
The specified ablateID is the fix ID of a fix ablate command which has been previously specified in the input script for use with the read_isurf command and (optionally) to perform ablation during a simulation. It stores the grid corner point values for each grid cell.
The output file is written in the same binary format as the read_isurf command reads in.
Because implicit-surface corner point values are not stored in restart files, a write_isurf followed by a later read_isurf is the way to save and restore the state of an ablating surface. To reproduce the saved state exactly, write the file with precision double (the default) so partially decremented corner point values are not truncated to integers, and read it back with the read_isurf precision double and push no options. The default read_isurf push yes re-normalizes corner point values as they are read and so will not reproduce a mid-ablation state exactly; see the push keyword discussion on the read_isurf doc page.
Restrictions: none
Related commands:
Default:
The optional keyword default is precision double.