** ====================================================================================== 
**
**             ___________                          _______
**     _______/          /_____ ________ __ _   ___/       \ ____
**    /      /___    ___/      /__   __/  /  \ /  /\__ _   /     \
**   /  //  /   /   /  /  //  /  /  / /  /    /  /___/_   /  //  /
**  /  ___ /   /   /  /  _   /  /  / /  /  /    //       /  //  /
** /__/       /___/  /__//__/  /__/ /__/__/ \__//_______/______/
**
** Authors:  Dave A. May          (dave.may@erdw.ethz.ch)           
**           Laetitia Le Pourhiet (laetitia.le_pourhiet@upmc.fr)    
**           Jed Brown            (jedbrown@mcs.anl.gov)            
**
** git url: https://bitbucket.org/jedbrown/ptatin3d.git 
** commit hash: [out-of-date] Execute "make releaseinfo" to update to the most recent revision 
** log: [out-of-date] Execute "make releaseinfo" to update to the most recent revision 
**                                                                       
** TATIN_CFLAGS = -std=gnu99 -O0 -g -Wall -Wno-switch -Wno-unused-function -Wno-incompatible-pointer-types-discards-qualifiers -Wno-tautological-constant-out-of-range-compare -Wno-enum-conversion -Wno-type-safety
**                                                                       
** WARNING pTatin3d appears to have been compiled with debug options 
** For significant performance improvements, please consult the file makefile.arch  
** Adjust TATIN_CFLAGS to include aggressive compiler optimizations 
**                                                                       
** ====================================================================================== 
RheologyConstantsInitialise: global viscosity cut-off, min= 1.000000e-100, max = 1.000000e+100  
DataBucketView(SEQ): ("Material Constants")
  L                  = 200 
  buffer             = 0 
  allocated          = 200 
  nfields registered = 17 
    [  0]: field name  ==>>     MaterialConst_MaterialType : Mem. usage = 3.20e-03 (MB) 
    [  1]: field name  ==>>   MaterialConst_ViscosityConst : Mem. usage = 1.60e-03 (MB) 
    [  2]: field name  ==>>       MaterialConst_ViscosityZ : Mem. usage = 4.80e-03 (MB) 
    [  3]: field name  ==>>    MaterialConst_ViscosityArrh : Mem. usage = 1.28e-02 (MB) 
    [  4]: field name  ==>>      MaterialConst_ViscosityFK : Mem. usage = 3.20e-03 (MB) 
    [  5]: field name  ==>>     MaterialConst_PlasticMises : Mem. usage = 3.20e-03 (MB) 
    [  6]: field name  ==>>        MaterialConst_PlasticDP : Mem. usage = 9.60e-03 (MB) 
    [  7]: field name  ==>>     MaterialConst_DensityConst : Mem. usage = 1.60e-03 (MB) 
    [  8]: field name  ==>> MaterialConst_DensityBoussinesq : Mem. usage = 4.80e-03 (MB) 
    [  9]: field name  ==>>          MaterialConst_SoftLin : Mem. usage = 3.20e-03 (MB) 
    [ 10]: field name  ==>>         MaterialConst_SoftExpo : Mem. usage = 3.20e-03 (MB) 
    [ 11]: field name  ==>>        EnergyMaterialConstants : Mem. usage = 1.44e-02 (MB) 
    [ 12]: field name  ==>>        EnergyConductivityConst : Mem. usage = 1.60e-03 (MB) 
    [ 13]: field name  ==>>    EnergyConductivityThreshold : Mem. usage = 6.40e-03 (MB) 
    [ 14]: field name  ==>>              EnergySourceConst : Mem. usage = 1.60e-03 (MB) 
    [ 15]: field name  ==>>              EnergySourceDecay : Mem. usage = 3.20e-03 (MB) 
    [ 16]: field name  ==>> EnergySourceAdiabaticAdvection : Mem. usage = 1.60e-03 (MB) 
  Total mem. usage                                                      = 8.00e-02 (MB) 
[pTatin] Created output directory: analytics_solcx 
[pTatin] Created log file: analytics_solcx/ptatin.log-2017.07.30_16:35:18 
[pTatin] Writing output to existing directory: analytics_solcx 
[pTatin] Created options file: analytics_solcx/ptatin.options-2017.07.30_16:35:18 
[pTatin] Created options file: analytics_solcx/ptatin.options 
  [pTatinModel]: Registering model [0] with name "template"
  [pTatinModel]: Registering model [1] with name "viscous_sinker"
  [pTatinModel]: Registering model [2] with name "Gene3D"
  [pTatinModel]: Registering model [3] with name "Gene3DNueve"
  [pTatinModel]: Registering model [4] with name "indentor"
  [pTatinModel]: Registering model [5] with name "rift3D"
  [pTatinModel]: Registering model [6] with name "rift3D_T"
  [pTatinModel]: Registering model [7] with name "sierra"
  [pTatinModel]: Registering model [8] with name "advdiff_example"
  [pTatinModel]: Registering model [9] with name "delamination"
  [pTatinModel]: Registering model [10] with name "Riftrh"
  [pTatinModel]: Registering model [11] with name "Rift_oblique3d"
  [pTatinModel]: Registering model [12] with name "geomod2008"
  [pTatinModel]: Registering model [13] with name "multilayer_folding"
  [pTatinModel]: Registering model [14] with name "submarinelavaflow"
  [pTatinModel]: Registering model [15] with name "ex_subduction"
  [pTatinModel]: Registering model [16] with name "iplus"
  [pTatinModel]: Registering model [17] with name "subduction_initiation2d"
  [pTatinModel]: Registering model [18] with name "convection2d"
  [pTatinModel]: Registering model [19] with name "thermal_sb"
  [pTatinModel]: Registering model [20] with name "sd3d"
  [pTatinModel]: Registering model [21] with name "pas"
  [pTatinModel]: Registering model [22] with name "pd"
  [pTatinModelDynamic]: Dynamically registering model with name "static_box"
  [pTatinModelDynamic]: Dynamically registering model with name "static_box_thermomech"
  [pTatinModelDynamic]: Dynamically registering model with name "analytics_vv"
  [pTatinModel]: Registering model [23] with name "analytics_vv"
  [pTatinModel]: -ptatin_model "analytics_vv" was detected
pTatin3dRestart: Required to specify a suffix for your restart file via -restart_prefix
pTatin3dRestart: Unable to restart job
BCList: Mem. usage (min,max) = 8.11e-02,8.11e-02 (MB) 
VolumeQuadratureCreate_GaussLegendreStokes:
	Using 3x3 pnt Gauss Legendre quadrature
DataBucketView(SEQ): ("GaussLegendre StokesCoefficients")
  L                  = 1944 
  buffer             = 1 
  allocated          = 1945 
  nfields registered = 1 
    [  0]: field name  ==>>              QPntVolCoefStokes : Mem. usage = 9.34e-02 (MB) 
  Total mem. usage                                                      = 9.34e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 0]")
  L                  = 108 
  buffer             = 1 
  allocated          = 109 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 1.22e-02 (MB) 
  Total mem. usage                                                      = 1.22e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 1]")
  L                  = 108 
  buffer             = 1 
  allocated          = 109 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 1.22e-02 (MB) 
  Total mem. usage                                                      = 1.22e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 2]")
  L                  = 108 
  buffer             = 1 
  allocated          = 109 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 1.22e-02 (MB) 
  Total mem. usage                                                      = 1.22e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 3]")
  L                  = 108 
  buffer             = 1 
  allocated          = 109 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 1.22e-02 (MB) 
  Total mem. usage                                                      = 1.22e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 4]")
  L                  = 324 
  buffer             = 1 
  allocated          = 325 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 3.64e-02 (MB) 
  Total mem. usage                                                      = 3.64e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 5]")
  L                  = 324 
  buffer             = 1 
  allocated          = 325 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 3.64e-02 (MB) 
  Total mem. usage                                                      = 3.64e-02 (MB) 
  MaterialPointsStokes: Using one2one projection
[[Swarm initialization: 0.0002 (sec)]]
SwarmMPntStd_AssignUniquePointIdentifiers : max_pid = 0 
[[Swarm->coordinate assignment: 1944 points : 0.0006 (sec)]]
************************** Starting _DataExCompleteCommunicationMap ************************** 
max_nnz = 0 
Mat Object: 1 MPI processes
  type: seqaij
  rows=1, cols=1
  total: nonzeros=0, allocated nonzeros=0
  total number of mallocs used during MatSetValues calls =0
    not using I-node routines
************************** Ending _DataExCompleteCommunicationMap [setup time: 1.7118e-04 (sec)] ************************** 
[[SwarmDMDA3dDataExchangerCreate: time = 1.9193e-04 (sec)]]
************************** Starting _DataExCompleteCommunicationMap ************************** 
max_nnz = 0 
Mat Object: 1 MPI processes
  type: seqaij
  rows=1, cols=1
  total: nonzeros=0, allocated nonzeros=0
  total number of mallocs used during MatSetValues calls =0
    not using I-node routines
************************** Ending _DataExCompleteCommunicationMap [setup time: 5.9128e-05 (sec)] ************************** 
[[SwarmDMDA3dDataExchangerCreate: time = 7.2002e-05 (sec)]]
[PSwarmSetUpCoords_FillBox]
Mesh size (6 x 6 x 2) : MG levels 2  
         level [ 0]: global Q2 elements (3 x 3 x 1) 
         level [ 1]: global Q2 elements (6 x 6 x 2) 
[r   0]: level [ 0]: local Q2 elements  (3 x 3 x 1) 
[r   0]: level [ 1]: local Q2 elements  (6 x 6 x 2) 
[r   0]: level [ 0]: element range [0 - 2] x [0 - 2] x [0 - 0] 
[r   0]: level [ 1]: element range [0 - 5] x [0 - 5] x [0 - 1] 
VolumeQuadratureCreate_GaussLegendreStokes:
	Using 3x3 pnt Gauss Legendre quadrature
DataBucketView(SEQ): ("GaussLegendre StokesCoefficients")
  L                  = 243 
  buffer             = 1 
  allocated          = 244 
  nfields registered = 1 
    [  0]: field name  ==>>              QPntVolCoefStokes : Mem. usage = 1.17e-02 (MB) 
  Total mem. usage                                                      = 1.17e-02 (MB) 
