Relevant printout of the input and output variables including the diffusivites and the growth rates and frequencies of the first and second roots from the dispersion relation are given below.
Inside the glf2d routine, diagnostic printout is provided if the variable lprint is greater than zero. Setting lprint to 6 prints the matrices and solution to the eigenvalue equation and setting it higher to 99 gives detailed output from the model. No diagnostic printout is given inside routine callglf2d.
c
c..printout
c
if ( lprint .gt. 0 ) then
write(6,100)
do j=0,jmaxm
write (6,110) rho(j), te_m(j), ti_m(j), ne_m(j), ni_m(j)
& , zeff_exp(j), q_exp(j), shat_exp(j)
enddo
c
write(6,120)
do j=0,jmaxm
write(6,110) rho(j), zpte_m(j), zpti_m(j)
& , zpne_m(j), zpni_m(j)
enddo
c
write(6,130)
do j=0,jmaxm
write (6,110) rho(j), diff_m(j), chie_m(j), chii_m(j)
& , etaphi_m(j), etapar_m(j), etaper_m(j), exch_m(j)
enddo
c
write(6,140)
do j=0,jmaxm
write(6,110) rho(j), egamma_m(j), gamma_p_m(j)
& , anrate_m(j),anrate2_m(j), anfreq_m(j), anfreq2_m(j)
enddo
endif
c
100 format(t5,'rho',t13,'Te',t21,'Ti',t29,'ne',t37,'ni'
& , t44,'Zeff',t53,'q',t60,'shear',t68,'#prof')
110 format (11(0pf8.4))
120 format(/,t5,'rho',t12,'zpte',t20,'zpti',t28,'zpne'
& , t36,'zpni',t68,'#log-grad')
130 format(/,t5,'rho',t12,'diff',t20,'chie',t28,'chii'
& , t35,'etaphi',t43,'etapar',t51,'etaper',t60
& , 'exch',t68,'#chi')
140 format(/,t5,'rho',t10,'egamma_m',t19, 'gamma_p',t27
& , 'anrate', t35,'anrate2',t43,'anfreq',t51,'anfreq2'
& , t68,'#gamma')
c
900 continue
c
stop
end