Brief GLF23 documentation as it relates to TRANSP ------------------------------------------------ K. Indireshkumar April 2, 2002 1) Invoking GLF23 in TRANSP -------------------------- Currently GLF23 model is used in Ion and Electron power balance equations. GLF23 model is invoked with the settings NKIMOD=13 and FIEFAC=1.0 (ion power balance) and NKEMOD=13 and NLXEIN=.TRUE. (electron power balance) in the namelist. GLF23 model in TRANSP currently requires that conductivity data be made available in Ufiles. The data is used where the conductivity values returned by GLF23 are zero. In order to invoke a coupled Te-Ti solver (new feature), in additon to the above switches, set NLPREDIC=.TRUE. in namelist. 2) Other optional namelist switches and variables ------------------------------------- In addition to the above switches needed to invoke GLF23, additional switches and variables may be set to control the application of the model. The real variable XRELFAC is a simple smoothing (in time) option that smoothes the rapid fluctuations in transport coefficients returned from GLF23. For example, electron conductivity chie is taken in the form, chie(t) = (1.0-XRELFAC)*chie(t-1)+XRELFAC*chie_glf23 where chie_glf23 is the conductivity returned from GLF23 and chie(t) is the value actually used in the transport equations. Hence set, XRELFAC=0.0 to 1.0 (0.1 default) The variables XIMAXGLF and XIMINGLF control the spatial extent over which the transport coefficients returned from GLF23 are used. The model is applied only in those locations where the normalized toroidal flux XI (radial coordinate) is such that XIMINGLF <= XI <= XIMAXGLF. XIMINGLF =0.0 to 1.0 (0.0 default) XIMAXGLF =0.0 to 1.0 (1.0 default) and XIMINGLF < XIMAXGLF. 3) GLF23 model control switches (variables) and their equivalents in TRANSP ------------------------------------------------------------ The following may also be set to control GLF23 TRANSP Switch Default Possible GLF23 (variable) values (switch) ------------- -------- ------- -------- LEIGEN 0 0 or 1 LEIGEN LJSHOOT 0 0 or 1 JSHOOT LIGLF 0 0 or 1 IGLF LIGRAD 1 0 or 1 IGRAD LDENGRAD 2 2 or 3 IDENGRAD ALPHA_E 0 0 or 0.6 ALPHA_E (See also Ref 1) LX_ALPHA 0 0 or 1 X_ALPHA L_DELAY 0 0 I_DELAY (See also Ref 1) LBT_FLAG 0 0 or 1 BT_FLAG LTPORT 1,1,1,1,0 0 or 1 (for all) ITPORT_PT(1:5) LNROOT 8 8 or 12 NROOT (See Ref 2) LROTSTAB 0 0 or 1 IROTSTAB (See Ref 2) 4) GLF23 related output (a) The following are rplot profile outputs XKIGLF23 -- Ion conductivity XKEGLF23 -- Electron Conductivity DIFFIGLF -- Ion plasma diffusivity ETPHIGLF -- Toroidal velocity diffusivity ETPARGLF -- Parallel velocity diffusivity ETPERGLF -- Perpendicular velocity diffusivity EXBGLF -- EXB shear rate GRAT1GLF -- Growth rate of leading mode GRAT2GLF -- Growth rate of 2nd leading mode FRAT1GLF -- Leading mode frequency FRAT2GLF -- 2nd leading mode frequency (Note: As of 12/10/02, density and momenta are not being predicted using GLF23. Hence, DIFFIGLF, ETPHIGLF, ETPARGLF, and ETPERGLF are not being used. They are for informational purposes only.) b) The multigraph CHI_GLF displays all diffusivities and conductivities. 5) References ---------- 1) Testglf.tex (Testglf.ps) documentation on GLF23 NTCC module. 2) Comments at the beginning of callglf2d.for in the GLF23 suite of programs from NTCC (in the library $CODESYSDIR/source/glf for Transp users) Notes: Some of the varialbes and switches are not in Ref. 1, but are in Ref. 2. Also, Ref. 2 contains good documentation on the switch ITPORT_PT(1:5).