I want to make a code for linear interpolation in fortran. From an input file there are 2 columns. I want to read the data first then apply linear interpolation to the data.

The input file includes below information.

3000 0.000 3050 0.019 3100 0.068 3150 0.167 3200 0.278 3250 0.398 3300 0.522 3350 0.636 3400 0.735 3450 0.813 3500 0.885 3550 0.940 3600 0.980 3650 1.000 3700 1.000 3750 0.974 3800 0.918 3850 0.802 3900 0.590 3950 0.355 4000 0.194 4050 0.107 4100 0.046 4150 0.003 4200 0.000

please make a code in FORTRAN.

I need to be done with fortran.