***************************************************************** * * * ****** * * ******* ******* * * *** * * * * * * * * * * * ** * * * * * * * * * * * * * * * * * * * * * ****** ******* * * ***** * * * * * * * * * * * * * * * * * * * * * * * * * * * * ** * * * * * * * * ******* ******* * * *** * * * * ======================================================== * * version 15.04.00E * * January 2008 * * * * nova, supernova, stellar & planetary model atmospheres * * * * by * * * * Peter H. Hauschildt * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * E. Baron * * Dept. of Physics and Astronomy * * University of Oklahoma * * Norman, OK 73019-0261 * * * * and * * * * France Allard * * CRAL/ENS * * 46, Allee d'Italie * * 69634 Lyon Cedex 07 * * France * * * * and * * * * Andreas Schweitzer * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Jason P. Aufdenberg * * National Optical Astronomy Observatory * * 950. N . Cherry Ave * * Tucson, AZ 85719 * * USA * * * * and * * * * C. Ian Short * * Department of Physics * * Florida Atlantic University * * Boca Raton, FL * * USA * * * * and * * * * Jason W. Ferguson * * Dept. of Physics * * Wichita State University * * Wichita, Kansas * * * * and * * * * Eric J. Lentz * * JIHIR/Physics Division * * Oak Ridge National Laboratory * * Oak Ridge, TN 37831-6354 USA * * and * * Department of Physics and Astronomy * * University of Tennessee * * * * and * * * * Travis S. Barman * * Planetary Research Center * * Lowell Observatory * * Flagstaff, Arizona * * * * and * * * * Darko Jevremovic * * Dept. of Physics and Astronomy * * University of Oklahoma * * Norman, OK 73019-0261 * * * * and * * * * Derek Homeier * * Universitaetssternwarte Goettingen * * Geismarlandstrasse 11 * * 37083 Goettingen * * Germany * * * * and * * * * Alexander Petz * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Birgit Fuhrmeister * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Sebastian Knop * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Michael Schuelke * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Matthias Dehn * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Alexander Christopher Wawrzyn * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Christine Johnas * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Daan van Rossum * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * * * * and * * * * Simon Huegelmeyer * * Universitaetssternwarte Goettingen * * Geismarlandstrasse 11 * * 37083 Goettingen * * Germany * * * * and * * * * Dennis Jack * * Hamburger Sternwarte * * Gojenbergsweg 112 * * 21029 Hamburg * * Germany * ***************************************************************** *-- notes: * this file contains the main program (phoenix) and some * administrative subroutines (mostly MPI related) * for a more complete description, get the 'manual' from PHH * or from http://dilbert.physast.uga.edu/~yeti * *-- personal note: * i'm trying to write this code as structured as possible, but * while reading or using it, never forget the following: ****************************************************************** * 'Real Programmers Don't Use Pascal!' (E. Post, 1982) * * and * * 'Scientific Progress Goes 'boink'?!' (Waterson, 1991) * * and * * 'Out Out!! You demons of stupidity!!' (S. Adams) * ****************************************************************** * and, of course, the following describes the spirit of this * * code very well: * ****************************************************************** * You've asked for it, you've GOT it! * ****************************************************************** * *-- last changes: * *15.04.00E: patches to allow compile with Xcode 2.4+, gfortran, * Linux/ifort/x86_64 (nathan, from Andy). * (11/Jan/2008, PHH) *15.04.00D: more 2D improvments, (rho,z) mode declared * operational, added Q_andD mode flag, added * variable_stepsize, fixed s3main compilation, * added Simon's latest disk stuff! * (19/Jun/2007, 25 days to the wedding party, PHH) * (14/Jun/2007, 1 month to the wedding party, PHH) * (13/Jun/2007, 1 month to the wedding, PHH) *15.04.00C: merged with Knop's latest GR updates, NMS3 now fully * operational incl. NLTE. Added Chris' voigt branch * fix. * (05/Jun/2007, PHH) *15.04.00B: added Eddie's debug and KAPCAL patches, added initial * 2D (rho,z) mode, added Dennis' patches to fix cio.c, * added Eddie's format fix. * (04/Jun/2007, sister's birthday, PHH) *15.04.00A: fixed 2D mode, added Eddie's patches for build and * Compiler.options, added Simon's disk mode!, added * Eddie's fixes to vcapcnt, cas_dissolved and cas_sort.., * added Daan's new ultc_ad.f, completed 2D FS. * (14/May/2007, PHH) *15.03.03C: added Daan's new ul_ad version and his new treatment * of UL terms in NLTE. This can change results. * Completed and tested port to PS3 (CBE in 64bit PPE * mode with gfortran) * (04/Apr/2007, PHH) *15.03.03B: added grep20.f utility, fixed smp mode compile errors * for phoenix (smp is experimental), added Eddie's PWR5 * fixes. Added Eddie's new -qhot to fix new IBM AIX * compiler bugs. initial support for CBE/PS3, fixed * sematics of read_chrom3d.f, added support for gfortran * (on CBE), fixed numerous little sematics problems * (26/Mar/2007, PHH) *15.03.03A: fixed machine.ifc-mpi.f, added Daan's new ultc_ad.f, * new reflective BCs in 3DRT, fixed intensity init bug * in 3DRT/LC/tracker, fixed nasty buggy call to ions * from ppress (mode was not passed) * (06/Mar/2007, PHH) *15.03.02B: updated default.mki * (22/Feb/2007, PHH) *15.03.02A: added use_restart_tstd option * (22/Feb/2007, PHH) *15.03.01A: added workaround for ifort/em64t bugs to 2ifc.pl, * updated OSX/Intel flags following Eddie's tip. * fixed verbose/debug/LTE bug in cas_voigt_lines. * added SP3=.true (64 bit clear) to ifc options. * added Costa's fixes for NLTE/external, fixed * misalignments, fixed /snmod/ common errors, * changed fipol to perform more stable interpolation * for constant y values, aligned gauss_lines and * voigt_lines to 8 bytes. * (13/Feb/2007, PHH) *15.03.00D: added Daan's latest active damping updates, 3DRT * additions, s3main.f changed. * (30/Jan/2007, PHH) *15.03.00C: added Eddie's Compiler.options (Xcode 2.4.1 support) * fixed openMP mode for 3DRT (it is still slower than serial) * fixed build bug for pw4-64bit-mpi mode, added Simon's * build and Compiler.options for G95 on Intel-Macs, * fixed absoft/OSX/ppc64 compiler options, fixed em64t * ifort compiler options, added Daan's active damping for UL, * added Eddie's makehydro.f, 2D fixes * (24/Jan/2007, PHH) *15.03.00B: minor cosmetics, 3DRT additions. fixed ancient UL * bug found by Daan! fixed s3main bug for 2lvl comparison * added Eddie's allocated-mode patches. fixed bug in * pve_grid%flux computation in MPI mode * (02/Jan/2007, PHH) *15.03.00A: restored behavior of original UL method for compat, * added Daan's latest improvements to active UL * (22/Dec/2006, PHH) *15.02.03B: fixed compat. problem with new Tcorr. * (21/Dec/2006, PHH) *15.02.03A: first periodic BC (PBC) version of 3DRT, Lstar support * for PBCs, OS_iteration support (Jordan and GS), changed * opacity averaging in 3DRT/LC to match S3R2T (tested with * PBCs), added Daan's active UL temperature correction, * update to Daan's active UL Tcorr * (20/Dec/2006, PHH) *15.02.02J: updates to s3main.f, updated Eddie's support * for Mac/Intel to autodect and optimizations. added * test_smp.f to test openMP stuff, updated SMP options, * added new 'tracking' LC solver. extended qtab to (0:..) * to prevent trouble if species are not available in kapcal. * added Cassandra code to check if an ion is in EOS to prevent * disasters. * (22/Nov/2006, PHH) *15.02.02I: fixed output in 3DRT, added Grotrian diagrams for all * atoms to the manual (file size increase!), added Eddie's * 3DRT/SC patches, added Derek's updates * (20/Oct/2006, PHH) *15.02.02H: fixed nthdealloc call returning error on iitac * (03/Oct/2006, PHH) *15.02.02G: added Eddie's patches for NERSC profiling, added * support for iitac with Pathscale/MPI * (03/Oct/2006, PHH) *15.02.02F: added LOADPTS to 3DRT/Makefile.in * (29/Sep/2006, PHH) *15.02.02E: changed allocatable to pointer in nms3_module.f * to work around old compilers, added scripts for * easy compile with pathscale compilers * (29/Sep/2006, PHH) *15.02.02D: updated Eddie's updates for non-thermal processes * to in_sn.150202 and nthmod.f * (28/Sep/2006, PHH) *15.02.02C: added Eddie's updates for non-thermal processes * to in_sn.150202 and nthmod.f * (27/Sep/2006, PHH) *15.02.02B: made new PhoenixAtoms.150202 with He I and He II * line profiles for external model atoms, moved * profile files to PhoenixAtoms.150202! * (27/Sep/2006, PHH+TSB) *15.02.02A: added external model atom file select options for * H, He and Ne to replace phoenix internal data if * so desired. set default value of v0 to 0.0d0. * added Eddie's latest patches for Intel-Macs * (25/Sep/2006, PHH) *15.02.01A: added Matthias patches and updates for use with DRIFT, * fixed Compiler.options for Columbia * (19/Sep/2006, Michael H's birthday!, PHH) *15.02.00D: s3main updates for comp. with 3D, fixed grelmolb*.f * so that processing will continue after unknown line * code detected * (04/Sep/2006, PHH) *15.02.00C: merged with Knop's new GR routine versions. NMS3 mode * in phoenix now MPI parallelized. fixed NMS3 Makefiles * added Andy's 3DRT/Makefile.in patch * (31/Aug/2006, PHH) *15.02.00B: fixed Makefile for s3main. updated Compiler.options * and build for Intel V9.1 compilers on x86_64, fixed * input.f bug found by Travis. added new 1502 settings * for FeH_version and (new) CrH_version that set all * correction factors to unity (use gf as in master database) * (28/Aug/2006, PHH) *15.02.00A: implemented driver in LC_solver.f, added conversion code * and support for Bernath line lists (FeH, CrH, TiH), * fixed machine.macosx.f and machine.macosx-mpi.f to * remove dgeev and second stubs, fixed Compiler.options * according to Eddie. Added options qCrH_B,qTiH_B,qFeH_B * to select Bernath's level data for Q's * (default is .false.), fixed bug in LC_solver_general * (21/Aug/2006, PHH) *15.01.00C: added re_open to aix-mpi version (IBM fucked up gpfs) * LC_solver re-design begins, added tridiag Lstar to * LC_solver_general, added Eddie's updates for alloc mode, * added patched gendealloc.f from Eddie, added cj's * pertuber code=4 patches, added Andy's FreeBSD6 patch for * build and Compiler.options, expanded load_guess to * include *all* ionpot data with g=1 * (28/Jul/2006, PHH) *15.01.00B: added re_open to aix-mpi version (IBM fucked up gpfs) * LC_solver re-design begins, added tridiag Lstar to * LC_solver_general, added Eddie's updates for alloc mode, * added patched gendealloc.f from Eddie, added cj's * pertuber code=4 patches * (19/Jul/2006, PHH) * (04/Jul/2006, semi-final day, whooo!, PHH) *15.01.00A: added Knop's GR RT and NMS3 updates, allows LTE NMS3 * use in phoenix. * (08/Jun/2006, PHH) *15.00.00F: added error checks to read_pve, fixed Makefile.in for * regression, added Lstar save/read to 3DRT * (08/Jun/2006, PHH) *15.00.00E: maintenance and sync updates * (25/May/2006, PHH) *15.00.00D: fixed numfrac alloc bug * (24/May/2006, PHH) *15.00.00C: fixed to*.pl bugs * (24/May/2006, PHH) *15.00.00B: added AGS2005 abundances to default.mki and to * abundance printout * (23/May/2006, PHH) *15.00.00A: release as 15.0.0 with layer=128 (!) as default. * (23/May/2006, PHH) *14.03.00A: added code for new atomic line list conversion. fixed * bug in merge_list.f, added code to distinguish old * and new (derived from Kurucz) O I model atom data, * added more elements (up to Rg, see MISC/blockdata.f) * added new zusum with totally new data structure etc. * fixed allocatable version compile. Added Eddie's * fixes to build, fixed iwhprof bug in fuzzlinb.f * fixed dummy continua & LTE line bug, added Derek's * cas_profile patches, manual updates, added V9 compiler * for altix, fixed bug in P&B H I coll. rate routine, * zusum compat. fixes, 3DRT SC fixes, improved pressure * ionization treatment, added Eddie's SC mods, fixed xlf * options to 'signedzero' and changed xlf90 to xlf95, * chanded module-eqv. alloc-on files to circumvent * intel f95 namespace bugs (ifort v9 stinks), added Andy's * FreeBSD fix for superlu, fixed pathscale options, * added G95-MPI target lam711, 64 bit for big runs) * (22/May/2006, PHH) * (29/Mar/2006, solar eclipse in the rain version! PHH) *14.02.00D: flux vector computation in 3DRT, G95-Linux and G95-XP * targets (both essentially untested), added vcapcnt.f * patch from Matthias (for drift), added several attempts * to c_read_block in order to avoid fatal I/O errors. * added line list checks to grelcmbp.f and grelmolbp.f * to catch deadlocks in case of different line databases * on different processors. fixed dozusu_general auto- * ionizing level bug in qas found by Alex Petz * (06/Mar/2006, PHH) * (27/Feb/2006, Helau! Alaaf! PHH) *14.02.00C: g95 support (MacOSX, 32bit and 64bit), * started reflective BCs in 3DRT, added Eddie's pwr5 * fixes, added Matthias' changes for DRIFT mode * (17/Feb/2006, PHH) *14.02.00B: mk_columbia fixed, fixed 3DRT MPI setup, fixed * bug in 3drt result *output* * (30/Jan/2006, PHH) *14.02.00A: added Eddie's NERSC patches, added zusum stored data * limits of Birgit (this can change the results!), * added Birgit's 3D stuff, added Eddie's nthmod.f, * added Eddie's fix for ltelinbgauss.f, plenty of fixes * for kapcal_3D, fixed He- vec. opac bug. fixed bug in * alalbf.f, updated kapcal, removed otab (didn't work), * fixed almgbf.f, made "al*.f" and "val*.f" consistent, * fixed h2olin,hcnlin,c2h2lin for HITRAN2004 in phoenix.f, * added Eddie's pwr5 fix (chianti.f),added output to * cas_alloc * (23/Jan/2006, new year, new bugs!, PHH) *14.01.03A: added Eddie's/Knop's fix for nms3_frml_snl.f, added * Travis' support for NAS Columbia (MPI mode) * added Eddie's NERSC patches * (21/Dec/2005, it's the winter solstice!, PHH) *14.01.02C: improving SC_solver.f, added arb. wide Lstar to * SC_solver.f, added Eddie's new input.f * (16/Nov/2005, PHH) *14.01.02B: added Eddie's phoenix.f fixes (tau=1 point insert), * deactivated Intel ifort bug fix in 3DRT/OS_iteration.f * as Intel ifort version 9 does NOT seem to have the * problem anymore, fixed ppress.f f bug that messed * with eheu but should not have. * (11/Nov/2005, Karneval starts! Yuk. PHH) *14.01.02A: added Ian's new model atom files for Ba I-II and Sr I-II * this requires in_sn.140102 * (04/Nov/2005, PHH) *14.01.01A: added Eddie's build updates, added Frances compos.f * and gradad.f, added Matthias' showPhxDiffs.sh and * dot20.pl utilities * (04/Nov/2005, PHH) *14.01.00A: added France's improvements to settling mode * (17/Oct/2005, PHH) *14.00.00D: made new "mk_ibm_cines", set linalg=4 as default, * edited 2aix*.pl to compensate for IBM (essl still default) * cleaned up 3DRT directory, adding info output * to molstat.f, fixed H2O line indicator printout bug * discovered by France. * (13/Oct/2005, PHH) *14.00.00B: phoenix/3D updates. fixed Compiler.options for * FreeBSD target, added cas_profiles.f fixes from France, * fixed the fixes, added France's rainout.f, added * France's version of Travis' rainout.f * (30/Sep/2005, PHH) *14.00.00A: adding support for simple 3D/2D mode (no lines, yet), * added France's TiO fel updates: ames_TiO_version=1400 * and tio_shomate=.true. by default., added France's * sqtio.f update, added irrad mode to 3D.f, added * Alex' 2 photon process updates, added Costa's updates * for HITRAN2004 etc., added eos_3D.f, added Derek's * patches for cas_profiles.f * (21/Sep/2005, PHH) * (19/Sep/2005, Happy Birthday Du Wicht!, PHH) *13.14.00B: added line transfer restart to 3DRT, some Makefile * and library maintenance * (07/Sep/2005, PHH) *13.14.00A: added Travis' profile type (3) for full NLTE profiles * with density interpolation to cas_profiles.f * cas_profiles.f: changed default (no data) Lorentz core * to pure H2 to avoid trouble with new perturbers, * added Derek's updates to cas_profiles.f, added 2lvl line * transport to 2D and 3D RT, added irradiated mode to RT2D * fixed build for HLRN (fools and idiots!), added initial * support for ifort on Opterons in MPI mode * (31/Aug/2005, PHH) *13.13.01A: removed adgrad from opacity tables, added 2D restricted * version of 3DRT (in 2D.f) and fixed it (seems to work) * (11/Aug/2005, PHH) *13.13.00E: changed pi in physconst.inc to 20 digits * using Maple 9, added Eddie's Compiler.options updates, * line RT updates for 3DRT, hierarchical MPI parallelization * of 3DRT, fixed plenty of stupid little bugs in 2lvl.f * (03/Aug/2005, PHH) *13.13.00D: 3DRT arb. Lstar width updates, final version for * paper I * (18/July/2005, PHH) *13.13.00C: more 3DRT updates, fixed Altix support, fixed * nms3_solve_for_I_superLU_DIST.f bug found by Knop, * fixed s3main.f bug found by Knop * (06/July/2005, PHH) *13.13.00B: updated 3DRT, added support for Absoft 9.2 on OSX, * updated itanic2 (altix) compiler options, fixed MPI * BC issues in 3DRT, added wide LStar framework to 3DRT * (01/July/2005, PHH) *13.13.00A: added Birgit's contrib. function add-on, added * Costa's VO-Plez, updated 3DRT * (17/June/2005, PHH) *13.12.00D: added Eddie's option changes, more 3DRT updates * (16/June/2005, PHH) *13.12.00C: added Costa's molmass for 48TiO fix (63.0429->63.9429) * (11/May/2005, PHH) *13.12.00B: added OSX Tiger Accelerate routines to machine.macosx.f * fixed dust_version flag * (11/May/2005, PHH) *13.12.00A: added Jason F's new dust opacity updates (activated * by default), 3DRT updates, OSX Tiger support (alpha) * (02/May/2005, PHH) *13.11.02A: added Costa's patches for TiO(Plez), codes 193-197. * (25/Apr/2005, PHH) *13.11.01B: fixed double declararion of pb_h_coll_rates, added * MJS's odeint update, more 3DRT updates, added MJS's * cio.c typo fix. * (25/Apr/2005, PHH) *13.11.01A: added Jason A's new H I bb collision rate routines * using Przybilla & Butler data (option pb_h_coll_rates, * default is .TRUE.) * (19/Apr/2005, PHH) *13.11.00B: fixed bug in Jordan and Gauss-Seidel solvers of 3DRT, * fixed Makefile.in for debug mode. * (07/Apr/2005, PHH) *13.11.00A: added Matthias' patches for DRIFT interfaces, * more SLU patches for 3DRT, added Jordan and * Gauss-Seidel solvers for 3DRT * (01/Apr/2005, April fool's Ed. with foolish bugs! PHH) *13.10.00B: added Eddie's _SuperLU for 3DRT * (28/Mar/2005, PHH) *13.10.00A: added Eddie's patches for adaptive grid version * (21/Mar/2005: it's spring!, PHH) *13.09.02D: more 3DRT updates, changed default fmt_id in * default.mki, added Derek's test_spectrum_output.f * to directory tree * (21/Mar/2005, PHH) *13.09.02C: 3DRT full OS version update, added Derek's patches. * (04/Mar/2005, PHH) *13.09.02B: added Matthias's updates to default,mki and manual.tex * (14/February/2005, it's the big St. Valentines Day * Massacre yet again! PHH) *13.09.02A: added Matthias's routines to use C. Hellings discons. * added Travis' rainout fixes to csppress.f, rainout.f * and kapdust.f * (14/February/2005, it's the big St. Valentines Day * Massacre yet again! PHH) *13.09.01A: added He1_C6_fix: * variable to deactivate He I C6 computation. The computed * value of C6 for He I is ludicrous, therefore vdW broadening * can be turned off (C6=0) for He I lines (this is the new default) * This variable is also used in Cassandra (if He I is in NLTE) * (28/January/2005, PHH) *13.09.00B: added JPA's fix for input.f * (19/January/2005, PHH) *13.09.00A: added JWF's latest dust changes (more species, defaults), * added dust_version flag. Default is NEW mode (1309), set * to 1308 to recover OLD setup. Warning: default CHANGES * results of models with dust opacities. You *must* set * dust_version=1308 in the input namelist to recover old * behavior, 1308 will *ignore* any input usedust settings! * Added Derek's modified thread_bind for GWDG machine * removed time wasting timing calls in grelmolbp*.f * added Ian's chiant.f patch fixing Sr symbols. * (13/January/2005, PHH) *13.08.05B: fixes for ANU Alphaserver SC, added Eddie's Compiler.options * fix for LAM 7.1, added Alex' cas_solve fix * (31/December/2004, year's end edition, PHH) *13.08.05A: added Birgit's changes * (22/November/2004, PHH) *13.08.04C: added 3DRT subdir * (22/November/2004, PHH) *13.08.04B: added NMS3/nms3_solve_for_I_superLU_DIST.f serial fix * from mjs. Fixed serial link of s3main. Added fixed * SuperLU_Dist, added support for ia32/ifort/MPI, * added NMS3 parallelization stuff * (11/November/2004, HELAU! ALAAF! PHH) *13.08.04A: added use_rkqc switch to allow trickier hydro * integrations * (31/Oct/2004, PHH, Halloween edition with extra * scary & evil bugs!) *13.08.03C: added Eddie's latest SuperLU_dist patches, more patches * from Eddie, added support for Intel Compilers for ia32em * on x86_64 machines (serial), streamlined Makefile.in * and Compiler.options to use CC and CC_QD throughout. * (29/October/2004, PHH) *13.08.03B: fixed c6corr declaration in compmoldampnlte, updated * build and Compiler.options for ia32 Intel 8.1 compilers, * fixed djext, dkext bug in s3moms (ugly, man!) * (26/October/2004, PHH) *13.08.03A: made c6corr an input variable, default is * backward compatible * (15/October/2004, PHH) *13.08.02G: added Eddie's latest SuperLU_dist patches * (11/October/2004, PHH) *13.08.02F: added Intel MKL support, added Eddie's Compiler.options * patches. *13.08.02E: fixed block2cache_mol allocation bugs (revert to * static arras, added checks for bounds) * (04/October/2004, PHH) *13.08.02D: added Eddie's patches for MISC * (01/October/2004, PHH) *13.08.02C: added MPI mode for altix Itanium2 systems, splitted * bench in 3 parts to circumvent spurious compiler tricks * (01/October/2004, PHH) *13.08.02B: fixed line_util routines etc. made blocking pointers * allocatable for more flexibility. Adapted to current * OSF1/AXP machines. Added experimental MPI mode for * OSF1/AXP with vendor MPI/elan libs. Added support * for Altix systems (Itanium2) using Intel compilers * version 8.1 (includes full opt support for QD) * (24/September/2004, PHH) *13.08.02A: added support for BT (Barber-Tennyson) water line list, * updated molstat, showmlines, linestat and showlines * for Altivec use. * (11/September/2004, PHH) *13.08.01A: added Knop's all layer rf_out option (set lay_rf <= 0) * added Eric's alloc-mode updates, fixed initialization * bug in avhem_Vardya.f * (05/September/2004, PHH) * (31/August/2004, Qingdao special edition, PHH) *13.08.00B: fixed Makefile.in for showlines, fixed use_opac_table * bug, added output of gf and chi for ID mode * (18/August/2004, PHH) *13.08.00A: increased maxblock_atm, added Alex' new code to allow * tstd = tau(Ross) by reading and using opacity tables. * (28/July/2004, PHH) *13.07.04D: fixed array-slice access bug for ngauss=0 and nlor=0 * in cas_insert*lines.f * (22/July/2004, PHH) *13.07.04C: added Jason F's changes to He- and H2- opacities, * default is 'on' (see kapcal.f and vcapcnt.f) * (21/July/2004, PHH) *13.07.04B: added debug version of Jimmy's EOS for testing, added * Eddie's SuperLU_DIST changes. cas_solve_stateq skipped * for iter=0 in order to save spectrum computing time. * fixed qjanaf data handling for T>6000K to avoid * disasters. fixed linestat and showlines. * (17/July/2004, PHH) *13.07.04A: modified dump level 0 behavior. fixed denormalized * IEEE 754 numbers leading to disaster on Opterons but * working fine on G5's problem. added use_log_bi_interpol, * set default of use_log_bi_interpol=.true. * (27/June/2004, PHH) *13.07.03G: fixed NLTE level =0 load * for continua, added Alex' bug fixes for opac table mode * (23/June/2004, PHH) *13.07.03F: fixed NAG detected compile time bug in interfaces.f, * timing output more readable * (16/June/2004, PHH) *13.07.03E: added latest compos.f from France, added Derek's * massive improvements of cas_profiles (near wings now * treated as Voigt profiles with correct damping constants) * added Derek's patches to build, Compiler.options etc. * (11/June/2004, Happy Birthday France!, PHH) *13.07.03D: fixed data type bug of zetmol and vdsqml in * mpi_support.f (changed from real*8 to real*4 for * optimization). * (08/June/2004, Venus transit version!, PHH) *13.07.03C: added Birgit's last chromosphere changes! Fixed * prefac_cache_mode_mol bug in grelmol*.f * added latest compos.f from France, * (03/June/2004, PHH) *13.07.03B: added mk_seneca, also for local I/O modes. added * vectorized exp's to atomic line selection, * note: small linmax (4096) dramatically increase speeds * for MPI versions with slow interconnect! fixed * Makefile.in for line_util/molecular_lines *13.07.03A: changed LTE line cache strategy for improved * performance, in particular in MPI modes. * recommendation: use many small caches in MPI mode * so that each process can hold all data in memory. * example: 100M lines -> linmaxm=1024, ncache_mol_g=23000 * for 8 MPI processes. Needs about 1GB/process! changed * cache replacement strategy to round robin. note: * for unknown reasons, holding all mol_v data in memory * is substantially slower than having small cache sizes * and reading data from disk. This also (indirectly) affects * speed in mol_g routines! added prefac_cache_mode_mol with * default zero to speed up mol_gauss (??) line computation. * released for testing. * (28/May/2004, on the way to Lyon, PHH) * (29/May/2004, Happy Birthday Mom!, PHH) * (30/May/2004, PHH) *13.07.02F: fixed re_open bug for AIX-MPI mode. * (18/May/2004, PHH) *13.07.02E: limit wl range to Lagrange continuum to make * nova/SN/wind spectrum generation faster, * fixed grelcmbp bug found by Birgit. Added * patches from Ian for linux-f95 target. added * fix for grelmolb*.f out-of-bounds error found * by Birgit and Andy * (12/May/2004, PHH) *13.07.02D: added workaround to fix NAG f95 (4.2 and 5.0) * compiler bug manifesting itself in KAPCAL/cssibf.f * (30/Apr/2004, PHH) *13.07.02C: fixed init bug in cas_profiles, updated compiler * options for x86_64 in MPI mode, added mk_nathan * (28/Apr/2004, PHH) *13.07.02B: fixed bugs in FPPRESS2 (still not operational), * more work on interpolation-fort.20-routine. * (25/Apr/2004, PHH) *13.07.02A: fixed uselin=2 for VO (no VO bands either!), added * initial version of read-fort.20-with-interpolation, * fixed write(7,) statements that ignored the input * format (if I find out who did that, I am going to * be very unpleasant!!), added new EOS to code tree in * FPPRESS2, 'make' will now produce phoenix2 and csppress2 * that are linked against the new EOS. It does not * work yet. * (23/Apr/2004, PHH) *13.07.01E: changed cas_solve to (a) allow local out-of-range * detection and (b) use ocnlte as failure backup. * (14/Apr/2004, PHH) *13.07.01D: added Peter Nugent's NAG f95 v5 interface etc fixes * (11/Apr/2004, Happy Easter!, PHH) *13.07.01C: added new uselin mode: 0: JOLA, 1: line from database, * 2: ignore both for easier testing. Initial * parallelization of NMS3 phase 1. * (07/Apr/2004, PHH) *13.07.01B: removed quick exit from dozusu_general to fix potential * bi change problem, fixed kapcal bug for NLTE Rayleigh * scattering, fixed structure mode * (27/Mar/2004, PHH+EAB) *13.07.01A: set default arithmetic precision in cas_solve to * 256bit, can be changed with the input switch of course. * (13/Mar/2004, PHH) *13.07.00B: added support for features requested by Jason F. * added Birgit's chromosphere mode changes * fixed zusum for pure LTE mode * (13/Mar/2004, PHH) *13.07.00A: changed treatment of NLTE ionization in EOSs to * improve accuracy. Added Eddie's ultimate improvements * to the ionization treatment. * (11/Mar/2004, PHH) *13.06.05A: fixed bugs in grelmolbp's. Added vectorization support * and optimized line selection routines (30% faster on * G5's!!), fixed deallocation of not allocated memory for test * versions of grel*.f (only relevant for empty line lists), * added Birgit's ichrom2 mode fix, added Travis' Itanic2 * options fix, added optimizations to Cassandra routines, * added AIX HPMLIB monitoring calls to a few Cassandra routines * (must be activated by hand) and options. changed AIX * options to *always* compile in 64bit mode. support for * intel fortran/C version 8, converted to release version * of xlf/OSX, switch default to global precision selection * for cas_solve_stateq, fixed minor problems to allow * running in xlf/OSX debug mode, added timing output, * added Eddie's Compiler.options fixes, fixes for Itanium2 * (03/Mar/2004, PHH) *13.06.04E: removed real loop vars in phoenix.f (opac table setup) * and ligre.f for better f95 compliance. * (20/Jan/2004, PHH) *13.06.04D: added Andy's nms3_os_setup patch (timing printout), * added Andy's patches to make ifc version 8.0 work * better (dreaded piece of junk). * (27/Jan/2004, PHH) *13.06.04C: fixed local scratch file mode bug found by Travis * (26/Jan/2004, PHH) *13.06.04B: added Eddie's fix to Rayleigh scattering in kapcal.f * and vcapcnt.f * (23/Jan/2004, PHH) *13.06.04A: added Birgit's patches for colmass as independent * variable rather that ststd, added Athlon64 to build * discovery for Linux, cleaned up QD calls fuer X86 * archs, need to fix options for QD with ifc (-O0, done) * or to switch to g++ to compile QD on X86 boxes! * (20/Jan/2004, PHH) *13.06.03E: fixed R* printing bug in print_radiative_rates * (15/Jan/2004, PHH) *13.06.03D: added Eddie's patches to phoenix.f, input.f and * cas_solve..f * (12/Jan/2004, Happy New Year! PHH) *13.06.03C: added MPI_Get_processor_name for more info output * in MPI mode, updated mpitest.f * (24/Dec/2003, HO. HO. HO. PHH) *13.06.03B: added nlte_neg_ion_approx flag to approximate correct * treatment of negative ions for the NLTE case. This also * corrects qas for the neutral species by the same factor * (21/Dec/2003, PHH) *13.06.03A: added Alex's patches for ff opacs and Cassandra, fixed * prfem reset bug in cas_gauss/voigt_lines. fixed myerr * declarations for flush (SGI) calls, added Eddie's patches * (19/Dec/2003, PHH) *13.06.02D: added Eddie's rholaw patch of 10/Dec/2003 * (15/Dec/2003, PHH) *13.06.02C: added Eddie's patches of 02/Dec/2003, added fixed * options for Itanic2-mpi and ifc-mpi targets, added * Eddie's patches, added Andy's showlines Makefile patch. * added SGI 'flush' fix from Dmitry Monin * (04/Dec/2003, PHH) *13.06.02B: added local MPI mode for atomic lines (grelcmbp_local) * made alo_line_factor,alo_cnt_factor inputs, * added kill_line_alo_thres, restored line coupling * to both upper and lower level (this helps dramatically) * (30/Nov/2003, PHH) *13.06.02A: Added Eddie's S3R2T_alloc/s3com.inc, added Eddie's * NMS3 fixes, added max_log_bi_change to cas_solve, * set default nlte_line_mode=3, new chromosphere mode!, * replaced MPI_REAL8 with MPI_DP, which is now defined * as MPI_DOUBLE_PRECISION in phx-mpi.inc. Commented * debug code left in grelcmbp.f (nobody but me ever reads * the output files, eh??). Fixed MPI options for IFC and * OSX-XLF, added FreeBSD-IFC target from Andy, added * X86_64-NAG target for Opterons. * (10/Nov/2003, PHH) *13.06.01C: Added Travis' Itanic support (alpha version, compiler * seems to have a wealth of bugs). Improved Itanic support, * changed cio.c to include (needed for Itanic), * added latest NMS3 changes to tree, set alo_line_factor * and alo_cnt_factor to 1 (affects qd_ solver!), removed * report options for ifc (printout takes more time than * compiling!). * (14/Oct/2003, PHH) *13.06.01B: Changed NG setup in nms3_solve_2lvl routine for NMS3, * removed SuperLU info output, fixed algamtn call sequence * bug found by EJL * (07/Oct/2003, PHH+EAB) *13.06.01A: new nms3_solve_rt driver routine for NMS3, * added 2lvl atom test solver to NMS3, changed to * use new optical.30Sept2003.d (updated refs) in new * in_sn.130601 directory * (04/Oct/2003, PHH+EAB) *13.06.00F: Changed compiler options for OSX/XLF to use * gcc 3.3 libraries. changed mpi_real4 to mpi_real in * csppress.f * (01/Oct/2003, PHH) *13.06.00E: call algamtn for mol=0 to avoid massive trouble. * (24/Sep/2003, PHH) *13.06.00D: added mk_ scripts for OSX, fixed 2macosx-xlf.pl to * work around csekund problems, fixed xlC options for QD. * (24/Sep/2003, PHH) *13.06.00C: added Derek's patch to machine.macosx.f * (19/Sep/2003, Happy Birthday MH! PHH) *13.06.00B: changed cas_line select call to work around OSX-xlf bug * (18/Sep/2003, PHH) *13.06.00A: added France's EOS patches, added Eddie's NMS3 updates, * added support for xlf on OSX, converted csppress to * namelist input (see csppress.mki for an example). * finished xlC options for QD on OSX/XLF. Updated bench.txt * with first results for G5. Added Derek's new profile * convolution for cassandra's special profiles. Fixed high-pressure * ionization reduction problem in ionsah and sahnlt. Added * ATLAS-BLAS libs to xlf for OSX * (17/Sep/2003, PHH) *13.05.02B: added Eddie's bug fix to NMS3. Updated IFC debug * compiler options for 7.1 * (24/Aug/2003, PHH) *13.05.02A: fixed unit conversion bug in dchi of dozusu_general.f * (20/Aug/2003, PHH) *13.05.01F: more NMS3 additions, added Andy's patch to phoenix.f, * added Derek's patch to cas_profiles.f * (18/Aug/2003, PHH) *13.05.01E: added Eric's pach for cas_line_select.f, added Andy's * updated Makefile.in, fixed cas_profile_read bug found * by Derek * (15/Aug/2003, PHH) *13.05.01D: added changes from Andy for new ifc, more s3main etc * changes * (14/Aug/2003, PHH) *13.05.01C: added changes from Eddie: sobid and binary_compositions * are now inputs, various s3main and NMS3 updates. * (12/Aug/2003, PHH) *13.05.01B: added alloc fixes from Eric, cas_profile bug fix from * Derek. * (10/Aug/2003, PHH) *13.05.01A: added hminus_meth to switch between H- * opacity setups, added dchi_factor to turn ionization * potential reduction on/off, added no_dchi_for_neg_ions * to turn on/off ionization potential reduction on/off * for negative ions (incl. H-), all in casscom.f, added * Eric's alloc patches, moved newcia to input and into * casscom.f to allow testing * (07/Aug/2003, PHH) *13.05.00C: added nms3_I_solver to NMS3, added some untility codes * (04/Aug/2003, PHH) *13.05.00B: added didl_discr_meth to NMS3 * (01/Aug/2003, PHH+EAB) *13.05.00A: added didl_discr_meth variable to smoothly switch * S3R2T dI/dlambda discretization methods. * (31/Jul/2003, PHH+EAB) *13.04.02J: fixed and reenabled AIX cprusage output * (29/Jul/2003, PHH) *13.04.02I: fixed SuperLU AIX-64bit compiler options * (28/Jul/2003, PHH) *13.04.02H: fixed internal to SuperLU blas call errors (learn Fortran, foools!) * (26/Jul/2003, PHH) *13.04.02G: fixed ifc compilation (24/Jul/2003, PHH) *13.04.02F: deactivate Derek's timing on AIX (they crash the code), * integrated SuperLU into the package (untested), * added '-debug' to build for easy access to * debug Compiler options, added more fixes for Absoft/OSX * compiler bugs (24/Jul/2003, PHH) *13.04.02E: added Derek's new timing routines for OSX, FreeBSD * and AIX (24/Jul/2003, PHH) *13.04.02D: fixed bug in select_lines calling NLTE line select * even if inlte=0 (23/Jul/2003, PHH) *13.04.02C: changed openMP setup in Cassandra (old setup too * complex for many systems!). dchi-correction reactivated * in cas_solve and retro-fitted into sahnlt.f so that France's * EOS uses it to (full internal consistency between ions, * ppress and Cassandra!). Added ions_4_pe call to newton * for perfect consistency with Cassandra even at very * low T's (for NLTE in chromospheres etc). zusum now uses * ionization data consistent with Cassandra. * !!!!!!! All this WILL change EOS results if PPRESS is used! * (23/Jul/2003, PHH) *13.04.02B: improved compiler options for ifc, smarter build, * AIX 5.2 xlf openMP options added (untested) * (23/Jul/2003, PHH) *13.04.02A: added Alex' changes (incl. new bremsstrahlung) * (22/Jul/2003, PHH) *13.04.01F: fixed openMP directives (21/Jul/2003, PHH) *13.04.01E: fixed various problems with cas_profile.f, added new * logical variable profile_extrap to enable/disable * profile extrapolation toward zero. * (17/Jul/2003, PHH) *13.04.01D: deactivated dchi-correction is cas_solve. to prevent * inconsistencies at high pressures. Warning printout added. * (03/Jul/2003, PHH) *13.04.01C: modified behavior of cas_greli: <=0 makes all NLTE lines * primary, 1d300 makes them all secondary. The previous setting * of doprf is retained. (03/Jul/2003, PHH) *13.04.01B: added Alex's Cassandra line select updates (01/Jul/2003, PHH) *13.04.01A: added France's EOS updates (30/Jun/2003, PHH) *13.04.00E: few updates to default.mki (30/Jun/2003, PHH) *13.04.00D: added new near-wing Lorentz profiles to cas_profiles * (11/June/2003, Happy Birthday France!!! I won't say * which, of course. I want to live. PHH) *13.04.00C: fixed bug in input.f * (07/June/2003, PHH) *13.04.00B: working on cas_solve to improve restart and error * handling. fixed case(3) bug in cas_voigt_lines.f, * added compat parameter to setup_grains to allow * comparing to old versions (set old_grain_sizes=.true. * in the namelist input to recover <130400 * results for default grain sizes!) * (05/June/2003, PHH) *13.04.00A: added Travis \& Dave's changes to the dust opacities * and treatment, including better size distributions. * this WILL change the results of dusty models! added * a few minor output fixes * (29/May/2003, Happy Birthday Mom!!!!!!!! PHH) *13.03.00H: fixed limit check bug in interp_new.f found by Travis * (28/May/2003, PHH) *13.03.00G: added new version of cas_solve with better error * handling, still experimental. Added Eddie's fixes to * He I model, NMS3 and sobidnlt.f, added fixes to cas_solve * added Eddie/Darko's more He I changes, added cas_solve * that uses Eddie's fix method * (26/May/2003, PHH) *13.03.00F: fixed layer>50 compile error in input.f, added * print_radiative_rates, added database check to bfcs * for safety (He I/II with CHIANTI!), switched cas_solve... * to global mode, added cas_solve reset in case of * failure, added Alexander's fixes to the Manual and * wrres/restart_file, added input parm check to kapcal.f * added upper T boundary check to interp_new if Mr. Freeze * is active (get him Batman!). * (21/May/2003, PHH) *13.03.00E: added Alex's fixes that I missed in 'D'! * (12/May/2003, PHH) *13.03.00D: added Andy's patch for make_chromos, added Alex's * lates patches for APED. * (09/May/2003, PHH) *13.03.00C: fixed openMP calls (not tested, no SMP box!), fixed * a few minor bugs. (05/May/2003, PHH) *13.03.00B: updated CHIANTI4 routines to ignore missing proton * rates, fixed some trivial bugs in CHIANTI4 and APED, * added Eddie's fix to init_irrad.f * (22/Apr/2003, PHH) *13.03.00A: added CHIANTI4 and APED routines of Alex. * (19/Apr/2003, PHH) *13.02.01B: reset allmwt=.false. as default for compatibility, * set nlte_line_mode=3 as default for best accuracy * (19/Apr/2003, PHH) *13.02.01A: added France's new convection treatmen for large * optical depths, introduces new 'new_mwt' input * (defaults to true!), changed default of allmwt to .true. * (18/Apr/2003, PHH) *13.02.00V: fixed 'partial NLTE' bug in cas_solve... * (19/Apr/2003, PHH) *13.02.00U: added Derek's AIX-PWR4-64bit (serial mode) compiler * options fix, added nlte_line_mode=3 [best option] * (18/Apr/2003, PHH) *13.02.00T: added nlte_line_mode input, added Eric's alloc * patches (02/Apr/2003, PHH) *13.02.00S: improved NLTE line profile code (26/Mar/2003, PHH) *13.02.00R: added special NLTE line profile code (23/Mar/2003, PHH) *13.02.00Q: changed gradad.f to return -999 in case of trouble * rather than 'stop' so that opacity tables can be * run at very low T's (21/Mar/2003, PHH) *13.02.00P: added France's patches (21/Mar/2003, PHH) *13.02.00O: more cas_profiles work (20/Mar/2003, PHH) *13.02.00N: added Travis SGI options fix, more cas_profiles * work in progress, added Eddie's OMP fixes, changed * errset call to prevent ESSL printout (16/Mar/2003, PHH) *13.02.00M: implemented prefac_cache_mode=0 (11/Mar/2003, PHH) *13.02.00L: finalized new line caching method (09/Mar/2003, PHH) *13.02.00K: updated new line caching method (07/Mar/2003, PHH) *13.02.00J: fixed rainout error handling, added global precision setting * as compile-time parameter to cas_solve (07/Mar/2003, PHH) *13.02.00I: fixed ID mode allocation bug (06/Mar/2003, PHH) *13.02.00H: adapted cas_solve for mixed NLTE/LTE cases (05/Mar/2003, PHH) *13.02.00G: added France's rainout update (04/Mar/2003, PHH) *13.02.00E: added Andy's patch fixing skipselect output * (01/Mar/2003, PHH) *13.02.00E: added Eddie's patch (23/Feb/2003, PHH) *13.02.00D: added Eric's IBM and alloc patches (19/Feb/2003, PHH) *13.02.00C: added France's patches for settling. * (14/Feb/2003: It's MY St. Valentine's Day Massacre!, PHH) *13.02.00B: added Eric's alloc-mode patches for irradiation. * (14/Feb/2003: It's the St. Valentine's Day Massacre!, PHH) *13.02.00A: changed NLTE mode to allow running LTE models (so that * the better NLTE profiles can be used!) * THIS CHANGES RESULTS IF NOT ALL AVAILABLE NLTE LEVELS * ARE USED. The old behavior can be restored by setting * compat=.true. in Cassandra/cas_setup*.f * fixed bug in resphot (ndim=0!) * added Alex' changes to build * (11/Feb/2003, PHH) *13.01.01K: fixed FreeBSD options, added Eddie's s3main.f (01/Feb/2003, PHH) *13.01.01J: fixed Linux SMP options (25/Jan/2003, PHH) *13.01.01I: added (untested) OpenMP support for ifc/Linux, made * ifc standard compiler for Linux, improved mk_* scrits * (23/Jan/2003, PHH) *13.01.01H: fixed Andy's Makefile.in improvements (22/Jan/2003, PHH) *13.01.01G: added Eric's alloc-mode patches and Andy's Makefile.in * improvements (22/Jan/2003, PHH) *13.01.01F: added patch from Derek to mollinb_voigt/gauss (stupid yeti!) * (21/Jan/2003, PHH) *13.01.01E: added patch from Derek to mollinb_voigt (stupid yeti!) * (21/Jan/2003, PHH) *13.01.01D: added prefac caching for LTE molecular Voigt lines * (19/Jan/2003, PHH) *13.01.01C: added prefac caching for LTE atomic Voigt lines * (19/Jan/2003, PHH) *13.01.01B: more detailed cache statistics output added * (18/Jan/2003, PHH) *13.01.01A: changed atomic/molecular LTE lines to cleanly * separate Voigt and Gauss line procedures, added * Voigt statistics output, added 2ifc.pl and updated * Compiler.options for INTEL ifc. * (18/Jan/2003, PHH) *13.01.00J: added Andy's mktar patch (17/Jan/2003, PHH) *13.01.00I: made NAG compiler default for OSX, removed etime.o linking * for OSX and FreeBSD (15/Jan/2003, PHH) *13.01.00H: added mk_linux_ifc, modified pptab2ascii.f to be more * general, use "export F_UFMTENDIAN=big:10,11" on Linux * to process BigEndian (the only real way to store data!) * on those pesky IA-32 machines (14/Jan/2003, PHH) *13.01.00G: adapted cio.c and Compiler.options so that Linux/ifc * can process 64bit offset files (13/Jan/2003, PHH) *13.01.00F: fixed machine.ifc.f to actually work as intended, * added P4 ifc benchmark result to Manual/bench.txt * (12/Jan/2003, PHH) *13.01.00E: update to compile with Intel ifc/icc/icpc under * Linux (09/Jan/2003, PHH) *13.01.00D: added Eddie's patches for phoenix.f (08/Jan/2003, PHH) *13.01.00C: added Eddie's patch for phoenix.f to make NLTE restart * work for all modes * added latest irradiation patches from Travis * (PHH, 06/Jan/2003) *13.01.00B: added Jason's updated wind routines and additional * patches (24/Dec/2002, Xmas Edition! HO-HO-HO!!!, PHH) *13.01.00A: fix NLTE mode bug found by Eddie (20/Dec/2002, PHH) * added new EOS fix and France's kapdust.f updates * (19/Dec/2002, PHH) * added Eric's new makehydro.f (10/Dec/2002, PHH) * Made lg_cutoff an input * added tio_shomate as input (05/Dec/2002, PHH) * Added these changes from France (05/Dec/2002, PHH): * qjanaf.f renamed qirwin.f (25/Nov/2002, FA) * generalized input/output (MISC/write_sab.f) and * use of settling stratified abundances SEHEU to * serve hydro as well as cloudy_day upon restart * included input/output of TOTN matrix (MISC/write_pi.f) * for restart from .18 file * Manual updated for description of settling paramaters * cloudy_day restart rethinked for EGPs (istartsolar) * csppress.f and csppressICE02.dqs updated for sch2, * shomates, eossc, settleos to specify upon start * also when running futur EOS tables (settleos=f). * 13.0.01L: added Eddie's Makefile.in (25/Nov/2002, PHH) * 13.0.01K: added Eddie's FUZZ/Makefile.in & Compiler.options * (16/Nov/2002, PHH) * 13.0.01J: added Eddie's NMS3 tree (14/Nov/2002, PHH) * 13.0.01I: fixed errset call in machine.aix-mpi.f to prevent * segfault on startup on pSeries 690's (PHH, 04/Nov/2002) * 13.0.01H: fixed bug found by Darko that hits NLTE wind models * that extrapolate b_i's due to grid changes (this is * now allowed for the time being), fixed build to * complain less (01/Nov/2002, PHH) * 13.0.01G: corrected target for IBM pSeries 690 machines, * added build support for Pwr4 MPI, added mk_hanni * to compile on the HLRN machines [great boxes!] * (31/Oct/2002, PHH, Halloween edition with extra * scary mega-bugs etc!) * 13.0.01F: added Eddie's changes to spultc.f, my sdot() fix * was cooler ... (21/Oct/2002, PHH) * 13.0.01E: Peter's 40th birthday Special Edition! Guaranteed * to be fortified with EXTRA bugs! (17/Oct/2002, PHH) * 13.0.01D: updated Compiler.options and Makefile.in for more * aggressive Mac OS X NAG f95 options and updated * Compiler.options for new NERSC IBM (Power3 running * AIX 5.1) (16/Oct/2002, PHH) * 13.0.01C: updated Compiler.options, added Eddie's new * cmflux_j, fctjcm and input 'use_cmflux_j' to turn * that stuff on * (15/Oct/2002, PHH) * 13.0.01B: fixed Eddie's spultc.f, added db2_new, fixed * zusum call in bi_2_ni.f (14/Oct/2002, PHH) * 13.0.01A: added Eddie's spultc.f (11/Oct/2002, PHH) * 13.0.00Y: update to Makefile for line_util (03/Oct/2002, PHH) * 13.0.00X: added Andy's filename fix to algamtn, * made cas_solve... QD switches an input to allow better * adaptation to various modeling conditions, * added Andy's Makefile.in changes * (19/Sep/2002, MH's birthday! Happy 35th, Michael!, PHH) * 13.0.00W: more switch fiddling for cas_solve, removed more stupid * write(25,'s etc. (17/Sep/2002, PHH) * 13.0.00V: tightened switch to 128/256bit arithmetic in cas_solve * to improve interoperability between machines * (06/Sep/2002, PHH) * 13.0.00U: fix for 2hp.pl to work around HP(S)UX bug, * fix bi_2_ni for dummy NLTE species, * (03/Sep/2002, PHH) * 13.0.00T: added Eddie's Compiler.options & build (30/Aug/2002, PHH) * 13.0.00S: fixed cas_dissolved bug, added Andy's Compiler.options * (29/Aug/2002, PHH) * 13.0.00R: changed cas_setup to use only specified number * of levels, mostly for testing (23/Aug/2002, PHH) * 13.0.00Q: reversed Travis's S3R2T changes due to problems * (23/Aug/2002, PHH) * 13.0.00P: added Ni VII to algam data. (22/Aug/2002, PHH) * 13.0.00O: added Eric's alloc patches. (21/Aug/2002, PHH) * 13.0.00N: added Darko's patches. (15/Aug/2002, PHH) * 13.0.00M: fixed col_rate_matrix_fast bug found by Darko that * hit chianti species without all levels (15/Aug/2002, PHH) * 13.0.00L: updated Compiler.options, added abundance_scaler.f * utility of JasonF. (12/Aug/2002, PHH) * 13.0.00K: added support for NAG f95 for Mac OS X, added * dozusu_general qas fix from Darko (09/Aug/2002, PHH) * 13.0.00J: fixed discon call to zusum for NLTE found by * Jason F. (06/Aug/2002, PHH) * 13.0.00I: fixed call stop_exit error in zusum found by Ian * tested 130000H add-on, seems to work (05/Aug/2002, PHH) * 13.0.00H: added H-like asymp Q for NLTE, changed treatment of * 1-level dummy NLTE atoms Q's dozusu_general.f) (04/Aug/2002, PHH) * 13.0.00G: working on the 64bit mode ... (31/July/2002, PHH) * 13.0.00F: added 64bit MPI version for IBM SP from Travis and Eddie * note that new flags are needed for the .dqs * (29/July/2002, PHH, first version created in Hamburg!) * 13.0.00E: fixed bfcs threshold bug discovered by Jason and killed * by Eddie (17/July/2002, PHH) * 13.0.00D: added special case for He II in bfcs (16/July/2002, PHH) * 13.0.00C: updated compiler options for solaris (16/July/2002, PHH) * 13.0.00B: added Eddie's mk_sp_nersc, Compiler.options and * NTH related stuff, * fixed He II NLTE in vcapcnt found by JPA (15/July/2002, PHH) * 13.0.00A: added code to ignore molecular lines with code > nmols * to grelmol*.f to allow old code version to read newer * line list versions. Declared the beta test over and put * version 13 as new development version for now. Last v12 * is 12.3.1A. * 13.0.00BETA-31: reverted to TiO algam data (zusumol.f) but * kept updated discon.f * (04/July/2002: Happy Birthday USA!, PHH) * 13.0.00BETA-30: changed d00 in FPPRESS/sqtio.f (29/June/2002, PHH) * 13.0.00BETA-29: updated gradad, ppress, interp_new and zusum.f, * added Jason's new windmail.f (27/June/2002, PHH) * 13.0.00BETA-28: updated regression.f (16/June/2002, PHH) * 13.0.00BETA-27: new discon.f, added Eric's updates (13/June/2002, PHH) * 13.0.00BETA-26: added Andy's and Eddie's fixes (13/June/2002, PHH) * 13.0.00BETA-25: cleaned up molecular b-f opacities, * added mol_bf_opac_on input (12/June/2002, PHH) * 13.0.00BETA-24: added most of Andy's patches (12/June/2002, PHH) * 13.0.00BETA-23: fixed compiler options for FreeBSD-NAG * (11/June/2002, Woah, France's birthday! PHH) * 13.0.00BETA-22: forced cas_solve to LTE if out of range * (11/June/2002, Woah, France's birthday! PHH) * 13.0.00BETA-20: fixed input.f problem on Mac OS X related to * common blocks used in private module (10/June/2002, PHH) * 13.0.00BETA-19: adopted c_line_list_info() to allow blocksizes * larger than 16bit (10/June/2002, PHH) * 13.0.00BETA-18: -17 didn't work, fixed problem with new internal * subroutine. Moved hydstop=1 stop after settling * setup. (08/June/2002, PHH) * 13.0.00BETA-17: moved 'get_start_info' call BEFORe "eos_init" call * to allow correct settling abundance restart * (07/June/2002, PHH) * 13.0.00BETA-16: added France's updates to compos.f, discon.f, sqcrh.f, * zusumol.f (07/June/2002, PHH) * 13.0.00BETA-15: added France's changes to FPPRESS (discon.f, * sqch4.f, sqcrh.f, sqtio.f, zusumol.f), saved * old routines, this will change the results! * (06/June/2002, PHH) * 13.0.00BETA-14: added Kurucz level data for P and Ni from *.atom * files to algam.05jun2002.data and modified zusum.f * to achieve 100% compatibility with v12 zusum and * v12 results for P and Ni (06/June/2002, PHH) * 13.0.00BETA-13: MPI mode compiles on Mac OSX (03/June/2002, PHH) * 13.0.00BETA-12: added Derek's new STDS CH4 lines, they are included * in the new molec_15_May2002.bin file (15/May/2002, PHH) * 13.0.00BETA-11: fixed memory leak in cas_solve (12/May/2002, PHH) * 13.0.00BETA-10: added Eric's alloc-mode patches (09/May/2002, PHH) * 13.0.00BETA-9: added Eddie's NLTE ID updates (08/May/2002, PHH) * 13.0.00BETA-8: added Eddie's updates (06/May/2002, PHH) * 13.0.00BETA-7: added output patches from France (03/May/2002, PHH) * 13.0.00BETA-6: added Andy's etime and nofpe patches (02/May/2002, PHH) * 13.0.00BETA-5: added France's latest settling codes (01/May/2002, PHH) * 13.0.00BETA-4: added alloc-mode fixes from Eric * (01/May/2002 [Whoa, Mayday, I should be @home], PHH) * 13.0.00BETA-3: added alloc-mode patch from Eric (30/Apr/2002, PHH) * 13.0.00BETA-2: added build diff from Andy (29/Apr/2002, PHH) * 13.0.00BETA-1: added fixes from Eddie and Eric (26/Apr/2002, PHH) * 13.0.00BETA: Beta test version number 1 (26/Apr/2002, PHH) * 13.0.00ALPHA: starting to include new Cassandra (25/Apr/2002, PHH) * 12.2.00C: added France's printout fixes (15/May/2002, PHH) * 12.2.00B: added France's printout fixes (03/May/2002, PHH) * 12.2.00A: added France's new settling codes (01/May/2002, PHH) * 12.1.00D: added France's new sqh2sc.f (26/Apr/2002, PHH) * 12.1.00C: added France's new zusumol.f and sqh2sc.f (24/Apr/2002, PHH) * 12.1.00B: more bugfixes from Darko's added (23/Apr/2002, PHH) * 12.1.00A: added Darko's chromosphere mode updates (22/Apr/2002, PHH) * 12.0.02F: fixed He I line bug in atomod found by Darko (19Apr/2002, PHH) * 12.0.02E: reversed discon.f for N2 again (15/Apr/2002, PHH) * 12.0.02D: reversed discon.f for N2 (05/Apr/2002, PHH) * 12.0.02C: changed discon data for N2, added insanity check * to gradad.f (04/Apr/2002, PHH) * 12.0.02B: fixed a time waster in 2gam.f (01/Apr/2002, * "FOOLS! I'll destroy them ALL!" [Gary Larson], PHH) * 12.0.02A: added "no_qd" switch (JasonF's request) (28/Mar/2002) * 12.0.01K: added new optical.18mar2002.d, gnsize.f, kapdust.f * to include NH3/c opacities (18/Mar/2002, PHH) * 12.0.01J: added 'reset' to logic determining solver_flag in newton.f * to catch EOS reset problem found by Jason F. (18/Mar/2002, PHH) * 12.0.01I: fixed uhext allocation bug found by Darko (06/Mar/2002, PHH) * 12.0.01H: added new optical.27feb2002.d to gnsize.f (27/Feb/2002, PHH) * 12.0.01G: changes to NMS3 for OSX and bug fixes (13/Feb/2002, PHH) * 12.0.01F: added Eddie's NMS3 (12/Feb/2002, PHH) * 12.0.01E: fixed initver bug for Na I (11/Feb/2002, PHH) * 12.0.01D: added alloc-mode fixes from Eric (07/Feb/2002, PHH) * 12.0.01C: fixed dgef call bug found by Eric (06/Feb/2002, PHH) * 12.0.01B: bug fixes and clean up in caseosb_n.f (31/Jan/2002, PHH) * 12.0.01A: added France's new AMES-TiO setup as 1200 (25/Jan/2002, PHH) * 12.0.00G: pptabl_new.f generalization fix (17/Jan/2002, PHH) * 12.0.00F: added Jason F's options for O2000, HPUX fixes * (that compiler is BAD!), minor updates in new * Cassandra (not active). (16/Jan/2002, PHH) * 12.0.00E: more fixes for new Cassandra (14/Jan/2002, PHH) * 12.0.00D: fixes for new Cassandra (10/Jan/2002, PHH) * 12.0.00C: EOS improvements, NMS3 fixes * (07/Jan/2002, PHH, long live the Euro!) * 12.0.00B: starting to implement DKL line treatment, * fixed csppress restart problem (24/Dec/2001, PHH, HO-HO-HO!) * 12.0.00ALPHA-12: added lg(T) csppress table mode (18/Dec/2001, PHH) * 12.0.00ALPHA-10: synced FPPRESS with 11.10 (sqh2sc) (17/Dec/2001, PHH) * 12.0.00ALPHA-7: synced FPPRESS with 11.10 (qjanaf) (30/Nov/2001, PHH) * 12.0.00ALPHA-3: added France's new EOS version (14/Nov/2001, PHH) * 12.0.00ALPHA: added QP 128bit and 256bit arithmetic support, * implemented use of QP in new Cassandra and retro-fitted * it into old Cassandra (fully automated). This is the * initial version of this, will compile only on Mac OSX * and FreeBSD. C++ compiler now required. (6/Nov/2001, PHH) * 11.11.0A: made sch2 an input parameter (19/Dec/2001, PHH) * 11.10.0D: added France's new H2 Q routine (17/Dec/2001, PHH) * 11.10.0C: added vo_version=1110 (10/Dec/2001, PHH) * 11.10.0B: fixed extrapolation problem in qjanaf.f (30/Nov/2001, PHH) * 11.10.0A: added new version of FPPRESS as default (old version * is FPPRESS.V11.9) (14/Nov/2001, PHH) * 11.9.02F: even more new Cassandra updates, added Eric's alloc * fix (31/Oct/2001, PHH: BOOOO HALLOWEEN!!!!! SCARY!) * 11.9.02E: added Travis's bug fixes, added new Cassandra version * updated (25/Oct/2001, PHH) * 11.9.02D: added Ian's update to machine.nag-mpi.f (24/Oct/2001, PHH) * 11.9.02C: added new cassandra parts to source tree for smooth * transition, no effect on default code (20/Oct/2001, PHH) * 11.9.02B: more chianti additions, manual patches (10/Oct/2001, PHH) * 11.9.02A: added Travis' stuff for irradiation. (24/Sep/2001, PHH) * 11.9.01A: changed casvgl and casgwl to use larger Doppler and * Lorentz width for a line rather that just deepest layer; * (21/Sep/2001, PHH) * 11.9.00H: reverted to old version EOS for backward compatibility * until new one is completely finished (21/Sep/2001, PHH) * 11.9.00G: bkdecod_f90 optimized away where statement * (15/Sep/2001, PHH) * 11.9.00F: fixed opacity table mode checks (30/Aug/2001, PHH) * 11.9.00E: added FA's latest EOS updates, added f-value multiplier * of 0.25 to SCAN-H2O in bkdecode_f90 to allow general use * of SCAN water (this MUST be undone if Jorgensen's QH2O is * used, see warning message! * This version changes the results of the EOS sligthly * due to improved thermodynamic data (17/Aug/2001, phh) * 11.9.00D: added Eddie's and Eric's add-ons and fixes (17/Aug/2001, phh) * 11.9.00C: added Andy's build V4, chianti progress (15/Aug/2001, phh) * 11.9.00B: added Andy's and Erics patches (13/Aug/2001, phh) * 11.9.00A: added France JANAF updates (09/Aug/2001, phh) * 11.8.01B: corrected/improved outer radius iteration for new_giant_mode * to faster reach convergence through better starting * guesses for R(out) and g(out) (01/Aug/2001, phh) * 11.8.01A: merge with EJL's latest 'allocate' changes (02/Jul/2001, phh) * 11.8.00B: secondary merge with FA's dust settling version * (28/May/2001, phh) * 11.8.00A: prelim merge with FA's dust settling version * (24/May/2001, phh) * 11.7.04E: added molecules to line database, added support * for France's SGI O3000 (24/May/2001, phh) * 11.7.04D: updated support (stage 1) for Sun's in MPI mode with * version 5 compiler (24/May/2001, phh) * 11.7.04C: added opac_out support for opacuty table mode by * moving capdust out of capcnt and added structure * output to binary file (21/May/2001, phh) * 11.7.04B: updated support to handle ALL atomic|molecular * lines as Voigt lines (opac tables) (19/May/2001, PHH) * 11.7.04A: support for SCAN-H2O added (09/May/2001, PHH) * 11.7.03D: various ALLOC-ON/OFF bug fixes (09/May/2001, PHH) * 11.7.03C: added work-around for MPI data structure bug on IBM SP3s * (02/may/2001, eab) * added fixes for Mac OS X (04/May/2001, PHH) * 11.7.03B: EJL's minor fixes, PHH's port to Mac OS X/MPI (02/May/2001, PHH) * 11.7.03A: added more alloc-version chanes from EJL (30/Apr/2001, PHH) * 11.7.02D: added doid for cmf models and updated windmain.f * (24/Apr/2001, JPA+PHH) * 11.7.02C: fixed fmt_id/stat switched bug (23/Apr/2001, PHH) * 11.7.02B: added FA's changes for SGI O3000 @CINES (23/Apr/2001, PHH) * 11.7.02A: added input fmt's to spectrum output (19/Apr/2001, PHH) * 11.7.01C: added vo_version to takeband.f (18/Apr/2001, PHH) * 11.7.01B: added Eddie's wdump.f update (17/Apr/2001, PHH) * 11.7.01A: added 10.5 _version's (15/Apr/2001, PHH) * 11.7.00A: new _version scheme to select fel version (11/Apr/2001, PHH) * 11.6.01A: added FA's and JWF's FPPRESS changes (09/Apr/2001, PHH) * 11.6.00A: added EJL's allocatable versions of LTELINES et al (08/Apr/2001, PHH) * 11.5.01A: fixes to Mac OS X mode, cio.c cleanup (06/Apr/2001, PHH) * 11.5.00B: fixed minor S3R2T_alloc problems (18/Mar/2001, PHH) * 11.5.00A: added Eric's S3-alloc speed improvements (14/Mar/2001, PHH) * 11.4.01A: added Eddie's fixed [rw]dump.f to MISC (02/Mar/2001,PHH) * 11.4.00A: added France's new takeband.f and bkdecode_f90.f to * code (replacing the old versions), this version also * uses a new molec6a.vocah.11.4.data file (01/Mar/2001,PHH) * 11.3.00A: added JPA's changes (windmain et al) (28/Feb/2001,PHH) * 11.2.01A: updated cshep.f anbd resphot.f routines (28/Feb/2001,PHH) * 11.2.00C: further additions to CHIANTI access module (23/Feb/2001,PHH) * 11.2.00B: further additions to CHIANTI access module (10/Feb/2001,PHH) * 11.2.00A: Mac OS X port COMPLETE, added EJL's changes: * removal of MEV in default setup, allocate mode * improvements (07/Feb/2001, PHH) * 11.1.01C: added initial Mac OS X port (05/Feb/2001, PHH) * 11.1.01B: added TSB's fixes for irradiation modes (02/Feb/2001, PHH) * 11.1.01A: added EAB's new sobidnlt.f, fixed wavcon call in * grelcmbp_local.f (17/jan/2001, PHH) * 11.1.00D: added EJL's fix to sobidnlt.f (12/jan/2001, PHH) * 11.1.00C: fixed NMS3 OP iterations (05/jan/2001, PHH, first edition of * the NEW millenium!) * 11.1.00B: made 'wlapp' input parameter (14/dec/2000, PHH) * 11.1.00A: added Jason's depth dependen vturb describtion (13/Dec/2000, PHH) * 11.0.00G: fixed bugs in cio.c and line_data_bcast (21/nov/2000 AS+PHH) * 11.0.00F: added Jason's PHXQ Perl batch queueing system (06/Nov/2000, phh) * 11.0.00E: minor fixes to NMS3 (Nov/2000, phh) * 11.0.00D: added Jason A's port for DEC AXP EV6 boxes (24/Oct/2000, phh) * 11.0.00C: fixed out-of-bounds bug in collwldata1 (19/Oct/2000, phh) * 11.0.00B: numerous bug fixes for merge, non-alloc mode is * now standard! (11/Oct/2000, phh) * 11.0.00A: initial merge with OU's add-on's (27/Sep/2000, phh) * 10.9.10B: split MPI and diffusion routines out of phoenix.f * into MISC subdirectory, added performance * enhancements for diffusion routines (21/Sep/2000, phh) * 10.9.10A: merged with diffusion add-on by Alain (14/Sep/2000, phh) * 10.9.9A: added Ian's Mn I+II NLTE atoms (12/Sep/2000, phh) * 10.9.8B: fixed table-direct mode of the EOS! see * ppress.f (23/Aug/2000, phh) * 10.9.8A: addition of Multi-Group mode (05/Jul/2000, phh) * 10.9.7K: merged with Jason's changed to wind modes, adapted * ratup*.f to better account for line profiles * (04/July/2000: Happy Birthday USA!, phh) * 10.9.7J: changed format to write in fort.7 file to allow for * more digits and easy reading * (04/July/2000: Happy Birthday USA!, phh) * 10.9.7I: fixed mode-bug in kapcal.f (27/June/2000, phh) * 10.9.7H: fixed csgrain0 bugs, removed no longer required * SGI fix. (25/Jun/2000, phh) * 10.9.7G: started to add support for NEC SX-5 (24/June/2000, phh) * 10.9.7F: boosted spectrum (unit 7) output to 15 digits for * high-res spectroscopy; fixed makefiles for * the @$%@! Sun (23/June/2000, phh) * 10.9.7E: added convection common to csgrain.f to "fix" a * SGI O2000 f90 problem (ntc set to zero). Funky ... * (23/June/2000, phh) * 10.9.7D: added support for Sun f90 version >=5.0 to build and * Compiler.options, untested. (21/June/2000, phh) * 10.9.7C: added ezl,ezlmold to line_data_bast to fix stupid * parallel machines (Hi SGI!) (20/June/2000, phh) * 10.9.7B: added Andy's freeBSD mode (with C I/O) (16/June/2000, phh) * 10.9.7A: merged with FA's dust updates. no comment (16/June/2000, phh) * 10.9.6D: added fixes to SGI Origin2000 MPI mode (15/June/2000, phh) * 10.9.6C: added fixes from CIS and AS (14/June/2000, phh) * 10.9.6B: added Ian's latest changes for chromosphere models, * passed initial verification tests (difference in H b_i * appear due to new H collision rates!). * (16/May/2000, phh) * 10.9.6A: merged with Darko's and Ian's chromosphere * additions (10/May/2000, phh) * 10.9.5B: added Eddie's new Makefile's! (10/May/2000, phh) * 10.9.5A: automatic distribution of core chars (see S3R2T/s3init.f), * automatic and adaptive setting of nwlnodes (if <= 0) * for the MPI modes (01/May/2000, phh) * 10.9.4B: added support for additional chars (19/Apr/2000, phh) * 10.9.4A: added latest patches from Andreas (19/Apr/2000, phh) * 10.9.3B: added Solaris C I/O mode for PGI (07/Apr/2000, phh) * 10.9.3A: added Linux C I/O mode for PGI (05/Apr/2000, phh) * 10.9.2B: added Linux/PGI/MPICH mode (25/Mar/2000, phh) * 10.9.2A: added France's NaCl/c opacity & setups (22/Mar/2000, phh) * 10.9.1D: updated NMS3 version (22/Mar/2000, phh) * 10.9.1C: updated lstar for additional characteristics * (16/Mar/2000, phh) * 10.9.1B: fixed default sapphire setup according to France * (07/Mar/2000, phh) * 10.9.1A: merged with CIS's chromosphere add-ons (04/Mar/2000, phh) * 10.9.0A: changed cloud treatment to scale height above * convection zone setup. (PHH, 03/Mar/2000) * 10.8.8E: fixed convection bugs for new_giant_mode=t * (01/Mar/2000, PHH) * 10.8.8D: increased default set of species for adgrad etc, * according to France, see mkinput.f (29/Feb/2000, PHH) * 10.8.8C: latest version of gradad.f and sqh2o.f from France * ((phh, 28/Feb/2000) * 10.8.8B: latest version of gradad.f, enlarged maxcor to 50 * for improved RF outputs (phh, 14/Feb/2000) * 10.8.8A: made ncore, rcout, rcin input parameters for better * control of s3r2t (phh, 11/Feb/2000) * 10.8.7D: added gradad changes from FA (PHH, 11/Feb/2000) * 10.8.7C: added vconv changes from FA (PHH, 04/Feb/2000) * 10.8.7B: andy found and fixed an evil bug in linecom.inc * nmols needed to be updated (PHH, 03/Feb/2000) * 10.8.7A: changed meaning of gnexp (FA), added updated fel * scaling factors for bkdecode_f90 (PHH, 31/Jan/2000) * 10.8.6D: opacity table bug fixes, implemented skip_select for * grelcmbp (phh, 29/jan/2000) * 10.8.6C: added H2-H2 and H2-He new CIA tables, fa: unit change in * entropy returned by gradad (phh, 25/jan/2000) * 10.8.6B: grain size distribution parameter input, extended * opacity table mode to include therm. data (phh, 19/Jan/2000) * 10.8.5C: initial support for MPI+openMP on AIX43 (phh, 11/Jan/2000) * 10.8.5B: various bug fixes and adjustments (phh, 10/Jan/2000) * 10.8.5A: added skip_select(default .false.) to allow overriding * multiple molecular line selections to save time during * iterations (phh, 03/Jan/2000) * 10.8.4B: updated bkdecode_f90 to adjust fel's (phh, 29/Dec/99) * 10.8.4A: added France's new CIA routines (KAPCAL/cia.f) to * replace the old ones (PHH, 20/Dec/99) * 10.8.3B: forced pscal=1.d0 in ppress.f to fix various probles, * set in_sn.1008 as standard input file location and * version (phh,26/Nov/99) * 10.8.3A: latest patches from Andy, added CH4, VO and CrH * lines (phh, 24/Nov/99) * 10.8.2E: added FA's fixed KAPCAL routines (phh, 15/Nov/99) * 10.8.2D: cleaned up opacity table mode, needs more testing before * going into production; added support for new Sun 64bit * compiler (phh,13/Nov/99) * 10.8.2C: added simple pressure ionization support to dozusu * casvoi*, and radratbl added auto-switch to simple NLTE * opacity formula for H I in opacity table mode (PHH, 10/Nov/99) * 10.8.2A: made dvmax, voigt_cutoff and dgmax (NLTE line windows) * inputs (phh, 07/Nov/99) * 10.8.1I: added latest patches from Andy (phh, 03/Nov/99) * 10.8.1H: added alpha-version SMP support to ltelin*.f, mollin*.f, * this might slow down the 'gauss' routines in SMP mode. * (phh, 30/oct/99) * 10.8.1G: added alpha-version SMP support to casvoi*.f, casabs*.f * and ratup*.f using extended openMP directives. May not * work correctly, experimental use only (PHH, 25/Oct/99) * 10.8.1F: numerous fixes, merges with Andy, cloud stuff polishing * (phh, 17/Oct/99, MY birthday!) * 10.8.1A: added clouds. (phh, 01/Oct/99) * 10.8.0C: merged with Andy, added Manual to source tree for * easier documentation updates. (phh, 28/sep/99) * 10.8.0B: fixed numerous dusty bugs brought to us by people who * should know better. You think I am annoyed? Heck, yes * indeed I am annoyed. That was the PG-13 version of what * I really thought when my NAG compiler found all this * stuff. (PHH, 26/Sep/99) * 10.8.0A: added 'usedust' parameter to switch/modify dust opacities * for individual species * 10.8.0A: added new dust opacities/treatment from WSU (phh, 13/Sep/99) * 10.7.4B: new options etc for AIX 4.3 on PWR3 boxes (phh, 25/Aug/99) * 10.7.4A: CIO for atomic and fuzz lines (phh, 19/Aug/99) * 10.7.3A: support for very large files (phh, 16/Aug/99) * 10.7.2C: additional support programs (cio.c, c_merge_list.f) * for use with line lists. (phh, 02/Aug/99) * 10.7.2B: support for FeH and H16OD (AMES) lines (phh, 23/Jul/99) * removed default for CH4 adgrad (doesn't work for low * Teff). * 10.7.2: added new CH4 Q's and support (zusumol,sqch4) and * updated default list of convection species to * include CH4, NOTE: this might crash some older * setups, either include CH4 in their EOS or overwrite * the set of species to remove CH4! (phh, 22/Jul/1999) * 10.7.1E: yet another fix for irrad models, changed output of * ingoing spectrum (phh, 02/jun/1999) * 10.7.1D: added ingenious fix to allow parallel line selection on * Sun SMP's, don't try this at home (phh, 31/may/1999) * 10.7.1C: added "mr_freeze" to FPPRESS/interp_new.f to freeze EOS * when table runs out, enables by default (phh, 14/may/1999) * 10.7.1B: increases max for NLTE level etc, Sun-MPI version does * NOT support parallel line selection (bug in library?) * it might work with serial line selec (phh, 12/may/99) * 10.7.1A: better support for Sun-MPI boxes (not fully debugged), * added 'rdng' switch to read NG models (phh, 10/may/99) * 10.7.0d: switched MEV to library (phh, 26/apr/99) * 10.7.0C: added fix for IO errors on IBM SP, blocked off unit 9 * () to save IO (phh, 21/apr/99) * 10.7.0B: blocks turned off in cass, mpi_reduce adapted for * small memory (phh, 15/apr/99) * 10.7.0A: added Cs I, Sr I-II, Rb I and Ba I-II NLTE (cis 09/Apr/99) * 10.6.0B: initial support for the PGI f90 compiler (Linux/Solaris), * use "build PGI-F90" to use it. Untested (phh, 04/Apr/99) * 10.6.0A: switched to new S3R2T version with improved dI/dlambda * discretization. (phh, 02/Apr/1999) * 10.5.4B: new local file version for molec. scratch, untested * on anything but T3E (phh, 22/Mar/99) * 10.5.4A: added support for prescribed irradiating spectrum (phh, 23/Feb/99) * 10.5.3A: updated with latest EOS from France, add the new * treatment for dust settling from France (phh, 15/feb/1999) * 10.5.2B: updated with latest EOS from France (phh, 11/feb/1999) * 10.5.2A: new H2O/c (ice) EOS from France added (phh, 04/feb/1999) * 10.5.1B: C IO implemented for molecular line lists on Sun * and T3E (phh, 29/Jan/1999) * 10.5.1A: new flag chosing to insert tstd=1 into the grid, * insert_tau_eq_one (default: 1) (phh, 25/jan/1999) * 10.5.0F: Cray T3E C-conversion routine (beta), opac_out=3 * mode (optical depth output) added (phh, 19/jan/99) * 10.5.0E: added cio.c routines for Solaris 64bit and faster * line list IO handling (phh, 15/Jan/99) * 10.5.0D: made level arrays allocatable, Sun port improved * split fppress.f (now FPPRESS/) (phh, 08/Jan/99) * 10.5.0C: added support for PGCC/Linux compiler (experimental!), * switched to EAB's new grel*.f routines that support * syncing to prevent T3E crashes (phh, 30/dec/98) * 10.5.0B: splitted nlte.a into 2 libs (Linux trouble) and * splitted misc.f into a library MISC/misc.a (phh, 29/dec/98) * 10.5.0A: a heck of a lot new NLTE species added by cis, * added the latest dust opacity updates and J-TiO fel * provided by FA, change of AMES-TiO fel (untested!), * merged version (phh, 10/dec/98) * 10.4.1D: turned dust opacity off during line selection to * increase performce (the dust opacity routines are * extremely primitive, phh, 01/dec/1998) * 10.4.1C: added collwldata1 to save memory.(phh, 25/nov/98) * 10.4.1B: reduction of memory usage (phh, 20/nov/98) * 10.4.1A: lots of bug fixes in caseos (phh, 13/nov/98) * 10.4.0B: added latest EOS changes by jwf (phh, 30/oct/98) * 10.4.0A: merged with cis's version, adds Fe IV-VI, Ni I-VI to NLTE. * note: Fe II sampling turned off, turn on sweeper to * reduce numer of lambda points! (phh, 28/oct/98) * 10.3.1A: fixed mk_* scripts (phh, 23/oct/98) * 10.3.0E: TiO fel updated to Langhoff 1997, sqh2o replaced by JANAF * for test in fppress/zusumol (phh, 15/oct/1998) * 10.3.0D: added support for Cray/Sun with IEEE line list format. * This might be VERY slow. I curse thee, stupid compiler * designer! May you end up working on Windoze crap! * (phh, 11/oct/1998) * 10.3.0C: molx now uses linecom, the Sun's seem to require this * No idea why (phh,10/oct/98) * 10.3.0B: updated convection detection criterion to be in sync * with options in fcthyd (phh,06/oct/98) * 10.3.0A: added David Schwenke's TiO lines (phh,21/sep/98) * 10.2.0B: added Ian's latest changes to caseosb_n (phh,16/sep/98) * 10.2.0: merged latest version of Jason's wind code (jpa,09/sep/98) * 10.1.0A: modified convection to allow for g_eff < 0 (phh,31/aug/98) * 10.1.0: merged with CIS's addition of Fe IV and improved * handling of underflows in caseosb_n, added JA's fix * of tau-grid logic (phh, 27/aug/98) * 10.0.2A: added interfaces for Cassandra, added support for use of * K-CN B-X lines with J-CN A-X lines (phh, 19/aug/98) * 10.0.1D: added interfaces for s3r2t (phh, 17/aug/98) * 10.0.1C: eliminated /algam/ block (phh, 21/jul/98) * 10.0.1B: added ramping of dust opacities (phh, 15/jul/98) * 10.0.1A: added output of R,Teff for each wavelength (phh, 14/jul/98) * 10.0.1 : added France's latest dust opacs! (phh, 13/jul/98) * 10.0.0C: added Andreas' latest changes! (phh, 18/jun/98) * 10.0.0B: improved flexible EOS table format! Set pptab_version=1 * to use it, 0 for the old format. (phh, 18/jun/98) * 10.0.0A: initial merge with Andy's Sybil version. Molecular * NLTE is not fully tested, usage is NOT recommended. * This merged version was created by PHH * (phh, 04/jun/98, woah, my sister's birthday!) * ran the whole crap through the debugger a few times, * it even seems to work correctly. Added an experimental * Solaris version, as of now not fully tested. * (phh, 11/june/98 woah, France's birthday!) * 9.1.5B: changes by Jason F. added (phh, 08/may/1998) * 9.1.5A: changes by Eddie added (phh, 07/may/1998) * 9.1.5: merged France's last updates to fppress.f and KAPCAL/, * fixed a number of bugs (cis+phh, 06/may/1998) * 9.1.4: added Jason's new He I 19-level model atom (phh, 16/mar/98) * added France's changes to TiO fel's (phh, 17/mar/98) * A: added Ian's OP cross-section update (phh, 01/apr/98, no joke!) * 9.1.3: modified gauss line search to minimize IO (phh, 14/feb/98) * added variable record size support for LTE line selection * routines, added new J-TiO fel's to bkdecode_f90 (phh, 17/feb/98) * made gausswin an input parameter (phh,19/feb/98) * A: splitted s3r2t.f into subdir S3R2T, new build (AS) (PHH, 13/mar/1998) * 9.1.2: updated to include C_p and (drho/dT)_P interpolation from * SC-EOS table; added latest fppress fixed from FA (phh,16/jan/98) * C: updated wlyman, new atomod from Eddie. (phh, 27/jan/98) * D: fixed bug inb phx_mpi_init (selec_fuzz) (phh, 30/jan/98) * E: changed capcen calc. in grel*.f to f90 array syntax to * circumvent SGI f90 bug (phh, 04/feb/98) * 9.1.1: added Isabelle Baraffe's routine to interpolate adgrad in * SC-EOS table as an option (use_eossc=.true.) (phh,30/dec/97) * 9.1.0: added K-shell photo-ionization for NLTE species * and 2-photon decay for H and He I (eab, 11/dec/97) * 9.0.1: capmff changed such that asi and amg set to zero if NLTE of Si * or Mg is used to prevent catastrophies, fixed bug in fuzzlinb * related to reading the first block in parallel mode (phh, 19/nov/97) * 9.0.0: new scratch file format and parallel line selection added (EAB+PHH) * 8.1.5: splitted fkapcal.f in KAPCAL/ for the T3E, added new * statistics output (phh, 05/sep/97) * 8.1.4: added Cs (phh, 29/jul/97) * 8.1.3B: set eosrho accuracy (1e-5), added Pg check to * grelcmb, grelmolb and fuzzb, added MPI support for * opacity table mode (phh,25/jul/97) * 8.1.3A: fixed MPI deadlock bug (blocking BCASTs) (phh,21/jul/97) * 8.1.3: added France's latest additions (phh,09/Jul/97) * 8.1.2: added France's latest additions (phh,04/Jun/97) * 8.1.1E: added level dissolution effects for hydrogen (phh,16/may/97) * 8.1.1D: changed units of flux printed to unit 7 to erg/cm^3/s (phh,13/may/97) * 8.1.1C: added input for tfac and taufac (phh,09/may/97) * 8.1.1B: added j_pipeline switch to enable j_pipeline (phh,07/may/97) * 8.1.1: added Eddie's latest improvements to MPI-S3 (eab+phh,06/may/97) * 8.1.0B: basic tests, new stuff from Jason incorporated (phh,30/apr/97) * 8.1.0A: better WL cluster parallelization, bug fixes (eab+phh,29/apr/97) * 8.0.1: new EOS, merged with Eddie's version for WL parallel (phh,02/apr/97) * 7.5.0A: added nsynch parameter to set wl-loop synch interval (phh,05/feb/97) * 7.5.0: beginning wavelength parallelization, cleanup (eab+phh,31/jan/97) * 7.4.2: line-array modularized, cleanup (phh, 30/jan/97) * 7.4.1: added nlte terms option, f90 cleanup (phh, 05/jan/97) * 7.4.0: merged with kaptab (phh, 11/dec/96) * 7.3.3: added ion group splitting to Cassandra (phh, 03/dec/96) * 7.3.2: merged with Jason's latest versions (phh+ja, 17/nov/96) * 7.3.1: MPI improvements, cassadmin changes (phh, 13/nov/96) * 7.3.0: merged with Eddie's version, NLTE for Fe+Co I-II (phh, 24/oct/96) * 7.2.2: merged with Jason's wind version (phh, 24/oct/96) * 7.2.1: allow convection only if T rises to the inside (phh, 17/oct/96) * 7.2.0b: added check to wlpkcm to avoid very close lambda * points (phh,10/oct/96) * 7.2.0a: fixed minor bugs, improved convection criterion (phh,30/sep/96) * 7.2.0: added Jason's wind code (ja+phh,09/aug/96) * 7.1.3: merged with Eddie's changes (phh,08/aug/96) * 7.1.2: added CH4 optical opacity table (fa+phh, 26/jul/96) * updated s3r2t, changed windvel interpretation, * added nmue_pp input for PP-ALI (phh, 07/aug/96) * 7.1.1: merged with Eddie's additions for non-uniform compositions and * radioactive heating (phh, 18/jun/96) * 7.1.0a: added tioeps flag for J-TiO epsilon band lines (phh, 06/jun/96) * 7.1.0: added F90 module and dynamic allocation mode for machines * that (a) support it, (b) can handle it, and (c) do it * efficiently. Can be switch with cas-on.pl and cas-off.pl. * works best in parallel mode (phh, 21/may/96) * 7.0.1: added better NLTE H I profiles, new threshold for NLTE Voigt * lines based on gf values, increased number of H I levels * possible to 30, implemented cache for secondary lines, * implemented better searchwindow alg. for fuzz lines, * corrected numerous minor annoying problems (phh, 03/may/96) * 7.0.0a: added output of radial grid to wrres/rdres (phh, 23/apr96) * 7.0.0 : MPI parallel version, initial release (phh, 13/apr/96) * 6.2.2 : more MPI additions, some are even tested! (phh, 04/mar/96) * 6.2.1.1 : added 'restricted' spectrum synthesis mode (phh, 20/feb/95) * 6.2.1 : added Andreas's changes to damping constants (now we can * use the data in the file directly), added output of the * angular radiation field for limb darkening/brightening * stuff (AS+PHH, 15/feb/96) * 6.2.0 : alpha support for MPI, * new molecular lines 140-146 added (phh, 09/jan/96) * 6.1.3 : minor bug fixes, adapted to NAG f90 (phh, 17/dec/95) * 6.1.2 : bug fixes in cassandra (phh, 16/dec/95) * 6.1.1 : minor bug fixes (phh, 15/dec/95) * 6.1.0 : added wavelength grid sweeping, b-f cross-section cutoff * in vectorized NLTE b-f cs routines in nlte.f, * cassandra optimization, streamlining (phh,03/dec/95) * 6.0.0 : added n_i solver for cassandra (phh, 27/nov/95) * 5.13.0: added YO and ZrO lines (PHH, 27/nov/95) * 5.12.0: added HCN, C2H2, CIA RT OS tables, * added H3+ lines from Neale + Tennyson (phh, 13/nov/95) * 5.11.0: direct access IO for PP tables, opac_out flag (08/nov/95,phh) * 5.10.1: added 'tau_conv_min' flag for giant models (07/nov/95, phh) * 5.10.0: added new EOS (FA+PHH, 25/oct/95, beta) * 5.9.00: added CNOSSI species (september 95, EAB+PHH) * 5.8.05: added support for nebular boundary conditions in s3 * (eab, 15/aug/95) * 5.8.04: added fallback for EOS to cure trouble (phh, 15/aug/95) * 5.8.03: changes casad2 to insert Ca II, Fe II, Ti I+II, Co II * continua only if these species are in NLTE. Otherwise, * user Verner&Yakovlev data (Mathisen for Ca II) to make * the code run faster (phh, 04/aug/95) * 5.8.02: opt. NLTE b-f cs computation (phh, 02/aug/95) * 5.8.01: fixed table mode for new EOS (phh, 29/jul/95) * 5.8.00: added Ti II and Co II model atoms (eab+phh, 28/jul/95) * 5.7.02: updated Ca II model atom to 87 level (eab+phh, 26/jul/95) * 5.7.01: updated part. fct. for H3+ (fa, 20/jul/95) * 5.7.00: latest EOS by France, 95 new molecules (phh, 11/jul/95) * 5.6.02: added updated EOS by France (phh, 04/jul/95) * 5.6.01: added NASA (D. Carbon) SiO lines (phh, 01/jul/95) * 5.6.00: added Ti I NLTE (phh, 27/jun/95) * 5.5.17: added improved vdW broadening to cassandra (phh, 26/jun/95) * 5.5.16: added support for non-uniform compositions and * gamma deposition function (eab+pen 17/may/95, * included phh 25/jun/95) * 5.5.15: added radiation pressure IO and eff. factor (phh) * 5.5.14: added tracking molecular line routines and * species rad. acceleration output (phh, 01/jun/95) * 5.5.13: updated line ID (Voigt only!) routines (phh, 27/may/95) * 5.5.12: replace elmu0 by elmass in compos (phh, 18/may/95) * 5.5.11: added printout of normalized abundances by number and ratio * relative to solar abundances (phh, 16/may/95) * 5.5.10: signal handler to stop job for HUP, will finish iteration * automatic scaling of temperature structure for dwarfs and * giants to improve initial guess for iter > 0 * (phh, 11/may/95) * 5.5.9: changed signal handler to immediately stop job for TERM * (phh, 09/may/95) * 5.5.8: changed behavior of locate to correct ltelines (phh) * 5.5.7: fixed minor bug in ltelinb and ltelinbgauss (phh, 26/apr/95) * 5.5.6: changed input of log(g) to allow O stars with winds (phh, 21/apr/95) * 5.5.5: use LTE in discon to have correct reaction rates (phh, 13/apr/95) * 5.5.4: fixed minor problems (phh, 12/apr/95) * 5.5.3: seperate search windows for atomic (dlilam) and molecular * lines (dlimol) (phh, 11/apr/95) * 5.5.2: optimized ltelinb and mollinb's, added stop for delta T<1K * (phh, 11/apr/95) * 5.5.1: changed bandwidth to 40 for s3r2t/ppali retry (phh, 09/apr/95) * 5.5.0: added Andy's iwhprof+molprof=3 (better vdW C6) plus * the splitted Gauss+Voigt arrays for better performance * (AS+PHH, 08/apr/95) * 5.4.5: corrected and updated calls to Verner etal b-f opacity * by using a 'zetel'-like approach, cf. absver1, initver1 and * zustver in fkapcal.f (phh, 23/mar/95) * 5.4.4: added NLTE for Li I (phh, 16/mar/95) * 5.4.3: added output of Pgas and Pe to wrres (phh, 16/mar/95) * 5.4.2: further optimized Voigt profiles (phh, 15/mar/95) * 5.4.1: cassandra blocking alg. includes caching of blocks (eab+phh, 13/mar/95) * 5.4.0: cassandra blocking alg. implemented (eab+phh, 02/mar/95) * 5.3.10: adapted to compile & run somewhat smoothly on HP+DEC boxes * (phh,28//feb/95) * 5.3.9: removed 30MB array by optimizing interp/pptabl (phh 16/feb/95) * 5.3.8: included Jorgensen H2O opacity table (from DRA) (phh 14/feb/95) * 5.3.7: increased Mg II levels (eab, pen, phh 06/feb/95) * 5.3.6: added O I NLTE (eab+phh, 18/jan/95) * 5.3.5: added automatic Gauss search window to ltelines.f (phh, 16/jan/95) * 5.3.4: replaced Yakovlev etal cs by Verner + Yakovlev edition * of the enlarged dataset, not activated yet (phh, 12/jan/95) * 5.3.3: added JOLA ID feature (no avtio, avh2o or avcn) (phh, 22/dec/94) * 5.3.2: made taulin available as input taulinrd (phh, 22/dec/94) * 5.3.1: merged with Andy's version, include more molecular lines, * splitted blocking alg. to use two files for atomic and * molecular lines, LTE lines not in separate file ltelines.f * to enhance modularity (phh, 17/dec/94) * 5.3.0: include grains a' la Alexander in the opacity (fa, phh, 16/dec/94) * 5.2.0: include grains a' la Alexander in the EOS (fa, phh, 10/dec/94) * 5.1.0: include high-T Yakovlev etal b-f data (phh, 21/nov/94) * 5.0.1: splitted machine dependent routines into machine.f (phh, 09/nov/94) * 5.0.0: Fe II bulk NLTe lines, ETLA code removed (phh+eab, 07/nov/94) * 4.9.3: added doid flag and code for lineID output (phh, 22/oct/94) * WARNING: doid = .true. will slow down the code significantly! * 4.9.2: added HITRAN92 lines (as, 05/oct/94) * 4.9.1: added convec_opacity_flag for different convection (phh, 07/sep/94) * 4.9.0: Fe II (yes, that's single ionized iron) NLTE (eab+phh, 01/sep/94) * 4.8.1: read in tables of extra points for NLTE lines (phh,30/aug/94) * 4.8.0: added Cassandra-Voigt and support for wlapp as subset (phh,24/aug/94) * 4.7.9: call ratupd only if inlte is non-zero, performance (phh,27/jul/94) * 4.7.8: call s3redf only if ultcor is false, performance (phh,25/jul/94) * 4.7.7: added molprof flag for molecular profile function (phh,26/may/94) * 4.7.6: minor bug fixes and improvements (phh,20/may/94) * 4.7.5: added nltewl routine to nlte.f, updated grelcm(b) (phh,09/may/94) * 4.7.4: added recover code for singular matrices in fppress (phh,09/may/94) * 4.7.3: added Jorgensen TiO lines to misc.f (phh,04/may/94) * 4.7.2: added 'blocked line algorithm' (misc.f) (phh,02/may/94) * 4.7.1: added 'docastri' switch (phh,25/apr/94) * 4.7.0: Ne I NLTE (EAB), computed line albedo (EAB), bug fixes (PHH+EAB), * dfdw and vdfdw changes (FA), merge of OKC and PHX versions (22/apr/94) * 4.6.0: added vectorized version of kapcal (recall mode only!) and * kapmol to the opacity module for better performance. (16/apr/94, phh) * 4.5.3: added allmwt input parameter to perform all convective iterations * using mwt=.true. and mit=.false., gives better convergence * for models that are unstable using the other method (07/apr/94, phh) * 4.5.2: changed grelcm to use new (BK's) line list format (04/apr/94, phh) * 4.5.1: changed partition function of TiO, no Janaf [fppress] (27/mar/94,phh) * 4.5.0: added more molecules and improved selection scheme (16/feb/94,phh) * 4.4.11: added otab feature to kapcal (14/feb/94, phh) * 4.4.10: added skipselect feature to grelcm and grelmol (09/feb/94, phh) * 4.4.9: added printout of radiation pressure and gradient (04/feb/94, phh) * 4.4.8: corrected problems for table=true & fastmode=1 operation (01/feb/94, phh) * 4.4.7: adapted phoenix, fcthyd & ablak for convection in giants (27/jan/94, phh) * 4.4.6: added windbeta parameter for special velocity laws (31/dec/93, phh) * 4.4.5: added windvel flag for special velocity laws (30/dec/93, phh) * 4.4.4: fixed bug in kapmol-calls to avtio and avcn (09/dec/93, phh) * 4.4.3: support for 3D radiative transfer (agnst,agnend) (phh, 03/dec/93) * 4.4.2: support for AGN models (modtyp = 6) (phh, 18/nov/93) * 4.4.1: Voigt profiles incorporated for molecular lines (phh, 01/nov/93) * 4.4.0: Voigt profiles incorporated for metal lines (phh, 01/nov/93) * 4.3.2: finalized molecular line treatment, merged list structures, * added parameters for molecular lines (ezlmol, gremlin, epsmol) * (phh, 26/oct/93) * 4.3.1: molecular lines added in version 1 (too be replaced) (phh, 06/oct/93) * 4.3.0: added D. Mihalas levels to algam data, added support for * higher ionization stages, changed zusu (france) to accommodate * changes and reduce memory reqs. set newalg to true in phoenix * to use new algam data, false will use old data (phh, 24/sep/93) * 4.2.7: added new Miller H2O opacities (phh) * 4.2.6: added additional species of new LTE line list (phh) * 4.2.5: skip Lagrange->Euler RF transformation if static mode (phh) * 4.2.4: minor modification to cassandra (reverse call to kapcal) (phh) * 4.2.3: corrected non-thermal initilzation bug (eab) * 4.2.2: corrected Bug found by Ed, incorporated his changes (eab+phh) * 4.2.1: line of nlte-candidates are thrown out, if INLTE .ne. 0. * before this version, these line were thrown out only if the species * was treated in NLTE. Cassandra always includes these transitions, * so we have to throw them out of the LTE list!! (phh) * 4.2.0: incorporated eb's He II and Na I treatment, minor cosmetic * changes. (phh, 13/aug/93, friday the 13th: Are You Ready For Freddy?) * 4.1.0: convection and stellar modes (giant and dwarf) implemented * and compared to test with France's code. Added ieos=2: use France' * pp tables instead of calculating (much faster). (phh, 11/aug/93) * 4.0.1: fixed problems with the EOS's, added new s3r2t.f version, * inculded new variable for ALO bandwidth. (phh, 23/jul/93) * 4.0.0: exchanged the old EOS with the molecular EOS of * France Allard. Added Brent-style solver for pure ion * eos, added molecules consistently. (phh, 21/jul/93) * 3.3.3: added Unsoeld-Lucy type temperature correction, * used if tsw .lt. 0 in unit 5 input file (phh, 15/may/93) * 3.3.2: fixed bug in bi2ni (he 2) (eb+phh, 24/mar/93) * 3.3.1: added modtyp=3 (v0 interpreted as vmax, linear velocity law!) * (phh+eb, 23/mar/93) * 3.3.0: added true 'arbitrary' density profiles: rdump,wdump,rholaw * for better communication with other codes. level 0 dumps implemented. * cleanup code for modifications in phoenix. (11/mar.1993, phh) * 3.2.3: oldsty stuff added for input file type detection, * s3r2t.f updated for plane-parallel transfer * 3.2.0: merged with ed barons modifications (density law, He I NLTE) * 3.1.1: introduced taulin variable for s3r2t interpolation selection: * if(dt+dt1 .le. taulin) -> linear interpol, parabolic otherwise * helps to prevent crashes in s3r2t to to bad interpolations (phh) * fixed Balmer line (H) bug in grelcm (h line now removed from list) * (22/nov/92, phh) * added cassandra support for badly approximated nlte lines, * see radrat for details (method like in ETLA). (11/oct/92, phh) * variable cleanup in phoenix, cmflux and fcthcm (block /cmprim/) * (29/sep/92, phh) * added rdump subroutine to read in dump level 1 files. (09/feb/92, phh) * changed definitions in mean opacities to be consistent with * mihalas**2 (kj,kb: absorption, kh,ross: extinction) (07/feb/92, phh) * remove bug in dinscm (mev.f) regarding planck function (18/jan/92, phh) * negative vtau1 gives mass-loss in msun/yr for nova-models (22/oct/91) * read-in of irtmth and ezl, (19/oct/91 phh) * temperature correction for s3r2t implemented (13/oct/91 phh) * approximative treatment of weak nlte line (lambda .gt. 1.2d4a) * implemented (phoenix.f,nlte.f,misc.f; common /nltapp/) * (14/jun/91 phh) * s3r2t radiative transfer added (21/may/91 phh), alpha version * line selection routine generalized (apr/91 phh) * optimized version, code for holger added (apr/91 phh) * generalized nlte version (oct/90 phh) * comoving frame version (28/nov/89 phh) c********************************************************************** cHPMLIB#include "f_hpm.h" /* !HPMLIB */ c c program phoenix c --------------- implicit none ...