Welcome! 

[Home] [PMWIN] [3D Visualization] [ASMWIN] [Downloads] [Links]

by Wen-Hsing Chiang Last Update: December 05, 2003
This page describes the internal cell-by-cell binary files of PMWIN.
A complete list of these files can be found in Appendix A of the manual
of PMWIN.

1. GRID GEOMETRY AND LOCATION:
The grid geometry and location are saved in the Grid Specification File
with the extension GRD. The format of this file is given below:

File Format
1. Data:  NROW   NCOL
2. Data:  X      Y     ANGLE
3. Data:  DELR(NCOL)
4. Data:  DELC(NROW)

Explanation of Fields Used in Input Instructions
NROW    is the number of model rows.
NCOL    is the number of model columns.
X       is the x-coordinate of the top-left corner of the model grid.
Y       is the y-coordinate of the top-left corner of the model grid.
ANGLE   is the rotation angle expressed in degrees and measured
        countercolckwise from the positive x- axis.
DELR    is the cell width along rows. Read one value for each
        of the NCOL columns. This is a single array with one value for
        each column.
DELC    is the cell width along columns. Read one value for each of the
        NROW rows. This is a single array with one value for each row.  

2. TIME-INDEPENDENT MODEL DATA FILES

Time-indenpendent model data (such as hydraulich conductivity or porosity)
are saved in time-independent model data file, which has a header of
1000*4 bytes. The header is not used here. Following the header, PMWIN
uses REAL (4 bytes) to save cell-by-cell model data seguentially.
The following example illustrates a model data file with NJ columns,
NI rows and 2 layers.

  HEADER (1000*4 bytes)

  f(1,1)  f(1,2) ... f(1,J) .... f(1,NJ)   ---
  f(2,1)  f(2,2) ... f(2,J) .... f(2,NJ)     |
    .       .           .           .        |     2D matrix for layer 1
    .       .           .           .        |     NI*NJ*4 bytes
  f(I,1)  f(I,2) ... f(I,J) .... f(I,NJ)     |
    .       .           .           .        |
    .       .           .           .        |
  f(NI,1) f(NI,2)... f(NI,J) ....f(NI,NJ)  ---

  f(1,1)  f(1,2) ... f(1,J) .... f(1,NJ)   ---
  f(2,1)  f(2,2) ... f(2,J) .... f(2,NJ)     |
    .       .           .           .        |     2D matrix for layer 2
    .       .           .           .        |     NI*NJ*4 bytes
  f(I,1)  f(I,2) ... f(I,J) .... f(I,NJ)     |
    .       .           .           .        |
    .       .           .           .        |
  f(NI,1) f(NI,2)... f(NI,J) ....f(NI,NJ)  ---


The file numbers and extensions of time-independent data files are
listed below:
3        TOP    top elevation of layers
4        BOT    Bottom elevation of layers
5        IBD    IBOUND array
6        POR    Effective Porosity
7        STO    specific storage
8        CON    horizontal hydraulic conductivity
9        LEA    vertical hydraulic conductivity
10       HEA    starting head
40       TIC    ICBUND Array (MT3D)
41       TSC    Starting Concentration (MT3D)
42       TOB    Observation Points (MT3D)
51       LKN    vertical leakance
52       DWA    Wetting threshold (MODFLOW BCF2)
53       TAL    longitudinal dispersivity (MT3D)
60       HTC    Transmissivity
61       SCC    Confined storage coefficient
62       YLD    Specific Yield
70       WAL    Direction (horizontal flow barrier)
71       WAC    Kf/Wall_thickness (horizontal flow barrier)
81       WBL    Subregions for the water budget calculation
83       _83    Recycle
85       C85    Location of reservoirs (PMWIN 4.1 or higher)
86       C86    Bottom Elevation of the reservoir (PMWIN 4.1 or higher)
87       C87    Bed conductivity of the reservoir (PMWIN 4.1 or higher)
88       C88    Bed Thickness of the reservoir (PMWIN 4.1 or higher)
89       C89    Layer indicator of the reservoir (PMWIN 4.1 or higher)
91       C91    Bulk density of the soil matrix (PMWIN 4.1 or higher)
92       C92    First sorption constant (PMWIN 4.1 or higher)
93       C93    Second sorption constant (PMWIN 4.1 or higher)
94       C94    First-order rate constant for the dissolved phase
                (PMWIN 4.1 or higher)
95       C95    First-order rate constant for the sorbed phase
                (PMWIN 4.1 or higher)
97       C97    digitizer (PMWIN 4.1 or higher)

** Do not use the file numbers 98, 99 and 100 for any array, these numbers
are reserved **

3. TIME-DEPENDENT MODEL DATA FILES

Each time-dependent file has a corresponding record of 1000*2 bytes
(INTEGER) in the file filename.RCD. Where filename is the name of the
model. Each integer of this record corresponds to a stress period
(max. 1000 stress periods are allowed). 
The beginning position (in bytes) of a record is
   (FN-1) * 2000 + 1
Where FN is the file number listed below. If the i-th integer has the
value j, the data of the j-th period will be used for the i-th stress period.

A time-dependent model data file has a header of 1000 * 4 bytes (REALs).

The value S in the n-th REAL indicates that the data of S-th stress
period is saved in the n-th position of the file.

In the following example, a model has 2 layers and 8 stress periods.
The data of the 1st, 3rd and 5th stress periods have been specified and
saved in the position 1, 3 and 2, respectively. In this case, the
specified data for the stress period 5 will be used for the stress
periods 6, 7 and 8, too. Similarly, the data of the first
stress period will be used for the 1st and 2nd periods. The data of the
stress period 3 will be used for the 3rd and 4th periods. This data
structure allows a user to jump between stress periods in the Data
Editor without having to save data of all stress periods.


RECORD in the RCD file:
Position (=Stress period) -> 1  2  3  4  5  6  7  8  ....
Values:                      1  1  3  3  5  5  5  5  ....
                                _________|________|
                                |
                                |
HEADER (1000*4 bytes):       1  5  3  0  0  0  0  0  0 ....
                                |
                                --------------------------------|
           position 1:       2D-matrix for layer 1, period 1    |
   (begin: 4001 bytes)       2D-matrix for layer 2, period 1    |
                                                                |
           position 2:       2D-matrix for layer 1, period 5 <--|
  (4001+(2-1)*NI*NJ*4)       2D-matrix for layer 2, period 5 <--|

           position 3:       2D-matrix for layer 1, period 3
  (4001+(3-1)*NI*NJ*4)       2D-matrix for layer 2, period 3