comments with suggestions and white space

source/f90/doeels.f90

@@ -113,10 +113,15 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, comment, comment_size,
   nout = 1 + nw / 20
   ener = 8065 * e0
   psia = phia / 180 * pi
+! *** kms
   psii = theta / 180 * pi
   cospsi = dcos(psii)
   sinpsi = dsin(psii)
   tanpsi = dtan(psii)
+! *** f2023 replacement suggestion
+!  cospsi = dcosd(phia)
+!  sinpsi = dsind(phia)
+!  tanpsi = dtand(phia)
   prefac = dsqrt(255500 / e0)/(137 * cospsi)
   facru = psia / cospsi * dsqrt(0.2624664d0 * e0)
   elleps = (1.0d0 - phia / phib) * (1.0d0 + phia / phib)

tests/seteps/setepsdriver-f90-new.f90

@@ -63,11 +63,11 @@ program setepsdriver
   wn = wnmin
   allocate (eps(layers))
 
-  open(unit = 12, file = 'build-f90/epsout')
+  open(unit = 12, file = 'build-f90-new/epsout')
   do i = 0, noPoints
-      call seteps(layers, nos, oscType, osc, epsinf, wn, eps, layers)
-      write (12, '(30g15.7)') wn,  (eps(j), j = 1, layers)
-      wn = wn + dwn
+    call seteps(layers, nos, oscType, osc, epsinf, wn, eps, layers)
+    write (12, '(30g15.7)') wn,  (eps(j), j = 1, layers)
+    wn = wn + dwn
   end do
   close (unit = 12)