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 /Lineabund.f | |
| parent | c396a4bc7a3282de9d2191de0e1076fec9df0783 (diff) | |
| download | moog-70fd14e98bc476a84aef352eaf063cfcba44ab83.tar.gz | |
Use modern equality operators
Diffstat (limited to 'Lineabund.f')
| -rwxr-xr-x | Lineabund.f | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/Lineabund.f b/Lineabund.f index 9c709ea..c7d807d 100755 --- a/Lineabund.f +++ b/Lineabund.f @@ -23,7 +23,7 @@ c*****find the c-o-g gf that matches the observed RW rwlgobs = dlog10(width(lim1)/wave1(lim1)) gfobs = gftab(ntabtot) do i=2,ntabtot - if (rwtab(i) .gt. rwlgobs) then + if (rwtab(i) > rwlgobs) then gfobs = gftab(i-1) + (gftab(i)-gftab(i-1))* . (rwlgobs-rwtab(i-1))/(rwtab(i)-rwtab(i-1)) go to 15 @@ -37,7 +37,7 @@ c*****the computed RW rwlgcal = dlog10(w(ncurve)/wave1(lim1)) gfcal = gftab(ntabtot) do i=2,ntabtot - if (rwtab(i) .gt. rwlgcal) then + if (rwtab(i) > rwlgcal) then gfcal = gftab(i-1) + (gftab(i)-gftab(i-1))* . (rwlgcal-rwtab(i-1))/(rwtab(i)-rwtab(i-1)) go to 20 @@ -52,9 +52,9 @@ c*****square root of the proposed gf shift. 20 error = (w(ncurve)-width(lim1))/width(lim1) ratio = 10.**(gfobs-gfcal) ncurve = ncurve + 1 - if (dabs(error) .ge. 0.0015 .and. ncurve .lt. 20) then + if (dabs(error) >= 0.0015 .and. ncurve < 20) then rwlcomp = dlog10(w(ncurve-1)/wave1(lim1)) - if (rwlcomp .gt. -4.7) then + if (rwlcomp > -4.7) then gf1(ncurve) = gf1(ncurve-1)*dsqrt(ratio) do i=1,ntau kapnu0(lim1,i) = kapnu0(lim1,i)*dsqrt(ratio) @@ -71,7 +71,7 @@ c*****square root of the proposed gf shift. c*****if the observed and computed RWs are close, do a final gf adjustment c*****and finish with one more line recomputation - if (ncurve .eq. 30) then + if (ncurve == 30) then write (nf1out,1001) write (nf2out,1001) endif @@ -84,7 +84,7 @@ c*****and finish with one more line recomputation wid1comp(lim1) = w(1) diff = dlog10(gf1(ncurve)/gf(lim1)) abundout(lim1) = abundin + diff - if (ncurve .ne. 1) then + if (ncurve /= 1) then write (nf1out,1002) ncurve endif |