xref: /phasta/phSolver/common/d2wall.c (revision 8746ab438bbda91291f8cdd62b94f8385f2d26f1)

#include <stdlib.h>
#include <FCMangle.h>
#include <new_interface.h>
#include <stdio.h>
#include <string.h> /*memset*/
#include <assert.h>
#include "common_c.h"
#include "phastaIO.h"
#include "phIO.h"
#include "phString.h"
#include "syncio.h"
#include "posixio.h"
#include "streamio.h"
#include "setsyncioparam.h"
#include "tmrc.h"

void
read_d2wall(  int* pid,
              int* numnp,
              double* array1,
              int* foundd2wall ) {
    int isize, nitems;
    int iarray[10];
    int j;
    double iotime = 0;
    for ( j = 0; j < 10; j++) {
       /*Initialize iarray to 0 so that we can assess the result of readheader*/
       iarray[j] = 0;
    }

    int nfields;
    int numparts;
    int irank;
    int nprocs;

    /*  Retrieve and compute the parameters required for SyncIO */
    numparts = workfc.numpe;
    irank = *pid; /* workfc.myrank; */
    nprocs = workfc.numpe;

    /* Calculate number of parts each proc deal with and where it start and end ... */
    int nppp = numparts/nprocs;/* nppp : Number of parts per proc ... */
    assert(nppp==1);
    int startpart = irank * nppp +1;/* Part id from which I (myrank) start ... */
    int endpart = startpart + nppp - 1;/* Part id to which I (myrank) end ... */

    phio_fp handle;
    char filename[255],path[255];
    memset((void*)filename,0,255);
    *foundd2wall = 0;
    /* First we try to read dwal from the restart files. */

    iotime = TMRC();
    if( outpar.input_mode == -1 )
      streamio_setup_read(&handle, streamio_get_gr());
    else if( outpar.input_mode == 0 )
      posixio_setup(&handle, 'r');
    else if( outpar.input_mode > 0 )
      syncio_setup_read(outpar.nsynciofiles, &handle);
    phio_constructName(handle,"restart",filename);
    phstr_appendInt(filename, timdat.lstep);
    phstr_appendStr(filename, ".");
    phio_openfile(filename, handle);

    int i;
    for ( i = 0; i < nppp; i++) { /*This loop is useful only if several parts per processor*/
      nitems = 2;
      phio_readheader(handle, "dwal", (void*)iarray, &nitems, "double", phasta_iotype);

      if (iarray[0] == (*numnp)) {
        if (irank==0) {
          printf("d2wall field found in %s\n",filename);
        }
        *foundd2wall = 1;
        isize = (*numnp);
        phio_readdatablock(handle, "dwal", (void*)(array1), &isize, "double", phasta_iotype );
      }
      else { /*d2wall fields was not found in the restart file*/
        *foundd2wall = 0;
        if (irank==0) {
          printf("d2wall field not found in %s - trying d2wall files now\n",filename);
        }
      }
    }
    phio_closefile(handle);
    iotime = TMRC() - iotime;
    if (irank==0) {
      printf("time to read d2wall (seconds) %f\n",iotime);
      printf("\n");
    }
}