Printout

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