Commit 05b5b145 authored by Michal Bursa's avatar Michal Bursa
Browse files

TLUSTY added

parents
COMMON/FIXALP/ABSO(MFREQ),EMIS(MFREQ),SCAT(MFREQ),
* REIT(MDEPTH),REIN(MDEPTH),REIM(MDEPTH),
* REIP(MLVEXP,MDEPTH),
* AREIT(MDEPTH),AREIN(MDEPTH),AREIM(MDEPTH),
* AREIP(MLVEXP,MDEPTH),
* CREIT(MDEPTH),CREIN(MDEPTH),CREIM(MDEPTH),
* CREIP(MLVEXP,MDEPTH),
* REIX(MDEPTH),CREIX(MDEPTH),
* REDX(MDEPTH),REDT(MDEPTH),REDN(MDEPTH),
* REDM(MDEPTH),REDP(MLVEXP,MDEPTH),
* REDXM(MDEPTH),REDTM(MDEPTH),REDNM(MDEPTH),
* REDMM(MDEPTH),REDPM(MLVEXP,MDEPTH),
* REDTP(MDEPTH),REDNP(MDEPTH),REDXP(MDEPTH),
* REDMP(MDEPTH),REDPP(MLVEXP,MDEPTH),
* HEIT(MDEPTH),HEIN(MDEPTH),HEIM(MDEPTH),
* HEIP(MLVEXP,MDEPTH),
* HEITM(MDEPTH),HEINM(MDEPTH),HEIMM(MDEPTH),
* HEIPM(MLVEXP,MDEPTH),
* HEITP(MDEPTH),HEINP(MDEPTH),HEIMP(MDEPTH),
* HEIPP(MLVEXP,MDEPTH),
* EHET(MDEPTH),EHEN(MDEPTH),ERET(MDEPTH),EREN(MDEPTH),
* EHEP(MLVEX3,MDEPTH),EREP(MLVEX3,MDEPTH),
* APT(MLVEXP,MDEPTH),APN(MLVEXP,MDEPTH),
* AAPT(MLVEX3,MDEPTH),AAPN(MLVEX3,MDEPTH),
* CAPT(MLVEX3,MDEPTH),CAPN(MLVEX3,MDEPTH),
* APP(MLVEXP,MLVEXP,MDEPTH),
* AAPP(MLVEX3,MLVEX3,MDEPTH),
* CAPP(MLVEX3,MLVEX3,MDEPTH),QTLAS,
* IFALI,IFPOPR,irprec,ifprec,itold1,itold2,itlas
COMMON A(MTOT,MTOT), B(MTOT,MTOT), C(MTOT,MTOT),
* E(MTOT,MTOT),
* VECL(MTOT), Y1(MTOT), Y2(MTOT),
* PSI0(MTOT), PSIM(MTOT), PSIP(MTOT),
* RAD0(MTOT), RADM(MTOT), RADP(MTOT),
* FKM(MFREX), FK0(MFREX), FKP(MFREX),
* ABSOM(MFREX), ABSO0(MFREX), ABSOP(MFREX),
* EMISM(MFREX), EMIS0(MFREX), EMISP(MFREX),
* SCATM(MFREX), SCAT0(MFREX), SCATP(MFREX),
* DABTM(MFREX), DABT0(MFREX), DABTP(MFREX),
* DEMTM(MFREX), DEMT0(MFREX), DEMTP(MFREX),
* DABNM(MFREX), DABN0(MFREX), DABNP(MFREX),
* DEMNM(MFREX), DEMN0(MFREX), DEMNP(MFREX),
* DABMM(MFREX), DABM0(MFREX), DABMP(MFREX),
* DEMMM(MFREX), DEMM0(MFREX), DEMMP(MFREX),
* WDEPM(MFREX), WDEP0(MFREX), WDEPP(MFREX),
* SBFM(MLEVEL), SBF0(MLEVEL), SBFP(MLEVEL),
* HEX(MLEVEL), REX(MLEVEL), REXA(MLEVEL),
* DSBFM(MLEVEL), DSBF0(MLEVEL), DSBFP(MLEVEL),
* SUMDCH(MLEVEL),
* DRCHM(MLEVEL,MFREX), DRETM(MLEVEL,MFREX),
* DRCH0(MLEVEL,MFREX), DRET0(MLEVEL,MFREX),
* DRCHP(MLEVEL,MFREX), DRETP(MLEVEL,MFREX)
COMMON/EXPRAD/
* ABSOEX(MFREX,MDEPTH),EMISEX(MFREX,MDEPTH),
* SCATEX(MFREX,MDEPTH),
* DABTEX(MFREX,MDEPTH),DEMTEX(MFREX,MDEPTH),
* DABNEX(MFREX,MDEPTH),DEMNEX(MFREX,MDEPTH),
* DABMEX(MFREX,MDEPTH),DEMMEX(MFREX,MDEPTH),
* DRCHEX(MLVEXP,MFREX,MDEPTH),
* DRETEX(MLVEXP,MFREX,MDEPTH)
COMMON/BPOCOM/ESEMAT(MLEVEL,MLEVEL),BESE(MLEVEL),
* ATT(MLEVEL),ANN(MLEVEL)
C
PARAMETER (MPAG=MLEVEL/6+1)
C
CHARACTER*40 FIDATA(MION),FIODF1(MION),FIODF2(MION),FIBFCS(MION)
CHARACTER*10 TYPLEV(MLEVEL)
CHARACTER*4 TYPION(MION)
C
COMMON/ATOPAR/AMASS(MATOM),ABUND(MATOM,MDEPTH),NUMAT(MATOM),
* N0A(MATOM),NKA(MATOM),nref(matom),iatex(30),
* nrefs(matom,mdepth),
* iifix(matom),iatref,modref
common/atomas/amas(100)
COMMON/IONPAR/FF(MION),CHARG2(MION),
* NFIRST(MION),NLAST(MION),NNEXT(MION),
* IZ(MION),IUPSUM(MION),ICUP(MION),ilte(mion),
* iltion(mion)
COMMON/LEVPAR/ENION(MLEVEL),G(MLEVEL),NQUANT(MLEVEL),
* IATM(MLEVEL),IEL(MLEVEL),ILK(MLEVEL),ilin(mlevel),
* iltlev(mlevel),indlev(mlevel),
* imodl(mlevel),iiexp(mlevel),iifor(mlevel),
* ipzert(mlevel),
* igzert(mlevel),indlgz(mlevel),iinonz(mlevel),
* NLVEXP,NLVFOR,NLVEXZ,LBPFX
COMMON/TRAPAR/FR0(MTRANS),OSC0(MTRANS),CPAR(MTRANS),
* FRQMX(MTRANS),FR0PC(MTRANS),OMECOL(MLEVEL,MLEVEL),
* XGRAD,STRL1,STRL2,STRLX,
* ILOW(MTRANS),IUP(MTRANS),INDEXP(MTRANS),
* KFR0(MTRANS),KFR1(MTRANS),ILUCTR(MTRANS),
* IFC0(MTRANS),IFC1(MTRANS),
* IFR0(MTRANS),IFR1(MTRANS),ITRA(MLEVEL,MLEVEL),
* IPROF(MTRANS),ICOL(MTRANS),INTMOD(MTRANS),
* ITRCON(MTRANS),IDIEL(MTRANS),
* IJTF(MTRANS),LCOMP(MTRANS),LINE(MTRANS)
COMMON/PHOSET/S0CS(MLEVEL),ALFCS(MLEVEL),BETCS(MLEVEL),
* GAMCS(MLEVEL),IBF(MLEVEL)
COMMON/TOPCS/CTOP(MFIT,MCROSS), !sigma=alog10(sigma/10^-18) of fit point
+ XTOP(MFIT,MCROSS) ! x = alog10(nu/nu0) of fit point
COMMON/VOIPAR/GAMAR(MVOIGT),STARK1(MVOIGT),STARK2(MVOIGT),
* STARK3(MVOIGT),VDWH(MVOIGT)
COMMON/HECRAT/COLHE1(19,19)
COMMON/TRACOR/LEXP(MTRANS),LALI(MTRANS)
COMMON/TRAALI/NFFIX,IFSUB,IFLEV
COMMON/AUXIND/IATH,IATHE,IELH,IELHM,IELHE1,IELHE2,IOPADD
COMMON/IONFIL/FIDATA,FIODF1,FIODF2,FIBFCS
COMMON/IONDAT/IATI(MION),IZI(MION),NLEVS(MION),NLLIM(MION)
COMMON/OSCHYD/OSH(20,20)
COMMON/PRINTP/NPGPOP,IIPR(6,MPAG),TYPLEV,TYPION
c
PARAMETER (MTRPRD=5)
COMMON/PRDPAR/DOPTR(MTRPRD,MDEPTH),COHER(MTRPRD,MDEPTH),
* PJBAR(MTRPRD,MDEPTH),RJBAR(MTRPRD,MDEPTH),
* XPDIV,
* IPRD(MTRANS),ITRTOT(MTRPRD),NTRPRD,IFPRD
C
C Parameters that specify dimensions of arrays
C
PARAMETER (MATOM = 99, ! maximum number of explicit atoms
* MION = 190, ! maximum number of explicit ions
* MLEVEL = 295, ! maximum number of explicit levels
* MLVEXP = 290, ! maximum number of explicit linearized levels
* MTRANS = 3000, ! maximum number of all transitions
* ! which are somehow taken into account
* MDEPTH = 70, ! maximum number of depth points
* MDEPTC = 70, ! maximum number of depth points (Compton)
* MFREQ = 1500, ! maximum number of frequency points
c * MFREQ = 1480, ! maximum number of frequency points
* MFREQC= 2000, ! maximum number of freq.points in continuum
* MFREX = 340, ! maximum number of linearized frequencies
* MMU = 80, ! maximum number of angle points
* MMUC = 80, ! maximum number of angle points (Compton)
* MANGEX = 80, ! maximum number of angle points (external rad.)
* MFRCOR = 1, ! maximum number of frequency points
* ! in one line for correction weights
* MFETA = 1, ! maximum number of frequency points
* ! in one line for ETA procedure (outdated)
* MCROSS= 800, ! maximum number of photoionization
* ! cross-sections + related quantities
* MBF = 800, ! maximum number of b-f transitions
* MFIT = 225, ! maximum number of fit points (OP b-f c.s)
* MITJ= 220, ! maximum number of overlapping transitions
* MMCDW = 14, ! maximum number of levels with pseudocont.
* MMER = 2, ! maximum number of merged levels
* MZZ = 10, ! maximum charge for occup.prob. ions
* MLINFR = 55, ! maximum number of lines at given freq.
* MFREQL = 2000, ! maximum number of frequencies per line
* MVOIGT= 175, ! maximum number of lines with Voigt
* ! profile
* NLMX = 80, ! highest hydrogenic level considered
* MSMX = 1, ! size of largest matrix kept in memory in SOLVE
* MTOT = 300, ! maximum number of model unknowns
* ! (ie. maximum dimension of vector PSI)
c * MFREQ1= 1, ! =1 for ISPLIN<5; =MFREQ otherwise
* MFREQ1=MFREQ,
* MFRTAB= 2, ! maximum number of freqeuncies in opac.table
* MTABT = 50, ! max.number of temperatures in opac.table
* MTABR = 50, ! max.number of densities in opac.table
* MCOND= 1, ! dummy
* MLEVE3= 1, ! =1 for diag.precond; =MLEVEL for tridia.prec.
* MLVEX3=MLVEXP, ! =1 for diag.oper.; =MLVEXP for tridiag. .
* MTRAN3=MTRANS) ! =1 for diag.oper; =MTRANS for tridiag.
C
PARAMETER (MAROWS = MFREX + 2, ! maximum number of rows of matrix
* ! A of complete linearization
* MCCOLS = MLEVEL+ 5) ! maximum number of columns of
* ! matrix C of compl. linearization
C
C Basic physical constants
C
PARAMETER (H = 6.6256D-27, ! Planck constant h
* BOLK = 1.38054D-16, ! Boltzmann constant k
* HK = 4.79928144D-11, ! h/k
* CAS = 2.997925D18, ! light speed c (A/s)
* EH = 2.17853041D-11, ! ionizaton energy of hydrogen
* BN = 1.4743D-2, ! 2*h/c**3, c -light speed
* SIGE = 6.6516D-25, ! Thomson scattering c-s
* SIG4P = 4.5114062D-6, ! Stefan-Boltzmann const/4pi
* PI4H = 1.8966D27, ! 4pi/h
* PCK = 4.19168946D-10, ! 4pi/c
* HMASS = 1.67333D-24) ! mass of hydrogen atom
C
C Basic mathematical constants
C
PARAMETER (UN = 1.0D0,
* HALF = 0.5D0,
* TWO = 2.0D0)
C
C Unit number
C
PARAMETER (IBUFF = 95)
C
C Basic parameters
C
COMMON/BASNUM/NATOM,NION,NLEVEL,NTRANS,ND,NFREQ,NFREQC,NFREQE,
* IOPTAB,IDISK,IZSCAL,IDMFIX,IHESO6,IFMOL,IFENTR,
* NFREQL,NLEV0,ICOLHN,IOSCOR,ILGDER,IFRYB,IFRSET,
* NFREAD,NELSC,NTRANC,IOVER,JALI,IBC,IUBC,
* IRDER,ILMCOR,IFDIEL,IFALIH,IFTENE,ITNDRE,
* ILPSCT,ILASCT,IRTE,IDLTE,IBFINT,INTRPL,ICHANG,
* NATOMS,IPSLTE,ISPODF,ITLUCY,NRETC,IFRAYL,IFPRAD
COMMON/INPPAR/TEFF,GRAV,
* YTOT(MDEPTH),WMM(MDEPTH),WMY(MDEPTH),
* xmstar,xmdot,rstar,alpha0,
* QGRAV,EDISC,DZETA,RELDST,
* visc,zeta0,zeta1,dmvisc,fractv,
* omeg32,wbarm,wbar,alphav,pgas0,
* bergfc,
* ISPLIN,IRSPLT,ivisc,ibche,LTE,LTGREY,LCHC,LRESC
COMMON/MATKEY/NN,NN0,INHE,INRE,INPC,INSE,INZD,INMP,NDRE,insel
COMMON/INVINT/XI2(NLMX),XI3(NLMX)
COMMON/RUNKEY/CHMAX,ITER,NITER,NITZER,INIT,LAC2,LFIN
COMMON/CONKEY/HMIX0,crflim,
* NCONIT,ICONV,INDL,IPRESS,ITEMP,ICBEG,
* itmcor,iconre,ideepc,ndcgap,IDCONZ
COMMON/OPCKEY/NCON,IRSCT,IOPHMI,IOPH2P,IOPHE1,IOPHE2,
* IOPHL1,IOPHL2,IPHE2C,IFMOFF
COMMON/PRINTS/IPRINT,IPRING,IPRIND,IPRINP,ICOOLP,ICHCKP,
* IPOPAC,IPRINI
COMMON/PSILIM/DPSILG,DPSILT,DPSILN,DPSILD
COMMON/CENTRL/ZND,IFZ0
c
COMMON/ANGLES/AMU(MMU),WTMU(MMU),FMU(MMU),NMU
COMMON/ANGEXT/AMUEX(MANGEX),WMUEX(MANGEX),NANGEX
c
common/comptn/amuc(mmuc),wtmuc(mmuc),
* amuc1(mmuc),amuc2(mmuc),amuc3(mmuc),
* amuj(mmuc),amuk(mmuc),amuh(mmuc),amun(mmuc),
* calph(mmuc,mmuc),cbeta(mmuc,mmuc),
* cgamm(mmuc,mmuc),RADZER,FRLCOM,
* SIGEC(MFREQ),ijorig(mfreq)
common/angnum/nmuc
common/compti/nedd,nsti,islab,ilbc,icompt,icomst,icomde,
* icombc,icmdra,knish,itcomp,icomve,icomrt,
* ichcoo,icomgr
common/comite/ncfor1,ncfor2,nccoup,ncitot,ncfull
common/mlcons/aconml,bconml,cconml
common/taursl/taurs(mdepth)
IMPLICIT REAL*8 (A-H,O-Z), LOGICAL*1 (L)
PARAMETER (MITER = 200,
* MLAMBD = 100)
COMMON/LAMBDA/NLAMBD,IFFIX(MITER),NETEXP(MITER),NETFIX(MITER),
* IETEXP(MITER,MLAMBD),IETFIX(MITER,MLAMBD),
* NITLAM(MITER+1),IELCOR
COMMON/CHNMAT/INHE0(MITER),INRE0(MITER),INPC0(MITER),
* INDL0(MITER),INSE0(MITER),INMP0(MITER),
* NN00(MITER),NDRE0(MITER),LCHMAT,LIROST
COMMON/ACCEL/ORELAX,ITEK,IACC,IACC0,IACD,KANT(MITER),KSNG,
* LSNG(MTOT),LASO,LRES2
COMMON/ACCLP/ILAM,IACPP,IACC0P,IACDP,LAC2P
COMMON/ACCLT/IACLT,IACLDT
COMMON/CHNAD/CHMAXT,NLAMT,ILDER,IBPOPE
C
c PARAMETER (MFREQP = 300000)
PARAMETER (MFREQP = MFREQ)
PARAMETER (MLINH = 78,
* MHT = 7,
* MHE = 20,
* MHWL = 54)
INTEGER*2 ITRLIN
REAL*4 PRFLIN,BFCS
real*8 absopac(mtabt,mtabr,mfrtab)
C
COMMON/MODPAR/DM(MDEPTH),TEMP(MDEPTH),ELEC(MDEPTH),DENS(MDEPTH),
* ANTO(MDEPTH),ANMA(MDEPTH),ANH1(MDEPTH),ZD(MDEPTH),
* HKT1(MDEPTH),TK1(MDEPTH),HKT21(MDEPTH),SQT1(MDEPTH),
* TEMP1(MDEPTH),ELEC1(MDEPTH),DENS1(MDEPTH),
* DENSI(MDEPTH),DENSIM(MDEPTH),
* ELSCAT(MDEPTH),ALAB(MDEPTH),
c * ENRG(MDEPTH),
* DELDM(MDEPTH),DEDM1,DELDMZ(MDEPTH),
* THETAV(MDEPTH),VISCD(MDEPTH),DALPMX,XHYD,
* DMTOT,RRDIL,TEMPBD,ALPTAV,ALPGAV,NALP,IBETA
COMMON/LEVPOP/POPUL(MLEVEL,MDEPTH),BFAC(MLEVEL,MDEPTH),
* POPINV(MLEVEL,MDEPTH),
* POPGRP(MLEVEL),
* POP(MLEVEL),SBF(MLEVEL),DSBF(MLEVEL),USUM(MION)
COMMON/POPZR0/POPZER,POPZR2,POPZCH,RPOP0(MLVEXP,MDEPTH),
* IPZERO(MLEVEL,MDEPTH),IGZERO(MLVEXP,MDEPTH)
COMMON/CRSWPS/CRSW(MDEPTH),SWPFAC,SWPLIM,SWPINC,ICRSW
COMMON/WMCOMP/WNHINT(NLMX,MDEPTH),WNHEII(NLMX,MDEPTH),
* wop(mlevel,mdepth),ifwop(mlevel)
COMMON/REPART/REINT(MDEPTH),REDIF(MDEPTH),TAUDIV,IDLST
COMMON/TURBUL/VTURB(MDEPTH),VTURBS(MDEPTH),VTB,IPTURB
C
COMMON/LEVREF/SBPSI(MLEVEL,MDEPTH),SBLPSI(MLEVEL,MDEPTH),
* DSBPST(MLEVEL,MDEPTH),DSBPSN(MLEVEL,MDEPTH),
* ILTREF(MLEVEL,MDEPTH),IGUIDE(MLEVEL),
* ILTERF(MLEVEL,MDEPTH)
COMMON/LEVFIX/PT(MLEVEL,MDEPTH),PN(MLEVEL,MDEPTH),
* PP(MLEVEL,MDEPTH)
COMMON/LEVADD/USUMS(MION,MDEPTH),
* DUSMT(MION,MDEPTH),
* DUSMN(MION,MDEPTH),
* DIESIG(MION,MDEPTH)
COMMON/MRGPAR/SGM0(MMER),
* FRCH(MMER),
* SGEXT1(MMER,MDEPTH),
* GMER(MMER,MDEPTH),
* SGMSUM(NLMX,MMER,MDEPTH),
* SGMSUD(NLMX,MMER,MDEPTH),
* SGMG(MMER,MDEPTH),
* IMRG(MLEVEL),
* IIMER(MMER)
COMMON/UPSUMS/DUSUMT(MION),DUSUMN(MION)
C
COMMON/DWNPAR/ELEC23(MDEPTH),
* ACOR(MDEPTH),
* Z3(MZZ),
* DWC1(MZZ,MDEPTH),
* DWC2(MDEPTH),
* DWF1(MMCDW,MDEPTH),
C * DWFL(MFREQL,MDEPTH,MMCDW),
* MCDW(MTRANS),ITRCDW(MMCDW),NCDW
C
COMMON/OBFPAR/ITRBF(MBF)
C
COMMON/GFFPAR/GF0(MDEPTH),GF1(MDEPTH),GF2(MDEPTH),
* GF3(MDEPTH),GF4(MDEPTH),GF5(MDEPTH),GF6(MDEPTH),
* GF0D(MDEPTH),GF1D(MDEPTH),GF2D(MDEPTH),
* GF3D(MDEPTH),GF4D(MDEPTH),GF5D(MDEPTH),
* GF6D(MDEPTH),DELTT(MDEPTH)
C
COMMON/OFFPAR/SFF3(MION,MDEPTH),
* SFF2(MION,MDEPTH),
* DSFF(MION,MDEPTH),
* CFFN(MDEPTH),CFFT(MDEPTH)
C
COMMON/OTRPAR/ABTRA(MTRANS,MDEPTH),EMTRA(MTRANS,MDEPTH),
* DEMLT(MTRANS,MDEPTH)
C
COMMON/CURRNT/XKF(MDEPTH),XKF1(MDEPTH),XKFB(MDEPTH)
COMMON/CUROPA/ABSO1(MDEPTH),EMIS1(MDEPTH),SCAT1(MDEPTH),
* ABSOT(MDEPTH),ABSOE1(MFREX),
* EMEL1(MDEPTH),ABSO1L(MDEPTH),EMIS1L(MDEPTH),
* ABSOPR(MDEPTH),EMISPR(MDEPTH)
COMMON/TOTRAD/RAD(MFREQ1,MDEPTH),FHD(MFREQ),
* FAK(MFREQ1,MDEPTH),RADK(MFREQ1,MDEPTH),
* EXTRAD(MFREQ),EXTINT(MFREQ,MMU),HEXTRD(MFREQ),
* TRAD,WDIL,EXTOT,TSTAR
COMMON/CURRAD/RAD1(MDEPTH),ALI1(MDEPTH),FAK1(MDEPTH),
* radcm(mfreq,mdepth),
* RADL(MFREQL,MDEPTH),ABSALI(MFREQL,MDEPTH),
* alih1(mdepth)
COMMON/CURDER/DABT1(MDEPTH),DEMT1(MDEPTH),
* DABN1(MDEPTH),DEMN1(MDEPTH),
* DABM1(MDEPTH),DEMM1(MDEPTH),
* DABX1(MDEPTH),DEMX1(MDEPTH),
* DABP1(MLEVEL,MDEPTH),DEMP1(MLEVEL,MDEPTH),
* DRCH1(MLEVEL,MDEPTH),DRET1(MLEVEL,MDEPTH),
* ABSFF(MDEPTH),DABFT(MDEPTH),DABFN(MDEPTH),
* DSFDT(MDEPTH),DSFDN(MDEPTH),DSFDM(MDEPTH),
* DSFDP(MLVEXP,MDEPTH),
* DSFDTM(MDEPTH),DSFDNM(MDEPTH),
* DSFDPM(MLVEXP,MDEPTH),
* DSFDTP(MDEPTH),DSFDNP(MDEPTH),
* DSFDPP(MLVEXP,MDEPTH)
COMMON/CURTRI/ALIM1(MDEPTH),ALIP1(MDEPTH)
COMMON/EXPRAF/RADEX(MFREX,MDEPTH),FAKEX(MFREX,MDEPTH)
C
COMMON/TOTPRF/PRFLIN(MDEPTH,MFREQP)
COMMON/TOTFLX/FLTOT(MDEPTH),FLFIX(MDEPTH),FLEXP(MDEPTH),
* FCOOL(MDEPTH),FCOOLI(MDEPTH),FLRD(MDEPTH),
* FPRAD(MDEPTH),GRAD(MDEPTH),FPRD(MDEPTH),
* GRADF(MDEPTH,MFREQ)
C
COMMON/RRATES/RRU(MTRANS,MDEPTH),RRD(MTRANS,MDEPTH),
* DRDT(MTRANS,MDEPTH)
COMMON/RRTOFF/RDDP(MTRAN3,MDEPTH),RDDM(MTRAN3,MDEPTH)
COMMON/CRATES/COLRAT(MTRANS,MDEPTH)
C
COMMON/FRQALL/FREQ(MFREQ),W(MFREQ),PROF(MFREQ),WCH(MFREQ),
* JIK(MFREQ),IJX(MFREQ),IJBF(MFREQ),IFS0,
* KIJ(MFREQ),LSKIP(MDEPTH,MFREQ)
COMMON/FRQINT/FRCMAX,FRCMIN,FRLMAX,FRLMIN,CFRMAX,DFTAIL,
* TSNU,VTNU,DDNU,CNU1,CNU2,IELNU,NFTAIL
COMMON/LINOVR/NITJ(MFREQ),IJLIN(MFREQ),ITRLIN(MITJ,MFREQ)
COMMON/LINFRQ/NLINES(MFREQ)
COMMON/PHOEXP/AIJBF(MFREQ),BFCS(MCROSS,MFREQC),IFREQB(MFREQC)
COMMON/FREAUX/W0E(MFREQ),BNUE(MFREQ),WC(MFREQ),
* IJTC(MTRANS),IJALI(MFREQ),IJEX(MFREQ),IJFR(MFREX)
COMMON/COMPIF/LINEXP(MTRANS)
C
COMMON/SURFAC/FLUX(MFREQ),FH(MFREQ),Q0(MFREX),UU0(MFREX)
C
COMMON/FILES/PSY0(MTOT,MDEPTH),PSY1(MTOT,MDEPTH),
* PSY2(MTOT,MDEPTH),PSY3(MTOT,MDEPTH)
C
COMMON/PRESSR/PTOTAL(MDEPTH),PGS(MDEPTH),PRADT(MDEPTH),
& PRADA(MDEPTH)
COMMON/HEQAUX/PRD0,IHECOR
COMMON/OPMEAN/ABROSD(MDEPTH),SUMDPL(MDEPTH),
& ABPLAD(MDEPTH),ABPMIN
COMMON/CHARFX/QFIX(MDEPTH)
COMMON/WINDBL/ALBE(MFREQ),IWINBL
COMMON/DFEALI/DJMAX,NTRALI
COMMON/OPACAD/ABAD,EMAD,SCAD,DAT,DAN,DET,DEN,DST,DSN,DDN(MLEVEL)
COMMON/MODCON/FLXC(MDEPTH),DELTA(MDEPTH)
COMMON/RESDER/RSAT,RSBT,RSAN,RSBN,RSAX(MLEVEL),RSBX(MLEVEL)
COMMON/GRAYTS/TAUROS(MDEPTH),TAUFLX(MDEPTH),
* TAUTHE(MDEPTH),THETA(MDEPTH)
COMMON/TAURSS/TROSS(MDEPTH)
COMMON/TCONST/TDISK,ITCONS
C
COMMON/HYDADD/PHMOL(MDEPTH),ANEREL,IHM,IH2,IH2P
COMMON/ELDNSP/ANP,AHTOT,AHMOL
C
COMMON/RRVALS/RR(30,30),ABNDD(30,MDEPTH),ENEV(30,30),IONIZ(30),
* IREF,IREFA,LGR(30),LRM(30)
COMMON/STATEP/Q,QM,DQT,DQN,DQM,ENER,QREF,DQTR,DQNR,PFHYD
COMMON/ODFCHT/CHANT(MDEPTH)
COMMON/STRAUX/XK0(MLINH),XK,DBETA,BETAD,ADH,DIVH
COMMON/STDPAR/ELSTD,IDSTD
c
c parameters for hydrogen Stark broadening tables
c
COMMON/HYDPRF/PRFHYD(MLINH,MHWL,MHT,MHE),
* WLHYD(MLINH,MHWL),
* WLH(MHWL,MLINH),
* XTLEM(MHT,MLINH),
* XNELEM(MHE,MLINH),
* NWLHYD(MLINH),
* NWLH(MLINH),
* NTH(MLINH),
* NEH(MLINH),
* ILINH(4,22),
* IHYDPR
C
COMMON/LTEGRP/TAUFIR,TAULAS,ABROS0,TSURF,ALBAVE,DION0,
* DM1,ABPLA0,NNEWD,
* NDGREY,IDGREY
COMMON/COMPTF/DLNFR(MFREQ),BNUS(MFREQ),
* CDER10(MFREQ),CDER1P(MFREQ),CDER1M(MFREQ),
* CDER20(MFREQ),CDER2P(MFREQ),CDER2M(MFREQ),
* DELJ(MFREQ,MDEPTH)
COMMON/VISPAR/TVISC(MDEPTH),DTVIST(MDEPTH),DTVISN(MDEPTH),
* DTVISR(MDEPTH)
C
c parameters for the opacity tables
c
COMMON/TABLOP/ FRTB1,FRTB2,RTAB1,RTAB2,TTAB1,TTAB2
common/numbopac/ numfreq,numrho,numtemp
common/vectors/ tempvec(mtabt),rhovec(mtabr)
common/opacities/frtab(mfrtab),absopac
common/raytbl/raytab(mtabt,mtabr),raysc(mdepth)
common/binopa/ibinop
PARAMETER (MFODF = 180,
* MDODF = 3,
* MFRO = MFREQL,
* MHOD = 3,
c * MKULEV= 2,
c * MLINE = 2,
c * MCFE = 2 )
* MKULEV= 7000,
* MLINE = 2300000,
* MCFE =28000000 )
C
REAL*4 SIGFE
C
COMMON/ODFION/INODF1(MION),INODF2(MION),INBFCS(MION),
* IKOBS(MION)
C
COMMON/ODFCTR/FRODF(MLEVEL),NFRODF(MHOD),INDODF(MLEVEL),
* JNDODF(MTRANS)
COMMON/ODFFRQ/FROS(MFRO,MHOD),WNUS(MFRO,MHOD),XDO(3,MHOD),
* KDO(4,MHOD)
COMMON/ODFMOD/I1ODF(MLEVEL),I2ODF(MLEVEL),NQLODF(MLEVEL)
COMMON/ODFSTK/XKIJ(MHOD,NLMX),WL0(MHOD,NLMX),FIJ(MHOD,NLMX)
COMMON/SPLCOM/SIGFE(MDODF,MCFE),
& FRS1,FRS2,DXNU,XJID(MDEPTH),JIDI(MDEPTH),JIDR(MDODF),
& JIDS,JIDN,NFRS1,NFTT
COMMON/OPALIM/M1FILE(NLMX,MHOD),M2FILE(NLMX,MHOD),IMERG
COMMON/OPLIMT/ALLIM1,ABLIM1,ABLIM2,ABLIM3
COMMON/LEVCOM/EMKU(MLEVEL,2),YMKU(MLEVEL,2),
& XEV(MLEVEL,MION),XOD(MLEVEL,MION),EU(2*MLEVEL),
& JEN(2*MLEVEL),
& EEV(MKULEV),AEV(MKULEV),SEV(MKULEV),WEV(MKULEV),
& EOD(MKULEV),AOD(MKULEV),SOD(MKULEV),WOD(MKULEV),
& KSEV(MKULEV),KSOD(MKULEV),NEVKU(MION),NODKU(MION),
& NLEVKU,NLINKU,KEVE,KODD
C
C Parameters that specify dimensions of arrays
C
PARAMETER (MATOM = 10, ! maximum number of explicit atoms
* MION = 190, ! maximum number of explicit ions
* MLEVEL = 295, ! maximum number of explicit levels
* MLVEXP = 290, ! maximum number of explicit linearized levels
* MTRANS = 9000, ! maximum number of all transitions
* ! which are somehow taken into account
* MDEPTH = 70, ! maximum number of depth points
* MDEPTC = 1, ! maximum number of depth points (Compton)
* MFREQ = 2000, ! maximum number of frequency points
c * MFREQ = 1480, ! maximum number of frequency points
* MFREQC= 2000, ! maximum number of freq.points in continuum
* MFREX = 400, ! maximum number of linearized frequencies
* MMU = 3, ! maximum number of angle points
* MMUC = 10, ! maximum number of angle points (Compton)
* MFRCOR = 1, ! maximum number of frequency points
* ! in one line for correction weights
* MFETA = 1, ! maximum number of frequency points
* ! in one line for ETA procedure (outdated)
* MCROSS=R 800, ! maximum number of photoionization
* ! cross-sections + related quantities
* MBF = 800, ! maximum number of b-f transitions
* MFIT = 225, ! maximum number of fit points (OP b-f c.s)
* MITJ= 220, ! maximum number of overlapping transitions
* MMCDW = 14, ! maximum number of levels with pseudocont.
* MMER = 2, ! maximum number of merged levels
* MZZ = 10, ! maximum charge for occup.prob. ions
* MLINFR = 55, ! maximum number of lines at given freq.
* MFREQL = 2000, ! maximum number of frequencies per line
* MVOIGT= 175, ! maximum number of lines with Voigt
* ! profile
* NLMX = 80, ! highest hydrogenic level considered
* MSMX = 1, ! size of largest matrix kept in memory in SOLVE
* MTOT = 300, ! maximum number of model unknowns
* ! (ie. maximum dimension of vector PSI)
* MFREQ1= 1, ! =1 for ISPLIN<5; =MFREQ otherwise
c * MFREQ1=MFREQ,
* MFRTAB= 2, ! maximum number of freqeuncies in opac.table
* MTABT = 50, ! max.number of temperatures in opac.table
* MTABR = 50, ! max.number of densities in opac.table
* MCOND= 1, ! dummy
* MLEVE3= 1, ! =1 for diag.precond; =MLEVEL for tridia.prec.
* MLVEX3=MLVEXP, ! =1 for diag.oper.; =MLVEXP for tridiag. .
* MTRAN3=MTRANS) ! =1 for diag.oper; =MTRANS for tridiag.
C
PARAMETER (MAROWS = MFREX + 2, ! maximum number of rows of matrix
* ! A of complete linearization
* MCCOLS = MLEVEL+ 5) ! maximum number of columns of
* ! matrix C of compl. linearization
C
C Basic physical constants
C
PARAMETER (H = 6.6256D-27, ! Planck constant h
* BOLK = 1.38054D-16, ! Boltzmann constant k
* HK = 4.79928144D-11, ! h/k
* CAS = 2.997925D18, ! light speed c (A/s)
* EH = 2.17853041D-11, ! ionizaton energy of hydrogen
* BN = 1.4743D-2, ! 2*h/c**3, c -light speed
* SIGE = 6.6516D-25, ! Thomson scattering c-s
* SIG4P = 4.5114062D-6, ! Stefan-Boltzmann const/4pi
* PI4H = 1.8966D27, ! 4pi/h
* PCK = 4.19168946D-10, ! 4pi/c
* HMASS = 1.67333D-24) ! mass of hydrogen atom
C
C Basic mathematical constants
C
PARAMETER (UN = 1.0D0,
* HALF = 0.5D0,
* TWO = 2.0D0)
C
C Unit number
C
PARAMETER (IBUFF = 95)
C
C Basic parameters
C
COMMON/BASNUM/NATOM,NION,NLEVEL,NTRANS,ND,NFREQ,NFREQC,NFREQE,
* IOPTAB,IDISK,IZSCAL,IDMFIX,IHESO6,IFMOL,IFENTR,
* NFREQL,NLEV0,ICOLHN,IOSCOR,ILGDER,IFRYB,IFRSET,
* NFREAD,NELSC,NTRANC,IOVER,JALI,IBC,IUBC,
* IRDER,ILMCOR,IFDIEL,IFALIH,IFTENE,ITNDRE,
* ILPSCT,ILASCT,IRTE,IDLTE,IBFINT,INTRPL,ICHANG,
* NATOMS,IPSLTE,ISPODF,ITLUCY,NRETC,IFRAYL,IFPRAD
COMMON/INPPAR/TEFF,GRAV,
* YTOT(MDEPTH),WMM(MDEPTH),WMY(MDEPTH),
* xmstar,xmdot,rstar,alpha0,
* QGRAV,EDISC,DZETA,RELDST,
* visc,zeta0,zeta1,dmvisc,fractv,
* omeg32,wbarm,wbar,alphav,pgas0,
* bergfc,
* ISPLIN,IRSPLT,ivisc,ibche,LTE,LTGREY,LCHC,LRESC
COMMON/MATKEY/NN,NN0,INHE,INRE,INPC,INSE,INZD,INMP,NDRE,insel
COMMON/INVINT/XI2(NLMX),XI3(NLMX)
COMMON/RUNKEY/CHMAX,ITER,NITER,NITZER,INIT,LAC2,LFIN
COMMON/CONKEY/HMIX0,crflim,
* NCONIT,ICONV,INDL,IPRESS,ITEMP,ICBEG,
* itmcor,iconre,ideepc,ndcgap,IDCONZ
COMMON/OPCKEY/NCON,IRSCT,IOPHMI,IOPH2P,IOPHE1,IOPHE2,
* IOPHL1,IOPHL2,IPHE2C,IFMOFF
COMMON/PRINTS/IPRINT,IPRING,IPRIND,IPRINP,ICOOLP,ICHCKP,
* IPOPAC,IPRINI
COMMON/PSILIM/DPSILG,DPSILT,DPSILN,DPSILD
COMMON/CENTRL/ZND,IFZ0
c
COMMON/ANGLES/AMU(MMU),WTMU(MMU),FMU(MMU),NMU
c
common/comptn/amuc(mmuc),wtmuc(mmuc),
* amuc1(mmuc),amuc2(mmuc),amuc3(mmuc),
* amuj(mmuc),amuk(mmuc),amuh(mmuc),amun(mmuc),
* calph(mmuc,mmuc),cbeta(mmuc,mmuc),
* cgamm(mmuc,mmuc),RADZER,FRLCOM,
* SIGEC(MFREQ),ijorig(mfreq)
common/angnum/nmuc
common/compti/nedd,nsti,islab,ilbc,icompt,icomst,icomde,
* icombc,icmdra,knish,itcomp,icomve,icomrt,
* ichcoo,icomgr
common/comite/ncfor1,ncfor2,nccoup,ncitot,ncfull
common/mlcons/aconml,bconml,cconml
common/taursl/taurs(mdepth)
C
INCLUDE 'fpvm3.h'
C
PARAMETER (MHOST = 30,
& MFRPAR = 2000)
REAL*4 TIMNOD
CHARACTER NODNAM*30
COMMON/VMNODE/TIMNOD(MHOST),NODNAM(MHOST),
& KTID(MHOST),KTIDS(MHOST),IKH(MHOST),
& NIHOS(MHOST),
& NTASK,NHOST,NARCH,KENCOD,IPVMON
pro allch
;
chan,'ca',19
chan,'fe',25
chan,'mg',11
chan,'ne',9
chan,'ni',27
chan,'s',15
chan,'si',13
;
return
end
****** Levels
3.98503e+15 1. 3 ' ar1 ' 0 0. 0
****** Continuum transitions
1 1002 1 15 99 0 0 0.000E+00 0.000E+00
2.892E+01
8 1 2
1.177
2
0.0047
0.08