Use modern equality operators

1 files changed, 41 insertions, 41 deletions

diff --git a/Specplot.f b/Specplot.f
index 03ea153..76992d1 100755
--- a/Specplot.f
+++ b/Specplot.f
@@ -18,7 +18,7 @@ c******************************************************************************
       real*8 inx1,inx2,iny1,iny2
 
 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, linitle, smitle
       endif
@@ -44,24 +44,24 @@ c*****write smoothing information at the top of the plot
       call sm_label (smitle)
 
 
-      if (isoitle(1:10) .eq. '          ') then
+      if (isoitle(1:10) == '          ') then
          isoitle(1:16) = 'no isotopic data'
       endif
       call sm_relocate (-0.120,1.075)
       call sm_label (isoitle(1:120))
-      if (numiso .gt. 3) then
+      if (numiso > 3) then
          call sm_relocate (-0.120,1.045)
          call sm_label (isoitle(121:240))
       endif
-      if (control .eq. 'gridplo' .or.
-     .    control .eq. 'gridsyn' .or.
-     .    control .eq. 'gridend') then
+      if (control == 'gridplo' .or.
+     .    control == 'gridsyn' .or.
+     .    control == 'gridend') then
          write (nf6out,3002) isoitle(1:120)
          write (nf6out,3002) isoitle(121:240)
       endif
 
 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
@@ -75,7 +75,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
@@ -89,7 +89,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'
@@ -103,15 +103,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
@@ -120,7 +120,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),
@@ -138,10 +138,10 @@ c*****plot the synthetic spectra
          noff = 80*(j-1)
          call sm_lweight (2.2)
          call sm_label (abitle(noff+1:noff+80))
-         if ((control .eq. 'gridplo' .or.
-     .        control .eq. 'gridsyn' .or.
-     .        control .eq. 'gridend') .and. 
-     .        whichwin.eq.'1of1') then
+         if ((control == 'gridplo' .or.
+     .        control == 'gridsyn' .or.
+     .        control == 'gridend') .and.
+     .        whichwin=='1of1') then
             write (nf6out,3002) abitle(noff+1:noff+80)
          endif
       enddo
@@ -150,10 +150,10 @@ 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' .or.
-     .       choice.eq.'g') then
+         if (choice=='h' .or. choice=='f' .or.
+     .       choice=='g') then
             call sm_lweight (4.0)
          else 
             call sm_lweight (2.2)
@@ -163,10 +163,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)
@@ -174,54 +174,54 @@ 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))
          endif
          call sm_label (obsitle)
       endif
-      if (iflip .eq. 1) then
+      if (iflip == 1) then
          call sm_relocate (xhi+0.05*(xlo-xhi),ylo+0.06*(yhi-ylo))
       else
          call sm_relocate (xlo+0.05*(xhi-xlo),ylo+0.06*(yhi-ylo))
       endif
       call sm_label (moditle)
-      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
@@ -235,7 +235,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
@@ -253,7 +253,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
@@ -283,12 +283,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
@@ -307,7 +307,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
@@ -317,7 +317,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