Use modern equality operators

1 files changed, 32 insertions, 32 deletions

diff --git a/Binplot.f b/Binplot.f
index edb727f..7a0f9e4 100755
--- a/Binplot.f
+++ b/Binplot.f
@@ -18,7 +18,7 @@ c******************************************************************************
 
 
 c*****for grid syntheses, dump out relevant information to a file
-      if (choice .eq. 'g') then
+      if (choice == 'g') then
          write (nf6out,3001) syncount
          write (nf6out,3002) obsitle, moditle, smitle
       endif
@@ -43,7 +43,7 @@ c*****write smoothing information at the top of the plot
 
 c*****write the non-varying isotopic information at the top of the plot
       do i=1,120
-         if (isoitle(i:i) .ne. ' ') go to 112
+         if (isoitle(i:i) /= ' ') go to 112
       enddo
       isoitle(1:16) = 'no isotopic data'
 112   call sm_relocate (-0.080,1.065)
@@ -51,7 +51,7 @@ c*****write the non-varying isotopic information at the top of the plot
  
 
 c*****define the real plot limits
-      if (xlo .lt. xhi) then
+      if (xlo < xhi) then
          call sm_limits (xlo,xhi,ylo,yhi)
          iflip = 0
       else
@@ -65,7 +65,7 @@ c*****define the real plot limits
 
 c*****draw and label the box for the spectra
       call defcolor (1)
-      if (whichwin .eq. '1of1') then
+      if (whichwin == '1of1') then
          idev = 1
          call sm_window (1,1,1,1,1,1)
       else
@@ -79,7 +79,7 @@ c*****draw and label the box for the spectra
       call sm_lweight (2.0)
       call sm_expand (0.8)
       call sm_box (1,2,4,4)
-      if (iflip .eq. 1) then
+      if (iflip == 1) then
          array = 'Wavenumber'
       else
          array = 'Wavelength'
@@ -93,15 +93,15 @@ c*****plot the synthetic spectra
       call sm_lweight (2.2)
       call sm_expand (0.7)
       do i=1,100
-         if (pec(i) .ne. 0) go to 111
+         if (pec(i) /= 0) go to 111
       enddo            
 111   do j=1,nsyn
-         if (choice.eq.'h' .or. choice.eq.'f' .or.
-     .       choice.eq.'g') then
+         if (choice=='h' .or. choice=='f' .or.
+     .       choice=='g') then
             call defcolor (8)
             call sm_ltype (j-1)
          else
-            if (smterm(1:3) .eq. 'x11') then
+            if (smterm(1:3) == 'x11') then
                call defcolor (j+1)
                call sm_ltype (0)
             else
@@ -110,7 +110,7 @@ c*****plot the synthetic spectra
             endif
          endif
          call sm_connect (xsyn,chunk(1,j),kount) 
