Using randomized random numbers stream:  Random()= 0.8973175465

Starting a CMS calculation for Neptune.

Process ID: 50745

 modelIndex= 3
 Increasing nl from 512 to 513 to allow acceleration.
 CMS: Init() called with  nl= 513 nInt= 16 np= 16 nq= 24
 R_eq= 24766 km.
 P_target= 5.374302587e-08 PU = 1 bar.
 Rotation period= 16:6:40 = 58000 seconds.
 Rotational parameter:  qRot= 0.02608061413
 Target:         J2=  3529.4000 J4=   -35.8000
 newtonStepSize= 1


 Setting initial densities to polytropic EOS.
 Density at center = 1.00192181

 Parameters: MC_parameters: r1= 0.8681910000 r2/r1= 0.5905829194 r3/r1= 0.2437657370 H1= 2.2450339941 H2= 0.4418363789 H3= 0.1055113790

 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 ^^^^^ Parsing file logN3_state.txt model: 1 ^^^^^^
 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

 Parameters: MC_parameters: r1= 0.8681910000 r2/r1= 0.5905829194 r3/r1= 0.2437657370 H1= 2.2450339941 H2= 0.4418363789 H3= 0.1055113790

 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Opening file "logN3_lambda.txt" to read 1 columns of floating point values: Columns= 1 
Read 513 lines from file "logN3_lambda.txt".
Opening file "logN3_rho.txt" to read 1 columns of floating point values: Columns= 1 
Read 513 lines from file "logN3_rho.txt".
 n= 0 jTotal[n]= -1
 n= 1 jTotal[n]= 0
 n= 2 jTotal[n]= 0.003529390779
 n= 3 jTotal[n]= 0
 n= 4 jTotal[n]= -3.5816301e-05
 n= 5 jTotal[n]= 0
 n= 6 jTotal[n]= 5.40320236e-07
 n= 7 jTotal[n]= 0
 n= 8 jTotal[n]= -7.019567568e-09
 n= 9 jTotal[n]= 0
 n= 10 jTotal[n]= -8.93500222e-10
 n= 11 jTotal[n]= 0
 n= 12 jTotal[n]= -2.275565793e-09
 n= 13 jTotal[n]= 0
 n= 14 jTotal[n]= 1.131461671e-09
 n= 15 jTotal[n]= 0
 n= 16 jTotal[n]= -1.91352051e-09
 Target:         J2=  3529.4000 J4=   -35.8000
 Final: Jn*10^6: J2=  3529.3908 J4=   -35.8163 J6=     0.5403 J8=    -0.0070 J10=    -0.0009 J12=    -0.0023 J14=   0.001131 J16=  -0.001914
 Accuracy of equipotentials= 9.363975294e-07
Opening file "logN3_pressure.txt" to read 1 columns of floating point values: Columns= 1 
Read 513 lines from file "logN3_pressure.txt".
 Pressure accuracy= 9.709974303e-08

 Volume integration check: total mass      = 0.9999999331
 Moment of inertia                         = 0.2382657597
 Moment of inertia (Hubbard's Eq.5, no DR) = 0.2382657599
 Angular momentum                          = 0.03847871346

 Rotation period= 16:6:40 = 58000 seconds.

 ================================================================== 

Program finished properly.
