subroutine om_water_init use basin_module use input_file_module use maximum_data_module use channel_data_module use hydrograph_module use sd_channel_module use constituent_mass_module implicit none character (len=80) :: titldum = "" ! |title of file character (len=80) :: header = "" ! |header of file integer :: eof = 0 ! |end of file integer :: imax = 0 !units |description logical :: i_exist ! |check to determine if file exists integer :: ichi = 0 !none |counter eof = 0 imax = 0 inquire (file=in_init%om_water, exist=i_exist) if (.not. i_exist .or. in_init%om_water == "null") then allocate (om_init_water(0:0)) allocate (om_init_name(0:0)) else do open (105,file=in_init%om_water) read (105,*,iostat=eof) titldum if (eof < 0) exit read (105,*,iostat=eof) header if (eof < 0) exit do while (eof == 0) read (105,*,iostat=eof) titldum if (eof < 0) exit imax = imax + 1 end do db_mx%om_water_init = imax allocate (om_init_water(0:imax)) allocate (om_init_name(0:imax)) rewind (105) read (105,*,iostat=eof) titldum if (eof < 0) exit read (105,*,iostat=eof) header if (eof < 0) exit do ichi = 1, db_mx%om_water_init read (105,*,iostat=eof) titldum if (eof < 0) exit backspace (105) read (105,*,iostat=eof) om_init_name(ichi), om_init_water(ichi) if (eof < 0) exit end do close (105) exit enddo endif return end subroutine om_water_init