diff options
Diffstat (limited to 'Partfn.f')
| -rwxr-xr-x | Partfn.f | 42 |
1 files changed, 42 insertions, 0 deletions
diff --git a/Partfn.f b/Partfn.f new file mode 100755 index 0000000..5463afb --- /dev/null +++ b/Partfn.f @@ -0,0 +1,42 @@ + + subroutine partfn (atom,jmark) +c****************************************************************************** +c This routine computes arrays or single partition functions; +c It uses the data and formulation of Kurucz's ATLAS9 program. +c However, the partition functions have been updated for several species, +c and those are calculated in subroutine *newpart*. +c****************************************************************************** + + implicit real*8 (a-h,o-z) + include 'Atmos.com' + include 'Linex.com' + include 'Quants.com' + + + iat = 10*nint(atom) - 1 + at = dfloat(iat)/10. + iatom = nint(atom) + + +c*****compute partition functions for 4 ionization states of an element. + do k=1,4 + iat = iat + 1 + at = dfloat(iat)/10. + if (partflag(iatom,k) .gt. 0) then + do i=1,ntau + u(jmark,k,i) = partnew(at,k,i) + enddo + else + do i=1,ntau + u(jmark,k,i) = ucalc(at,i) + enddo + endif + enddo + + + return + end + + + + |