Code: MAS Version: 0.66d_ccmc_v01 Updated on: 04/30/2003 Source file: mas47d_ccmc_v01.F Run ID: mas Run started on: 2009/04/07 Run started at: 14:22:00 Ran on machine: zmac.predsci.com Machine type: Darwin ### Parameters: NR = 61 NT = 71 NP = 64 NHISTMAX = 2002 ### Input file contents: &invars option='streamer' fldtype='potential' bingauss=.true. bnfile='br_m9_filt3.dat' eqtype='parker' np1d=130 onedfile='parker1.8.pw' rhor0=2.00 bcr0type='1dchar' bcr1type='1dchar' tmax=120 ntmax=40000 dtmax=0.050 dtmin=0.005 ifideal=0 slund=1.00e3 visc=0.01 rsifile=' ' rl=29. g0=.823 ifrho=1 iftemp=1 ifvdgv=1 ifpc=1 rfrac=.4,1. drratio=60.,4. nfrmesh=5 tfrac=.5,1. dtratio=.5,2. nftmesh=5 mmodes=0,1,2,3,4,5,6,7,8,9 tpltxint=120. ihistint=1 ifprec=1 trsdump=20. upwindv=1. cfl=.4 isitype=1 dformat='hdf' plotlist='br','bt','bp','jr','jt','jp','rho','t','p', 'vr','vt','vp' tnode=0.,10.0 vnode=0.,0. ishearprof=2 arotate=1. brotate=0. crotate=0. ihst=23 jhst=33 khst=1 parchar=.false. ubzero=.true. nfiltub=3 he_frac=0. radloss=0. tcond=0. ifaw=0 tnode_ch=0.,500. q0phys_ch=0.,0. tbc0=1.80e6 tbc1=0. upwinda=0. upwindar=1. upwindat=1. upwindap=1. upwindphi=1. strmax=.4 fdphi=.false. si_aggressive=.true. / ### Parameter values: &INVARS NTMAX = 40000, TMAX = 120.000000000000 , DTMAX = 5.000000000000000E-002, DTMIN = 5.000000000000000E-003, RL = 29.0000000000000 , G0 = 0.823000000000000 , IFVDGV = 1, SIMULT = 1.10000000000000 , DTMULT = 1.00000000000000 , STRMAX = 0.400000000000000 , ISITYPE = 1, CFL = 0.400000000000000 , OMEGAETA = 0.500000000000000 , MMODES = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 23*-1, UPWINDA = 0.000000000000000E+000, UPWINDV = 1.00000000000000 , PDMA = 1.00000000000000 , IFIDEAL = 0, IETATYPE = 1, SLUND = 1000.00000000000 , RMCELL = 0.000000000000000E+000, RMCELL_FACR = 1.00000000000000 , RMCELL_FACT = 1.00000000000000 , RMCELL_FACP = 1.00000000000000 , RMCELL_ETAMAX = 1.000000000000000E+020, RMCELL_UNIFT = F, NRETAP = 1, RETAP = 20*1.00000000000000 , ETAPR = 20*1.00000000000000 , NTETAP = 1, TETAP = 20*0.000000000000000E+000 , ETAPT = 20*1.00000000000000 , IVISTYPE = 1, VISC = 1.000000000000000E-002, RECELL = 0.000000000000000E+000, RECELL_FACR = 1.00000000000000 , RECELL_FACT = 1.00000000000000 , RECELL_FACP = 1.00000000000000 , RECELL_VISCMAX = 1.000000000000000E+020, RECELL_UNIFT = F, NRVISP = 1, RVISP = 20*1.00000000000000 , VISPR = 20*1.00000000000000 , NTVISP = 1, TVISP = 20*0.000000000000000E+000 , VISPT = 20*1.00000000000000 , IFRSOUT = 1, RSIFILE = , OPTION = streamer , EQTYPE = parker , FLDTYPE = potential , IFPREC = 1, ICGTYPE = 0, EPSCG = 1.000000000000000E-008, NCGMAX = 500, SIGMA = 1.00000000000000 , GAMMA = 1.05000000000000 , NCGHIST = 0, NWAVE = 1, LWAVE = 1, EPSWAVE = 1.000000000000000E-006, EPSCGA = 1.000000000000000E-008, DRRATIO = 60.0000000000000 , 4.00000000000000 , 8*1.00000000000000 , RFRAC = 0.400000000000000 , 9*1.00000000000000 , DTRATIO = 0.500000000000000 , 2.00000000000000 , 8*1.00000000000000 , TFRAC = 0.500000000000000 , 9*1.00000000000000 , NFRMESH = 5, NFTMESH = 5, IFRHO = 1, IFTEMP = 1, IFPC = 1, PREDFLOW = F, BETAPC_P = 0.500000000000000 , BETAPC_V = 0.500000000000000 , SI_AGGRESSIVE = T, FAC_CFLV = 1.00000000000000 , BETAPC_SI = 1.00000000000000 , IHISTINT = 1, THISTINT = 0.000000000000000E+000, IPLTXINT = 0, TPLTXINT = 120.000000000000 , ONEDFILE = parker1.8.pw , NP1D = 130, DFORMAT = hdf , PLOTLIST = br bt bp jr jt jp rho t p vr vt vp , B0 = 1.00000000000000 , BINGAUSS = T, RHOR0 = 2.00000000000000 , PR0 = 0.100000000000000 , DTHMAX = 0.698100000000000 , TH0 = 1.57079632700000 , TNODE = 0.000000000000000E+000, 10.0000000000000 , 98*-1.00000000000000 , VNODE = 100*0.000000000000000E+000 , RHOFAC1 = 4.000000000000000E-002, RHOFAC2 = 0.104000000000000 , IHST = 23, JHST = 33, KHST = 1, IRSDUMP = 0, TRSDUMP = 20.0000000000000 , IFOHMIC = 0, IFHVISC = 0, BCR1TYPE = 1dchar , BCR0TYPE = 1dchar , ISETRHO = 0, OMEGARHO = 0.100000000000000 , UPWINDAR = 1.00000000000000 , UPWINDAT = 1.00000000000000 , UPWINDAP = 1.00000000000000 , UPWINDPHI = 1.00000000000000 , IFEQCS = 0, DIPANGLE = 0.000000000000000E+000, BNFILE = br_m9_filt3.dat , GBC0ERR = 1.000000000000000E-005, PARCHAR = F, UBZERO = T, NFILTUB = 3, NFILTNE0 = 0, TCOND = 0.000000000000000E+000, RADLOSS = 0.000000000000000E+000, TBC0 = 1800000.00000000 , TBC1 = 0.000000000000000E+000, TNE0 = 1.00000000000000 , TNODE_CH = 0.000000000000000E+000, 500.000000000000 , 98*-1.00000000000000 , Q0_POLE_CH = 100*1.00000000000000 , Q0_EQ_CH = 100*1.00000000000000 , HLEN_POLE_CH = 100*0.500000000000000 , HLEN_EQ_CH = 100*0.500000000000000 , Q0PHYS_CH = 100*0.000000000000000E+000 , QALFVEN_CH = 100*0.000000000000000E+000 , TC_R = 10.0000000000000 , TC_DR = 5.00000000000000 , TC_FAC0 = 1.00000000000000 , TC_FAC1 = 1.00000000000000 , HE_FRAC = 0.000000000000000E+000, BR00 = 0.000000000000000E+000, IFAW = 0, RHO_AW = 1.00000000000000 , ISHEARPROF = 2, AROTATE = 1.00000000000000 , BROTATE = 0.000000000000000E+000, CROTATE = 0.000000000000000E+000, ALPHA_NOCOLL = 0.000000000000000E+000, TC_NOCOLL_R = 10.0000000000000 , TC_NOCOLL_DR = 5.00000000000000 , NERADMIN = -10000000000.0000 , THPROF_TH = 0.700000000000000 , THPROF_DTH = 8.000000000000000E-002, AWTHPROF_TH = 0.700000000000000 , AWTHPROF_DTH = 8.000000000000000E-002, EMGFLUX = F, BMIN = 5.000000000000000E-002, TNODE_EF = 100*-1.00000000000000 , BRFILE_EF = NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEF AULT NO_DEFAULT NO_DEFAULT NO_DEFAULT NO_DEFAULT , AWTHPROF2_TH = 2.44159265358979 , VRBCR0EF = exb , VTBCR0EF = zero , VPBCR0EF = zero , IFRADBC = 1, IFBHEAT = 0, BHEATFILE = brbyb2.dat , FDPHI = F, IFIMPLRAD = 1, IFETA_DYN = 0, ETA_MIN = 0.000000000000000E+000, ETA_MAX = 0.000000000000000E+000, ETA_JMIN = 30.0000000000000 , ETA_DJ = 5.00000000000000 , ETA_JSLOPE = 1.00000000000000 , BSQFLOOR = 0.000000000000000E+000, FJFRAC = 0.500000000000000 , T0_TC = 250000.000000000 , ALPHA_TC = 2.50000000000000 , DT0_TC = 20000.0000000000 , NFILT_KAPPA = 0, SPLIT_VISC = T, GAMMA1 = 0.000000000000000E+000, GAMMA2 = 0.000000000000000E+000, DRGAMMA = 0.000000000000000E+000, RGAMMA = 0.000000000000000E+000 / ### COMMENT from START: ### Upwind coefficients for the A equation: UPWINDAR = 1.00000000000000 UPWINDAT = 1.00000000000000 UPWINDAP = 1.00000000000000 ### COMMENT from GENMESH: ### Mesh information for coordinate r: Number of mesh segments = 2 Segment chi r 0.0000 1.0000 1 0.8573 12.6000 2 1.0000 30.0000 Mesh point location: i r dr dr(i)/dr(i-1) 1 1.00000 2 1.01782 0.01782 3 1.03633 0.01852 1.03920 4 1.05609 0.01975 1.06671 5 1.07741 0.02132 1.07961 6 1.10051 0.02311 1.08354 7 1.12557 0.02505 1.08426 8 1.15273 0.02716 1.08431 9 1.18219 0.02945 1.08431 10 1.21412 0.03194 1.08431 11 1.24875 0.03463 1.08431 12 1.28630 0.03755 1.08431 13 1.32702 0.04072 1.08431 14 1.37117 0.04415 1.08431 15 1.41904 0.04787 1.08431 16 1.47095 0.05191 1.08431 17 1.52724 0.05629 1.08431 18 1.58827 0.06103 1.08431 19 1.65445 0.06618 1.08431 20 1.72620 0.07176 1.08431 21 1.80401 0.07781 1.08431 22 1.88838 0.08437 1.08431 23 1.97986 0.09148 1.08431 24 2.07905 0.09919 1.08431 25 2.18661 0.10756 1.08431 26 2.30323 0.11663 1.08431 27 2.42969 0.12646 1.08431 28 2.56681 0.13712 1.08431 29 2.71549 0.14868 1.08431 30 2.87671 0.16122 1.08431 31 3.05152 0.17481 1.08431 32 3.24107 0.18955 1.08431 33 3.44661 0.20553 1.08431 34 3.66947 0.22286 1.08431 35 3.91112 0.24165 1.08431 36 4.17315 0.26203 1.08431 37 4.45727 0.28412 1.08431 38 4.76534 0.30808 1.08431 39 5.09939 0.33405 1.08431 40 5.46161 0.36222 1.08431 41 5.85436 0.39276 1.08431 42 6.28023 0.42587 1.08431 43 6.74201 0.46178 1.08431 44 7.24273 0.50071 1.08431 45 7.78566 0.54293 1.08431 46 8.37436 0.58871 1.08431 47 9.01271 0.63835 1.08433 48 9.70503 0.69232 1.08454 49 10.45688 0.75185 1.08599 50 11.27745 0.82056 1.09139 51 12.18351 0.90606 1.10420 52 13.20262 1.01911 1.12477 53 14.37217 1.16955 1.14761 54 15.73497 1.36280 1.16524 55 17.33577 1.60079 1.17464 56 19.22096 1.88520 1.17766 57 21.43603 2.21507 1.17498 58 24.00748 2.57145 1.16089 59 26.90362 2.89614 1.12627 60 30.00000 3.09638 1.06914 ### COMMENT from GENMESH: ### Mesh information for coordinate t: Number of mesh segments = 2 Segment chi t 0.0000 0.0000 1 0.5000 1.5708 2 1.0000 3.1416 Mesh point location: i t dt dt(i)/dt(i-1) 1 0.00000 2 0.06144 0.06144 3 0.12233 0.06089 0.99109 4 0.18229 0.05996 0.98470 5 0.24113 0.05885 0.98146 6 0.29882 0.05769 0.98036 7 0.35537 0.05654 0.98013 8 0.41079 0.05542 0.98011 9 0.46511 0.05432 0.98011 10 0.51834 0.05324 0.98011 11 0.57052 0.05218 0.98011 12 0.62166 0.05114 0.98011 13 0.67179 0.05012 0.98011 14 0.72091 0.04913 0.98011 15 0.76906 0.04815 0.98011 16 0.81625 0.04719 0.98011 17 0.86251 0.04625 0.98011 18 0.90784 0.04533 0.98011 19 0.95227 0.04443 0.98011 20 0.99582 0.04355 0.98011 21 1.03850 0.04268 0.98011 22 1.08033 0.04183 0.98011 23 1.12133 0.04100 0.98011 24 1.16151 0.04018 0.98011 25 1.20090 0.03939 0.98011 26 1.23950 0.03860 0.98011 27 1.27734 0.03783 0.98011 28 1.31442 0.03708 0.98011 29 1.35076 0.03634 0.98011 30 1.38638 0.03562 0.98011 31 1.42129 0.03491 0.98011 32 1.45552 0.03422 0.98018 33 1.48907 0.03356 0.98067 34 1.52205 0.03297 0.98254 35 1.55460 0.03255 0.98721 36 1.58699 0.03240 0.99521 37 1.61955 0.03255 1.00482 38 1.65252 0.03297 1.01296 39 1.68608 0.03356 1.01778 40 1.72030 0.03422 1.01971 41 1.75521 0.03491 1.02022 42 1.79083 0.03562 1.02029 43 1.82718 0.03634 1.02029 44 1.86426 0.03708 1.02029 45 1.90209 0.03783 1.02029 46 1.94069 0.03860 1.02029 47 1.98008 0.03939 1.02029 48 2.02026 0.04018 1.02029 49 2.06126 0.04100 1.02029 50 2.10309 0.04183 1.02029 51 2.14578 0.04268 1.02029 52 2.18932 0.04355 1.02029 53 2.23375 0.04443 1.02029 54 2.27909 0.04533 1.02029 55 2.32534 0.04625 1.02029 56 2.37253 0.04719 1.02029 57 2.42068 0.04815 1.02029 58 2.46981 0.04913 1.02029 59 2.51993 0.05012 1.02029 60 2.57107 0.05114 1.02029 61 2.62325 0.05218 1.02029 62 2.67649 0.05324 1.02029 63 2.73080 0.05432 1.02029 64 2.78622 0.05542 1.02029 65 2.84277 0.05654 1.02029 66 2.90046 0.05769 1.02027 67 2.95931 0.05885 1.02004 68 3.01927 0.05996 1.01889 69 3.08016 0.06089 1.01554 70 3.14159 0.06144 1.00899 ### COMMENT from LOAD1D: ### A 1D equilibrium was read in. File name: parker1.8.pw Base density, RHOR0 = 2.00000000000000 ### COMMENT from READBR: ### Reading Br file: br_m9_filt3.dat NTIME = 0 TIME = 0.000000000000000E+000 ### COMMENT from READBR: ### Normalization factor for the normal magnetic field: Maximum field magnitude before normalization = 52.3813305169217 Maximum field magnitude after normalization = 23.7353458213313 Scaling factor for the magnetic field, BSCALE = 2.20689139780074 [Original field = BSCALE * normalized field.] ### COMMENT from READBR: ### Computed flux balance: Positive flux = 7.32550866510453 Negative flux = -7.34559559562492 ### COMMENT from READBR: ### Flux balance correction: BR00 (monopole Br field magnitude) = -1.597967639341078E-003 ### COMMENT from POTFLD: The (interpolated) photospheric Br field has been written to file: br_photo.hdf ### COMMENT from LOADETA: ### Uniform resistivity profile used. IETATYPE = 1 SLUND = 1000.00000000000 ### COMMENT from LOADVISC: ### Uniform viscosity profile used. IVISTYPE = 1 VISC = 1.000000000000000E-002 ### COMMENT from START: The point for collection of history diagnostics is located at: IHST = 23 , R = 1.97985566111912 JHST = 33 , T = 1.48907440491756 KHST = 1 , P = 0.000000000000000E+000 ### COMMENT from PDUMP: ### Wrote out the requested fields. File sequence number = 1 NTIME = 0 TIME = 0.000000000000000E+000 ### COMMENT from WRRSTRT: Writing to restart file: rsmas001 NTIME = 451 TIME = 20.0269095499109 ### COMMENT from WRRSTRT: Writing to restart file: rsmas002 NTIME = 980 TIME = 40.0237302126526 ### COMMENT from WRRSTRT: Writing to restart file: rsmas003 NTIME = 1513 TIME = 60.0360734909224 ### COMMENT from DUMPHIST: Wrote time histories to file: hmas01 Wrote v histories to file: vmas01 NTIME = 2002 TIME = 78.3827682446249 IHIST = 2002 ### COMMENT from WRRSTRT: Writing to restart file: rsmas004 NTIME = 2046 TIME = 80.0346981073220 ### COMMENT from WRRSTRT: Writing to restart file: rsmas005 NTIME = 2578 TIME = 100.026698274375 ### COMMENT from PDUMP: ### Wrote out the requested fields. File sequence number = 2 NTIME = 3109 TIME = 120.012477343323 ### COMMENT from WRRSTRT: Writing to restart file: rsmas006 NTIME = 3109 TIME = 120.012477343323 ### COMMENT from DUMPHIST: Wrote time histories to file: hmas02 Wrote v histories to file: vmas02 NTIME = 3109 TIME = 120.012477343323 IHIST = 1108 ### COMMENT from WRRSTRT: Writing to restart file: rsmas NTIME = 3109 TIME = 120.012477343323 Code: MAS Version: 0.66d_ccmc_v01 Updated on: 04/30/2003 Source file: mas47d_ccmc_v01.F Run ID: mas Run ended on: 2009/04/07 Run ended at: 16:06:18 Ran on machine: zmac.predsci.com Machine type: Darwin ### End of run ... NTIME = 3109 TIME = 1.2001248E+02 ### Time advance information: Average time step used = 3.860163311139354E-002 ----------------------------------------- ### Field solver convergence information: ----------------------------------------- ### Semi-implicit advance: Total number of solves = 410388 Total number of iterations = 14010095 Average number of iterations per solve = 34 ### Implicit viscosity advance: Total number of solves = 205194 Total number of iterations = 3452759 Average number of iterations per solve = 16 ### Implicit resistivity advance: Total number of solves = 202085 Total number of iterations = 376072 Average number of iterations per solve = 1 ### Implicit thermal conductivity advance: Total number of solves = 0 Total number of iterations = 0 Average number of iterations per solve = 0 ### Potential field solution: Total number of solves = 63 Total number of iterations = 4140 Average number of iterations per solve = 65 ### CPU time used = 6243.33 seconds.