Removed bug: G = 0 was not defined in the subroutine called input_field

src/field_io_mod

src/io/field_io_mod.f

@@ -62,6 +62,7 @@ subroutine input_field(field,field_unit)
  real(long) :: swap(N_monomer)
  type(version_type) :: version
 
+ G = 0
  ! Input file format version (e.g., `format 1 0')
  call input_version(version, field_unit)
 
@@ -80,6 +81,7 @@ subroutine input_field(field,field_unit)
  end if
  end do
 
+
  end subroutine input_field
  !==============================================================
 
@@ -122,6 +124,7 @@ subroutine output_field(field,field_unit,group_name,N_basis_out)
  call output(N_monomer,'N_monomer',o=field_unit)
  call output(N_output,'N_star',o=field_unit)
 
+
  fmt = '('//trim(int_string(N_monomer))//'ES20.12,4X,'
  fmt = trim(fmt)//trim(int_string(dim))//'I4,I6'//')'
  do i = 1, N_output

src/pscf_pd.fp.f

@@ -957,7 +957,7 @@ program pscf_pd
  open(unit=field_unit,file=trim(input_file),status='old')
  call input_field(rho,field_unit)
  close(field_unit) 
-  
+
  do alpha=1,N_monomer
  do i=1,N_star
  omega_basis(alpha,i) = sum(chi(:,alpha)*rho(:,i))