included from readmajorcube.f

INTERP.FOR (text/plain), 7.26 KB, created by Tom Fogal on 2005-08-19 11:26:42 UTC

c *********************************************************************
c
c Purpose :: Define the named COMMON blocks for interpolated fields.
c
c Programmer :: Dr. Douglas J. Larson
c Creation Date : Thursday, March 20, 1997
c
c The general idea behind the layout of this interpolation gridding
c scheme can be expressed in a C like way:
c
c struct GridParameters {
c        char filename[80];    /* The name of the data file on disk */
c        enum boolean TimeDep=FALSE; /* TRUE-Time dependent */
c        int dims;             /* Is this 1D, 2D, or 3D ? */
c        int coorsys;          /* 1-Cartesian, 2-Cylindrical, 3-Spherical */
c        int rangetotal[3];    /* Number of boxes in the range */
c        double voxelsize;     /* Volume of cube cell */
c        double time[3];       /* 1:min, 2:max, 3:current */
c        double range[6];      /* Range start and end for each axes */
c        double rangedelta[3]; /* Range delta for each axes */
c        struct GridParameters *nextgrid_ptr;
c };
c
c where we have an array of GridParameter structures for each set of
c precomputed data. Because the computer has limited memory resources,
c we have chosen to keep only the parameters of the available datasets 
c in memory and have a separate set of buffers for the actual data. 
c This entails creating a way of keeping track of what has been read 
c into RAM and what is on disk. The solution is a another structure for 
c bookkeeping:
c
c struct GridBookkeeping {
c        int filecount;         /* Number of data files available */
c        int LastMajorCubeUsed;
c        int MajorCubeNumber;
c        struct GridParameters *GridParameter_ptr;
c };
c 
c struct GridBuffer {
c        double Field_x[MAXGRID];
c        double Field_y[MAXGRID];
c        double Field_z[MAXGRID];
c };
c
c struct FieldGrid {
c        struct GridBookkeeping *GridBookkeeping_ptr;
c        int ActiveBuffer[MAXBUFFERS];
c        int BufferCount[MAXBUFFERS];  /* Number of buffer hits */
c        struct GridBuffer[MAXBUFFERS];
c };
c
c Anyway, that is what ran through my head as I wrote this originally.
c
c struct FieldGrid *ElectricField_ptr;
c struct FieldGrid *EPotential_ptr;
c struct FieldGrid *MagneticField_ptr;
c
c Interpolation inputs & variables
c Key:  BI  B-field Interpolation
c       EI  E-field Interpolation
c       PI  Electric Potential Interpolation
c
      integer maxIsets  ! Maximum number of Interpolation data sets
      integer maxvalues ! ReadMajorCube uses this
      integer maxactive ! Limit of data sets read in at one time
      integer maxxcount,maxycount,maxzcount