BCList: Mem. usage (min,max) = 1.41e-02,1.41e-02 (MB) 
*** Rheology update for RHEOLOGY_VISCOUS selected ***
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-02,1.00e+02; log10(max/min) 4.00e+00; cpu time 3.91e-05 (sec)
Level [1]: Coarse grid type :: Re-discretisation :: assembled operator 
Level [0]: Coarse grid type :: Re-discretisation :: assembled operator 
   [[ COMPUTING FLOW FIELD FOR STEP : 0 ]]
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-02,1.00e+02; log10(max/min) 4.00e+00; cpu time 6.79e-05 (sec)
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-02,1.00e+02; log10(max/min) 4.00e+00; cpu time 6.60e-05 (sec)
  0 KSP Component U,V,W,P residual norm [ 0.000000000000e+00, 2.304766918639e-02, 0.000000000000e+00, 0.000000000000e+00 ]
  1 KSP Component U,V,W,P residual norm [ 3.959593366510e-14, 2.288427628901e-02, 1.782174455335e-14, 1.933682550788e-03 ]
  2 KSP Component U,V,W,P residual norm [ 2.904052760644e-12, 1.823811781865e-02, 1.257465193425e-12, 9.365744204264e-03 ]
  3 KSP Component U,V,W,P residual norm [ 2.371173998711e-11, 2.448885217989e-05, 9.943277106528e-12, 7.508736644163e-04 ]
  4 KSP Component U,V,W,P residual norm [ 2.896658994316e-11, 2.517274620152e-08, 1.211078779337e-11, 2.408676494360e-05 ]
  5 KSP Component U,V,W,P residual norm [ 4.145930903039e-11, 1.631222660687e-10, 1.732892317483e-11, 1.918547641301e-06 ]
  6 KSP Component U,V,W,P residual norm [ 2.107356308714e-11, 2.843835510181e-11, 8.805943481707e-12, 7.262614527762e-07 ]
  7 KSP Component U,V,W,P residual norm [ 4.412963620197e-13, 3.811516472117e-13, 1.824221631894e-13, 5.782873544888e-08 ]
  8 KSP Component U,V,W,P residual norm [ 3.849387534180e-13, 3.070895497752e-13, 1.585312725285e-13, 4.426225256866e-10 ]
  9 KSP Component U,V,W,P residual norm [ 3.879607760933e-13, 3.095014704450e-13, 1.597768585384e-13, 9.604546639610e-13 ]
 10 KSP Component U,V,W,P residual norm [ 1.783346817241e-13, 1.422677728558e-13, 7.344438305988e-14, 2.689522058759e-13 ]
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-02,1.00e+02; log10(max/min) 4.00e+00; cpu time 4.39e-05 (sec)
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-02,1.00e+02; log10(max/min) 4.00e+00; cpu time 6.41e-05 (sec)
  writing pvdfilename analytics_solcx/timeseries_passive_pswarm.pvd 
  writing pvdfilename analytics_solcx/timeseries_vp.pvd 
[[DESIGN FLAW]] pTatin3d_ModelOutput_VelocityPressure_Stokes: require better physics modularity to extract (u,p) <---| (X) 
[[DESIGN FLAW]] pTatinOutputMeshVelocityPressureVTS_v0_binary: only printing P0 component of pressure field 
pTatin3d_ModelOutput_VelocityPressure_Stokes() -> step000000_vp.(pvd,pvts,vts): CPU time 1.16e-03 (sec) 
  writing pvdfilename analytics_solcx/timeseries_mpoints_std.pvd 
pTatin3d_ModelOutput_MPntStd() -> step000000_mpoints_std.(pvd,pvtu,vtu): CPU time 1.55e-03 (sec) 
[inorms] +1.666666666667e-01,+2.677293357054e-03,+9.530743238464e-02,+9.534502906535e-02,+4.886039454018e-03
  TimeStep control(StkCourant): | current = 1.0000e+30 : trial = 1.3084e+00 [accepted] | ==>> dt used = 1.3084e+00 |
  TimeStep control(StkSurfaceCourant): | current = 1.3084e+00 : trial = 1.0000e+32 | ==>> dt used = 1.3084e+00 |
  timestep[] dt_courant = 1.3084e+00 
[[ModelDestroy_Template]]
[[ModelDestroy_ViscousSinker]]
[[ModelDestroy_Gene3D]]
[[ModelDestroy_Gene3DNueve]]
[[ModelDestroy_Indentor]]
[[ModelDestroy_Rift3D]]
[[ModelDestroy_Rift3D_T]]
[[ModelDestroy_Sierra]]
[[ModelDestroy_AdvDiffExample]]
[[ModelDestroy_Delamination]]
[[ModelDestroy_Riftrh]]
[[ModelDestroy_Rift_oblique3d]]
[[ModelDestroy_MultilayerFolding]]
[[ModelDestroy_SubmarineLavaFlow]]
[[ModelDestroy_ExSubduction]]
[[ModelDestroy_iPLUS]]
[[ModelDestroy_Subduction_Initiation2d]]
[[ModelDestroy_Thermal_Convection2d]]
[[ModelDestroy_AnlVV]]
