read_infile Subroutine

public subroutine read_infile(filename)

Open and read chease file named "filename". Populates module level variables with the output.

Arguments

Type IntentOptional Attributes Name
character(len=*), intent(in) :: filename

Contents

Source Code


Source Code

subroutine read_infile(filename)
  implicit none
  character(len=*), intent(in) :: filename
  integer :: infile
  open(newunit = infile, file = filename)
  read(infile, *) ; read(infile, *) npsi_chease ; read(infile, *) ; read(infile, *) nchi_chease
  if (debug) write(*, *) npsi_chease, "<---npsi_chease"
  if (debug) write(*, *) nchi_chease, "<---nchi_chease"

  read(infile, *) ; read(infile, *) r0exp_chease
  read(infile, *) ; read(infile, *) b0exp_chease

  allocate( psi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) psi_chease
  allocate( chi_chease(nchi_chease)) ; read(infile, *) ; read(infile, *) chi_chease
  allocate( rgeom_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) rgeom_chease
  allocate( ageom_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) ageom_chease
  allocate( q_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) q_chease
  allocate( dqdpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dqdpsi_chease
  allocate( d2qdpsi2_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) d2qdpsi2_chease
  allocate( p_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) p_chease
  allocate( dpdpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dpdpsi_chease
  allocate( f_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) f_chease
  allocate( fdfdpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) fdfdpsi_chease
  allocate( v_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) v_chease
  allocate( rho_t_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) rho_t_chease
  allocate( shear_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) shear_chease
  allocate( dsheardpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dsheardpsi_chease
  allocate( kappa_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) kappa_chease
  allocate( delta_lower_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) delta_lower_chease
  allocate( delta_upper_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) delta_upper_chease
  allocate( dvdpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dvdpsi_chease
  allocate( dpsidrhotor_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dpsidrhotor_chease
  allocate( gdpsi_av_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) gdpsi_av_chease
  allocate( radius_av_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) radius_av_chease
  allocate( r_av_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) r_av_chease
  allocate( te_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) te_chease
  allocate( dtedpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dtedpsi_chease
  allocate( ne_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) ne_chease
  allocate( dnedpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dnedpsi_chease
  allocate( ti_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) ti_chease
  allocate( dtidpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dtidpsi_chease
  allocate( ni_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) ni_chease
  allocate( dnidpsi_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) dnidpsi_chease
  allocate( zeff_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) zeff_chease
  allocate( signeo_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) signeo_chease
  allocate( jbsbav_chease(npsi_chease)) ; read(infile, *) ; read(infile, *) jbsbav_chease

  allocate( g11_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) g11_chease
  allocate( g12_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) g12_chease
  allocate( g22_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) g22_chease
  allocate( g33_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) g33_chease
  allocate( b_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) b_chease
  allocate( dbdpsi_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) dbdpsi_chease
  allocate( dbdchi_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) dbdchi_chease
  allocate( dpsidr_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) dpsidr_chease
  allocate( dpsidz_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) dpsidz_chease
  allocate( dchidr_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) dchidr_chease
  allocate( dchidz_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) dchidz_chease
  allocate( jacobian_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) jacobian_chease
  allocate( r_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) r_chease
  allocate( z_chease(npsi_chease, nchi_chease)) ; read(infile, *) ; read(infile, *) z_chease

  close(unit = infile)
end subroutine read_infile