Skip to content
Snippets Groups Projects
Commit 8e6f9bec authored by kamischi's avatar kamischi
Browse files

fix compilation of eels and eels-boson.

also refactor implicit function zero in surlos.
parent 3f46f6fc
No related branches found
No related tags found
No related merge requests found
module mod_doeels
module doeels_mod
contains
subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, comment, comment_size, &
layers, neps, nper, name, name_size, thick, epsinf, nos, osc, osc_size,&
......@@ -235,4 +235,4 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, comment, comment_size,
120 format(40x, f10.2, d12.4, f10.2)
130 format(2x, f5.1, 3x, f11.3, d14.5)
end subroutine doeels
end module mod_doeels
end module doeels_mod
......@@ -9,6 +9,7 @@ program eels_boson
! ******************************************************************
use get_commandline_options_mod
use doeels_mod
implicit none
......
......@@ -8,6 +8,8 @@ program eels
! * *
! ******************************************************************
use doeels_mod
implicit none
double precision :: e0, theta, phia, phib, wmin, wmax, dw
......
module surlos_mod
contains
logical function zero(z)
double complex, intent(in) :: z
zero = (dble(z) == 0.0d0) .and. (dimag(z) == 0.0d0)
end function zero
double precision function surlos(dk, eps, thick, layers, nper)
! ******************************************************************
......@@ -22,7 +28,7 @@ double precision function surlos(dk, eps, thick, layers, nper)
integer, intent(in) :: layers, nper
integer :: lmax
logical :: static, zero, skip
logical :: static, skip
integer :: lstart, n
double precision, allocatable :: arg(:)
double precision :: argmin, argmax, cn, cnm1, epsmac, sn, snm1, tn
......@@ -33,7 +39,7 @@ double precision function surlos(dk, eps, thick, layers, nper)
epsmac = epsmac / 2
argmax = 0.5d0 * dlog(2 / epsmac)
zero(z) = (dble(z) == 0.0d0) .and. (dimag(z) == 0.0d0)
! zero(z) = (dble(z) == 0.0d0) .and. (dimag(z) == 0.0d0)
! write (*,*) 'surlos:'
! write (*,*) 'thick: ', size(thick)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment