

   Setting up GPlates boundary velocity plugin.

   Input point 1 spherical coordinates: 1.571 4.870
   Input point 1 normalized cartesian coordinates: 0.157 -0.988 -0.000
   Input point 1 rotated model coordinates: 0.000 1.000 0.000
   Input point 2 spherical coordinates: 1.571 5.240
   Input point 2 normalized cartesian coordinates: 0.503 -0.864 -0.000
   Input point 2 rotated model coordinates: -0.362 0.932 -0.000

   Model will be rotated by -170.97 degrees around axis -0.00 0.00 -1.00
   The ParaView rotation angles are: -0.00 -0.00 -170.97
   The inverse ParaView rotation angles are: -0.00 -0.00 170.97

   Loading GPlates data boundary file ASPECT_DIR/data/boundary-velocity/gplates/current_day.gpml.


   Loading new gplates velocity file did not succeed.
   Assuming constant boundary conditions for rest of model run.

Number of active cells: 192 (on 3 levels)
Number of degrees of freedom: 2,832 (1,728+240+864)

*** Timestep 0:  t=0 years, dt=0 years
   Solving temperature system... 0 iterations.
   Solving Stokes system (GMG)... 65+0 iterations.

   Postprocessing:
     Writing graphical output:           output-burnman_seismic_test/solution/solution-00000
     RMS, max velocity:                  0.504 m/year, 3.08 m/year
     Temperature min/avg/max:            273 K, 2055 K, 3500 K
     Heat fluxes through boundary parts: -4.725e+07 W, 3.896e+07 W

Termination requested by criterion: end time



