如何在fortran中找到给定范围内的孪生素数对

I want to make a program in fortran that lists all the pairs of twin prime numbers in a range, n1 to n2, given by the user. I did this, but it doesn't work and I don't know how to fix it:I want to make a program in fortran that lists