diff options
Diffstat (limited to 'Damping.f')
| -rwxr-xr-x | Damping.f | 62 |
1 files changed, 31 insertions, 31 deletions
@@ -13,43 +13,43 @@ c****************************************************************************** j = linnumber iwave = int(wave1(j)) iatom10 = nint(10.*atom1(j)) - if (dampnum(j) .lt. 0.) dampnum(j) = 10.**dampnum(j) + if (dampnum(j) < 0.) dampnum(j) = 10.**dampnum(j) c*****for a few lines, explicit detailed broadening terms have c appeared in the literature, and so do these lines with a c sepaarate subroutine - if (itru .eq. 0) then + if (itru == 0) then c Ca II - if (iatom10 .eq. 201) then - if (iwave .eq. 8498 .or. - . iwave .eq. 8542 .or. - . iwave .eq. 8662 .or. - . iwave .eq. 3933) then + if (iatom10 == 201) then + if (iwave == 8498 .or. + . iwave == 8542 .or. + . iwave == 8662 .or. + . iwave == 3933) then call trudamp (j) damptype(j) = 'TRUEgam' return endif c CH - elseif(iatom10 .eq. 1060) then - if (iwave .eq. 3693) then + elseif(iatom10 == 1060) then + if (iwave == 3693) then call trudamp (j) damptype(j) = 'TRUEgam' return endif c Ca I - elseif (iatom10 .eq. 200) then - if (iwave.eq.6717 .or. iwave.eq.6318 - . .or. iwave.eq.6343 .or. iwave.eq.6361) then + elseif (iatom10 == 200) then + if (iwave==6717 .or. iwave==6318 + . .or. iwave==6343 .or. iwave==6361) then call trudamp (j) damptype(j) = 'TRUEgam' return endif c Ca I autoionization - elseif (iatom10 .eq. 200) then - if (iwave.eq.6318 .or. - . iwave.eq.6343 .or. - . iwave.eq.6361) then + elseif (iatom10 == 200) then + if (iwave==6318 .or. + . iwave==6343 .or. + . iwave==6361) then call trudamp (j) damptype(j) = 'TRUEgam' return @@ -82,19 +82,19 @@ c c6 done as in dampingopt = 0 c*****these damping calculations are done at each atmosphere level - if (linprintopt .gt. 2) write (nf1out,1001) j, wave1(j) + if (linprintopt > 2) write (nf1out,1001) j, wave1(j) do i=1,ntau ich = nint(charge(j)) v1 = dsqrt(2.1175d8*t(i)*(1.0/amass(j)+1.008)) c*****first calculate an Unsold approximation to gamma_VanderWaals - if (atom1(j) .gt. 100.) then + if (atom1(j) > 100.) then ebreakup = 7.0 else ebreakup = chi(j,ich) endif - if (e(j,1).ge.ebreakup .or. e(j,2).ge.ebreakup) then + if (e(j,1)>=ebreakup .or. e(j,2)>=ebreakup) then unsold = 1.0e-33 else unsold = dabs(1.61d-33*(13.598*charge(j)/(ebreakup - @@ -105,15 +105,15 @@ c*****first calculate an Unsold approximation to gamma_VanderWaals c*****dampingopt = 0 or c*****dampingopt = 1 and no Barklem data - if (dampingopt .eq. 0 .or. - . (dampingopt.eq.1 .and. gambark(j).lt.0)) then - if (dampnum(j) .eq. 0.0) then + if (dampingopt == 0 .or. + . (dampingopt==1 .and. gambark(j)<0)) then + if (dampnum(j) == 0.0) then damptype(j) = 'UNSLDc6' gammav = 17.0*unsold**0.4*v1**0.6*numdens(1,1,i) - elseif (dampnum(j) .lt. 1.0d-15) then + elseif (dampnum(j) < 1.0d-15) then damptype(j) = ' MYc6' gammav = 17.0*dampnum(j)**0.4*v1**0.6*numdens(1,1,i) - elseif (dampnum(j) .lt. 1.0d-04) then + elseif (dampnum(j) < 1.0d-04) then damptype(j) = 'MYgamma' gammav = dampnum(j)*(t(i)/10000.)**0.3*numdens(1,1,i) else @@ -124,23 +124,23 @@ c*****dampingopt = 1 and no Barklem data c*****dampingopt = 1 with extant Barklem data - elseif (dampingopt.eq.1 .and. gambark(j).gt.0.) then + elseif (dampingopt==1 .and. gambark(j)>0.) then damptype(j) = 'BKgamma' gammav = . gambark(j)*(t(i)/10000.)**alpbark(j)*numdens(1,1,i) c*****dampingopt = 2 - elseif (dampingopt .eq. 2) then + elseif (dampingopt == 2) then damptype(j) = 'BLKWLc6' gammav = 17.0*((1.0 + 0.67*e(j,1))*unsold)**0.4* . v1**0.6*numdens(1,1,i) c*****dampingopt = 3 - elseif (dampingopt .eq. 3) then + elseif (dampingopt == 3) then damptype(j) = 'NXTGNc6' - if (dampnum(j) .le. 1.0d-10) dampnum(j) = 1.0 + if (dampnum(j) <= 1.0d-10) dampnum(j) = 1.0 c6h = dabs(1.01d-32*(charge(j)**2)* . (13.598/(ebreakup - e(j,1)))**2 - 1.61d-33* . (13.598/(ebreakup-e(j,2)))**2) @@ -160,7 +160,7 @@ c*****dampingopt = 3 c*****compute radiative broadening either by an approximate formula or c*****the value in Barklem.dat) - if (gamrad(j).ne.0.0 .and. dampingopt .eq. 1) then + if (gamrad(j)/=0.0 .and. dampingopt == 1) then gammar = gamrad(j) else gammar = 2.223d15/wave1(j)**2 @@ -169,7 +169,7 @@ c*****the value in Barklem.dat) c*****now Stark broadening (approximate formulae) excdiff = chi(j,nint(charge(j))) - e(j,2) - if (excdiff .gt. 0.0 .and. atom1(j).lt.100.) then + if (excdiff > 0.0 .and. atom1(j)<100.) then effn2 = 13.6*charge(j)**2/excdiff else effn2 = 25. @@ -180,7 +180,7 @@ c*****now Stark broadening (approximate formulae) c*****now finish by summing the gammas and computing the Voigt *a* values gammatot = gammar + gammas + gammav a(j,i) = gammatot*wave1(j)*1.0d-8/(12.56636*dopp(j,i)) - if (linprintopt .gt. 2) write (nf1out,1002) i, gammar, + if (linprintopt > 2) write (nf1out,1002) i, gammar, . gammas, gammav, gammatot, a(j,i) enddo return |