diff options
Diffstat (limited to 'Gammabark.f')
| -rwxr-xr-x | Gammabark.f | 22 |
1 files changed, 11 insertions, 11 deletions
diff --git a/Gammabark.f b/Gammabark.f index 94a697e..40794d4 100755 --- a/Gammabark.f +++ b/Gammabark.f @@ -23,8 +23,8 @@ c*****on first entry to this routine, read damping data from either c 'Barklem.dat' or 'BarklemUV.dat', depending on the wavelength region c of the linelist; read Van der Waals params, if radiative damping c data are present - if (firstread .eq. 0) then - if (wave1(nlines) .gt. 3000.) then + if (firstread == 0) then + if (wave1(nlines) > 3000.) then nwant = 35 else nwant = 36 @@ -33,7 +33,7 @@ c data are present call blankstring (line) read (nwant,1001,end=10) line read (line,*) wavebk(k), idbk(k), gammabk(k), alphabk(k) - if (line(34:) .ne. ' ') then + if (line(34:) /= ' ') then read(line(34:),*) gammarad(k) else gammarad(k) = 0.0 @@ -48,20 +48,20 @@ c*****identify the Barklem list positions of the wavelength limits of c the input line list wavemin = 10000000. do j=1,nlines+nstrong - if (wave1(j) .lt. wavemin) wavemin = wave1(j) + if (wave1(j) < wavemin) wavemin = wave1(j) enddo wavemax = 0. do j=1,nlines+nstrong - if (wave1(j) .gt. wavemax) wavemax = wave1(j) + if (wave1(j) > wavemax) wavemax = wave1(j) enddo do k=1,numbark - if (wavemin-wavebk(k) .lt. 1.0) then + if (wavemin-wavebk(k) < 1.0) then nummin = k exit endif enddo do k=nummin,numbark - if (wavebk(k)-wavemax .gt. 1.0) then + if (wavebk(k)-wavemax > 1.0) then nummax = k exit endif @@ -73,19 +73,19 @@ c*****search for Barklem data gambark(j) = -1. alpbark(j) = -1. gamrad(j) = -1. - if (atom1(j) .gt. 100.) cycle + if (atom1(j) > 100.) cycle iatom10 = nint(10.*atom1(j)) do k=nummin,nummax waveerror = -(wave1(j) - wavebk(k))/wavebk(k) iii = nint(10.*idbk(k)) - if (dabs(waveerror).lt.5.0d-06 .and. - . iii .eq. iatom10) then + if (dabs(waveerror)<5.0d-06 .and. + . iii == iatom10) then gamrad(j) = gammarad(k) gambark(j) = 10.**gammabk(k) alpbark(j) = (1.-alphabk(k))/2. exit endif - if (waveerror .gt. 5.0d-06) exit + if (waveerror > 5.0d-06) exit enddo enddo |