Merge branch 'develop' into develop

.github/ISSUE_TEMPLATE/bug_report.md

@@ -2,7 +2,7 @@
 name: Bug report
 about: Create a report to help us improve
 title: ''
-labels: ''
+labels: bug
 assignees: ''
 
 ---

.github/ISSUE_TEMPLATE/code-enhancement.md

@@ -0,0 +1,20 @@
+---
+name: Code enhancement
+about: Suggest an improvement for this project
+title: ''
+labels: enhancement
+assignees: ''
+
+---
+
+**Is your feature request related to a problem? Please describe.**
+A clear and concise description of what the problem is. Ex. I'm always frustrated when [...]
+
+**Describe the solution you'd like**
+A clear and concise description of what you want to happen.
+
+**Describe alternatives you've considered**
+A clear and concise description of any alternative solutions or features you've considered.
+
+**Additional context**
+Add any other context or screenshots about the feature request here.

.gitignore

@@ -15,5 +15,4 @@ uedge.egg-info/
 dist/uedge-8.0.0-py3.7-macosx-10.9-x86_64.egg
 uedge.egg-info
 *.c
-*.f
-
+*.f

aph/aph.v

@@ -1,7 +1,7 @@
 aph   # Atomic physics (hydrogenic)
 
 ***** Physical_constants:
-ev	real	/1.6022e-19/   # 1 electron volt in Joules
+ev_aph	real	/1.6022e-19/   # 1 electron volt in Joules
 m_prot  real    /1.67e-27/     # proton mass
 
 ***** Data_input:
@@ -15,7 +15,7 @@ aphdir	character*120 # name of directory containing data files
 data_directory	character*120 # another dirname containing data files. This is to be be passed in
 
 ***** Ionization_energy:
-erad	real [eV] /25./ # tot elec engy loss/ioniz (rad+binding) if istabon=0
+erad	real [eV] /25./ +input # tot elec engy loss/ioniz (rad+binding) if istabon=0
 
 ***** Rtdata:
 # hydrogenic rate table data from ADPAK via Braams' rate code
@@ -145,23 +145,23 @@ readehr2(fname:string)	subroutine
 
 ***** Aphwrk:
 # working arrays for 2-d spline interpolation
-nxdata  integer
-nydata  integer
-xdata(1:nxdata)             _real
-ydata(1:nydata)             _real
-fdata(1:nxdata,1:nydata)    _real
-ldf                          integer
-iflag                        integer
-kxords  integer /4/               # order of spline fit versus log(te)
-                                  # kxords=4 (default) is cubic interpolation
-kyords  integer /4/               # order of spline fit versus log10(ne)
-                                  # kyords=4 (default) is cubic interpolation
-xknots(1:nxdata+kxords)     _real
-yknots(1:nydata+kyords)     _real
-workh(1:nxdata*nydata+2*kxords*(nxdata+1)) _real # work array
-rsacoef(1:nxdata,1:nydata)  _real # spline coeff's for ionization
-rracoef(1:nxdata,1:nydata)  _real # spline coeff's for recombination
-rqacoef(1:nxdata,1:nydata)  _real # spline coeff's for line emission
+nxdata_aph  integer
+nydata_aph  integer
+xdata_aph(1:nxdata_aph)             _real
+ydata_aph(1:nydata_aph)             _real
+fdata_aph(1:nxdata_aph,1:nydata_aph)    _real
+ldf_aph                          integer
+iflag_aph                        integer
+kxords_aph  integer /4/               # order of spline fit versus log(te)
+                                  # kxords_aph=4 (default) is cubic interpolation
+kyords_aph  integer /4/               # order of spline fit versus log10(ne)
+                                  # kyords_aph=4 (default) is cubic interpolation
+xknots_aph(1:nxdata_aph+kxords_aph)     _real
+yknots_aph(1:nydata_aph+kyords_aph)     _real
+workh(1:nxdata_aph*nydata_aph+2*kxords_aph*(nxdata_aph+1)) _real # work array
+rsacoef(1:nxdata_aph,1:nydata_aph)  _real # spline coeff's for ionization
+rracoef(1:nxdata_aph,1:nydata_aph)  _real # spline coeff's for recombination
+rqacoef(1:nxdata_aph,1:nydata_aph)  _real # spline coeff's for line emission
 
 ***** Subs:
 # Subroutines that can be called from the parser