diff options
Diffstat (limited to 'Oneline.f')
| -rwxr-xr-x | Oneline.f | 44 |
1 files changed, 22 insertions, 22 deletions
@@ -19,9 +19,9 @@ c the line profile c*****get started; calculate an initial step size; wavestep only is c*****used in synpop - if (imode .eq. 0) gf1(ncurve) = gf(lim1) + if (imode == 0) gf1(ncurve) = gf(lim1) dellam(1) = 0. - if (wavestep .eq. 0.) then + if (wavestep == 0.) then st1 = wave1(lim1)*dopp(lim1,jtau5)/2.997929e10/5. st1 = dfloat(ifix(10000.*sngl(st1)))/10000. else @@ -32,16 +32,16 @@ c*****used in synpop c*****calculate continuous opacity and intensity/flux at line wavelength wave = wave1(lim1) - if (abs(wave-waveold) .gt. 30.) then + if (abs(wave-waveold) > 30.) then waveold = wave call opacit(2,wave) - if (imode.ne.2 .and. modprintopt.ge.2) + if (imode/=2 .and. modprintopt>=2) . write(nf1out,1002) wave,(kaplam(i),i=1,ntau) call cdcalc(1) first = 0.4343*cd(1) flux = rinteg(xref,cd,dummy1,ntau,first) - if (imode .ne. 2) then - if (iunits .eq. 1) then + if (imode /= 2) then + if (iunits == 1) then write (nf1out,1003) 1.d-4*wave,flux else write (nf1out,1004) wave,flux @@ -51,14 +51,14 @@ c*****calculate continuous opacity and intensity/flux at line wavelength c*****check the wavelength step size; expand/contract as necessary - if (wavestep .eq. 0.) then + if (wavestep == 0.) then wave = wave1(lim1) call taukap call cdcalc(2) first = 0.4343*cd(1) d(1) = rinteg(xref,cd,dummy1,ntau,first) do k=1,30 - if (k .eq. 30) then + if (k == 30) then write (*,1010) wave stop endif @@ -68,17 +68,17 @@ c*****check the wavelength step size; expand/contract as necessary first = 0.4343*cd(1) d(2) = rinteg(xref,cd,dummy1,ntau,first) d2d1 = d(2)/d(1) - if (d2d1 .le. 0.2) then + if (d2d1 <= 0.2) then st1 = st1/1.5 - elseif (d2d1 .le. 0.6) then + elseif (d2d1 <= 0.6) then st1 = st1/1.2 - elseif (d2d1.gt.0.60 .and. d2d1.lt.0.80) then - if (imode .ne. 2) write (nf1out,1001) lim1, + elseif (d2d1>0.60 .and. d2d1<0.80) then + if (imode /= 2) write (nf1out,1001) lim1, . 1000.*width(lim1), storig, st1, k exit - elseif (d2d1 .ge. 0.80) then + elseif (d2d1 >= 0.80) then st1 = st1*1.6 - elseif (d2d1 .ge. 0.90) then + elseif (d2d1 >= 0.90) then st1 = st1*2.1 endif st1 = dble(idnint(10000.*st1))/10000. @@ -96,7 +96,7 @@ c until the depth is very small in the line wing call cdcalc(2) first = 0.4343*cd(1) d(n) = rinteg(xref,cd,dummy1,ntau,first) - if (linprintopt.ge.3 .and. n.eq.1 .and. imode.eq.2) then + if (linprintopt>=3 .and. n==1 .and. imode==2) then do i=1,ntau dummy1(i) = xref(i)*cd(i) enddo @@ -104,21 +104,21 @@ c until the depth is very small in the line wing cdmean = rinteg(xref,dummy1,dummy2,ntau,first)/ . rinteg(xref,cd,dummy2,ntau,first) do i=1,ntau - if (cdmean .lt. cd(i)) exit + if (cdmean < cd(i)) exit enddo write (nf1out,1005) lim1, cdmean, i, xref(i) do i=1,ntau - if (taunu(i)+taulam(i) .ge. 1.) exit + if (taunu(i)+taulam(i) >= 1.) exit enddo write (nf1out,1006) lim1, i, dlog10(tauref(i)), . dlog10(taulam(i)), dlog10(taunu(i)) endif - if (d(n)/d(1) .lt. 0.0050) then + if (d(n)/d(1) < 0.0050) then ndepths = n exit endif - if (n .eq. maxsteps) then - if (d(n).gt.0.001 .and. imode.ne.2) then + if (n == maxsteps) then + if (d(n)>0.001 .and. imode/=2) then write (nf1out,1007) ndepths = maxsteps endif @@ -127,7 +127,7 @@ c until the depth is very small in the line wing c*****finish the calculation - if (imode .ne. 2) write (nf1out,1008) (d(j),j=1,ndepths) + if (imode /= 2) write (nf1out,1008) (d(j),j=1,ndepths) do n=2,ndepths d(ndepths+n-1) = d(n) enddo @@ -142,7 +142,7 @@ c*****finish the calculation enddo first = 2*dellam(ndep)*d(ndep) w(ncurve) = rinteg(dellam,d,dinteg,ndep,first) - if (imode .ne. 2) then + if (imode /= 2) then ew = 1000.0*w(ncurve) gflog = dlog10(gf1(ncurve)) rwlog = dlog10(w(ncurve)/wave1(lim1)) |