subroutine cli_smeas use climate_module use input_file_module use time_module use maximum_data_module implicit none integer :: i = 0 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 integer :: iyr = 0 !none |number of years logical :: i_exist !none |check to determine if file exists integer :: istep = 0 ! | integer :: iyr_prev = 0 !none |previous year integer :: iyrs = 0 ! | real :: solrad = 0. !MJ/m^2 |solar radiation for the day in HRU eof = 0 imax = 0 !! read all measured daily solar radiation data inquire (file=in_cli%slr_cli, exist=i_exist) if (.not. i_exist .or. in_cli%slr_cli == "null") then allocate (slr(0:0)) allocate (slr_n(0)) else do open (107,file=in_cli%slr_cli) 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 allocate (slr(0:imax)) allocate (slr_n(imax)) rewind (107) read (107,*,iostat=eof) titldum if (eof < 0) exit read (107,*,iostat=eof) header if (eof < 0) exit do i = 1, imax read (107,*, iostat=eof) slr_n(i) if (eof < 0) exit end do rewind (107) read (107,*,iostat=eof) titldum if (eof < 0) exit read (107,*,iostat=eof) header if (eof < 0) exit do i = 1, imax read (107,*,iostat = eof) slr(i)%filename if (eof < 0) exit !!!!!weather path code if (in_path_slr%slr == "null") then open (108,file = slr(i)%filename) else open (108,file = TRIM(ADJUSTL(in_path_slr%slr))//slr(i)%filename) endif !!!!!weather path code read (108,*,iostat=eof) titldum if (eof < 0) exit read (108,*,iostat=eof) header if (eof < 0) exit read (108,*,iostat=eof) slr(i)%nbyr, slr(i)%tstep, slr(i)%lat, slr(i)%long, & slr(i)%elev if (eof < 0) exit ! the precip time step has to be the same as time%step allocate (slr(i)%ts(366,slr(i)%nbyr), source = 0.) ! read and save start jd and yr read (108,*,iostat=eof) iyr, istep if (eof < 0) exit slr(i)%start_day = istep slr(i)%start_yr = iyr backspace (108) if (iyr > time%yrc) then slr(i)%yrs_start = iyr - time%yrc else ! read and store entire year slr(i)%yrs_start = 0 end if do read (108,*,iostat=eof) iyr, istep, solrad if (eof < 0) exit if (iyr >= time%yrc .and. istep >= time%day_start) exit end do backspace (108) iyr_prev = iyr iyrs = 1 do read (108,*,iostat=eof) iyr, istep, slr(i)%ts(istep,iyrs) if (eof < 0) exit if (istep == 365 .or. istep == 366) then read (108,*,iostat=eof) iyr, istep if (eof < 0) exit backspace (108) if (iyr /= iyr_prev) then iyr_prev = iyr iyrs = iyrs + 1 end if end if end do close (108) ! save end jd and year slr(i)%end_day = istep slr(i)%end_yr = iyr end do close (107) exit end do endif db_mx%slrfiles = imax return end subroutine cli_smeas