TRANSP Grid Analysis 200741M35 D3D tr.tail
Reading in namelist from toray.in ...
LinLiu ECCD / <B>
Ray pattern = .gauss. Number of rays = 30
running GAFIT version 4.0 (12 July 2012)
WARNING: Ignoring unrecognized command-line argument "200741M35".
Only "-trace" and "-verbose" are valid options,
though they may be abbreviated to two characters.
reading user-provided GAFIT namelist file
TorGA_read_psiin:: nr = 129 ( 65 )
TorGA_read_psiin:: nz = 129 ( 65 )
Ray= 1 ISTOP= 7 IFLAG= 2 AT 0.163799E+03 0.747732E+01 0.276633E+02
Ray= 2 ISTOP= 7 IFLAG= 2 AT 0.163879E+03 0.857299E+01 0.276550E+02
Ray= 3 ISTOP= 7 IFLAG= 2 AT 0.164145E+03 0.780786E+01 0.266109E+02
Ray= 4 ISTOP= 7 IFLAG= 2 AT 0.163851E+03 0.659390E+01 0.269761E+02
Ray= 5 ISTOP= 7 IFLAG= 2 AT 0.163858E+03 0.656051E+01 0.284700E+02
Ray= 6 ISTOP= 7 IFLAG= 2 AT 0.163952E+03 0.777392E+01 0.289234E+02
Ray= 7 ISTOP= 7 IFLAG= 2 AT 0.163912E+03 0.882438E+01 0.296347E+02
Ray= 8 ISTOP= 7 IFLAG= 2 AT 0.164218E+03 0.948575E+01 0.286914E+02
Ray= 9 ISTOP= 7 IFLAG= 2 AT 0.163932E+03 0.967982E+01 0.272241E+02
Ray= 10 ISTOP= 7 IFLAG= 2 AT 0.163786E+03 0.926096E+01 0.259383E+02
Ray= 11 ISTOP= 7 IFLAG= 2 AT 0.163852E+03 0.833384E+01 0.252136E+02
Ray= 12 ISTOP= 7 IFLAG= 2 AT 0.163832E+03 0.718417E+01 0.251110E+02
Ray= 13 ISTOP= 7 IFLAG= 2 AT 0.163732E+03 0.611863E+01 0.256586E+02
Ray= 14 ISTOP= 7 IFLAG= 2 AT 0.163862E+03 0.540114E+01 0.268472E+02
Ray= 15 ISTOP= 7 IFLAG= 2 AT 0.164151E+03 0.523104E+01 0.283198E+02
Ray= 16 ISTOP= 7 IFLAG= 2 AT 0.164312E+03 0.566216E+01 0.295796E+02
Ray= 17 ISTOP= 7 IFLAG= 2 AT 0.164251E+03 0.658101E+01 0.302818E+02
Ray= 18 ISTOP= 7 IFLAG= 2 AT 0.164268E+03 0.771780E+01 0.303812E+02
Ray= 19 ISTOP= 7 IFLAG= 2 AT 0.164308E+03 0.872498E+01 0.319771E+02
Ray= 20 ISTOP= 7 IFLAG= 2 AT 0.163916E+03 0.104421E+02 0.303813E+02
Ray= 21 ISTOP= 7 IFLAG= 2 AT 0.163820E+03 0.113362E+02 0.281827E+02
Ray= 22 ISTOP= 7 IFLAG= 2 AT 0.163749E+03 0.111999E+02 0.258102E+02
Ray= 23 ISTOP= 7 IFLAG= 2 AT 0.163866E+03 0.100388E+02 0.239744E+02
Ray= 24 ISTOP= 7 IFLAG= 2 AT 0.163729E+03 0.821829E+01 0.229892E+02
Ray= 25 ISTOP= 7 IFLAG= 2 AT 0.163688E+03 0.617765E+01 0.233161E+02
Ray= 26 ISTOP= 7 IFLAG= 2 AT 0.164041E+03 0.446618E+01 0.250083E+02
Ray= 27 ISTOP= 7 IFLAG= 2 AT 0.164130E+03 0.359366E+01 0.272521E+02
Ray= 28 ISTOP= 7 IFLAG= 2 AT 0.164220E+03 0.375944E+01 0.295875E+02
Ray= 29 ISTOP= 7 IFLAG= 2 AT 0.164362E+03 0.490043E+01 0.314617E+02
Ray= 30 ISTOP= 7 IFLAG= 2 AT 0.164286E+03 0.672465E+01 0.322828E+02
TorGa_wrtnetcdf, kopt= 1
TORAY-GA total CPU time 0.50 sec
Generic ECH post-processing: ps2trcom_ech...
antenna # iant = 1
powecha_use(iant) = 0.0000000000000000
ps_ts2 prof sum = 0.0000000000000000
trcom prof sum = 0.0000000000000000
ps_ts2 cur sum = 0.0000000000000000
trcom cur sum = 0.0000000000000000
antenna # iant = 2
powecha_use(iant) = 426614.14973532391
ps_ts2 prof sum = 426536.18117919972
trcom prof sum = 426536.18117919989
ps_ts2 cur sum = -2672.4171946986044
trcom cur sum = -2672.4171946986048
antenna # iant = 3
powecha_use(iant) = 590094.56821166223
ps_ts2 prof sum = 590013.86500476184
trcom prof sum = 590013.86500476161
ps_ts2 cur sum = 4193.8482720257407
trcom cur sum = 4193.8482720257380
antenna # iant = 4
powecha_use(iant) = 534395.75228690112
ps_ts2 prof sum = 534318.45006627066
trcom prof sum = 534318.45006627066
ps_ts2 cur sum = -4540.7259463062510
trcom cur sum = -4540.7259463062492
Total ECH heating: 1550868.4962502320
Total Current Drive: -3019.2948689791156
TORAYMAIN Bottom: Current directory= /scratch/shared/tr_jlestz/transp_compute/D3D/200741M35
*** FULL NEUTRAL SOURCE CALCULATION PERFORMED *** NSTEP= 1737
TA= 7.63000E-01 CPU TIME= 6.90180E-02 SECONDS. DT= 1.04410E-04
TEC1= 0.76300000000000034
TEC2= 0.76400000000000035
TECON= 0.50069999694824219
TECOFF= 2.4198999404907227
fe_outtim(1)= 9.9999999999999995E+033
fe_outtim(2)= 9.9999999999999995E+033
%check_save_state: nbflag = T
%check_save_state: iwrite_now = F
%check_save_state: check at wall_hours = 21.650403055555444
--> plasma_hash("gframe"): TA= 7.630000E-01 NSTEP= 1737 Hash code: 58124874
->PRGCHK: bdy curvature ratio at t= 7.6350E-01 seconds is: 7.7002E-02
% MHDEQ: TG1= 0.763000 ; TG2= 0.763500 ; DTG= 5.000E-04
* MHD EQUILIBRIUM CALCULATED, CPU TIME = 3.3440E-03 SECONDS
DATA R*BT AT EDGE: 3.2388E+02
->EQBDY_CHECK: last flux surface curvature ratio is: 7.7002E-02
%get_rygrid: no grids for Psi(R,Y)
% GFRAME - GEOMETRY TIMESTEP TG1= 0.763000 TO TG2= 0.763500 @ NSTEP 1737
GFRAME TG2 MOMENTS CHECKSUM: 1.2517604750806D+04
%POHCAL warning: q input AND Te prediction, POH < 0 suppressed!
?pt_dvs_solver: itmax= 400
?pt_dvs_solver: too many iterations.
?pt_dvs_solver: debug files written:
state files: pt_ta.debug.cdf pt_tb.debug.cdf
control obj: ctrl.debug.cdf
%POHCAL warning: q input AND Te prediction, POH < 0 suppressed!
%check_save_state: nbflag = F
%check_save_state: iwrite_now = T
%check_save_state: check at wall_hours = 21.879728888888849
%wrstf: start call wrstf.
%wrstf: open new restart file:200741M35RS.DAT
%wrstf: open200741M35RS.DAT_NEW, status= 0
%wrstf: record, ier= 1 (ier=1, OK)
%wrstf: psload_save
%ps_rezone_chk: no rezone operator available (no MHD equilibrium).
%wrstf: write marker file.
%wrstf: transp_rplot_state_backup
% RESTART RECORD WRITTEN AT TA= 7.6316967E-01
arrays allocated in pt_diff_data_init
arrays allocated in ptsolver_data_init
%WRSTF: virtual memory size = 1.407E+03 MB.
?pt_dvs_solver: itmax= 400
?pt_dvs_solver: too many iterations.
?pt_dvs_solver: debug files written:
state files: pt_ta.debug.cdf pt_tb.debug.cdf
control obj: ctrl.debug.cdf
?pt_dvs_solver: itmax= 400
?pt_dvs_solver: too many iterations.
?pt_dvs_solver: debug files written:
state files: pt_ta.debug.cdf pt_tb.debug.cdf
control obj: ctrl.debug.cdf
%POHCAL warning: q input AND Te prediction, POH < 0 suppressed!
%POHCAL warning: q input AND Te prediction, POH < 0 suppressed!
?pt_dvs_solver: itmax= 400
?pt_dvs_solver: too many iterations.
?pt_dvs_solver: debug files written:
state files: pt_ta.debug.cdf pt_tb.debug.cdf
control obj: ctrl.debug.cdf
%POHCAL warning: q input AND Te prediction, POH < 0 suppressed!
%check_save_state: nbflag = F
%check_save_state: iwrite_now = T
%check_save_state: check at wall_hours = 22.560566388888677
%wrstf: start call wrstf.
%wrstf: open new restart file:200741M35RS.DAT
%wrstf: open200741M35RS.DAT_NEW, status= 0
%wrstf: record, ier= 1 (ier=1, OK)
%wrstf: psload_save
%ps_rezone_chk: no rezone operator available (no MHD equilibrium).
%wrstf: write marker file.
%wrstf: transp_rplot_state_backup
% RESTART RECORD WRITTEN AT TA= 7.6323562E-01
arrays allocated in pt_diff_data_init
arrays allocated in ptsolver_data_init
%WRSTF: virtual memory size = 1.407E+03 MB.
%POHCAL warning: q input AND Te prediction, POH < 0 suppressed!
?pt_dvs_solver: itmax= 400
?pt_dvs_solver: too many iterations.
?pt_dvs_solver: debug files written:
state files: pt_ta.debug.cdf pt_tb.debug.cdf
control obj: ctrl.debug.cdf
%POHCAL warning: q input AND Te prediction, POH < 0 suppressed!
?pt_dvs_solver: itmax= 400
?pt_dvs_solver: too many iterations.
?pt_dvs_solver: debug files written:
state files: pt_ta.debug.cdf pt_tb.debug.cdf
control obj: ctrl.debug.cdf
?pt_dvs_solver: itmax= 400
?pt_dvs_solver: too many iterations.
?pt_dvs_solver: debug files written:
state files: pt_ta.debug.cdf pt_tb.debug.cdf
control obj: ctrl.debug.cdf
%check_save_state: nbflag = F
%check_save_state: iwrite_now = T
%check_save_state: check at wall_hours = 23.264946944444432
%wrstf: start call wrstf.
%wrstf: open new restart file:200741M35RS.DAT
%wrstf: open200741M35RS.DAT_NEW, status= 0
%wrstf: record, ier= 1 (ier=1, OK)
%wrstf: psload_save
%ps_rezone_chk: no rezone operator available (no MHD equilibrium).
%wrstf: write marker file.
%wrstf: transp_rplot_state_backup
% RESTART RECORD WRITTEN AT TA= 7.6340040E-01