*** VMEC6 Fixed Boundary Up-Down Asymmetric Equilibrium Code *** *** S.P. Hirshman (ORNL) & R.M.Wieland (PPPL) *** COMPUTATION PARAMETERS: (u = theta, v = zeta) --------------------------------------------- ns nu nv mu mv 201 22 1 9 0 CONFIGURATION PARAMETERS: ------------------------- nfp idata gamma phiedge(wb) curtor(A) betafac 0 1 0.000E+00 -2.555E+00 1.230E+00 1.000E+00 RUN CONTROL PARAMETERS: ----------------------- ncurr niter nsin nstep nvacskip ftol 0 2000 11 100 0 1.00E-09 R-Z FOURIER BOUNDARY COEFFICIENTS R = RBC*cos(m*u - n*v) + RBS*sin(m*u - n*v), Z = ZBC*cos(m*u - n*v) + ZBS*sin(m*u-n*v) -------------------------------------------------------------------------------------- mb nb rbc rbs zbc zbs raxis zaxis 0 0 1.6165E+00 0.0000E+00 -4.3074E-02 0.0000E+00 1.6839E+00 5.7198E-02 1 0 4.5385E-01 5.5993E-02 5.5993E-02 8.6838E-01 2 0 3.1969E-02 4.3972E-02 4.8108E-03 -2.0702E-02 3 0 2.0171E-02 -8.4784E-03 -2.9994E-03 -4.5277E-03 4 0 -4.0099E-03 -1.4940E-03 -3.3675E-04 3.5956E-03 5 0 -8.3469E-04 2.3719E-03 1.6493E-03 1.3282E-03 LEGEND ------ FSQR, FSQZ = Normalized Physical Force Residuals fsqr, fsqz = Preconditioned Force Residuals ----------------------- BEGIN FORCE ITERATIONS ----------------------- VMEC6 > betafac= 1. NS = 11 NO. FOURIER MODES = 9 ITER FSQR FSQZ FSQL fsqr fsqz DELT R00(0) Z00(0) WMHD BETA PHIFAC CURTOR CTOR_ ITER FSQR FSQZ FSQL fsqr fsqz DELT R00(0) Z00(0) WMHD BETA PHIFAC CURTOR CTOR_ 1 1.52E-01 2.69E-02 3.76E-02 3.17E-03 3.85E-04 1.10E+00 1.684E+00 5.720E-02 2.88310E+01 4.62E-03 1.000 1.23 1.00 100 2.19E-06 8.80E-07 8.82E-07 2.46E-10 1.83E-09 1.10E+00 1.679E+00 5.488E-02 2.71084E+01 4.97E-03 1.000 1.23 1.21 108 9.87E-07 4.59E-07 2.62E-07 8.83E-11 4.69E-10 1.10E+00 1.680E+00 5.539E-02 2.71083E+01 4.97E-03 1.000 1.23 1.21 VMEC6> COMPUTATIONAL TIME FOR NS= 11 IS 7.13E-02 SECONDS VMEC6> JACOBIAN RESETS = 0 VMEC6> ns= 11 :d(ln W)/dt = 3.172E-07 d(ln R0)/dt = 2.605E-05 NS = 31 NO. FOURIER MODES = 9 VMEC6> Turning ON Ip-edge Matching by varying Phi-edge. ITER FSQR FSQZ FSQL fsqr fsqz DELT R00(0) Z00(0) WMHD BETA PHIFAC CURTOR CTOR_ ITER FSQR FSQZ FSQL fsqr fsqz DELT R00(0) Z00(0) WMHD BETA PHIFAC CURTOR CTOR_ 1 7.29E-02 4.95E-02 1.38E-04 1.23E-08 3.34E-08 1.10E+00 1.680E+00 5.545E-02 2.71079E+01 5.02E-03 1.000 1.23 1.22 100 1.05E-06 4.49E-07 1.62E-08 3.85E-12 1.94E-11 1.10E+00 1.683E+00 5.756E-02 2.66414E+01 5.12E-03 0.991 1.23 1.23 200 3.92E-08 9.62E-09 1.87E-09 1.06E-13 2.26E-13 1.10E+00 1.684E+00 5.758E-02 2.57867E+01 5.29E-03 0.975 1.23 1.23 300 1.06E-08 2.57E-09 1.86E-10 4.32E-15 1.07E-14 1.10E+00 1.684E+00 5.770E-02 2.54901E+01 5.35E-03 0.970 1.23 1.23 400 3.61E-09 7.63E-10 2.73E-11 5.50E-16 1.08E-15 1.10E+00 1.684E+00 5.775E-02 2.53357E+01 5.38E-03 0.967 1.23 1.23 500 1.27E-09 2.56E-10 6.97E-12 1.13E-16 2.11E-16 1.10E+00 1.684E+00 5.776E-02 2.52598E+01 5.40E-03 0.966 1.23 1.23 519 9.92E-10 2.02E-10 5.01E-12 8.38E-17 1.53E-16 1.10E+00 1.684E+00 5.777E-02 2.52471E+01 5.40E-03 0.965 1.23 1.23 VMEC6> COMPUTATIONAL TIME FOR NS= 31 IS 5.41E-01 SECONDS VMEC6> JACOBIAN RESETS = 0 VMEC6> ns= 31 :d(ln W)/dt = 2.622E-05 d(ln R0)/dt = 2.462E-07 VMEC6> NUMBER OF JACOBIAN RESETS = 0 VMEC6> TOTAL COMPUTATIONAL TIME : 0.61 SECONDS VMEC6> TIME IN VACUUM LOOP : 0.00 SECONDS VMEC6> TIME TO READ IN DATA: 0.00 SECONDS VMEC6> TIME TO OUTPUT DATA : 0.00 SECONDS VMEC6> TIME IN PROFILE RECONSTRUCTION: 0.00 SECONDS %VMEC6 Rel. Grad-Shaf. error (excl. edge): 0.138148 * MHD EQUILIBRIUM CALCULATED, CPU TIME = 6.3672E-01 SECONDS DATA R*BT AT EDGE: 3.3093E+02 ->EQBDY_CHECK: last flux surface curvature ratio is: 1.1015E-01 %get_rygrid: no grids for Psi(R,Y) % GFRAME - GEOMETRY TIMESTEP TG1= 2.195000 TO TG2= 2.200000 @ NSTEP 1953 GFRAME TG2 MOMENTS CHECKSUM: 9.0982601503216E+03 *** FULL NEUTRAL SOURCE CALCULATION PERFORMED *** NSTEP= 1953 TA= 2.19500E+00 CPU TIME= 2.40234E+00 SECONDS. DT= 4.26185E-04 *** FULL NEUTRAL SOURCE CALCULATION PERFORMED *** NSTEP= 1954 TA= 2.19543E+00 CPU TIME= 2.27246E+00 SECONDS. DT= 5.32731E-04 *** FULL NEUTRAL SOURCE CALCULATION PERFORMED *** NSTEP= 1960 TA= 2.19746E+00 CPU TIME= 2.24512E+00 SECONDS. DT= 4.22985E-04 *** FULL NEUTRAL SOURCE CALCULATION PERFORMED *** NSTEP= 1966 TA= 2.19957E+00 CPU TIME= 2.34375E+00 SECONDS. DT= 4.31988E-04 %MFRCHK - LABEL "BALE0_SGF", # 3= -1.01021E-39 RESET TO ZERO %MFRCHK - LABEL "TQBA0_SGF", # 1= -2.01268E-41 RESET TO ZERO %MFRCHK - LABEL "BALE0_GF_D", # 3= -1.01021E-39 RESET TO ZERO %MFRCHK - LABEL "TQBA0_GF_D", # 1= -2.01268E-41 RESET TO ZERO %MFRCHK - LABEL "BALE0_SRC", # 3= -1.83675E-40 RESET TO ZERO %MFRCHK - LABEL "TQBA0_SRC", # 1= -2.20873E-41 RESET TO ZERO %MFRCHK - LABEL "BALE0_RC_D", # 3= -1.83675E-40 RESET TO ZERO %MFRCHK - LABEL "TQBA0_RC_D", # 1= -2.20873E-41 RESET TO ZERO %(TRANSP)elvis_request: requesting ElVis data transfer. !elvis_request: NLEND=.true. so sending elvis CLOSE command 4.2 Call trmpi_end (NORMAL EXIT) trmpi_end2 -I- __MPI undefined 4.2 TERMINATE THE RUN (NORMAL EXIT) CPU TIME USED (hours): 3.40441E+00 %kill_nubeam_server: no server READY file found. ------------ stderr ---------------- (mpi_share_env) process myid= 0 cwd: /local/tr_mcavedon/transp_compute/AUGD/18935A02 ==>runtrx: TRANSP run successful ==========(runtrx)====================== ==========TRANSP output conversion====== date: Wed May 23 14:06:26 EDT 2012 ( swift001.pppl.gov ) ==========(runtrx)====================== ---------------> starting: plotcon 18935A02 2012/05/23:14:06:26 %initcpl: MDS_CACHE enabled. %PoPlot -- reading .PLN files %POPLT2-- PROCESSING RUN 18935A02 SHOT NO. 18935 EXPECT 416 SCALAR FCNS, 1020 PROFILE FCNS OF TIME "MF" FILE RECORD SIZE = 40 WORDS (FLOATING PT) dmg_datbuf_expand call from dmgini_sized: isize= 0 18935A02MF.PLN size = 18M %initcpl: MDS_CACHE enabled. [mds_cache_disable: MDS+ cache disabled.] dmg_datbuf_expand call from dmgini_sized: isize= 0 (retry folding filename to lowercase) ...reading TF.PLN header data... cdfcon: NETcdf file datestamp : 23-May-2012 14:06:36 build_date: call sget_env build_date: call ufopen xshare_build.dat cdfcon: Transp common build date : Wed Mar 21 05:53:07 EDT Define Dimensions 8 define Scalar Fct 416 Define Multi Graphs 433 Write Profiles 1020 X 1 1 40 XB 2 2 40 THETA 3 7 80 RMJSYM 4 20 165 RMAJM 5 22 81 MCINDX 6 23 220 IBOLO 7 25 25 ILIM 8 31 15 read NF File : 417 416 Write Multigraph: 433 ...readback test of .CDF file... 1871 variables, 10 dimensions 9 att ...header check SUCCESSFUL; now check data. ...check profile data... plotcon: CDF files: /local/tr_mcavedon/transp_compute/AUGD/18935A02/18935A02.CDF /local/tr_mcavedon/transp_compute/AUGD/18935A02/18935A02PH.CDF %targz_pseq: no directory: 18935A02_replay (normal exit) --------------->plotcon: normal exit. 2012/05/23:14:06:36 ==>runtrx: TRANSP postprocessing OK ==========(runtrx)====================== ==========TRANSP finish and cleanup===== date: Wed May 23 14:06:36 EDT 2012 ( swift001.pppl.gov ) ==========(runtrx)====================== %finishup: terminating runtime transp with send_tr_elvis CLOSE %finishup -I- pppl.gov production run %finishup: copying TRANSP permanent output files to /u/tr_mcavedon/transp/result/AUGD.04 mv 18935A02_AUGD.REQUEST /u/tr_mcavedon/transp/result/AUGD.04/18935A02_AUGD.REQUEST mv 18935A02CC.TMP /u/tr_mcavedon/transp/result/AUGD.04/18935A02CC.TMP mv 18935A02.CDF /u/tr_mcavedon/transp/result/AUGD.04/18935A02.CDF mv 18935A02ex.for /u/tr_mcavedon/transp/result/AUGD.04/18935A02ex.for mv 18935A02_nubeam_init.dat /u/tr_mcavedon/transp/result/AUGD.04/18935A02_nubeam_init.dat mv 18935A02PH.CDF /u/tr_mcavedon/transp/result/AUGD.04/18935A02PH.CDF mv 18935A02.target /u/tr_mcavedon/transp/result/AUGD.04/18935A02.target mv 18935A02TR.DAT /u/tr_mcavedon/transp/result/AUGD.04/18935A02TR.DAT mv 18935A02.trexe /u/tr_mcavedon/transp/result/AUGD.04/18935A02.trexe mv 18935A02TR.INF /u/tr_mcavedon/transp/result/AUGD.04/18935A02TR.INF %finishup: retaining 18935A02tr.log mv 18935A02TR.MSG /u/tr_mcavedon/transp/result/AUGD.04/18935A02TR.MSG rm: No match. %finishup: cp -f /local/tr_mcavedon/transp_tmp/AUGD.04_18935A02.tar.gz /u/tr_mcavedon/transp/result/AUGD.04/AUGD.04_18935A02.tar.gz %finishup: wrote /u/tr_mcavedon/transp/result/AUGD.04/AUGD.04_18935A02.FILESREADY ==========(runtrx)====================== ==========>runtrx normal exit<========== date: Wed May 23 14:06:39 EDT 2012 ( swift001.pppl.gov ) ==========>runtrx runsite = pppl.gov<======