diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2021-08-04 15:18:08 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2021-08-04 15:18:08 -0400 |
| commit | 70fd14e98bc476a84aef352eaf063cfcba44ab83 (patch) | |
| tree | 2737f5c1c21c76824f49c6a1b9c7061952b646f5 /Opaccouls.f | |
| parent | c396a4bc7a3282de9d2191de0e1076fec9df0783 (diff) | |
| download | moog-70fd14e98bc476a84aef352eaf063cfcba44ab83.tar.gz | |
Use modern equality operators
Diffstat (limited to 'Opaccouls.f')
| -rwxr-xr-x | Opaccouls.f | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/Opaccouls.f b/Opaccouls.f index 3c4f59d..2410af8 100755 --- a/Opaccouls.f +++ b/Opaccouls.f @@ -19,12 +19,12 @@ c******************************************************************************` data c/ -2.268d30, 4.077d28, 1.035d28, 4.593d27, 2.371d27, . 1.229d27/ - if (freq .lt. z*z*3.28805d15/dfloat(n*n)) then + if (freq < z*z*3.28805d15/dfloat(n*n)) then coulx = 0. else coulx = 2.815d29/freq**3/dfloat(n**5)*z**4 - if (n .gt. 6) return - if (n .eq. 1) then + if (n > 6) return + if (n == 1) then coulx = coulx*coulbf1s(freq,z) else coulx = coulx*(a(n)+(b(n)+c(n)*(z*z/freq))*(z*z/freq)) @@ -64,7 +64,7 @@ c******************************************************************************` . 0.2701,0.2648,0.2595,0.2544,0.2493,0.2443,0.2394,0.2345,0.2298, . 0.2251,0.2205,0.2160,0.2115,0.2072,0.2029,0.1987/ - if (freq/z**2 .lt. 3.28805d15) then + if (freq/z**2 < 3.28805d15) then coulbf1s = 0. else elog = log10(freq/z**2/3.28805d15) |