1*59599516SKenneth E. Jansen subroutine input() 2*59599516SKenneth E. Jansenc 3*59599516SKenneth E. Jansenc---------------------------------------------------------------------- 4*59599516SKenneth E. Jansenc This routine inputs all the necessary data, allocates required array 5*59599516SKenneth E. Jansenc storage, and sets up the appropriate parameters for the processing. 6*59599516SKenneth E. Jansenc 7*59599516SKenneth E. Jansenc 8*59599516SKenneth E. Jansenc Farzin Shakib, Winter 1987. 9*59599516SKenneth E. Jansenc Zdenek Johan, Winter 1991. (Fortran 90) 10*59599516SKenneth E. Jansenc---------------------------------------------------------------------- 11*59599516SKenneth E. Jansenc 12*59599516SKenneth E. Jansen include "commonM2NFixBnd.h" 13*59599516SKenneth E. Jansen include "mpif.h" 14*59599516SKenneth E. Jansen 15*59599516SKenneth E. Jansen external endata 16*59599516SKenneth E. Jansen 17*59599516SKenneth E. Jansen integer, allocatable :: nsons(:) 18*59599516SKenneth E. Jansenc 19*59599516SKenneth E. Jansen character*8 date 20*59599516SKenneth E. Jansen character*80 card 21*59599516SKenneth E. Jansen 22*59599516SKenneth E. Jansenc assigned in phasta.cc 23*59599516SKenneth E. Jansenc numpe=npe 24*59599516SKenneth E. Jansenc myrank=mrank 25*59599516SKenneth E. Jansen 26*59599516SKenneth E. Jansen call MPI_BARRIER(MPI_COMM_WORLD, ierr) 27*59599516SKenneth E. Jansen 28*59599516SKenneth E. Jansen epsM = sqrt(epsilon(one)) 29*59599516SKenneth E. Jansenc 30*59599516SKenneth E. Jansenc.... read in and block all data 31*59599516SKenneth E. Jansenc 32*59599516SKenneth E. Jansen call readnblk() 33*59599516SKenneth E. Jansenc 34*59599516SKenneth E. Jansenc....return 35*59599516SKenneth E. Jansenc 36*59599516SKenneth E. Jansen return 37*59599516SKenneth E. Jansen end 38