subroutine field_read use input_file_module use maximum_data_module use topography_data_module implicit none integer :: ith = 0 !none |counter character (len=80) :: titldum = ""! |title of file character (len=80) :: header = "" ! |header of file integer :: eof = 0 ! |end of file integer :: imax = 0 !none |determine max number for array (imax) and total number in file logical :: i_exist !none |check to determine if file exists eof = 0 imax = 0 !! read all data from topo.dat inquire (file=in_hyd%field_fld, exist=i_exist) if (.not. i_exist .or. in_hyd%field_fld == "null") then allocate (field_db(0:0)) else do open (107,file=in_hyd%field_fld) read (107,*,iostat=eof) titldum if (eof < 0) exit read (107,*,iostat=eof) header if (eof < 0) exit do while (eof == 0) read (107,*,iostat=eof) titldum if (eof < 0) exit imax = imax + 1 end do db_mx%field = imax allocate (field_db(0:imax)) rewind (107) read (107,*,iostat=eof) titldum if (eof < 0) exit read (107,*,iostat=eof) header if (eof < 0) exit do ith = 1, db_mx%field read (107,*,iostat=eof) titldum if (eof < 0) exit backspace (107) read (107,*,iostat=eof) field_db(ith) if (eof < 0) exit end do exit enddo endif close (107) return end subroutine field_read