-         if (iflip .eq. 1) then
+         if (iflip == 1) then
             call sm_relocate (xhi+0.045*(xlo-xhi),
      .                     ylo+(0.12+0.06*j)*(yhi-ylo))
             call sm_draw (xhi+0.005*(xlo-xhi),
@@ -134,9 +134,9 @@ c*****plot the synthetic spectra
 
    
 c*****plot the observed spectrum
-      if (plotopt .eq. 2) then
+      if (plotopt == 2) then
          call defcolor (1)
-         if (choice.eq.'h' .or. choice.eq.'f') then
+         if (choice=='h' .or. choice=='f') then
             call sm_lweight (4.0)
          else 
             call sm_lweight (2.2)
@@ -146,10 +146,10 @@ c*****plot the observed spectrum
          style(1) = 43.5
          call sm_ptype (style,1)
          mount = lim2obs - lim1obs + 1
-         if (mount .lt. 500) then
+         if (mount < 500) then
             call sm_points (xobs(lim1obs),yobs(lim1obs),mount)
          else
-            if (histoyes .eq. 1) then
+            if (histoyes == 1) then
                call sm_histogram (xobs(lim1obs),yobs(lim1obs),mount)
             else
                call sm_connect (xobs(lim1obs),yobs(lim1obs),mount)
@@ -157,7 +157,7 @@ c*****plot the observed spectrum
          endif
          call sm_lweight (2.2)
          call sm_expand (0.7)
-         if (iflip .eq. 1) then
+         if (iflip == 1) then
             call sm_relocate (xhi+0.05*(xlo-xhi),ylo+0.12*(yhi-ylo))
          else
             call sm_relocate (xlo+0.05*(xhi-xlo),ylo+0.12*(yhi-ylo))
@@ -165,7 +165,7 @@ c*****plot the observed spectrum
          call sm_label (obsitle)
       endif
       do i=1,2
-         if (iflip .eq. 1) then
+         if (iflip == 1) then
             call sm_relocate (xhi+0.05*(xlo-xhi),
      .                        ylo+0.01+0.06*(2-i)*(yhi-ylo))
          else
@@ -174,41 +174,41 @@ c*****plot the observed spectrum
          endif
          call sm_label (modbin(i))
       enddo
-      if (whichwin.eq.'1of1' .or. plotopt.ne.2) then
+      if (whichwin=='1of1' .or. plotopt/=2) then
          return
       endif
 
 
 c*****this section of code is executed only if a deviations plot is desired;
 c     find the starting and stopping points in the arrays for the deviations
-      if (xsyn(kount) .le. xobs(lim1obs)) go to 1000
-      if (xsyn(1) .gt. xobs(lim2obs)) go to 1000
-      if (xsyn(1) .gt. xobs(lim1obs)) go to 150
+      if (xsyn(kount) <= xobs(lim1obs)) go to 1000
+      if (xsyn(1) > xobs(lim2obs)) go to 1000
+      if (xsyn(1) > xobs(lim1obs)) go to 150
       lim3obs = lim1obs
       do k=2,kount
-         if (xsyn(k) .gt. xobs(lim3obs)) then
+         if (xsyn(k) > xobs(lim3obs)) then
             lim1syn = k - 1
             go to 155
          endif
       enddo
 150   lim1syn = 1
       do l=lim1obs,lim2obs
-         if (xsyn(lim1syn) .le. xobs(l)) then
+         if (xsyn(lim1syn) <= xobs(l)) then
             lim3obs = l 
             go to 155
          endif
       enddo
-155   if (xsyn(kount) .lt. xobs(lim2obs)) go to 160
+155   if (xsyn(kount) < xobs(lim2obs)) go to 160
       lim4obs = lim2obs
       do k=lim1syn,kount
-         if (xsyn(k) .gt. xobs(lim4obs)) then
+         if (xsyn(k) > xobs(lim4obs)) then
             lim2syn = k
             go to 165
          endif
       enddo
 160   lim2syn = kount
       do l=lim3obs,lim2obs
-         if (xsyn(lim2syn) .lt. xobs(l)) then
+         if (xsyn(lim2syn) < xobs(l)) then
             lim4obs = l - 1
             go to 165
          endif
@@ -221,7 +221,7 @@ c  of the synthetic spectra is considered sufficient
          lpoint = lim1syn
          devsigma = 0.
          do i=lim3obs,lim4obs
-170         if (xsyn(lpoint+1) .lt. xobs(i)) then
+170         if (xsyn(lpoint+1) < xobs(i)) then
                lpoint = lpoint + 1
                go to 170
             endif
@@ -235,7 +235,7 @@ c  of the synthetic spectra is considered sufficient
          
       
 c  from first set of deviations, define the plot limits, draw and label box
-         if (j .eq. 1) then
+         if (j == 1) then
             yup = -1000.
             ydown = +1000.
             do i=lim3obs,lim4obs
@@ -265,12 +265,12 @@ c  from first set of deviations, define the plot limits, draw and label box
 
 
 c  plot the array of deviations
-         if (choice.eq.'h' .or. choice.eq.'f' .or.
-     .       choice.eq.'g') then
+         if (choice=='h' .or. choice=='f' .or.
+     .       choice=='g') then
             call defcolor (8)
             call sm_ltype (j-1)
          else
-            if (smterm(1:3) .eq. 'x11') then
+            if (smterm(1:3) == 'x11') then
                call defcolor (j+1)
                call sm_ltype (0)
             else
@@ -289,7 +289,7 @@ c  plot the array of deviations
          call sm_draw(xhi-0.215*(xhi-xlo),
      .                 ydown+(0.10+0.06*j)*(yup-ydown))
          call sm_label (array)
-         if (choice .eq. 'g') then
+         if (choice == 'g') then
             noff = 80*(j-1)
             write (nf6out,3002) abitle(noff+1:noff+80)
             write (nf6out,3003) devsigma, velsh
@@ -298,7 +298,7 @@ c  plot the array of deviations
 
 
 c  reset the spectrum plot boundaries before exiting
-      if(xlo .lt. xhi) then
+      if(xlo < xhi) then
          call sm_limits (xlo,xhi,ylo,yhi)
          iflip = 0
       else