output_daily_final.f90

program output_daily_final

use parametersmod
use newsplinemod,   only : newspline

implicit none

character(6) :: header
real(sp),    dimension(53) :: temp_F, temp_C
integer(i4), dimension(53) :: date

real(sp),    dimension(36) :: temp
integer(i4), dimension(36) :: nk
real(sp),    dimension(2)  :: bcond
real(sp),    dimension(:), allocatable :: daydata

integer(i4)  :: i
integer(i4)  :: j
integer(i4)  :: n
integer(i4)  :: summ

!------
! Read temperature data file
open(10, file = "captaincook_monthly_temp_NOAA.txt", status = 'old')

read(10,*) header

do i = 1, 53

  read(10,*) date(i), temp_F(i), temp_C(i)

end do

!---
! Define number of days in each month
nk = [31,28,31,30,31,30,31,31,30,31,30,31,  &
      31,29,31,30,31,30,31,31,30,31,30,31,  &
      31,28,31,30,31,30,31,31,30,31,30,31]

!---
! Initiate temperature array. Copy first and last month to start and end.
summ = sum(nk)

allocate(daydata(summ))

daydata = 0.0

temp(1:36) = temp_C(1:36)

! Apply circular boundary conditions by copy last interval to the beginning and vice versa
bcond(1) = temp_C(36)
bcond(2) = temp_C(1)

!---
! Interpolation without any limit options
call newspline(temp,nk,bcond,daydata)

!---
! Print newspline output
n = 1

do i = 1, 36

  do j = 1, nk(i)

    write(*,*) temp(i), daydata(n)

    n = n + 1

  end do

end do


end program output_daily_final