The coding for the GLF23 driver begins as follows:
c@testglf.f
c 18-feb-03 version 1.6
c stand-alone driver for the GLF23 model
c written by Jon Kinsey, General Atomics
c---:----1----:----2----:----3----:----4----:----5----:----6----:----7-c
c
implicit none
c
c...declare variables
c
character cdate*24, line*132
c character ctime*8 ! used if calling routine clock w/ +U77 option
integer jpd
double precision epsilon
parameter ( jpd=50 )
c
double precision te_m(0:jpd), ti_m(0:jpd)
& , ne_m(0:jpd), ni_m(0:jpd), ns_m(0:jpd)
& , zpte_m(0:jpd), zpti_m(0:jpd), zpne_m(0:jpd), zpni_m(0:jpd)
& , angrotp_exp(0:jpd), egamma_exp(0:jpd), gamma_p_exp(0:jpd)
& , vphi_m(0:jpd), vpar_m(0:jpd), vper_m(0:jpd)
& , zeff_exp(0:jpd), bt_exp, bteff_exp(0:jpd), rho(0:jpd), arho_exp
& , gradrho_exp(0:jpd), gradrhosq_exp(0:jpd)
& , rmin_exp(0:jpd), rmaj_exp(0:jpd), rmajor_exp
& , q_exp(0:jpd), shat_exp(0:jpd), alpha_exp(0:jpd)
& , elong_exp(0:jpd), zimp_exp, amassimp_exp, amassgas_exp
& , alpha_e, x_alpha
c
double precision zpte_in, zpti_in, zpne_in, zpni_in, drho
c
double precision diffnem, chietem, chiitim
& , etaphim, etaparm, etaperm, exchm
& , diff_m(0:jpd), chie_m(0:jpd), chii_m(0:jpd), etaphi_m(0:jpd)
& , etapar_m(0:jpd), etaper_m(0:jpd), exch_m(0:jpd)
& , egamma_m(0:jpd), egamma_d(0:jpd,10), gamma_p_m(0:jpd)
& , anrate_m(0:jpd), anrate2_m(0:jpd)
& , anfreq_m(0:jpd), anfreq2_m(0:jpd)
c
integer lprint, nroot, jshoot, jmm, jmaxm, itport_pt(1:5)
& , igrad, idengrad, i_delay, j, k, leigen, irotstab, bt_flag, iglf
The variables in the input file in are read through
the namelist nlglf:
c
namelist /nlglf/ leigen, lprint, nroot, iglf, jshoot, jmm, jmaxm
& , itport_pt, irotstab, te_m, ti_m, ne_m, ni_m, ns_m
& , igrad, idengrad, zpte_in, zpti_in, zpne_in, zpni_in
& , angrotp_exp, egamma_exp, gamma_p_exp, vphi_m, vpar_m, vper_m
& , zeff_exp, bt_exp, bt_flag, rho, arho_exp
& , gradrho_exp, gradrhosq_exp
& , rmin_exp, rmaj_exp, rmajor_exp, q_exp, shat_exp
& , alpha_exp, elong_exp, zimp_exp, amassimp_exp, amassgas_exp
& , alpha_e, x_alpha
c