cc      parameter (maxxcount=130,maxycount=130,maxzcount=100)
cc      parameter (maxxcount=100,maxycount=100,maxzcount=100)
      parameter (maxxcount=75,maxycount=75,maxzcount=75)
      parameter (maxactive=2,maxIsets=50)
      parameter (maxvalues=maxxcount*maxycount*maxzcount)

      logical BI_TimeDep  ! TRUE-Time dependence ON, FALSE-static
      logical EI_TimeDep  ! TRUE-Time dependence ON, FALSE-static
      logical PI_TimeDep  ! TRUE-Time dependence ON, FALSE-static

      integer AI_dims(maxIsets)    ! Is this 1D, 2D, or 3D ?
      integer AI_coorsys(maxIsets) ! 1-Cartesian, 2-Cylindrical, 3-Spherical
      integer AI_filecount         ! Number of A vector potential files
      integer AI_rangetotal(3,maxIsets) ! Number of boxes in the range

      integer BI_dims(maxIsets)    ! Is this 1D, 2D, or 3D ?
      integer BI_coorsys(maxIsets) ! 1-Cartesian, 2-Cylindrical, 3-Spherical
      integer BI_filecount         ! Number of B-field files
      integer BI_rangetotal(3,maxIsets) ! Number of boxes in the range
      integer BI_read(maxactive,2) ! list of data sets now in RAM, history
      integer BI_LastMajorCubeUsed
      integer BI_MajorCubeNumber

      integer EI_dims(maxIsets)    ! Is this 1D, 2D, or 3D ?
      integer EI_coorsys(maxIsets) ! 1-Cartesian, 2-Cylindrical, 3-Spherical
      integer EI_filecount         ! Number of E-field files
      integer EI_rangetotal(3,maxIsets) ! Number of boxes in the range
      integer EI_read(maxactive,2) ! list of data sets now in RAM, history
      integer EI_LastMajorCubeUsed
      integer EI_MajorCubeNumber

      integer PI_InterpPotSegment  ! Segment number where potential defined
      integer PI_dims(maxIsets)    ! Is this 1D, 2D, or 3D ?
      integer PI_coorsys(maxIsets) ! 1-Cartesian, 2-Cylindrical, 3-Spherical
      integer PI_filecount         ! Number of Electric potential files
      integer PI_rangetotal(3,maxIsets) ! Number of boxes in the range
      integer PI_read(maxactive,2) ! list of data sets now in RAM, history
      integer PI_LastMajorCubeUsed
      integer PI_MajorCubeNumber

      character*128 BI_filenames(maxIsets) ! Array of file names, B
      character*128 EI_filenames(maxIsets) ! Array of file names, E
      character*128 PI_filenames(maxIsets) ! Array of file names, Potential

      DOUBLE PRECISION BI_voxelsize(maxIsets) ! Volume of cube cell
      DOUBLE PRECISION BI_time(3,maxIsets)    ! 1:min, 2:max, 3:current
      DOUBLE PRECISION BI_range(6,maxIsets)   ! Range start and end for each axes
      DOUBLE PRECISION BI_rangedelta(3,maxIsets) ! Range delta for each axes
      DOUBLE PRECISION BI_bx(maxxcount,maxycount,maxzcount,maxactive)
      DOUBLE PRECISION BI_by(maxxcount,maxycount,maxzcount,maxactive)
      DOUBLE PRECISION BI_bz(maxxcount,maxycount,maxzcount,maxactive)

      DOUBLE PRECISION EI_voxelsize(maxIsets) ! Volume of cube cell
      DOUBLE PRECISION EI_time(3,maxIsets)    ! 1:min, 2:max, 3:current
      DOUBLE PRECISION EI_range(6,maxIsets)   ! Range start and end for each axes
      DOUBLE PRECISION EI_rangedelta(3,maxIsets) ! Range delta for each axes
      DOUBLE PRECISION EI_ex(maxxcount,maxycount,maxzcount,maxactive)
      DOUBLE PRECISION EI_ey(maxxcount,maxycount,maxzcount,maxactive)
      DOUBLE PRECISION EI_ez(maxxcount,maxycount,maxzcount,maxactive)

      DOUBLE PRECISION PI_voxelsize(maxIsets) ! Volume of cube cell
      DOUBLE PRECISION PI_time(3,maxIsets)    ! 1:min, 2:max, 3:current
      DOUBLE PRECISION PI_range(6,maxIsets)   ! Range start and end for each axes
      DOUBLE PRECISION PI_rangedelta(3,maxIsets) ! Range delta for each axes
      DOUBLE PRECISION PI_data(maxxcount,maxycount,maxzcount,maxactive)
c ======================================================================
      common/AI_integer/AI_dims,AI_coorsys,AI_filecount,AI_rangetotal
c
      common/BI_chars/BI_filenames
      common/BI_logicals/BI_TimeDep
      common/BI_control/BI_read,BI_LastMajorCubeUsed,BI_MajorCubeNumber
      common/BI_integer/BI_dims,BI_coorsys,BI_filecount,BI_rangetotal
      common/BI_reals/BI_voxelsize,BI_time,BI_range,BI_rangedelta
      common/BI_values/BI_bx,BI_by,BI_bz
c
      common/EI_chars/EI_filenames
      common/EI_logicals/EI_TimeDep
      common/EI_control/EI_read,EI_LastMajorCubeUsed,EI_MajorCubeNumber
      common/EI_integer/EI_dims,EI_coorsys,EI_filecount,EI_rangetotal
      common/EI_reals/EI_voxelsize,EI_time,EI_range,EI_rangedelta
      common/EI_values/EI_ex,EI_ey,EI_ez
c
      common/PI_chars/PI_filenames
      common/PI_logicals/PI_TimeDep
      common/PI_control/PI_read,PI_LastMajorCubeUsed,PI_MajorCubeNumber,
     %                  PI_InterpPotSegment
      common/PI_integer/PI_dims,PI_coorsys,PI_filecount,PI_rangetotal
      common/PI_reals/PI_voxelsize,PI_time,PI_range,PI_rangedelta
      common/PI_values/PI_data
c