diff options
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 |