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 /Correl.f | |
| parent | c396a4bc7a3282de9d2191de0e1076fec9df0783 (diff) | |
| download | moog-70fd14e98bc476a84aef352eaf063cfcba44ab83.tar.gz | |
Use modern equality operators
Diffstat (limited to 'Correl.f')
| -rwxr-xr-x | Correl.f | 14 |
1 files changed, 7 insertions, 7 deletions
@@ -26,13 +26,13 @@ c****************************************************************************** c*****find the array subscripts for the minimum and maximum c wavelengths in the synthetic spectrum. do i=1,kount - if (xsyn(i) .ge. wavemin) then + if (xsyn(i) >= wavemin) then ilosyn = i go to 5 endif enddo 5 do i=ilosyn,kount - if (xsyn(i) .gt. wavemax) then + if (xsyn(i) > wavemax) then ihisyn = i - 1 go to 10 endif @@ -51,7 +51,7 @@ c*****dump the synthetic spectrum wavelengths and fluxes into working arrays c*****find the array subscript for the minimum wavelength in the c observed spectrum do j=1,lount - if (xobs(j) .ge. xgood(1)) then + if (xobs(j) >= xgood(1)) then jloobs = j go to 15 endif @@ -63,7 +63,7 @@ c spectrum that matches the wavelength step size of the c synthetic spectrum 15 j = jloobs do i=1,itotsyn - if (xgood(i) .gt. (xobs(j+1)+xobs(j))/2.) j = j + 1 + if (xgood(i) > (xobs(j+1)+xobs(j))/2.) j = j + 1 q = xgood(i) zgood(i) = yobs(j-1)*(q-xobs(j))*(q-xobs(j+1))/ . ((xobs(j-1)-xobs(j))*(xobs(j-1)-xobs(j+1))) + @@ -89,12 +89,12 @@ c*****interpolate to find the maximum of the correlation function corrmax = spy(1) mc = 1 do i=2,imax - if (spy(i) .gt. corrmax) then + if (spy(i) > corrmax) then corrmax = spy(i) mc = i endif enddo - if (mc.eq.1 .or. mc.eq.imax) then + if (mc==1 .or. mc==imax) then deltawave = 0. return endif @@ -113,7 +113,7 @@ c*****interpolate to find the maximum of the correlation function jmax = 1 corrmax = yinterp(1) do j=2,21 - if (yinterp(j) .gt. corrmax) then + if (yinterp(j) > corrmax) then jmax = j corrmax = yinterp(j) shiftmax = xinterp(j) |