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 /Synspec.f | |
| parent | c396a4bc7a3282de9d2191de0e1076fec9df0783 (diff) | |
| download | moog-70fd14e98bc476a84aef352eaf063cfcba44ab83.tar.gz | |
Use modern equality operators
Diffstat (limited to 'Synspec.f')
| -rwxr-xr-x | Synspec.f | 40 |
1 files changed, 20 insertions, 20 deletions
@@ -15,19 +15,19 @@ c****************************************************************************** c*****initialize the synthesis
write (nf1out,1101)
write (nf2out,1002) moditle(1:73)
- if (iunits .eq. 1) then
+ if (iunits == 1) then
write (nf2out,1103) start/1.d4, sstop/1.d4,
. step/1.d4, delta/1.d4
else
write (nf2out,1102) start,sstop,step,delta
endif
- if (iraf .eq. 1) then
+ if (iraf == 1) then
npoint = (sstop-start)/step
write (nf4out,1104) npoint,wave,wave,step,step
write (nf4out,1105)
write (nf4out,1106) moditle
do j=1,93
- if (pec(j) .gt. 0 ) then
+ if (pec(j) > 0 ) then
dummy1(j) = dlog10(xabund(j)) + 12.0
write (nf4out,1107) names(j),dummy1(j)
endif
@@ -35,7 +35,7 @@ c*****initialize the synthesis write (nf4out,1108) vturb(1)
write (nf4out,1109)
endif
- if (mode .ne. 4) then
+ if (mode /= 4) then
lim1line = 0
lim2line = 0
lim1obs = 0
@@ -56,15 +56,15 @@ c spectrum wavelength, if needed do n=1,kount
num = num + 1
wave = oldstart + (n-1)*step
- if (dabs(wave-wavl)/wave .ge. 0.001) then
+ if (dabs(wave-wavl)/wave >= 0.001) then
wavl = wave
call opacit (2,wave)
- if (modprintopt .ge. 2)
+ if (modprintopt >= 2)
. write (nf1out,1001) wave,(kaplam(i),i=1,ntau)
call cdcalc (1)
first = 0.4343*cd(1)
flux = rinteg(xref,cd,dummy1,ntau,first)
- if (iunits .eq. 1) then
+ if (iunits == 1) then
write (nf1out,1003) 1.d-4*wave,flux
else
write (nf1out,1004) wave,flux
@@ -75,9 +75,9 @@ c spectrum wavelength, if needed c*****find the appropriate set of lines for this wavelength, reading
c in a new set if this is the initial depth calculation or if
c needed because the line list end has been reached
- if (mode .eq. 3) then
+ if (mode == 3) then
20 call linlimit
- if (lim2line .lt. 0) then
+ if (lim2line < 0) then
call inlines (2)
call nearly (1)
go to 20
@@ -90,7 +90,7 @@ c needed because the line list end has been reached c*****compute a spectrum depth at this point; if there are no absorption
c lines in the interval then just set the depth to zero without
c extensive line calculations
- if (lineflag .lt. 0) then
+ if (lineflag < 0) then
d(num) = 0.
else
call taukap
@@ -98,45 +98,45 @@ c extensive line calculations first = 0.4343*cd(1)
d(num) = rinteg(xref,cd,dummy1,ntau,first)
endif
- if (mod(n,10) .eq. 0) then
- if (iraf .eq. 1) then
+ if (mod(n,10) == 0) then
+ if (iraf == 1) then
do j = 1,10
dd(num-10+j) = 1. - d(num-10+j)
enddo
write (nf4out,1110) (dd(num-10+j),j=1,10)
endif
- if (iunits .eq. 1) then
+ if (iunits == 1) then
wave3 = 1.d-4*(wave - 9.0*step)
write (nf1out,1112) wave3,(d(num-10+j),j=1,10)
else
wave3 = wave - 9.0*step
write (nf1out,1111) wave3,(d(num-10+j),j=1,10)
endif
- if (nf2out .gt. 0) write (nf2out,1110) (d(num-10+j),j=1,10)
+ if (nf2out > 0) write (nf2out,1110) (d(num-10+j),j=1,10)
endif
- if (num .ge. 5000) num = 0
+ if (num >= 5000) num = 0
enddo
c*****finish the synthesis
nn = mod(num,10)
- if (nn .ne. 0) then
- if (iraf .eq. 1) then
+ if (nn /= 0) then
+ if (iraf == 1) then
do j=1,nn
dd(num-nn+j) = 1. - d(num-nn+j)
enddo
write (nf4out,1110) (dd(num-nn+j),j=1,nn)
endif
- if (iunits .eq. 1) then
+ if (iunits == 1) then
wave3 = 1.d-4*(wave - 9.0*step)
write (nf1out,1112) wave3,(d(num-nn+j),j=1,nn)
else
wave3 = wave - 9.0*step
write (nf1out,1111) wave3,(d(num-nn+j),j=1,nn)
endif
- if (nf2out .gt. 0) write (nf2out,1110) (d(num-nn+j),j=1,nn)
+ if (nf2out > 0) write (nf2out,1110) (d(num-nn+j),j=1,nn)
endif
- if (iunits .eq. 1) then
+ if (iunits == 1) then
write (nf1out,1113) 1.d-4*wave
else
write (nf1out,1114) wave
|