program test

use COESA_module
use iso_fortran_env, wp => real64

implicit none

real(wp) :: h,d
integer :: i

do i = 1,101
    h=(i-1)*1000
    d = COESA_density(h)
    write(*,*) h/1000,d
end do

end program test