From 304c032e65f565589eacfedcca0e7a07c29832d7 Mon Sep 17 00:00:00 2001
From: Douglas Guptill <douglas.guptill@dal.ca>
Date: Tue, 9 Jun 2009 22:18:41 +0000
Subject: [PATCH] input file and stdout from grace

---
 test/input_vM_lvl5.txt  |  623 ++++++
 test/stdout_vM_lvl5.txt | 4229 +++++++++++++++++++++++++++++++++++++++
 2 files changed, 4852 insertions(+)
 create mode 100644 test/input_vM_lvl5.txt
 create mode 100644 test/stdout_vM_lvl5.txt

diff --git a/test/input_vM_lvl5.txt b/test/input_vM_lvl5.txt
new file mode 100644
index 00000000..2cb5cbdd
--- /dev/null
+++ b/test/input_vM_lvl5.txt
@@ -0,0 +1,623 @@
+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+!                                                                              |
+!              ||===\\                                                         | 
+!              ||    \\                                                        |
+!              ||     ||   //==\\   ||  ||   //==||  ||/==\\                   |
+!              ||     ||  ||    ||  ||  ||  ||   ||  ||    ||                  |
+!              ||    //   ||    ||  ||  ||  ||   ||  ||                        |
+!              ||===//     \\==//    \\==\\  \\==\\  ||                        |
+!                                                                              |
+!              Input File                                                      |
+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+
+This file is read by two subroutines in the code:
+- read_controlling_parameters
+- read_input_file
+The indications between [] indicate the type of the read parameter. It can be
+an integer [int], a real*8 [dp], a character chain [char] or a boolean [bool].
+
+
+CONTROLLING PARAMETERS
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]debug is a parameter that switches on/off various prints and outputs
+(the level of printing for error, warning and messages of the solver id%cntl(4)
+is set to the debug value).
+if debug is zero, no debugging 
+if debug is equal to one, this triggers the terminal display of some key parameters 
+if debug is equal to two, same as debug equal to one, and the code produces various 
+vtk files in the DEBUG subdirectories (surfaces, olsf, ...), as well as cross sections 
+for all nonlinear iterations. Careful, memory consuming !
+
+
+      debug = 1
+
+[bool]doDoRuRe is a flag that triggers the production of output files needed to produce
+the DoRuRes. DoRuRe stands for 'Douar Run Report'. 
+
+       doDoRuRe = F
+
+[bool]compute_qp_gram triggers the production of qpgrams for every grid.
+
+      compute_qpgram = F
+
+[bool]compute_reaction_forces toggles on/off the reaction forces computations.
+
+      compute_reaction_forces = F
+
+RESTART
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]irestart is a restart flag; if irestart is not 0, the run will restart from
+an output file given by [char]restartfile and at step.
+
+      irestart = 0
+
+      restartfile=OUT/time_0187.bin
+
+
+TIMESTEPPING
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[dp]dt is the time step length (if dt is negative, courant conidition is used and
+automatic time stepping is turned on)
+
+      dt=1.d-3
+
+[int]nstep is the number of time steps
+
+      nstep = 200
+
+[dp]courant is only used when dt is negative; it determines the size of the time
+step from the maximum value of the velocity field amplitude. The time step s
+the product of courant by the ratio of the smallest leaf size by the maximum
+velocity.
+
+      courant=.5d0
+
+[bool]normaladvect is a flag used to determine which algorithm to use to calculate
+the new geoletry of the normals to the surfaces at the nodes on the surfaces
+if normaladvect is T, the normals are advected using the velocity gradient
+if normaladvect is F, the normals are re-computed from the geometry of the
+surface
+
+      normaladvect = T
+
+GRID ITERATIONS
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]griditer is a flag that allows for nonlinear iterations; when positive, a fixed number
+(griditer) of iterations is permitted; when negative, the number of
+nonlinear iterations is determined by a convergence criterion.
+
+      griditer = 1
+
+[dp]octree_refine_ratio is the threshold value used to determine whether the octree 
+has converged or not. the larger the value, the less stringent the test.
+
+      octree_refine_ratio=1.025d0
+
+
+NONLINEAR ITERATIONS
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+
+[int]nonlinear_iterations is the maximum number of nonlinear iterations (i.e. the 
+iterations on a given constant grid)
+if nonlinear_iterations is positive, it simply is the number of nonlinear
+iterations performed for each grid. When negative it indicates an upper bound
+of nonlinear iterations, but the actual number of nonlinear iterations is
+determined by a convergence criterion (see the 'tol' parameter)
+
+      nonlinear_iterations = 1
+
+[dp]tol is the relative tolerance used to estimate convergence on the computed
+velocity field
+
+      tol=0.01d0
+
+[bool]adaptive_tol is a flag that toggles on/off the evolution of the tol parameter
+with the grid level: when velocity convergence is reached on a grid, the
+latest meaure of the velocity difference between the two last obtained
+solutions is put in tol, so that on the following generated grid, the solution
+reaches at least the same level of convergence. It allows to start with a not
+too stringent value of tol at uniform octree level that evolves with every
+grid, assuming that increasing the level of refinement of the octree allows to
+better capture the solution, hence allowing a tighter convergence. 
+
+      adaptive_tol = F
+
+OCTREES
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]leveluniform_oct is the level of uniform discretization of space; note that a
+level is a power of two used to divide the unit cube
+
+      leveluniform_oct = 5
+
+[int]levelmax_oct is maximum level of octree discretization
+
+      levelmax_oct = 6
+
+[bool]ismooth is a flag to impose an additional level of smoothing after refinement
+for the surfaces and strain rate. It ensures that no leaf is flanked by
+other leaves differing by more than 1 level of refinement
+
+      ismooth = F
+
+[int]noctreemax is the maximum size of any octree used in all computations
+
+      noctreemax=10000000
+
+[dp]refine_ratio is used to determine octree refinement based on a given criterion. 
+All leaves where the criterion is larger than refine_ratio times the maximum 
+of this criterion are refined 
+
+      refine_ratio=-10000
+
+[int]refine_criterion determines which refinement algorithm is to be used.
+Several criteria exist for the refinement of the osolve octree. 1 is the second 
+invariant of the deviatoric strain-rate tensor; 2 is the sum of the squares of the
+diagonal terms of the deviatoric strain-rate tensor; 3 is the second invariant
+of the deviatoric strain rate tensor timses the leaf size. any other value sets
+the criterion to zero and leads to no refinement. 
+
+      refine_criterion=1
+
+[int]initial_refine_level is the initial level at which the refinement of the
+octree will be performed. it has to be smaller than levelmax_oct
+this should be used (different from levelmax_oct) in case the flow is very 
+localized (nonlinear/plastic analysis)
+
+      initial_refine_level = 5
+
+[bool]renumber_nodes is a flag that can toggle on/off the renumbering of the nodes
+by mean of Sloan's algorithm (T/F)
+
+      renumber_nodes = T
+
+PRESSURE
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+[int]smoothing_type is a parameter allows to choose which type of smoothing is to 
+be applied to the pressure field: 0 is none, 1 is center->nodes->center, 2 is the 
+same, but weighted by neighbouring elemental volumes, 3 is regular grid+SPH, and 4 is SPH.  
+
+      smoothing_type = 1
+
+CLOUD
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]npmin and [int]npmax are used to update the 3D volumetric cloud. npmin corresponds
+to the minimum number of particles in any leaf; npmax is the maximum allowable
+number in any leaf
+
+      npmin=8
+      npmax=27
+
+
+FEM + DIVFEM + MUMPS
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]levelcut is the number of levels used to estimate the volume integrals in the
+divFEM algorithmi for cut cells; by testing, it has been estimated that a value
+of 2 is very accurate
+
+      levelcut=2
+
+[int]levelapprox is the number of levels used to estimate the remaining integrals
+using an improved version of Marthijn's clever algorithm, usually 3 is plenty
+
+      levelapprox=3
+
+[dp]penalty is a global penalty parameter used to impose the bad faces or
+incompatible faces linear constraints
+
+      penalty=1.d8
+
+[bool]excl_vol is a parameter that toggles off the assumption that lsf's are built
+on top of one another
+
+      excl_vol = F 
+
+
+TEMPERATURE
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+[dp] ztemp is the height interval between which a linear temperature gradient 
+is set: the temperature is 1 at the bottom, and 0 at ztemp.
+[dp] tempscale
+
+      calculate_temp = F
+
+      ztemp = .1d0
+
+      tempscale=250.d0
+
+MATERIALS
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]nmat is number of materials
+
+      nmat=1
+
+[int]material0 determines what is the material above the 1st surface (the free
+surface) if material0 is 0 then it is the void and the properties "0" are used
+for that part of the model; otherwise the material is one of the material,
+comprised between 1 and nmat
+
+      material0 = 0
+
+[dp]densityi, [dp]viscosityi and [dp]penaltyi are the density, viscosity and
+incompressibility used for material i; there should be nmat sets of material
+properties; there should also be a nil material if material0 has been set to 0
+[dp]expon is the nonlinear viscosity exponent
+[dp]diffusivity is the heat diffusivity
+[dp]heat is the heat production
+
+[char]plasticity_type is the type of plasticity
+- No    -> no plasticity, purely (nonlinear) viscous material
+- vM    -> von Mises yield criterion
+        -> [dp]plasticity_1st_param is the yield value
+- DPI   -> Drucker-Prager type of yield criterion
+        -> the yield locus passes through the outer apices of the Mohr-Coulomb hexagon
+        -> [dp]plasticity_1st_param is the angle phi
+        -> [dp]plasticity_2nd_param is the cohesion c
+- DPII  -> Drucker-Prager type of yield criterion
+        -> the yield locus passes through the inner apices of the Mohr-Coulomb hexagon
+        -> [dp]plasticity_1st_param is the angle phi
+        -> [dp]plasticity_2nd_param is the cohesion c
+- DPIII -> Drucker-Prager type of yield criterion
+        -> [dp]plasticity_1st_param is alpha
+        -> [dp]plasticity_2nd_param is k
+- MC    -> Mohr-Coulomb type of yield criterion
+        -> [dp]plasticity_1st_param is the angle phi
+        -> [dp]plasticity_2nd_param is the cohesion c
+
+      density0              = 0.d0
+      viscosity0            = 1.d-5
+      penalty0              = 1.d8
+      expon0                = 1.d0
+      diffusivity0          = 1.d0
+      heat0                 = 0.d0
+      activationenergy0     = 0.d0
+      plasticity_type0      = MC
+      plasticity_1st_param0 = 15.d0
+      plasticity_2nd_param0 = 1.72d-3
+      plasticity_3rd_param0 = 0.5d0
+      plasticity_4th_param0 = 1.5d0
+      plasticity_5th_param0 = 5.d0
+
+      density1              =-1.d0
+      viscosity1            = 1.3088d3
+      penalty1              = 1.d8
+      expon1                = 1.d0
+      diffusivity1          = 0.d0
+      heat1                 = 0.d0
+      activationenergy1     = 0.d0
+      plasticity_type1      = vM
+      plasticity_1st_param1 = 1.73362d-2
+      plasticity_2nd_param1 = 1.73362d-3
+      plasticity_3rd_param1 = 0.25d0
+      plasticity_4th_param1 = 1.25d0
+
+[dp]viscositymin and viscositymax are bounds on the viscosity (if negative bound is not imposed)
+These bounds are introduced to prevent the viscosity to reach unrealistic values, especially when
+using non-linear (power-law or brittle) rheologies
+
+      viscositymin=1.d-8
+      viscositymax=1.d8
+
+SURFACES
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]ns is number of surfaces to track
+
+      ns=1
+
+for each surface, one needs to define a levelt, itype, material, surface_type, activation_time 
+and surface_params .
+- [int]levelt is the inital level for the particles on the surface; to be accurate
+and avoid wholes in the surface during definition of the lsf, one should use
+levelt eq to levelmax_oct+1 for all surfaces as a minimum value; 
+- [int]itype should be 1 for foldable surfaces or 0 for nonfoldable surfaces; 
+- [int]material is the material type refering to the table of material available (max nmat); 
+- [dp]activation_time is the time the surface becomes active (before that time, it
+  is glued to the 0th surface). this parameter is useful when defining stratigraphic horizons; 
+  default is -1, ie the surface is not glued to the free surface
+- [int]surface_type is comprised between 1 and 8. 
+  1 corresponds to  a flat surface, 
+    -> surface_param_01 is the z level
+  2 to a rectangular emboss, 
+    -> surface_param_01 is the z level
+    -> surface_param_02 and 03 are x1,x2 
+    -> surface_param_04 and 05 are y1,y2 
+    -> surface_param_06 is the thickness 
+  3 to a convex spherical emboss, 
+    -> surface_param_01 is the z level
+    -> surface_param_02 and 03 are x0,y0 
+    -> surface_param_04 is the radius 
+  4 to concave spherical emboss,
+    -> surface_param_01 is the z level
+    -> surface_param_02 and 03 are x0,y0 
+    -> surface_param_04 is the radius 
+  5 to a double rectangular emboss, 
+    -> surface_param_01 is the z level
+    -> surface_param_02 and 03 are x1,x2 
+    -> surface_param_04 and 05 are x3,x4 
+    -> surface_param_06 and 07 are y1,y2 
+    -> surface_param_08 and 09 are y3,y4 
+    -> surface_param_10 is the thickness 
+  6 to a sinus, 
+    -> surface_param_01 is the z level
+    -> surface_param_02 is the wavelength 
+    -> surface_param_03 is the amplitude 
+  7 to a noisy surface, 
+    -> surface_param_01 is the z level
+    -> surface_param_02 is the noise amplitude
+  8 to a double sinus.
+    -> surface_param_01 is the z level
+    -> surface_param_02 is the x-wavelength 
+    -> surface_param_03 is the x-amplitude 
+    -> surface_param_04 is the y-wavelength 
+    -> surface_param_05 is the y-amplitude 
+- [int]leveloct is the level at which the octree will be refined in the vicinity of the surface.
+
+      levelt1             = 6
+      itype1              = 0
+      surface_type_1      = 1
+      rand1               = T
+      surface_param_01_1  = .1d0
+      material1           = 1
+      activation_time_1   = -1. 
+      leveloct1           = 6
+      stretch1            = 1.5d0
+      anglemax1           = 180.d0
+      criterion1          = 1
+      anglemaxoctree1     = 180.d0
+      spread_surface_points1 = 1
+
+[int]niter_move is the number of iterations used to update particle positions
+using an implicit, mid-point algorithm (default is 10)
+
+      niter_move = 10
+
+[dp]stretch is the maximum allowed increase in linear length between two initially
+adjacent particles on any surface; when this stretch is achieved, a new
+particle is inserted on the surface, half-way along the stretched edge
+
+      stretch = 1.5d0
+
+[dp]anglemax is the maximum allowed angle between two adjacent normals
+when the angle is reached a new point is inserted beteen the two points to
+reduce the angle between the two normals
+
+      anglemax=180.d0
+
+
+[int]criterion is criterion used to define the octree in the vicinity of the
+sufaces; criterion 1 corresponds to imposing that all leaves that are cut by any of the
+surfaces must be at level levelmax_oct; criterion 2 corresponds
+to imposing that discretization is proportional to the curvature of the
+surface; curvature is calculated from the local divergence of the normals.
+criterion 3 corresponds to imposing that all leaves that contain at
+least one particle of any surface is at levelmax_oct; 
+
+      criterion = 2
+
+[dp]anglemaxoctree is only defined for criterion 2; t is the maximum allowable angle
+between two adjacent normals; if the angle is greater than anglemaxoctree, the local
+octree leaves are forced to be at level levelmax_oct; otherwise they are
+proportionally larger (smaller levels) (default is 10)
+
+     anglemaxoctree = 10.d0
+
+
+REFINEMENT IN BOXES
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]nboxes is the number of boxes in which the user imposes a set level of
+discretization
+
+      nboxes = 0
+
+for each box we need to specify the two end corners and the level
+the syntax is [dp]boxnx0, [dp]boxnx1, [dp]boxny0, [dp]boxny1, [dp]boxnz0, [dp]boxnz1
+and boxnlevel where n is the box number
+      box1x0=0.d0
+      box1x1=1.d0
+      box1y0=0.d0
+      box1y1=1.d0
+      box1z0=0.0d0
+      box1z1=0.2d0
+      box1level=7
+
+      box2x0=0.1d0
+      box2x1=0.4d0
+      box2y0=0.0d0
+      box2y1=0.5d0
+      box2z0=0.0d0
+      box2z1=0.2d0
+      box2level=7
+
+      box3x0=0.6d0
+      box3x1=0.9d0
+      box3y0=0.5d0
+      box3y1=1.0d0
+      box3z0=0.0d0
+      box3z1=0.2d0
+      box3level=7
+
+REFINEMENT ON CUBE FACES
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+[bool]ref_on_faces toggles on/off the user imposed refinement on faces.
+for each of the six faces of the cube, on defines the level at which the
+desired area is to be refined. this area is given by bottom, top, left and
+right coordinates that are comprised between 0 and 1
+faces 1,2,3,4,5,6 respectively correspond to x=0,x=1,y=0,y=1,z=0,z=1
+
+      ref_on_faces = F
+
+      level_face1=5
+      b1=.02
+      t1=.51
+      l1=.11
+      r1=.81
+      level_face2=5
+      b2=.45
+      t2=.55
+      l2=.46
+      r2=.56
+      level_face3=5
+      b3=.0
+      t3=.2
+      l3=.3
+      r3=.7
+      level_face4=5
+      b4=.4
+      t4=.5
+      l4=.14
+      r4=.4
+      level_face5=5
+      b5=0.
+      t5=1.
+      l5=0.
+      r5=0.75
+      level_face6=5
+      b6=.26
+      t6=.56
+      l6=.16
+      r6=.86
+
+
+EROSION
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[bool]erosion is a flag that toggles on/off the erosion. 
+
+      erosion = T
+
+if erosion is on, one also needs to set the erosion level/height, [dp]zerosion
+Note that this is a first attempt at erosion; in future versions, DOUAR
+should be easily linked to a surface processes model like CASCADE
+
+      zerosion=.1d0
+
+If erosion is on, one also needs to define a length scale and a velocity scale to
+properly translate the uplift rate produced by DOUAR into something that is adequate
+for CASCADE; 
+[dp]length_scale is the scale of the unit model in km 
+[dp]velocity_scale is the velocity scale in km/Myr. 
+if length_scale is negative, erosion is assumed to be perfect (no call to CASCADE is needed)
+
+      length_scale=200.d0
+      velocity_scale=10.d0
+
+One then needs to define the erosion constants in CASCADE
+[dp]fluvial_erosion is the fluvial erosion constant in 1/m^2
+[dp]diffusion_erosion is the diffusion erosion constant in m^2/yr
+(4d-2 32d-2)
+
+      fluvial_erosion=2.d-2
+      diffusion_erosion=16.d-2
+
+One finally needs to specify the boundary conditions for CASCADE
+if [int]baselevelx0 is set to 1 then the boundary at x=0 is set at baselevel (water and sediment exit)
+if [int]baselevelx1 is set to 1 then the boundary at x=1 is set at baselevel (water and sediment exit)
+if [int]baselevely0 is set to 1 then the boundary at y=0 is set at baselevel (water and sediment exit)
+if [int]baselevely1 is set to 1 then the boundary at y=1 is set at baselevel (water and sediment exit)
+
+      baselevelx0=1
+      baselevelx1=1
+      baselevely0=0
+      baselevely1=0
+
+MATRIX VISUALISATION
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[bool]visualise_matrix allows the user to turn on the visual representation of the 
+matrices used in the code. Be careful, the generated postscript files are huge !
+
+      visualise_matrix = F
+
+CROSS SECTIONS
+_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
+
+[int]nsections is the number of cross-sections to be output. [dp]scale is the scale used to produce the postscript file.
+[int]xyz takes values 1, 2 or 3, and corresponds to planes defined by x=constant, y=constant, and z=constant resp.
+The [bool]flags are self explanatory. [char]colormap is the chosen colormap (jet or hot). [int]ncolours is the number of
+colours used to produce the plot. 
+
+nsections = 0
+
+xyz_1          = 1 
+slice_1        = 0.5001
+flag_press_1   = T
+flag_spress_1  = T 
+flag_e2d_1     = T
+flag_e3d_1     = F
+flag_lode_1    = F 
+flag_crit_1    = F 
+flag_grid_1    = T 
+flag_mu_1      = F   
+flag_u_1       = F
+flag_v_1       = F
+flag_w_1       = F
+flag_q_1       = F 
+flag_uvw_1     = F 
+flag_lsf_1     = F 
+flag_vfield_1  = F 
+flag_colour_1  = T
+flag_plastic_1 = F
+flag_velvect_1 = T
+scale_1        = 500.
+colormap_1     = jet
+ncolours_1     = 256
+
+xyz_2          = 2
+slice_2        = 0.901
+flag_press_2   = T
+flag_e2d_2     = T
+flag_e3d_2     = F 
+flag_lode_2    = F 
+flag_crit_2    = F 
+flag_grid_2    = F 
+flag_mu_2      = F   
+flag_u_2       = F
+flag_v_2       = F
+flag_w_2       = F
+flag_q_2       = F 
+flag_uvw_2     = F 
+flag_lsf_2     = F 
+flag_vfield_2  = F 
+flag_colour_2  = T
+flag_plastic_2 = F
+flag_velvect_2 = F
+scale_2        = 800.
+colormap_2     = jet 
+ncolours_2     = 256
+
+
+xyz_3          = 3
+slice_3        = 0.0010
+flag_press_3   = F 
+flag_e2d_3     = F
+flag_e3d_3     = F
+flag_lode_3    = F 
+flag_crit_3    = F 
+flag_grid_3    = F 
+flag_mu_3      = F   
+flag_u_3       = F
+flag_v_3       = F
+flag_w_3       = F
+flag_q_3       = F
+flag_uvw_3     = F 
+flag_lsf_3     = F 
+flag_vfield_3  = F 
+flag_colour_3  = T
+flag_plastic_3 = F
+flag_velvect_3 = F
+scale_3        = 400.
+colormap_3     = jet
+ncolours_3     = 256
+
diff --git a/test/stdout_vM_lvl5.txt b/test/stdout_vM_lvl5.txt
new file mode 100644
index 00000000..dbacc197
--- /dev/null
+++ b/test/stdout_vM_lvl5.txt
@@ -0,0 +1,4229 @@
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+ program called with no argument
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Computations  Total time   Step time   Increment
+                       Reading Input      0.0009      0.0009      0.0009
+                       Problem Setup      6.9290      6.9290      6.9280
+                     define surfaces      6.9290      6.9290      0.0000
+                          surface 01      6.9290      6.9290      0.0000
+                        define cloud      6.9421      6.9421      0.0130
+              define velocity octree      6.9422      6.9422      0.0001
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       1      7.0615
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces      7.0624      0.0009      0.0009
+                          surface 01      7.0624      0.0010      0.0000
+                                                                        S. 1:     4356points
+                      refine surface      7.0624      0.0010      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay      7.0673      0.0058      0.0048
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1      7.0714
+ ----------------------------------------------------------------------- 
+                 Create octree solve      7.0714      0.0000      0.0000
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov      7.0872      0.0158      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion      7.0963      0.0250      0.0092
+                                                                        max(crit)        =   0.000000E+00
+                                                                        max(crit,w_l_i_f)=      -Infinity
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement      7.0967      0.0253      0.0004
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon      7.1677      0.0963      0.0710
+             Imbed surface in osolve      7.2725      0.2011      0.1049
+                embedding surface  1      7.2726      0.2012      0.0001
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability      8.9183      1.8469      1.6457
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=        0 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=        0 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  F
+                     Change 3D cloud      8.9213      1.8499      0.0030
+        Interpolate velo onto osolve     10.3643      3.2929      1.4430
+                                                                        ov%pressure:      0.0000E+00 0.0000E+00
+                                                                        osolve%pressure:  0.0000E+00 0.0000E+00
+                      Find bad faces     10.4800      3.4087      0.1157
+                           Find void     10.6134      3.5420      0.1333
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions     10.6217      3.5503      0.0083
+                         wsmp setup1     10.6248      3.5535      0.0032
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2     10.8553      3.7839      0.2304
+ -----------------------------------
+ start of non-linear iteratio      1     10.8648
+ -----------------------------------
+                        build system     10.8649      0.0001      0.0001
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   1.22294E-04  1.30880E+03
+                                                                        viscosity capped in      7168
+                          wsmp solve     13.2342      2.3694      2.3693
+  ************************************************************************
+  
+       CONTAINS Runtime Modules of Watson Sparse Matrix Package V 9.1.7.
+        (C) Copyright IBM Corporation 1997, 2008.  All Rights Reserved.
+
+
+WARNING**: The environment variable MALLOC_TRIM_THRESHOLD_
+is not set. This may adversely affect the performance of
+WSMP. Please refer to WSMP documentation for more details.
+
+WARNING**: The environment variable MALLOC_MMAP_MAX_ is
+not set. This may adversely affect the performance of WSMP
+Please refer to WSMP documentation for more details.
+  ************************************************************************
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.4264500 s
+                                                                        wsmp: ordering time:               0.6810880 s
+                                                                        wsmp: symbolic fact. time:         0.1536419 s
+                                                                        wsmp: Cholesky fact. time:         0.4089160 s
+                                                                        wsmp: Factorisation            11735.2406520 Mflops
+                                                                        wsmp: back substitution time:      0.0150051 s
+                                                                        wsmp: from b to rhs vector:        0.0019119 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.6873171 s
+                                                                        =================================
+                                                                        u : -0.15175E+00  0.13100E+01
+                                                                        v : -0.12328E+00  0.65164E+01
+                                                                        w : -0.55999E-01  0.16432E+01
+                                                                        =================================
+                 Calculate pressures     14.9269      4.0621      1.6927
+                                                                        raw    pressures : -0.90782E+05  0.75862E+05
+                 Smoothing pressures     14.9875      4.1227      0.0606
+                                                                        smooth pressures : -0.13110E+05  0.48234E+04
+                do leaf measurements     14.9935      4.1286      0.0060
+       compute convergence criterion     15.0151      4.1503      0.0217
+                                                                        velocity_diff_norm = 3.2582059 (tol = 0.0100000)
+                                                                        increase_current_level -> F
+                                                                        current_level -> 5
+                    slicing the cube     15.0168      4.1520      0.0017
+Compute isostasy and adjust vertical     15.0169      4.1520      0.0001
+                  Compute divergence     15.0169      4.1520      0.0000
+                        wsmp_cleanup     15.0258      4.1610      0.0089
+              Reset surface geometry     15.0278      4.1629      0.0020
+ -----------------------------------------------------------------------
+           Temperature calculations      15.0284      4.1636      0.0006
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces     15.1829      4.3181      0.1545
+                   Advect surface  1     15.1830      4.3181      0.0000
+                                                                        S. 1:   11 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1     15.2303      4.3655      0.0473
+                 Update cloud fields     46.0719     35.2071     30.8416
+                    Advect the cloud     46.4870     35.6221      0.4150
+                        Write output     48.8053     37.9405      2.3184
+                    End of time step     49.8742     39.0093      1.0688
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       2     49.8742
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces     49.8742      0.0000      0.0000
+                          surface 01     49.8742      0.0000      0.0000
+                                                                        S. 1:     4356points
+                      refine surface     49.8743      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay     49.8788      0.0046      0.0046
+                                                                        S. 1:   36 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1     49.8835
+ ----------------------------------------------------------------------- 
+                 Create octree solve     49.8836      0.0000      0.0000
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov     49.8994      0.0158      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion     49.9100      0.0265      0.0107
+                                                                        max(crit)        =   9.026395E+01
+                                                                        max(crit,w_l_i_f)=   9.026395E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement     49.9147      0.0312      0.0046
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon     49.9857      0.1022      0.0710
+             Imbed surface in osolve     50.0881      0.2046      0.1023
+                embedding surface  1     50.0881      0.2046      0.0001
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability     51.7601      1.8766      1.6720
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud     51.7627      1.8792      0.0026
+        Interpolate velo onto osolve     53.0325      3.1490      1.2698
+                                                                        ov%pressure:     -0.1850E+05 0.1333E+05
+                                                                        osolve%pressure: -0.1311E+05 0.4823E+04
+                      Find bad faces     53.2420      3.3584      0.2094
+                           Find void     53.3753      3.4918      0.1333
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions     53.3832      3.4997      0.0080
+                         wsmp setup1     53.3865      3.5030      0.0032
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2     53.6183      3.7348      0.2318
+ -----------------------------------
+ start of non-linear iteratio      1     53.6285
+ -----------------------------------
+                        build system     53.6286      0.0000      0.0000
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   9.63581E-05  6.27944E-01
+                                                                        viscosity capped in      7168
+                          wsmp solve     56.0127      2.3842      2.3842
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0730410 s
+                                                                        wsmp: ordering time:               0.6769719 s
+                                                                        wsmp: symbolic fact. time:         0.1534011 s
+                                                                        wsmp: Cholesky fact. time:         0.4040279 s
+                                                                        wsmp: Factorisation            11877.2173027 Mflops
+                                                                        wsmp: back substitution time:      0.0147259 s
+                                                                        wsmp: from b to rhs vector:        0.0019500 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3243999 s
+                                                                        =================================
+                                                                        u : -0.18242E-04  0.10942E+01
+                                                                        v : -0.11415E+00  0.15268E+00
+                                                                        w : -0.52912E-01  0.76303E+00
+                                                                        =================================
+                 Calculate pressures     57.3426      3.7141      1.3298
+                                                                        raw    pressures : -0.38678E+00  0.14668E+00
+                 Smoothing pressures     57.4039      3.7754      0.0613
+                                                                        smooth pressures : -0.11488E+00  0.90245E-03
+                do leaf measurements     57.4099      3.7814      0.0060
+       compute convergence criterion     57.4323      3.8037      0.0224
+                                                                        velocity_diff_norm = 3.1911487 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 5
+                    slicing the cube     57.4339      3.8053      0.0016
+Compute isostasy and adjust vertical     57.4339      3.8054      0.0001
+                  Compute divergence     57.4339      3.8054      0.0000
+                        wsmp_cleanup     57.4429      3.8143      0.0089
+              Reset surface geometry     57.4450      3.8165      0.0022
+ -----------------------------------------------------------------------
+           Temperature calculations      57.4458      3.8173      0.0008
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces     57.6060      3.9775      0.1602
+                   Advect surface  1     57.6061      3.9775      0.0000
+                                                                        S. 1:    3 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1     57.6504      4.0219      0.0443
+                 Update cloud fields     83.4390     29.8105     25.7886
+                    Advect the cloud     83.8666     30.2381      0.4276
+                        Write output     86.2290     32.6004      2.3624
+                    End of time step     87.3195     33.6910      1.0906
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       3     87.3196
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces     87.3196      0.0000      0.0000
+                          surface 01     87.3197      0.0000      0.0000
+                                                                        S. 1:     4356points
+                      refine surface     87.3197      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay     87.3245      0.0049      0.0048
+                                                                        S. 1:    5 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1     87.3295
+ ----------------------------------------------------------------------- 
+                 Create octree solve     87.3296      0.0001      0.0001
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov     87.3454      0.0159      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion     87.3565      0.0269      0.0111
+                                                                        max(crit)        =   4.351508E+01
+                                                                        max(crit,w_l_i_f)=   1.302893E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement     87.3610      0.0315      0.0046
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon     87.4320      0.1024      0.0710
+             Imbed surface in osolve     87.5325      0.2030      0.1005
+                embedding surface  1     87.5326      0.2030      0.0001
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability     89.2035      1.8740      1.6710
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud     89.2065      1.8770      0.0030
+        Interpolate velo onto osolve     90.4735      3.1439      1.2670
+                                                                        ov%pressure:     -0.1343E+00 0.6504E-02
+                                                                        osolve%pressure: -0.1149E+00 0.9024E-03
+                      Find bad faces     90.6757      3.3461      0.2022
+                           Find void     90.8090      3.4795      0.1333
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions     90.8169      3.4874      0.0079
+                         wsmp setup1     90.8203      3.4908      0.0034
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2     91.0515      3.7220      0.2312
+ -----------------------------------
+ start of non-linear iteratio      1     91.0617
+ -----------------------------------
+                        build system     91.0617      0.0000      0.0000
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   6.43125E-04  9.59614E+00
+                                                                        viscosity capped in      7168
+                          wsmp solve     93.4024      2.3407      2.3407
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0795159 s
+                                                                        wsmp: ordering time:               0.6805451 s
+                                                                        wsmp: symbolic fact. time:         0.1520360 s
+                                                                        wsmp: Cholesky fact. time:         0.4030600 s
+                                                                        wsmp: Factorisation            11905.7413517 Mflops
+                                                                        wsmp: back substitution time:      0.0148430 s
+                                                                        wsmp: from b to rhs vector:        0.0020211 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3323109 s
+                                                                        =================================
+                                                                        u : -0.63225E-06  0.10127E+01
+                                                                        v : -0.13552E+00  0.16205E+00
+                                                                        w : -0.26949E-01  0.58584E+00
+                                                                        =================================
+                 Calculate pressures     94.7403      3.6786      1.3379
+                                                                        raw    pressures : -0.29642E+00  0.79373E-01
+                 Smoothing pressures     94.8019      3.7402      0.0616
+                                                                        smooth pressures : -0.10843E+00  0.00000E+00
+                do leaf measurements     94.8079      3.7463      0.0060
+       compute convergence criterion     94.8301      3.7684      0.0222
+                                                                        velocity_diff_norm = 0.2976647 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 5
+                    slicing the cube     94.8318      3.7701      0.0017
+Compute isostasy and adjust vertical     94.8318      3.7702      0.0001
+                  Compute divergence     94.8319      3.7702      0.0000
+                        wsmp_cleanup     94.8409      3.7792      0.0090
+              Reset surface geometry     94.8431      3.7814      0.0023
+ -----------------------------------------------------------------------
+           Temperature calculations      94.8439      3.7822      0.0008
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces     95.0041      3.9424      0.1602
+                   Advect surface  1     95.0041      3.9424      0.0000
+                                                                        S. 1:    3 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1     95.0485      3.9868      0.0444
+                 Update cloud fields    115.7792     24.7175     20.7307
+                    Advect the cloud    116.2104     25.1487      0.4312
+                        Write output    118.5814     27.5197      2.3710
+                    End of time step    119.6718     28.6101      1.0904
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       4    119.6719
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    119.6719      0.0000      0.0000
+                          surface 01    119.6719      0.0000      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    119.6719      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    119.6768      0.0050      0.0049
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    119.6818
+ ----------------------------------------------------------------------- 
+                 Create octree solve    119.6818      0.0001      0.0001
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    119.6978      0.0160      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    119.7089      0.0271      0.0111
+                                                                        max(crit)        =   3.303445E+01
+                                                                        max(crit,w_l_i_f)=   1.362159E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement    119.7137      0.0319      0.0048
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon    119.7854      0.1036      0.0717
+             Imbed surface in osolve    119.8895      0.2078      0.1041
+                embedding surface  1    119.8896      0.2078      0.0001
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability    121.5961      1.9143      1.7065
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    121.5992      1.9175      0.0031
+        Interpolate velo onto osolve    122.9167      3.2349      1.3174
+                                                                        ov%pressure:     -0.1247E+00 0.8834E-04
+                                                                        osolve%pressure: -0.1084E+00 0.0000E+00
+                      Find bad faces    123.0414      3.3596      0.1247
+                           Find void    123.1762      3.4944      0.1348
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions    123.1847      3.5029      0.0085
+                         wsmp setup1    123.1883      3.5066      0.0036
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2    123.4230      3.7412      0.2346
+ -----------------------------------
+ start of non-linear iteratio      1    123.4342
+ -----------------------------------
+                        build system    123.4342      0.0000      0.0000
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   5.59495E-04  1.02031E+02
+                                                                        viscosity capped in      7168
+                          wsmp solve    125.7785      2.3443      2.3442
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0782461 s
+                                                                        wsmp: ordering time:               0.6756001 s
+                                                                        wsmp: symbolic fact. time:         0.1535821 s
+                                                                        wsmp: Cholesky fact. time:         0.4049351 s
+                                                                        wsmp: Factorisation            11850.6085809 Mflops
+                                                                        wsmp: back substitution time:      0.0147769 s
+                                                                        wsmp: from b to rhs vector:        0.0019999 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3294241 s
+                                                                        =================================
+                                                                        u : -0.63763E-07  0.10012E+01
+                                                                        v : -0.14477E+00  0.16266E+00
+                                                                        w : -0.95903E-02  0.60784E+00
+                                                                        =================================
+                 Calculate pressures    127.1134      3.6792      1.3349
+                                                                        raw    pressures : -0.28159E+00  0.63161E-01
+                 Smoothing pressures    127.1751      3.7409      0.0618
+                                                                        smooth pressures : -0.10688E+00  0.00000E+00
+                do leaf measurements    127.1819      3.7477      0.0068
+       compute convergence criterion    127.2047      3.7705      0.0228
+                                                                        velocity_diff_norm = 0.0643483 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 5
+                    slicing the cube    127.2064      3.7722      0.0017
+Compute isostasy and adjust vertical    127.2065      3.7723      0.0001
+                  Compute divergence    127.2065      3.7723      0.0000
+                        wsmp_cleanup    127.2155      3.7814      0.0091
+              Reset surface geometry    127.2177      3.7835      0.0022
+ -----------------------------------------------------------------------
+           Temperature calculations     127.2185      3.7843      0.0008
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    127.3790      3.9448      0.1605
+                   Advect surface  1    127.3790      3.9448      0.0000
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    127.4236      3.9894      0.0445
+                 Update cloud fields    145.6892     22.2550     18.2657
+                    Advect the cloud    146.1248     22.6906      0.4355
+                        Write output    148.5280     25.0938      2.4032
+                    End of time step    149.6117     26.1775      1.0837
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       5    149.6117
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    149.6117      0.0000      0.0000
+                          surface 01    149.6117      0.0000      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    149.6118      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    149.6178      0.0061      0.0061
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    149.6229
+ ----------------------------------------------------------------------- 
+                 Create octree solve    149.6230      0.0001      0.0001
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    149.6388      0.0159      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    149.6500      0.0271      0.0112
+                                                                        max(crit)        =   3.209192E+01
+                                                                        max(crit,w_l_i_f)=   1.404212E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement    149.6548      0.0318      0.0048
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon    149.7257      0.1028      0.0709
+             Imbed surface in osolve    149.8298      0.2069      0.1041
+                embedding surface  1    149.8299      0.2069      0.0001
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability    151.5003      1.8774      1.6704
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    151.5033      1.8804      0.0030
+        Interpolate velo onto osolve    152.8151      3.1921      1.3118
+                                                                        ov%pressure:     -0.1219E+00 0.4505E-04
+                                                                        osolve%pressure: -0.1069E+00 0.0000E+00
+                      Find bad faces    152.9984      3.3755      0.1834
+                           Find void    153.1318      3.5089      0.1334
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions    153.1398      3.5169      0.0080
+                         wsmp setup1    153.1433      3.5203      0.0034
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2    153.3752      3.7522      0.2319
+ -----------------------------------
+ start of non-linear iteratio      1    153.3861
+ -----------------------------------
+                        build system    153.3861      0.0000      0.0000
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   4.71723E-04  7.81910E+02
+                                                                        viscosity capped in      7168
+                          wsmp solve    155.7868      2.4007      2.4007
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0740390 s
+                                                                        wsmp: ordering time:               0.6782150 s
+                                                                        wsmp: symbolic fact. time:         0.1556020 s
+                                                                        wsmp: Cholesky fact. time:         0.4042411 s
+                                                                        wsmp: Factorisation            11870.9547432 Mflops
+                                                                        wsmp: back substitution time:      0.0148351 s
+                                                                        wsmp: from b to rhs vector:        0.0019600 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3291750 s
+                                                                        =================================
+                                                                        u : -0.12700E-07  0.10000E+01
+                                                                        v : -0.14855E+00  0.16059E+00
+                                                                        w : -0.40837E-02  0.60952E+00
+                                                                        =================================
+                 Calculate pressures    157.1214      3.7354      1.3347
+                                                                        raw    pressures : -0.28456E+00  0.64289E-01
+                 Smoothing pressures    157.1845      3.7985      0.0631
+                                                                        smooth pressures : -0.10632E+00  0.00000E+00
+                do leaf measurements    157.1906      3.8046      0.0061
+       compute convergence criterion    157.2135      3.8274      0.0229
+                                                                        velocity_diff_norm = 0.0252260 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 5
+                    slicing the cube    157.2152      3.8291      0.0017
+Compute isostasy and adjust vertical    157.2153      3.8292      0.0001
+                  Compute divergence    157.2153      3.8292      0.0000
+                        wsmp_cleanup    157.2244      3.8384      0.0092
+              Reset surface geometry    157.2269      3.8408      0.0024
+ -----------------------------------------------------------------------
+           Temperature calculations     157.2277      3.8417      0.0008
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    157.3890      4.0030      0.1613
+                   Advect surface  1    157.3890      4.0030      0.0000
+                                                                        S. 1:    3 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    157.4338      4.0478      0.0448
+                 Update cloud fields    174.2075     20.8215     16.7737
+                    Advect the cloud    174.6450     21.2589      0.4375
+                        Write output    177.0699     23.6839      2.4249
+                    End of time step    178.1885     24.8024      1.1186
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       6    178.1887
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    178.1888      0.0001      0.0001
+                          surface 01    178.1888      0.0001      0.0001
+                                                                        S. 1:     4356points
+                      refine surface    178.1889      0.0002      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    178.1948      0.0061      0.0059
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    178.1996
+ ----------------------------------------------------------------------- 
+                 Create octree solve    178.1998      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    178.2156      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    178.2268      0.0271      0.0111
+                                                                        max(crit)        =   3.199963E+01
+                                                                        max(crit,w_l_i_f)=   1.437570E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement    178.2316      0.0320      0.0048
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon    178.3027      0.1030      0.0711
+             Imbed surface in osolve    178.4074      0.2078      0.1047
+                embedding surface  1    178.4076      0.2080      0.0002
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability    180.0790      1.8794      1.6714
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    180.0822      1.8825      0.0032
+        Interpolate velo onto osolve    181.4342      3.2346      1.3520
+                                                                        ov%pressure:     -0.1205E+00 0.1892E-02
+                                                                        osolve%pressure: -0.1063E+00 0.0000E+00
+                      Find bad faces    181.6239      3.4243      0.1897
+                           Find void    181.7573      3.5577      0.1335
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions    181.7656      3.5660      0.0083
+                         wsmp setup1    181.7699      3.5703      0.0043
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2    182.0026      3.8029      0.2326
+ -----------------------------------
+ start of non-linear iteratio      1    182.0134
+ -----------------------------------
+                        build system    182.0137      0.0002      0.0002
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   4.19760E-04  1.30880E+03
+                                                                        viscosity capped in      7168
+                          wsmp solve    184.4364      2.4230      2.4228
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0818019 s
+                                                                        wsmp: ordering time:               0.6810110 s
+                                                                        wsmp: symbolic fact. time:         0.1542459 s
+                                                                        wsmp: Cholesky fact. time:         0.4080091 s
+                                                                        wsmp: Factorisation            11761.3263706 Mflops
+                                                                        wsmp: back substitution time:      0.0148439 s
+                                                                        wsmp: from b to rhs vector:        0.0019732 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3421690 s
+                                                                        =================================
+                                                                        u : -0.31677E-12  0.10000E+01
+                                                                        v : -0.14996E+00  0.15814E+00
+                                                                        w : -0.17733E-02  0.60421E+00
+                                                                        =================================
+                 Calculate pressures    185.7838      3.7704      1.3474
+                                                                        raw    pressures : -0.29174E+00  0.70691E-01
+                 Smoothing pressures    185.8466      3.8331      0.0627
+                                                                        smooth pressures : -0.10626E+00  0.00000E+00
+                do leaf measurements    185.8529      3.8395      0.0064
+       compute convergence criterion    185.8755      3.8621      0.0225
+                                                                        velocity_diff_norm = 0.0170324 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 5
+                    slicing the cube    185.8773      3.8639      0.0018
+Compute isostasy and adjust vertical    185.8775      3.8641      0.0002
+                  Compute divergence    185.8775      3.8641      0.0000
+                        wsmp_cleanup    185.8866      3.8732      0.0091
+              Reset surface geometry    185.8892      3.8758      0.0026
+ -----------------------------------------------------------------------
+           Temperature calculations     185.8903      3.8768      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    186.0502      4.0368      0.1599
+                   Advect surface  1    186.0503      4.0369      0.0002
+                                                                        S. 1:    3 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    186.0947      4.0812      0.0443
+                 Update cloud fields    201.9391     19.9256     15.8444
+                    Advect the cloud    202.3814     20.3680      0.4423
+                        Write output    204.7923     22.7789      2.4110
+                    End of time step    205.9055     23.8921      1.1132
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       7    205.9057
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    205.9057      0.0001      0.0001
+                          surface 01    205.9058      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    205.9058      0.0002      0.0001
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    205.9122      0.0065      0.0064
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    205.9164
+ ----------------------------------------------------------------------- 
+                 Create octree solve    205.9166      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    205.9324      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    205.9435      0.0271      0.0111
+                                                                        max(crit)        =   3.199995E+01
+                                                                        max(crit,w_l_i_f)=   1.464236E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement    205.9482      0.0319      0.0047
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon    206.0193      0.1029      0.0711
+             Imbed surface in osolve    206.1224      0.2060      0.1031
+                embedding surface  1    206.1226      0.2062      0.0002
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability    207.7870      1.8706      1.6644
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    207.7902      1.8738      0.0031
+        Interpolate velo onto osolve    209.0930      3.1766      1.3028
+                                                                        ov%pressure:     -0.1200E+00 0.1606E-02
+                                                                        osolve%pressure: -0.1063E+00 0.0000E+00
+                      Find bad faces    209.3216      3.4053      0.2287
+                           Find void    209.4552      3.5388      0.1335
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions    209.4630      3.5467      0.0079
+                         wsmp setup1    209.4665      3.5501      0.0034
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2    209.6967      3.7803      0.2302
+ -----------------------------------
+ start of non-linear iteratio      1    209.7063
+ -----------------------------------
+                        build system    209.7065      0.0002      0.0002
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   3.80269E-04  1.30880E+03
+                                                                        viscosity capped in      7168
+                          wsmp solve    212.1489      2.4426      2.4424
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0692410 s
+                                                                        wsmp: ordering time:               0.6784201 s
+                                                                        wsmp: symbolic fact. time:         0.1543071 s
+                                                                        wsmp: Cholesky fact. time:         0.4050360 s
+                                                                        wsmp: Factorisation            11847.6578642 Mflops
+                                                                        wsmp: back substitution time:      0.0148780 s
+                                                                        wsmp: from b to rhs vector:        0.0020611 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3242190 s
+                                                                        =================================
+                                                                        u : -0.55090E-12  0.10000E+01
+                                                                        v : -0.15231E+00  0.15263E+00
+                                                                        w : -0.74072E-03  0.59509E+00
+                                                                        =================================
+                 Calculate pressures    213.4788      3.7725      1.3300
+                                                                        raw    pressures : -0.29933E+00  0.77822E-01
+                 Smoothing pressures    213.5429      3.8366      0.0641
+                                                                        smooth pressures : -0.10634E+00  0.00000E+00
+                do leaf measurements    213.5499      3.8436      0.0070
+       compute convergence criterion    213.5725      3.8662      0.0226
+                                                                        velocity_diff_norm = 0.0132731 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 5
+                    slicing the cube    213.5743      3.8680      0.0017
+Compute isostasy and adjust vertical    213.5744      3.8681      0.0001
+                  Compute divergence    213.5745      3.8682      0.0000
+                        wsmp_cleanup    213.5835      3.8772      0.0090
+              Reset surface geometry    213.5860      3.8797      0.0025
+ -----------------------------------------------------------------------
+           Temperature calculations     213.5870      3.8807      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    213.7466      4.0402      0.1596
+                   Advect surface  1    213.7467      4.0404      0.0002
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    213.7912      4.0849      0.0445
+                 Update cloud fields    231.0251     21.3188     17.2339
+                    Advect the cloud    231.4662     21.7599      0.4411
+                        Write output    233.8975     24.1912      2.4313
+                    End of time step    235.0059     25.2996      1.1084
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       8    235.0061
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    235.0062      0.0001      0.0001
+                          surface 01    235.0062      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    235.0063      0.0002      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    235.0119      0.0058      0.0056
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    235.0174
+ ----------------------------------------------------------------------- 
+                 Create octree solve    235.0176      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    235.0336      0.0162      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    235.0446      0.0273      0.0111
+                                                                        max(crit)        =   3.199999E+01
+                                                                        max(crit,w_l_i_f)=   1.484641E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement    235.0495      0.0321      0.0049
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon    235.1214      0.1040      0.0718
+             Imbed surface in osolve    235.2266      0.2092      0.1053
+                embedding surface  1    235.2268      0.2094      0.0002
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability    236.9232      1.9059      1.6964
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    236.9263      1.9089      0.0031
+        Interpolate velo onto osolve    238.3742      3.3568      1.4479
+                                                                        ov%pressure:     -0.1196E+00 0.9774E-03
+                                                                        osolve%pressure: -0.1063E+00 0.0000E+00
+                      Find bad faces    238.5092      3.4918      0.1350
+                           Find void    238.6441      3.6268      0.1349
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions    238.6525      3.6351      0.0084
+                         wsmp setup1    238.6566      3.6392      0.0041
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2    238.8919      3.8746      0.2353
+ -----------------------------------
+ start of non-linear iteratio      1    238.9033
+ -----------------------------------
+                        build system    238.9035      0.0002      0.0002
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   3.78953E-04  1.30880E+03
+                                                                        viscosity capped in      7168
+                          wsmp solve    241.3326      2.4293      2.4291
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0791759 s
+                                                                        wsmp: ordering time:               0.6803219 s
+                                                                        wsmp: symbolic fact. time:         0.1552930 s
+                                                                        wsmp: Cholesky fact. time:         0.4047692 s
+                                                                        wsmp: Factorisation            11855.4668610 Mflops
+                                                                        wsmp: back substitution time:      0.0149009 s
+                                                                        wsmp: from b to rhs vector:        0.0019600 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3367021 s
+                                                                        =================================
+                                                                        u : -0.59607E-12  0.10000E+01
+                                                                        v : -0.14717E+00  0.15143E+00
+                                                                        w : -0.48921E-03  0.58533E+00
+                                                                        =================================
+                 Calculate pressures    242.6751      3.7718      1.3425
+                                                                        raw    pressures : -0.30831E+00  0.85813E-01
+                 Smoothing pressures    242.7393      3.8360      0.0642
+                                                                        smooth pressures : -0.10651E+00  0.00000E+00
+                do leaf measurements    242.7457      3.8423      0.0064
+       compute convergence criterion    242.7689      3.8655      0.0232
+                                                                        velocity_diff_norm = 0.0118346 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 5
+                    slicing the cube    242.7706      3.8673      0.0017
+Compute isostasy and adjust vertical    242.7708      3.8675      0.0002
+                  Compute divergence    242.7708      3.8675      0.0001
+                        wsmp_cleanup    242.7798      3.8765      0.0089
+              Reset surface geometry    242.7824      3.8791      0.0026
+ -----------------------------------------------------------------------
+           Temperature calculations     242.7834      3.8801      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    242.9443      4.0410      0.1609
+                   Advect surface  1    242.9445      4.0411      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    242.9893      4.0860      0.0449
+                 Update cloud fields    258.9886     20.0852     15.9992
+                    Advect the cloud    259.4362     20.5328      0.4476
+                        Write output    261.9316     23.0283      2.4954
+                    End of time step    263.0471     24.1438      1.1155
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step       9    263.0473
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    263.0474      0.0001      0.0001
+                          surface 01    263.0474      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    263.0475      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    263.0531      0.0058      0.0056
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    263.0584
+ ----------------------------------------------------------------------- 
+                 Create octree solve    263.0587      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    263.0745      0.0161      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    263.0858      0.0274      0.0113
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   1.499447E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement    263.0910      0.0325      0.0051
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon    263.1621      0.1036      0.0711
+             Imbed surface in osolve    263.2663      0.2078      0.1042
+                embedding surface  1    263.2664      0.2080      0.0002
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability    264.9355      1.8771      1.6691
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    264.9385      1.8800      0.0030
+        Interpolate velo onto osolve    266.2787      3.2203      1.3402
+                                                                        ov%pressure:     -0.1194E+00 0.5359E-03
+                                                                        osolve%pressure: -0.1065E+00 0.0000E+00
+                      Find bad faces    266.4773      3.4188      0.1985
+                           Find void    266.6108      3.5523      0.1335
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions    266.6190      3.5605      0.0082
+                         wsmp setup1    266.6227      3.5643      0.0038
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2    266.8559      3.7974      0.2331
+ -----------------------------------
+ start of non-linear iteratio      1    266.8667
+ -----------------------------------
+                        build system    266.8669      0.0002      0.0002
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   3.55086E-04  1.30880E+03
+                                                                        viscosity capped in      7168
+                          wsmp solve    269.2786      2.4119      2.4118
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0754819 s
+                                                                        wsmp: ordering time:               0.6816471 s
+                                                                        wsmp: symbolic fact. time:         0.1535859 s
+                                                                        wsmp: Cholesky fact. time:         0.4066930 s
+                                                                        wsmp: Factorisation            11799.3863718 Mflops
+                                                                        wsmp: back substitution time:      0.0148981 s
+                                                                        wsmp: from b to rhs vector:        0.0019259 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3345132 s
+                                                                        =================================
+                                                                        u : -0.55656E-12  0.10000E+01
+                                                                        v : -0.14947E+00  0.15290E+00
+                                                                        w : -0.52167E-03  0.57760E+00
+                                                                        =================================
+                 Calculate pressures    270.6193      3.7526      1.3406
+                                                                        raw    pressures : -0.31625E+00  0.93330E-01
+                 Smoothing pressures    270.6830      3.8163      0.0637
+                                                                        smooth pressures : -0.10672E+00  0.00000E+00
+                do leaf measurements    270.6893      3.8226      0.0064
+       compute convergence criterion    270.7124      3.8457      0.0231
+                                                                        velocity_diff_norm = 0.0100831 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 5
+                    slicing the cube    270.7141      3.8474      0.0017
+Compute isostasy and adjust vertical    270.7143      3.8476      0.0002
+                  Compute divergence    270.7144      3.8476      0.0000
+                        wsmp_cleanup    270.7232      3.8565      0.0089
+              Reset surface geometry    270.7257      3.8590      0.0025
+ -----------------------------------------------------------------------
+           Temperature calculations     270.7267      3.8600      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    270.8864      4.0197      0.1597
+                   Advect surface  1    270.8865      4.0198      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    270.9312      4.0644      0.0446
+                 Update cloud fields    286.7809     19.9142     15.8498
+                    Advect the cloud    287.2280     20.3613      0.4471
+                        Write output    289.6612     22.7945      2.4332
+                    End of time step    290.7736     23.9069      1.1123
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      10    290.7737
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    290.7738      0.0001      0.0001
+                          surface 01    290.7738      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    290.7739      0.0002      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    290.7799      0.0061      0.0060
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    290.7838
+ ----------------------------------------------------------------------- 
+                 Create octree solve    290.7840      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct F
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    290.7999      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    290.8108      0.0270      0.0110
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   1.511215E+01
+                                                                        osolve%nleaves=   32768
+    Smoothen octree after refinement    290.8159      0.0321      0.0050
+                                                                        osolve%nleaves=   32768
+                   Build osolve icon    290.8870      0.1031      0.0711
+             Imbed surface in osolve    290.9907      0.2068      0.1037
+                embedding surface  1    290.9909      0.2071      0.0002
+                                                                        osolve%nleaves=   39936
+          Assessing octree stability    292.6510      1.8671      1.6601
+                                                                        current refine level     5
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    32768 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    39936 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    292.6538      1.8700      0.0029
+        Interpolate velo onto osolve    294.0016      3.2177      1.3477
+                                                                        ov%pressure:     -0.1193E+00 0.2277E-03
+                                                                        osolve%pressure: -0.1067E+00 0.0000E+00
+                      Find bad faces    294.2041      3.4202      0.2025
+                           Find void    294.3376      3.5538      0.1335
+                                                                        whole leaf in fluid     3072
+          Define boundary conditions    294.3452      3.5614      0.0076
+                         wsmp setup1    294.3489      3.5650      0.0036
+                                                                        n =    35151
+                                                                        nz=  1198722
+                         wsmp setup2    294.5799      3.7960      0.2310
+ -----------------------------------
+ start of non-linear iteratio      1    294.5901
+ -----------------------------------
+                        build system    294.5903      0.0002      0.0002
+                                                                        nelem per proc (min/max)         917        1083
+                                                                        viscosity range   3.76056E-04  1.30880E+03
+                                                                        viscosity capped in      7168
+                          wsmp solve    296.9896      2.3995      2.3993
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0699971 s
+                                                                        wsmp: ordering time:               0.6794989 s
+                                                                        wsmp: symbolic fact. time:         0.1538131 s
+                                                                        wsmp: Cholesky fact. time:         0.4043159 s
+                                                                        wsmp: Factorisation            11868.7567086 Mflops
+                                                                        wsmp: back substitution time:      0.0149090 s
+                                                                        wsmp: from b to rhs vector:        0.0019891 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            1.3248041 s
+                                                                        =================================
+                                                                        u : -0.47392E-12  0.10000E+01
+                                                                        v : -0.15038E+00  0.15322E+00
+                                                                        w : -0.64845E-03  0.57307E+00
+                                                                        =================================
+                 Calculate pressures    298.3199      3.7298      1.3303
+                                                                        raw    pressures : -0.32336E+00  0.99317E-01
+                 Smoothing pressures    298.3840      3.7939      0.0641
+                                                                        smooth pressures : -0.10693E+00  0.00000E+00
+                do leaf measurements    298.3906      3.8005      0.0066
+       compute convergence criterion    298.4134      3.8233      0.0227
+                                                                        velocity_diff_norm = 0.0081639 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    298.4152      3.8251      0.0018
+Compute isostasy and adjust vertical    298.4153      3.8252      0.0002
+                  Compute divergence    298.4154      3.8253      0.0000
+                        wsmp_cleanup    298.4244      3.8343      0.0090
+              Reset surface geometry    298.4270      3.8369      0.0025
+ -----------------------------------------------------------------------
+           Temperature calculations     298.4280      3.8379      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    298.5854      3.9953      0.1574
+                   Advect surface  1    298.5855      3.9954      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    298.6301      4.0400      0.0446
+                 Update cloud fields    316.4626     21.8725     17.8325
+                    Advect the cloud    316.9164     22.3263      0.4538
+                        Write output    319.3705     24.7804      2.4541
+                    End of time step    320.4993     25.9092      1.1288
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      11    320.4995
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    320.4996      0.0001      0.0001
+                          surface 01    320.4996      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    320.4996      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    320.5056      0.0061      0.0060
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    320.5096
+ ----------------------------------------------------------------------- 
+                 Create octree solve    320.5098      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    320.5257      0.0161      0.0159
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    320.5365      0.0269      0.0108
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   1.520419E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    320.5430      0.0334      0.0065
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    320.6775      0.1679      0.1345
+             Imbed surface in osolve    320.8657      0.3561      0.1881
+                embedding surface  1    320.8658      0.3562      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    323.4876      2.9780      2.6217
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    32768 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    39936 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  F
+                     Change 3D cloud    323.4907      2.9811      0.0032
+        Interpolate velo onto osolve    325.7133      5.2037      2.2226
+                                                                        ov%pressure:     -0.1193E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1127E+00 0.0000E+00
+                      Find bad faces    325.8697      5.3601      0.1564
+                           Find void    326.0870      5.5774      0.2173
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    326.1008      5.5912      0.0138
+                         wsmp setup1    326.1063      5.5967      0.0055
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    326.8381      6.3285      0.7318
+ -----------------------------------
+ start of non-linear iteratio      1    326.8688
+ -----------------------------------
+                        build system    326.8689      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   3.12793E-04  1.30880E+03
+                                                                        viscosity capped in     24576
+                          wsmp solve    329.6970      2.8282      2.8280
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0887301 s
+                                                                        wsmp: ordering time:               2.4339421 s
+                                                                        wsmp: symbolic fact. time:         0.5012250 s
+                                                                        wsmp: Cholesky fact. time:         3.7778420 s
+                                                                        wsmp: Factorisation            15788.4785550 Mflops
+                                                                        wsmp: back substitution time:      0.0778730 s
+                                                                        wsmp: from b to rhs vector:        0.0050821 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.8850181 s
+                                                                        =================================
+                                                                        u : -0.18699E-03  0.10000E+01
+                                                                        v : -0.15756E+00  0.14756E+00
+                                                                        w : -0.36089E-01  0.55778E+00
+                                                                        =================================
+                 Calculate pressures    336.5959      9.7271      6.8989
+                                                                        raw    pressures : -0.10275E+01  0.79132E+00
+                 Smoothing pressures    336.6703      9.8015      0.0745
+                                                                        smooth pressures : -0.12947E+00  0.00000E+00
+                do leaf measurements    336.6824      9.8136      0.0120
+       compute convergence criterion    336.7206      9.8518      0.0382
+                                                                        velocity_diff_norm = 0.1457759 (tol = 0.0100000)
+                                                                        increase_current_level -> F
+                                                                        current_level -> 6
+                    slicing the cube    336.7234      9.8547      0.0029
+Compute isostasy and adjust vertical    336.7236      9.8548      0.0002
+                  Compute divergence    336.7237      9.8549      0.0000
+                        wsmp_cleanup    336.7525      9.8837      0.0288
+              Reset surface geometry    336.7577      9.8889      0.0052
+ -----------------------------------------------------------------------
+           Temperature calculations     336.7588      9.8900      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    336.9229     10.0541      0.1641
+                   Advect surface  1    336.9231     10.0543      0.0002
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    336.9673     10.0986      0.0443
+                 Update cloud fields    354.3455     27.4767     17.3781
+                    Advect the cloud    355.0606     28.1918      0.7151
+                        Write output    358.9913     32.1226      3.9307
+                    End of time step    360.6920     33.8232      1.7007
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      12    360.6922
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    360.6922      0.0001      0.0001
+                          surface 01    360.6923      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    360.6923      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    360.6979      0.0057      0.0056
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    360.7016
+ ----------------------------------------------------------------------- 
+                 Create octree solve    360.7017      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             F
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    360.7177      0.0161      0.0159
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    360.7360      0.0345      0.0184
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   2.915290E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    360.7584      0.0569      0.0224
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    360.8930      0.1914      0.1345
+             Imbed surface in osolve    361.0827      0.3811      0.1897
+                embedding surface  1    361.0829      0.3813      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    363.7070      3.0055      2.6242
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    363.7100      3.0084      0.0029
+        Interpolate velo onto osolve    365.8571      5.1556      2.1471
+                                                                        ov%pressure:     -0.1434E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1295E+00 0.0000E+00
+                      Find bad faces    366.0194      5.3178      0.1622
+                           Find void    366.2367      5.5352      0.2174
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    366.2501      5.5486      0.0134
+                         wsmp setup1    366.2558      5.5542      0.0056
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    366.9886      6.2870      0.7328
+ -----------------------------------
+ start of non-linear iteratio      1    367.0190
+ -----------------------------------
+                        build system    367.0192      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   3.00250E-04  1.30880E+03
+                                                                        viscosity capped in     28672
+                          wsmp solve    369.8267      2.8077      2.8075
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0891190 s
+                                                                        wsmp: ordering time:               2.4336569 s
+                                                                        wsmp: symbolic fact. time:         0.5020540 s
+                                                                        wsmp: Cholesky fact. time:         3.7854180 s
+                                                                        wsmp: Factorisation            15756.8801067 Mflops
+                                                                        wsmp: back substitution time:      0.0790250 s
+                                                                        wsmp: from b to rhs vector:        0.0054049 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.8950250 s
+                                                                        =================================
+                                                                        u : -0.36832E-04  0.10000E+01
+                                                                        v : -0.15431E+00  0.14290E+00
+                                                                        w : -0.17754E-01  0.55937E+00
+                                                                        =================================
+                 Calculate pressures    376.7362      9.7172      6.9095
+                                                                        raw    pressures : -0.76122E+00  0.53638E+00
+                 Smoothing pressures    376.8107      9.7917      0.0745
+                                                                        smooth pressures : -0.12744E+00  0.00000E+00
+                do leaf measurements    376.8261      9.8071      0.0154
+       compute convergence criterion    376.8655      9.8465      0.0394
+                                                                        velocity_diff_norm = 0.0234088 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    376.8683      9.8493      0.0028
+Compute isostasy and adjust vertical    376.8685      9.8495      0.0002
+                  Compute divergence    376.8685      9.8495      0.0000
+                        wsmp_cleanup    376.8972      9.8782      0.0287
+              Reset surface geometry    376.9027      9.8836      0.0055
+ -----------------------------------------------------------------------
+           Temperature calculations     376.9038      9.8847      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    377.0676     10.0486      0.1638
+                   Advect surface  1    377.0677     10.0487      0.0002
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    377.1120     10.0930      0.0443
+                 Update cloud fields    394.6390     27.6199     17.5269
+                    Advect the cloud    395.3657     28.3467      0.7268
+                        Write output    399.3244     32.3053      3.9586
+                    End of time step    401.0014     33.9824      1.6770
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      13    401.0015
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    401.0016      0.0001      0.0001
+                          surface 01    401.0016      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    401.0017      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    401.0076      0.0061      0.0059
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    401.0119
+ ----------------------------------------------------------------------- 
+                 Create octree solve    401.0120      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    401.0280      0.0161      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    401.0468      0.0350      0.0188
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   2.928906E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    401.0691      0.0573      0.0223
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    401.2037      0.1918      0.1346
+             Imbed surface in osolve    401.3920      0.3801      0.1883
+                embedding surface  1    401.3922      0.3804      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    404.0241      3.0122      2.6319
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    404.0270      3.0151      0.0029
+        Interpolate velo onto osolve    406.1863      5.1744      2.1593
+                                                                        ov%pressure:     -0.1387E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1274E+00 0.0000E+00
+                      Find bad faces    406.3555      5.3437      0.1693
+                           Find void    406.5728      5.5610      0.2173
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    406.5865      5.5746      0.0137
+                         wsmp setup1    406.5921      5.5802      0.0056
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    407.3247      6.3128      0.7326
+ -----------------------------------
+ start of non-linear iteratio      1    407.3584
+ -----------------------------------
+                        build system    407.3586      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   2.98579E-04  1.30880E+03
+                                                                        viscosity capped in     28672
+                          wsmp solve    410.2218      2.8634      2.8632
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0976439 s
+                                                                        wsmp: ordering time:               2.4295881 s
+                                                                        wsmp: symbolic fact. time:         0.5020862 s
+                                                                        wsmp: Cholesky fact. time:         3.7778170 s
+                                                                        wsmp: Factorisation            15788.5831784 Mflops
+                                                                        wsmp: back substitution time:      0.0786929 s
+                                                                        wsmp: from b to rhs vector:        0.0055070 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.8916721 s
+                                                                        =================================
+                                                                        u : -0.37596E-05  0.10000E+01
+                                                                        v : -0.15112E+00  0.13944E+00
+                                                                        w : -0.10320E-01  0.55174E+00
+                                                                        =================================
+                 Calculate pressures    417.1276      9.7692      6.9058
+                                                                        raw    pressures : -0.77542E+00  0.55106E+00
+                 Smoothing pressures    417.2018      9.8434      0.0742
+                                                                        smooth pressures : -0.12755E+00  0.00000E+00
+                do leaf measurements    417.2145      9.8561      0.0127
+       compute convergence criterion    417.2548      9.8964      0.0403
+                                                                        velocity_diff_norm = 0.0118898 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    417.2576      9.8992      0.0028
+Compute isostasy and adjust vertical    417.2578      9.8994      0.0002
+                  Compute divergence    417.2578      9.8995      0.0000
+                        wsmp_cleanup    417.2870      9.9286      0.0292
+              Reset surface geometry    417.2923      9.9339      0.0053
+ -----------------------------------------------------------------------
+           Temperature calculations     417.2934      9.9350      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    417.4575     10.0991      0.1641
+                   Advect surface  1    417.4577     10.0993      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    417.5020     10.1436      0.0443
+                 Update cloud fields    435.9409     28.5825     18.4389
+                    Advect the cloud    436.6613     29.3030      0.7205
+                        Write output    440.7111     33.3527      4.0498
+                    End of time step    442.3959     35.0375      1.6848
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      14    442.3961
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    442.3962      0.0001      0.0001
+                          surface 01    442.3962      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    442.3963      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    442.4019      0.0058      0.0057
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    442.4053
+ ----------------------------------------------------------------------- 
+                 Create octree solve    442.4055      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    442.4215      0.0162      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    442.4400      0.0347      0.0185
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   2.948623E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    442.4623      0.0570      0.0223
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    442.5969      0.1916      0.1346
+             Imbed surface in osolve    442.7860      0.3807      0.1891
+                embedding surface  1    442.7862      0.3809      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    445.3975      2.9922      2.6113
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    445.4004      2.9951      0.0029
+        Interpolate velo onto osolve    447.5714      5.1661      2.1710
+                                                                        ov%pressure:     -0.1382E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1276E+00 0.0000E+00
+                      Find bad faces    447.7872      5.3819      0.2158
+                           Find void    448.0045      5.5992      0.2173
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    448.0177      5.6124      0.0133
+                         wsmp setup1    448.0231      5.6178      0.0054
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    448.7552      6.3499      0.7321
+ -----------------------------------
+ start of non-linear iteratio      1    448.7865
+ -----------------------------------
+                        build system    448.7866      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   2.94276E-04  1.30880E+03
+                                                                        viscosity capped in     28672
+                          wsmp solve    451.6416      2.8552      2.8550
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0904090 s
+                                                                        wsmp: ordering time:               2.4321890 s
+                                                                        wsmp: symbolic fact. time:         0.5038891 s
+                                                                        wsmp: Cholesky fact. time:         3.7824950 s
+                                                                        wsmp: Factorisation            15769.0566055 Mflops
+                                                                        wsmp: back substitution time:      0.0789418 s
+                                                                        wsmp: from b to rhs vector:        0.0054200 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.8936601 s
+                                                                        =================================
+                                                                        u : -0.76317E-06  0.10000E+01
+                                                                        v : -0.14878E+00  0.13828E+00
+                                                                        w : -0.10632E-01  0.54416E+00
+                                                                        =================================
+                 Calculate pressures    458.5494      9.7630      6.9078
+                                                                        raw    pressures : -0.79145E+00  0.56793E+00
+                 Smoothing pressures    458.6240      9.8376      0.0746
+                                                                        smooth pressures : -0.12765E+00  0.00000E+00
+                do leaf measurements    458.6401      9.8536      0.0160
+       compute convergence criterion    458.6793      9.8928      0.0392
+                                                                        velocity_diff_norm = 0.0078505 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    458.6821      9.8956      0.0028
+Compute isostasy and adjust vertical    458.6823      9.8958      0.0002
+                  Compute divergence    458.6823      9.8958      0.0000
+                        wsmp_cleanup    458.7118      9.9253      0.0295
+              Reset surface geometry    458.7171      9.9306      0.0053
+ -----------------------------------------------------------------------
+           Temperature calculations     458.7182      9.9317      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    458.8828     10.0964      0.1646
+                   Advect surface  1    458.8830     10.0965      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    458.9273     10.1409      0.0444
+                 Update cloud fields    475.7126     26.9262     16.7853
+                    Advect the cloud    476.4529     27.6664      0.7402
+                        Write output    480.5325     31.7460      4.0796
+                    End of time step    482.2407     33.4543      1.7082
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      15    482.2409
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    482.2409      0.0001      0.0001
+                          surface 01    482.2410      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    482.2410      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    482.2477      0.0068      0.0067
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    482.2511
+ ----------------------------------------------------------------------- 
+                 Create octree solve    482.2512      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    482.2672      0.0161      0.0159
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    482.2851      0.0340      0.0179
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   2.967086E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    482.3074      0.0563      0.0223
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    482.4420      0.1909      0.1345
+             Imbed surface in osolve    482.6330      0.3819      0.1910
+                embedding surface  1    482.6331      0.3820      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    485.2482      2.9972      2.6151
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    485.2509      2.9998      0.0027
+        Interpolate velo onto osolve    487.4471      5.1961      2.1962
+                                                                        ov%pressure:     -0.1379E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1277E+00 0.0000E+00
+                      Find bad faces    487.6492      5.3981      0.2020
+                           Find void    487.8665      5.6154      0.2173
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    487.8794      5.6283      0.0129
+                         wsmp setup1    487.8849      5.6338      0.0055
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    488.6111      6.3600      0.7262
+ -----------------------------------
+ start of non-linear iteratio      1    488.6442
+ -----------------------------------
+                        build system    488.6444      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   2.88088E-04  1.30880E+03
+                                                                        viscosity capped in     28672
+                          wsmp solve    491.4410      2.7968      2.7966
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0867100 s
+                                                                        wsmp: ordering time:               2.4317250 s
+                                                                        wsmp: symbolic fact. time:         0.5033760 s
+                                                                        wsmp: Cholesky fact. time:         3.7859449 s
+                                                                        wsmp: Factorisation            15754.6871588 Mflops
+                                                                        wsmp: back substitution time:      0.0784519 s
+                                                                        wsmp: from b to rhs vector:        0.0054822 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.8920090 s
+                                                                        =================================
+                                                                        u : -0.15777E-07  0.10000E+01
+                                                                        v : -0.14698E+00  0.13631E+00
+                                                                        w : -0.10792E-01  0.53755E+00
+                                                                        =================================
+                 Calculate pressures    498.3468      9.7025      6.9057
+                                                                        raw    pressures : -0.80449E+00  0.58380E+00
+                 Smoothing pressures    498.4207      9.7765      0.0740
+                                                                        smooth pressures : -0.12778E+00  0.00000E+00
+                do leaf measurements    498.4357      9.7915      0.0149
+       compute convergence criterion    498.4739      9.8297      0.0382
+                                                                        velocity_diff_norm = 0.0067848 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    498.4767      9.8325      0.0028
+Compute isostasy and adjust vertical    498.4769      9.8327      0.0002
+                  Compute divergence    498.4769      9.8327      0.0000
+                        wsmp_cleanup    498.5060      9.8618      0.0290
+              Reset surface geometry    498.5112      9.8670      0.0053
+ -----------------------------------------------------------------------
+           Temperature calculations     498.5123      9.8681      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    498.6766     10.0324      0.1643
+                   Advect surface  1    498.6767     10.0325      0.0002
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    498.7201     10.0759      0.0434
+                 Update cloud fields    514.6166     25.9723     15.8964
+                    Advect the cloud    515.3563     26.7121      0.7397
+                        Write output    519.4697     30.8255      4.1134
+                    End of time step    521.1991     32.5549      1.7295
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      16    521.1993
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    521.1994      0.0001      0.0001
+                          surface 01    521.1994      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    521.1994      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    521.2058      0.0065      0.0064
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    521.2101
+ ----------------------------------------------------------------------- 
+                 Create octree solve    521.2103      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    521.2262      0.0161      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    521.2443      0.0343      0.0181
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   2.983220E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    521.2667      0.0566      0.0224
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    521.4014      0.1913      0.1346
+             Imbed surface in osolve    521.5917      0.3816      0.1904
+                embedding surface  1    521.5920      0.3819      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    524.2117      3.0016      2.6198
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    524.2139      3.0039      0.0022
+        Interpolate velo onto osolve    526.4142      5.2041      2.2002
+                                                                        ov%pressure:     -0.1378E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1278E+00 0.0000E+00
+                      Find bad faces    526.6171      5.4070      0.2030
+                           Find void    526.8346      5.6245      0.2175
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    526.8473      5.6372      0.0127
+                         wsmp setup1    526.8527      5.6426      0.0054
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    527.5828      6.3727      0.7301
+ -----------------------------------
+ start of non-linear iteratio      1    527.6167
+ -----------------------------------
+                        build system    527.6169      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   2.81409E-04  1.30880E+03
+                                                                        viscosity capped in     28672
+                          wsmp solve    530.4102      2.7935      2.7933
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0885749 s
+                                                                        wsmp: ordering time:               2.4330220 s
+                                                                        wsmp: symbolic fact. time:         0.5014760 s
+                                                                        wsmp: Cholesky fact. time:         3.7836859 s
+                                                                        wsmp: Factorisation            15764.0933547 Mflops
+                                                                        wsmp: back substitution time:      0.0784702 s
+                                                                        wsmp: from b to rhs vector:        0.0054588 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.8910060 s
+                                                                        =================================
+                                                                        u :  0.00000E+00  0.10000E+01
+                                                                        v : -0.14574E+00  0.13456E+00
+                                                                        w : -0.10627E-01  0.53222E+00
+                                                                        =================================
+                 Calculate pressures    537.3155      9.6988      6.9053
+                                                                        raw    pressures : -0.81606E+00  0.59721E+00
+                 Smoothing pressures    537.3902      9.7735      0.0747
+                                                                        smooth pressures : -0.12792E+00  0.00000E+00
+                do leaf measurements    537.4034      9.7867      0.0132
+       compute convergence criterion    537.4434      9.8267      0.0400
+                                                                        velocity_diff_norm = 0.0060967 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    537.4462      9.8295      0.0028
+Compute isostasy and adjust vertical    537.4464      9.8297      0.0002
+                  Compute divergence    537.4464      9.8297      0.0000
+                        wsmp_cleanup    537.4756      9.8589      0.0292
+              Reset surface geometry    537.4808      9.8641      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations     537.4819      9.8652      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    537.6436     10.0269      0.1617
+                   Advect surface  1    537.6437     10.0270      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    537.6881     10.0714      0.0444
+                 Update cloud fields    554.1161     26.4994     16.4279
+                    Advect the cloud    554.8653     27.2486      0.7492
+                        Write output    559.0129     31.3962      4.1476
+                    End of time step    560.7327     33.1160      1.7198
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      17    560.7328
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    560.7329      0.0001      0.0001
+                          surface 01    560.7329      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    560.7330      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    560.7393      0.0065      0.0064
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    560.7426
+ ----------------------------------------------------------------------- 
+                 Create octree solve    560.7427      0.0001      0.0001
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    560.7587      0.0161      0.0159
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    560.7768      0.0342      0.0181
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   2.996925E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    560.7992      0.0566      0.0224
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    560.9337      0.1911      0.1345
+             Imbed surface in osolve    561.1227      0.3801      0.1890
+                embedding surface  1    561.1229      0.3803      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    563.7501      3.0076      2.6273
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    563.7524      3.0098      0.0023
+        Interpolate velo onto osolve    565.9709      5.2283      2.2185
+                                                                        ov%pressure:     -0.1377E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1279E+00 0.0000E+00
+                      Find bad faces    566.1782      5.4356      0.2072
+                           Find void    566.3955      5.6530      0.2174
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    566.4082      5.6657      0.0127
+                         wsmp setup1    566.4137      5.6711      0.0055
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    567.1428      6.4002      0.7291
+ -----------------------------------
+ start of non-linear iteratio      1    567.1764
+ -----------------------------------
+                        build system    567.1766      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   2.75089E-04  1.30880E+03
+                                                                        viscosity capped in     28672
+                          wsmp solve    569.9681      2.7917      2.7916
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0944800 s
+                                                                        wsmp: ordering time:               2.4319880 s
+                                                                        wsmp: symbolic fact. time:         0.5021811 s
+                                                                        wsmp: Cholesky fact. time:         3.7940481 s
+                                                                        wsmp: Factorisation            15721.0391103 Mflops
+                                                                        wsmp: back substitution time:      0.0784080 s
+                                                                        wsmp: from b to rhs vector:        0.0052350 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.9066861 s
+                                                                        =================================
+                                                                        u : -0.11323E-06  0.10000E+01
+                                                                        v : -0.14421E+00  0.13390E+00
+                                                                        w : -0.10438E-01  0.52778E+00
+                                                                        =================================
+                 Calculate pressures    576.8892      9.7128      6.9211
+                                                                        raw    pressures : -0.82527E+00  0.60808E+00
+                 Smoothing pressures    576.9634      9.7870      0.0741
+                                                                        smooth pressures : -0.12804E+00  0.00000E+00
+                do leaf measurements    576.9762      9.7999      0.0129
+       compute convergence criterion    577.0165      9.8401      0.0402
+                                                                        velocity_diff_norm = 0.0056287 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    577.0194      9.8430      0.0029
+Compute isostasy and adjust vertical    577.0195      9.8432      0.0002
+                  Compute divergence    577.0196      9.8432      0.0000
+                        wsmp_cleanup    577.0487      9.8723      0.0291
+              Reset surface geometry    577.0540      9.8776      0.0053
+ -----------------------------------------------------------------------
+           Temperature calculations     577.0550      9.8786      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    577.2190     10.0426      0.1640
+                   Advect surface  1    577.2192     10.0428      0.0002
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    577.2625     10.0861      0.0434
+                 Update cloud fields    594.6389     27.4625     17.3764
+                    Advect the cloud    595.3967     28.2203      0.7578
+                        Write output    599.5246     32.3483      4.1280
+                    End of time step    601.2704     34.0940      1.7457
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      18    601.2705
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    601.2706      0.0001      0.0001
+                          surface 01    601.2706      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    601.2707      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    601.2766      0.0061      0.0060
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    601.2799
+ ----------------------------------------------------------------------- 
+                 Create octree solve    601.2800      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    601.2960      0.0161      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    601.3139      0.0340      0.0179
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.008536E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    601.3355      0.0556      0.0216
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    601.4701      0.1902      0.1345
+             Imbed surface in osolve    601.6598      0.3799      0.1897
+                embedding surface  1    601.6600      0.3801      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    604.2778      2.9979      2.6178
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    604.2800      3.0001      0.0023
+        Interpolate velo onto osolve    606.6433      5.3634      2.3633
+                                                                        ov%pressure:     -0.1376E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1280E+00 0.0000E+00
+                      Find bad faces    606.8321      5.5522      0.1888
+                           Find void    607.0494      5.7695      0.2172
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    607.0620      5.7821      0.0126
+                         wsmp setup1    607.0674      5.7875      0.0055
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    607.7977      6.5178      0.7302
+ -----------------------------------
+ start of non-linear iteratio      1    607.8312
+ -----------------------------------
+                        build system    607.8314      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   2.69029E-04  1.30880E+03
+                                                                        viscosity capped in     28672
+                          wsmp solve    610.6499      2.8187      2.8185
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0901132 s
+                                                                        wsmp: ordering time:               2.4335821 s
+                                                                        wsmp: symbolic fact. time:         0.5037811 s
+                                                                        wsmp: Cholesky fact. time:         3.7978470 s
+                                                                        wsmp: Factorisation            15705.3134573 Mflops
+                                                                        wsmp: back substitution time:      0.0787120 s
+                                                                        wsmp: from b to rhs vector:        0.0054719 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.9098411 s
+                                                                        =================================
+                                                                        u : -0.50588E-06  0.10000E+01
+                                                                        v : -0.14134E+00  0.13352E+00
+                                                                        w : -0.10299E-01  0.52375E+00
+                                                                        =================================
+                 Calculate pressures    617.5736      9.7424      6.9236
+                                                                        raw    pressures : -0.83809E+00  0.62122E+00
+                 Smoothing pressures    617.6483      9.8170      0.0747
+                                                                        smooth pressures : -0.12817E+00  0.00000E+00
+                do leaf measurements    617.6614      9.8302      0.0131
+       compute convergence criterion    617.7007      9.8695      0.0394
+                                                                        velocity_diff_norm = 0.0056752 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    617.7037      9.8725      0.0029
+Compute isostasy and adjust vertical    617.7039      9.8726      0.0002
+                  Compute divergence    617.7039      9.8727      0.0000
+                        wsmp_cleanup    617.7329      9.9017      0.0290
+              Reset surface geometry    617.7384      9.9072      0.0055
+ -----------------------------------------------------------------------
+           Temperature calculations     617.7394      9.9082      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    617.9031     10.0719      0.1637
+                   Advect surface  1    617.9032     10.0720      0.0002
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    617.9468     10.1156      0.0436
+                 Update cloud fields    634.3657     26.5345     16.4189
+                    Advect the cloud    635.1310     27.2998      0.7653
+                        Write output    639.3206     31.4893      4.1895
+                    End of time step    641.0921     33.2609      1.7715
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      19    641.0922
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    641.0923      0.0001      0.0001
+                          surface 01    641.0923      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    641.0924      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    641.0992      0.0070      0.0068
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    641.1024
+ ----------------------------------------------------------------------- 
+                 Create octree solve    641.1026      0.0001      0.0001
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    641.1185      0.0161      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    641.1368      0.0344      0.0183
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.018356E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    641.1583      0.0559      0.0215
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    641.2928      0.1904      0.1345
+             Imbed surface in osolve    641.4795      0.3771      0.1867
+                embedding surface  1    641.4797      0.3773      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    644.0961      2.9937      2.6164
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    644.0984      2.9959      0.0023
+        Interpolate velo onto osolve    646.3426      5.2402      2.2443
+                                                                        ov%pressure:     -0.1375E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1282E+00 0.0000E+00
+                      Find bad faces    646.5541      5.4516      0.2114
+                           Find void    646.7713      5.6689      0.2173
+                                                                        whole leaf in fluid    24576
+          Define boundary conditions    646.7838      5.6814      0.0125
+                         wsmp setup1    646.7891      5.6867      0.0053
+                                                                        n =   101400
+                                                                        nz=  3738351
+                         wsmp setup2    647.5147      6.4123      0.7256
+ -----------------------------------
+ start of non-linear iteratio      1    647.5480
+ -----------------------------------
+                        build system    647.5481      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3641        4097
+                                                                        viscosity range   2.63387E-04  1.30880E+03
+                                                                        viscosity capped in     28672
+                          wsmp solve    650.3795      2.8315      2.8313
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0897181 s
+                                                                        wsmp: ordering time:               2.4317980 s
+                                                                        wsmp: symbolic fact. time:         0.5004609 s
+                                                                        wsmp: Cholesky fact. time:         3.7887712 s
+                                                                        wsmp: Factorisation            15742.9350331 Mflops
+                                                                        wsmp: back substitution time:      0.0785191 s
+                                                                        wsmp: from b to rhs vector:        0.0055020 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            6.8950889 s
+                                                                        =================================
+                                                                        u : -0.18988E-05  0.10000E+01
+                                                                        v : -0.13996E+00  0.13389E+00
+                                                                        w : -0.10105E-01  0.52009E+00
+                                                                        =================================
+                 Calculate pressures    657.2885      9.7405      6.9090
+                                                                        raw    pressures : -0.85186E+00  0.63770E+00
+                 Smoothing pressures    657.3632      9.8152      0.0748
+                                                                        smooth pressures : -0.12831E+00  0.00000E+00
+                do leaf measurements    657.3755      9.8276      0.0123
+       compute convergence criterion    657.4150      9.8670      0.0395
+                                                                        velocity_diff_norm = 0.0055571 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    657.4180      9.8700      0.0030
+Compute isostasy and adjust vertical    657.4182      9.8702      0.0002
+                  Compute divergence    657.4182      9.8703      0.0000
+                        wsmp_cleanup    657.4471      9.8991      0.0288
+              Reset surface geometry    657.4522      9.9042      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations     657.4533      9.9053      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    657.6182     10.0703      0.1650
+                   Advect surface  1    657.6184     10.0704      0.0002
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    657.6618     10.1138      0.0434
+                 Update cloud fields    674.9177     27.3697     17.2559
+                    Advect the cloud    675.6797     28.1317      0.7620
+                        Write output    679.8450     32.2970      4.1653
+                    End of time step    681.5807     34.0328      1.7358
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      20    681.5809
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    681.5810      0.0001      0.0001
+                          surface 01    681.5810      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    681.5811      0.0002      0.0001
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    681.5859      0.0049      0.0048
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    681.5892
+ ----------------------------------------------------------------------- 
+                 Create octree solve    681.5893      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    681.6053      0.0161      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    681.6233      0.0341      0.0180
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.026898E+01
+                                                                        osolve%nleaves=   54272
+    Smoothen octree after refinement    681.6449      0.0557      0.0216
+                                                                        osolve%nleaves=   54272
+                   Build osolve icon    681.7794      0.1902      0.1345
+             Imbed surface in osolve    681.9674      0.3782      0.1880
+                embedding surface  1    681.9676      0.3784      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    684.5859      2.9967      2.6183
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54272 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    684.5882      2.9990      0.0024
+        Interpolate velo onto osolve    686.8307      5.2415      2.2425
+                                                                        ov%pressure:     -0.1375E+00 0.0000E+00
+                                                                        osolve%pressure: -0.1283E+00 0.0000E+00
+                      Find bad faces    687.0289      5.4398      0.1983
+                           Find void    687.2462      5.6571      0.2173
+                                                                        whole leaf in fluid    24621
+          Define boundary conditions    687.2586      5.6695      0.0124
+                         wsmp setup1    687.2640      5.6748      0.0054
+                                                                        n =   101889
+                                                                        nz=  3757095
+                         wsmp setup2    687.9993      6.4101      0.7353
+ -----------------------------------
+ start of non-linear iteratio      1    688.0306
+ -----------------------------------
+                        build system    688.0307      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3658        4119
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     28791
+                          wsmp solve    690.8922      2.8616      2.8614
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0870979 s
+                                                                        wsmp: ordering time:               2.4473710 s
+                                                                        wsmp: symbolic fact. time:         0.5054629 s
+                                                                        wsmp: Cholesky fact. time:         3.9979291 s
+                                                                        wsmp: Factorisation            14940.9905232 Mflops
+                                                                        wsmp: back substitution time:      0.0763080 s
+                                                                        wsmp: from b to rhs vector:        0.0057020 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.1202030 s
+                                                                        =================================
+                                                                        u : -0.63375E-05  0.10000E+01
+                                                                        v : -0.13970E+00  0.13361E+00
+                                                                        w : -0.98687E-02  0.56087E+00
+                                                                        =================================
+                 Calculate pressures    698.0272      9.9966      7.1351
+                                                                        raw    pressures : -0.86359E+00  0.65168E+00
+                 Smoothing pressures    698.1014     10.0708      0.0741
+                                                                        smooth pressures : -0.12849E+00  0.43763E-04
+                do leaf measurements    698.1161     10.0855      0.0147
+       compute convergence criterion    698.1543     10.1237      0.0382
+                                                                        velocity_diff_norm = 0.0055830 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    698.1572     10.1266      0.0029
+Compute isostasy and adjust vertical    698.1574     10.1268      0.0002
+                  Compute divergence    698.1574     10.1269      0.0000
+                        wsmp_cleanup    698.1868     10.1563      0.0294
+              Reset surface geometry    698.1921     10.1615      0.0053
+ -----------------------------------------------------------------------
+           Temperature calculations     698.1931     10.1625      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    698.3561     10.3255      0.1630
+                   Advect surface  1    698.3562     10.3256      0.0001
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    698.4005     10.3699      0.0443
+                 Update cloud fields    715.8961     27.8655     17.4956
+                    Advect the cloud    716.6584     28.6279      0.7623
+                        Write output    720.8574     32.8269      4.1990
+                    End of time step    722.5987     34.5681      1.7412
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      21    722.5988
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    722.5989      0.0001      0.0001
+                          surface 01    722.5989      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    722.5990      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    722.6054      0.0066      0.0065
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    722.6087
+ ----------------------------------------------------------------------- 
+                 Create octree solve    722.6089      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    722.6249      0.0161      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    722.6427      0.0340      0.0179
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.034684E+01
+                                                                        osolve%nleaves=   54377
+    Smoothen octree after refinement    722.6644      0.0557      0.0217
+                                                                        osolve%nleaves=   54377
+                   Build osolve icon    722.7992      0.1905      0.1348
+             Imbed surface in osolve    722.9885      0.3798      0.1893
+                embedding surface  1    722.9887      0.3800      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    725.6015      2.9928      2.6128
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54272 leaves
+                                                                        nleaves_new=    54377 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    725.6040      2.9953      0.0025
+        Interpolate velo onto osolve    727.8505      5.2418      2.2464
+                                                                        ov%pressure:     -0.1375E+00 0.8459E-03
+                                                                        osolve%pressure: -0.1285E+00 0.4376E-04
+                      Find bad faces    728.0706      5.4619      0.2201
+                           Find void    728.2879      5.6792      0.2173
+                                                                        whole leaf in fluid    24659
+          Define boundary conditions    728.3004      5.6917      0.0125
+                         wsmp setup1    728.3058      5.6971      0.0054
+                                                                        n =   102042
+                                                                        nz=  3765411
+                         wsmp setup2    729.0382      6.4295      0.7323
+ -----------------------------------
+ start of non-linear iteratio      1    729.0719
+ -----------------------------------
+                        build system    729.0720      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3663        4125
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     28839
+                          wsmp solve    731.9534      2.8815      2.8813
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0918279 s
+                                                                        wsmp: ordering time:               2.5783439 s
+                                                                        wsmp: symbolic fact. time:         0.5032260 s
+                                                                        wsmp: Cholesky fact. time:         3.7912669 s
+                                                                        wsmp: Factorisation            15630.0435648 Mflops
+                                                                        wsmp: back substitution time:      0.0764549 s
+                                                                        wsmp: from b to rhs vector:        0.0055189 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.0470071 s
+                                                                        =================================
+                                                                        u : -0.20467E-04  0.10000E+01
+                                                                        v : -0.13990E+00  0.13351E+00
+                                                                        w : -0.10454E-01  0.53754E+00
+                                                                        =================================
+                 Calculate pressures    739.0143      9.9424      7.0609
+                                                                        raw    pressures : -0.87665E+00  0.66824E+00
+                 Smoothing pressures    739.0900     10.0182      0.0758
+                                                                        smooth pressures : -0.12860E+00  0.92608E-03
+                do leaf measurements    739.1033     10.0314      0.0133
+       compute convergence criterion    739.1424     10.0705      0.0390
+                                                                        velocity_diff_norm = 0.0059488 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    739.1452     10.0734      0.0029
+Compute isostasy and adjust vertical    739.1454     10.0735      0.0002
+                  Compute divergence    739.1455     10.0736      0.0000
+                        wsmp_cleanup    739.1751     10.1033      0.0297
+              Reset surface geometry    739.1803     10.1084      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations     739.1813     10.1094      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    739.3460     10.2741      0.1647
+                   Advect surface  1    739.3462     10.2743      0.0002
+                                                                        S. 1:    3 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    739.3905     10.3186      0.0444
+                 Update cloud fields    756.6507     27.5788     17.2602
+                    Advect the cloud    757.4169     28.3450      0.7662
+                        Write output    761.6245     32.5526      4.2076
+                    End of time step    763.3996     34.3277      1.7751
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      22    763.3998
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    763.3998      0.0001      0.0001
+                          surface 01    763.3999      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    763.3999      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    763.4061      0.0064      0.0062
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    763.4095
+ ----------------------------------------------------------------------- 
+                 Create octree solve    763.4096      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    763.4256      0.0161      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    763.4437      0.0342      0.0181
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.042006E+01
+                                                                        osolve%nleaves=   54496
+    Smoothen octree after refinement    763.4651      0.0557      0.0214
+                                                                        osolve%nleaves=   54496
+                   Build osolve icon    763.6003      0.1909      0.1352
+             Imbed surface in osolve    763.7970      0.3875      0.1967
+                embedding surface  1    763.7972      0.3877      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    766.3933      2.9839      2.5962
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54377 leaves
+                                                                        nleaves_new=    54496 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    766.3955      2.9861      0.0022
+        Interpolate velo onto osolve    768.6162      5.2067      2.2206
+                                                                        ov%pressure:     -0.1375E+00 0.5102E-02
+                                                                        osolve%pressure: -0.1286E+00 0.9261E-03
+                      Find bad faces    768.9049      5.4955      0.2888
+                           Find void    769.1222      5.7127      0.2172
+                                                                        whole leaf in fluid    24704
+          Define boundary conditions    769.1347      5.7252      0.0125
+                         wsmp setup1    769.1401      5.7306      0.0054
+                                                                        n =   102243
+                                                                        nz=  3774345
+                         wsmp setup2    769.8758      6.4663      0.7357
+ -----------------------------------
+ start of non-linear iteratio      1    769.9094
+ -----------------------------------
+                        build system    769.9095      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3671        4135
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     28897
+                          wsmp solve    772.7621      2.8527      2.8526
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0917790 s
+                                                                        wsmp: ordering time:               2.4338601 s
+                                                                        wsmp: symbolic fact. time:         0.5063221 s
+                                                                        wsmp: Cholesky fact. time:         3.9091449 s
+                                                                        wsmp: Factorisation            15408.5172312 Mflops
+                                                                        wsmp: back substitution time:      0.0751650 s
+                                                                        wsmp: from b to rhs vector:        0.0057039 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.0223429 s
+                                                                        =================================
+                                                                        u : -0.48524E-04  0.10000E+01
+                                                                        v : -0.14008E+00  0.13429E+00
+                                                                        w : -0.11414E-01  0.52596E+00
+                                                                        =================================
+                 Calculate pressures    779.7982      9.8888      7.0361
+                                                                        raw    pressures : -0.89737E+00  0.68898E+00
+                 Smoothing pressures    779.8745      9.9652      0.0763
+                                                                        smooth pressures : -0.12870E+00  0.12999E-02
+                do leaf measurements    779.8924      9.9831      0.0179
+       compute convergence criterion    779.9327     10.0233      0.0403
+                                                                        velocity_diff_norm = 0.0058029 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    779.9354     10.0261      0.0027
+Compute isostasy and adjust vertical    779.9356     10.0262      0.0002
+                  Compute divergence    779.9356     10.0263      0.0000
+                        wsmp_cleanup    779.9653     10.0559      0.0297
+              Reset surface geometry    779.9705     10.0611      0.0052
+ -----------------------------------------------------------------------
+           Temperature calculations     779.9716     10.0622      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    780.1364     10.2270      0.1648
+                   Advect surface  1    780.1365     10.2272      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    780.1807     10.2713      0.0442
+                 Update cloud fields    797.4067     27.4974     17.2260
+                    Advect the cloud    798.1765     28.2671      0.7697
+                        Write output    802.3895     32.4801      4.2130
+                    End of time step    804.1444     34.2350      1.7549
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      23    804.1445
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    804.1446      0.0001      0.0001
+                          surface 01    804.1446      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    804.1447      0.0002      0.0001
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    804.1512      0.0067      0.0065
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    804.1554
+ ----------------------------------------------------------------------- 
+                 Create octree solve    804.1556      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    804.1716      0.0162      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    804.1899      0.0345      0.0184
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.049019E+01
+                                                                        osolve%nleaves=   54587
+    Smoothen octree after refinement    804.2116      0.0562      0.0217
+                                                                        osolve%nleaves=   54587
+                   Build osolve icon    804.3471      0.1917      0.1355
+             Imbed surface in osolve    804.5383      0.3828      0.1911
+                embedding surface  1    804.5385      0.3831      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    807.1377      2.9823      2.5992
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54496 leaves
+                                                                        nleaves_new=    54587 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    807.1400      2.9845      0.0022
+        Interpolate velo onto osolve    809.5285      5.3731      2.3886
+                                                                        ov%pressure:     -0.1376E+00 0.7395E-02
+                                                                        osolve%pressure: -0.1287E+00 0.1300E-02
+                      Find bad faces    809.8064      5.6509      0.2779
+                           Find void    810.0237      5.8682      0.2173
+                                                                        whole leaf in fluid    24738
+          Define boundary conditions    810.0365      5.8811      0.0128
+                         wsmp setup1    810.0418      5.8864      0.0053
+                                                                        n =   102444
+                                                                        nz=  3782766
+                         wsmp setup2    810.7813      6.6258      0.7395
+ -----------------------------------
+ start of non-linear iteratio      1    810.8149
+ -----------------------------------
+                        build system    810.8151      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3679        4142
+                                                                        viscosity range   1.55919E-05  1.30880E+03
+                                                                        viscosity capped in     28953
+                          wsmp solve    813.6950      2.8800      2.8799
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0970750 s
+                                                                        wsmp: ordering time:               2.5262070 s
+                                                                        wsmp: symbolic fact. time:         0.5073020 s
+                                                                        wsmp: Cholesky fact. time:         4.1235340 s
+                                                                        wsmp: Factorisation            14664.1480719 Mflops
+                                                                        wsmp: back substitution time:      0.0815229 s
+                                                                        wsmp: from b to rhs vector:        0.0054460 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.3414631 s
+                                                                        =================================
+                                                                        u : -0.87301E-04  0.10000E+01
+                                                                        v : -0.13856E+00  0.13424E+00
+                                                                        w : -0.12690E-01  0.53177E+00
+                                                                        =================================
+                 Calculate pressures    821.0505     10.2356      7.3556
+                                                                        raw    pressures : -0.91974E+00  0.70703E+00
+                 Smoothing pressures    821.1275     10.3125      0.0769
+                                                                        smooth pressures : -0.12888E+00  0.15522E-02
+                do leaf measurements    821.1399     10.3250      0.0124
+       compute convergence criterion    821.1801     10.3652      0.0402
+                                                                        velocity_diff_norm = 0.0079786 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    821.1829     10.3680      0.0028
+Compute isostasy and adjust vertical    821.1831     10.3682      0.0002
+                  Compute divergence    821.1832     10.3682      0.0000
+                        wsmp_cleanup    821.2132     10.3983      0.0301
+              Reset surface geometry    821.2183     10.4033      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations     821.2193     10.4043      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    821.3843     10.5693      0.1650
+                   Advect surface  1    821.3844     10.5695      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    821.4287     10.6137      0.0442
+                 Update cloud fields    838.7706     27.9556     17.3419
+                    Advect the cloud    839.5540     28.7391      0.7834
+                        Write output    843.7910     32.9761      4.2370
+                    End of time step    845.5710     34.7561      1.7800
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      24    845.5713
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    845.5713      0.0001      0.0001
+                          surface 01    845.5714      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    845.5714      0.0002      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    845.5770      0.0057      0.0055
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    845.5812
+ ----------------------------------------------------------------------- 
+                 Create octree solve    845.5813      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    845.5972      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    845.6155      0.0344      0.0184
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.055746E+01
+                                                                        osolve%nleaves=   54622
+    Smoothen octree after refinement    845.6368      0.0556      0.0213
+                                                                        osolve%nleaves=   54622
+                   Build osolve icon    845.7709      0.1898      0.1341
+             Imbed surface in osolve    845.9553      0.3741      0.1843
+                embedding surface  1    845.9555      0.3743      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    848.5482      2.9671      2.5927
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54587 leaves
+                                                                        nleaves_new=    54622 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    848.5504      2.9693      0.0022
+        Interpolate velo onto osolve    850.7923      5.2112      2.2419
+                                                                        ov%pressure:     -0.1377E+00 0.8507E-02
+                                                                        osolve%pressure: -0.1289E+00 0.1552E-02
+                      Find bad faces    851.0893      5.5082      0.2970
+                           Find void    851.3066      5.7255      0.2173
+                                                                        whole leaf in fluid    24781
+          Define boundary conditions    851.3194      5.7383      0.0128
+                         wsmp setup1    851.3247      5.7436      0.0053
+                                                                        n =   102903
+                                                                        nz=  3801126
+                         wsmp setup2    852.0608      6.4796      0.7361
+ -----------------------------------
+ start of non-linear iteratio      1    852.0940
+ -----------------------------------
+                        build system    852.0942      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3696        4160
+                                                                        viscosity range   1.57325E-05  1.30880E+03
+                                                                        viscosity capped in     29069
+                          wsmp solve    855.0549      2.9609      2.9607
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0938649 s
+                                                                        wsmp: ordering time:               2.6364760 s
+                                                                        wsmp: symbolic fact. time:         0.5152869 s
+                                                                        wsmp: Cholesky fact. time:         4.1968420 s
+                                                                        wsmp: Factorisation            14741.1497875 Mflops
+                                                                        wsmp: back substitution time:      0.0823469 s
+                                                                        wsmp: from b to rhs vector:        0.0056560 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.5308700 s
+                                                                        =================================
+                                                                        u : -0.18641E-03  0.10001E+01
+                                                                        v : -0.31056E+00  0.13631E+00
+                                                                        w : -0.13930E-01  0.61753E+00
+                                                                        =================================
+                 Calculate pressures    862.5999     10.5059      7.5450
+                                                                        raw    pressures : -0.95656E+00  0.74841E+00
+                 Smoothing pressures    862.6776     10.5836      0.0776
+                                                                        smooth pressures : -0.12903E+00  0.16688E-02
+                do leaf measurements    862.6901     10.5961      0.0125
+       compute convergence criterion    862.7315     10.6375      0.0414
+                                                                        velocity_diff_norm = 0.0058599 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    862.7343     10.6403      0.0028
+Compute isostasy and adjust vertical    862.7345     10.6405      0.0002
+                  Compute divergence    862.7345     10.6405      0.0000
+                        wsmp_cleanup    862.7639     10.6699      0.0294
+              Reset surface geometry    862.7690     10.6750      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations     862.7700     10.6760      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    862.9354     10.8414      0.1653
+                   Advect surface  1    862.9355     10.8415      0.0002
+                                                                        S. 1:    4 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    862.9796     10.8856      0.0441
+                 Update cloud fields    881.3039     29.2099     18.3243
+                    Advect the cloud    882.0790     29.9850      0.7751
+                        Write output    886.3283     34.2343      4.2492
+                    End of time step    888.0732     35.9792      1.7449
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      25    888.0733
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    888.0734      0.0001      0.0001
+                          surface 01    888.0734      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    888.0735      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    888.0802      0.0069      0.0068
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    888.0833
+ ----------------------------------------------------------------------- 
+                 Create octree solve    888.0835      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    888.0992      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    888.1174      0.0341      0.0182
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.062182E+01
+                                                                        osolve%nleaves=   54755
+    Smoothen octree after refinement    888.1390      0.0557      0.0215
+                                                                        osolve%nleaves=   54755
+                   Build osolve icon    888.2736      0.1903      0.1346
+             Imbed surface in osolve    888.4631      0.3798      0.1895
+                embedding surface  1    888.4633      0.3801      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    891.0696      2.9863      2.6063
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54622 leaves
+                                                                        nleaves_new=    54755 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    891.0718      2.9885      0.0022
+        Interpolate velo onto osolve    893.3054      5.2221      2.2336
+                                                                        ov%pressure:     -0.1379E+00 0.1313E-01
+                                                                        osolve%pressure: -0.1290E+00 0.1669E-02
+                      Find bad faces    893.5907      5.5074      0.2853
+                           Find void    893.8080      5.7247      0.2173
+                                                                        whole leaf in fluid    24836
+          Define boundary conditions    893.8206      5.7373      0.0126
+                         wsmp setup1    893.8258      5.7425      0.0052
+                                                                        n =   103098
+                                                                        nz=  3810903
+                         wsmp setup2    894.5686      6.4853      0.7428
+ -----------------------------------
+ start of non-linear iteratio      1    894.5996
+ -----------------------------------
+                        build system    894.5998      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3702        4169
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29131
+                          wsmp solve    897.5516      2.9519      2.9518
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0944951 s
+                                                                        wsmp: ordering time:               2.5877302 s
+                                                                        wsmp: symbolic fact. time:         0.5145919 s
+                                                                        wsmp: Cholesky fact. time:         4.1559970 s
+                                                                        wsmp: Factorisation            14820.3742528 Mflops
+                                                                        wsmp: back substitution time:      0.0820019 s
+                                                                        wsmp: from b to rhs vector:        0.0055420 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.4407630 s
+                                                                        =================================
+                                                                        u : -0.27193E-03  0.10004E+01
+                                                                        v : -0.13541E+00  0.13233E+00
+                                                                        w : -0.13883E-01  0.52902E+00
+                                                                        =================================
+                 Calculate pressures    905.0059     10.4063      7.4543
+                                                                        raw    pressures : -0.97895E+00  0.76811E+00
+                 Smoothing pressures    905.0855     10.4859      0.0796
+                                                                        smooth pressures : -0.12899E+00  0.18574E-02
+                do leaf measurements    905.0979     10.4982      0.0124
+       compute convergence criterion    905.1393     10.5397      0.0415
+                                                                        velocity_diff_norm = 0.0056222 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    905.1423     10.5426      0.0029
+Compute isostasy and adjust vertical    905.1424     10.5428      0.0002
+                  Compute divergence    905.1425     10.5429      0.0000
+                        wsmp_cleanup    905.1720     10.5723      0.0295
+              Reset surface geometry    905.1770     10.5774      0.0050
+ -----------------------------------------------------------------------
+           Temperature calculations     905.1781     10.5784      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    905.3438     10.7442      0.1658
+                   Advect surface  1    905.3440     10.7444      0.0002
+                                                                        S. 1:    4 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    905.3884     10.7888      0.0444
+                 Update cloud fields    923.3833     28.7836     17.9949
+                    Advect the cloud    924.1596     29.5599      0.7763
+                        Write output    928.4043     33.8047      4.2448
+                    End of time step    930.1597     35.5601      1.7554
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      26    930.1599
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    930.1600      0.0001      0.0001
+                          surface 01    930.1600      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    930.1600      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    930.1664      0.0065      0.0063
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    930.1694
+ ----------------------------------------------------------------------- 
+                 Create octree solve    930.1696      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    930.1854      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    930.2037      0.0343      0.0183
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.068374E+01
+                                                                        osolve%nleaves=   54874
+    Smoothen octree after refinement    930.2252      0.0558      0.0215
+                                                                        osolve%nleaves=   54874
+                   Build osolve icon    930.3601      0.1907      0.1349
+             Imbed surface in osolve    930.5497      0.3803      0.1896
+                embedding surface  1    930.5499      0.3805      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    933.1667      2.9973      2.6168
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54755 leaves
+                                                                        nleaves_new=    54874 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    933.1690      2.9996      0.0022
+        Interpolate velo onto osolve    935.4066      5.2372      2.2376
+                                                                        ov%pressure:     -0.1377E+00 0.1166E-01
+                                                                        osolve%pressure: -0.1290E+00 0.1857E-02
+                      Find bad faces    935.6750      5.5056      0.2684
+                           Find void    935.8923      5.7229      0.2173
+                                                                        whole leaf in fluid    24881
+          Define boundary conditions    935.9050      5.7356      0.0127
+                         wsmp setup1    935.9101      5.7407      0.0051
+                                                                        n =   103245
+                                                                        nz=  3818451
+                         wsmp setup2    936.6583      6.4889      0.7482
+ -----------------------------------
+ start of non-linear iteratio      1    936.6903
+ -----------------------------------
+                        build system    936.6905      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3708        4176
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29180
+                          wsmp solve    939.6534      2.9630      2.9629
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0889661 s
+                                                                        wsmp: ordering time:               2.5114930 s
+                                                                        wsmp: symbolic fact. time:         0.5129402 s
+                                                                        wsmp: Cholesky fact. time:         3.8328819 s
+                                                                        wsmp: Factorisation            15893.9700519 Mflops
+                                                                        wsmp: back substitution time:      0.0745361 s
+                                                                        wsmp: from b to rhs vector:        0.0056670 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.0268240 s
+                                                                        =================================
+                                                                        u : -0.31850E-03  0.10010E+01
+                                                                        v : -0.13447E+00  0.13257E+00
+                                                                        w : -0.14663E-01  0.52529E+00
+                                                                        =================================
+                 Calculate pressures    946.6939     10.0036      7.0405
+                                                                        raw    pressures : -0.99610E+00  0.79182E+00
+                 Smoothing pressures    946.7722     10.0818      0.0782
+                                                                        smooth pressures : -0.12915E+00  0.20237E-02
+                do leaf measurements    946.7876     10.0972      0.0154
+       compute convergence criterion    946.8288     10.1384      0.0412
+                                                                        velocity_diff_norm = 0.0061200 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    946.8316     10.1412      0.0028
+Compute isostasy and adjust vertical    946.8318     10.1414      0.0002
+                  Compute divergence    946.8318     10.1415      0.0000
+                        wsmp_cleanup    946.8616     10.1713      0.0298
+              Reset surface geometry    946.8670     10.1767      0.0054
+ -----------------------------------------------------------------------
+           Temperature calculations     946.8680     10.1777      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    947.0328     10.3425      0.1648
+                   Advect surface  1    947.0330     10.3426      0.0002
+                                                                        S. 1:    5 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    947.0772     10.3869      0.0443
+                 Update cloud fields    964.3846     27.6942     17.3073
+                    Advect the cloud    965.1705     28.4802      0.7859
+                        Write output    969.4308     32.7404      4.2603
+                    End of time step    971.1881     34.4978      1.7573
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      27    971.1884
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces    971.1884      0.0001      0.0001
+                          surface 01    971.1885      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface    971.1885      0.0002      0.0001
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay    971.1945      0.0061      0.0060
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1    971.1975
+ ----------------------------------------------------------------------- 
+                 Create octree solve    971.1977      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov    971.2135      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion    971.2312      0.0337      0.0177
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.074196E+01
+                                                                        osolve%nleaves=   54951
+    Smoothen octree after refinement    971.2527      0.0552      0.0215
+                                                                        osolve%nleaves=   54951
+                   Build osolve icon    971.3878      0.1903      0.1351
+             Imbed surface in osolve    971.5772      0.3797      0.1894
+                embedding surface  1    971.5774      0.3799      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability    974.1633      2.9657      2.5858
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54874 leaves
+                                                                        nleaves_new=    54951 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud    974.1655      2.9679      0.0022
+        Interpolate velo onto osolve    976.4395      5.2420      2.2740
+                                                                        ov%pressure:     -0.1379E+00 0.1214E-01
+                                                                        osolve%pressure: -0.1292E+00 0.2024E-02
+                      Find bad faces    976.7819      5.5844      0.3424
+                           Find void    976.9993      5.8018      0.2174
+                                                                        whole leaf in fluid    24912
+          Define boundary conditions    977.0125      5.8150      0.0132
+                         wsmp setup1    977.0188      5.8212      0.0063
+                                                                        n =   103377
+                                                                        nz=  3825897
+                         wsmp setup2    977.7639      6.5664      0.7452
+ -----------------------------------
+ start of non-linear iteratio      1    977.7994
+ -----------------------------------
+                        build system    977.7996      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3712        4183
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29222
+                          wsmp solve    980.7481      2.9487      2.9485
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0923102 s
+                                                                        wsmp: ordering time:               2.7005229 s
+                                                                        wsmp: symbolic fact. time:         0.5153160 s
+                                                                        wsmp: Cholesky fact. time:         4.5112300 s
+                                                                        wsmp: Factorisation            13850.2066060 Mflops
+                                                                        wsmp: back substitution time:      0.0862119 s
+                                                                        wsmp: from b to rhs vector:        0.0055830 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.9115460 s
+                                                                        =================================
+                                                                        u : -0.33602E-03  0.10014E+01
+                                                                        v : -0.13241E+00  0.13195E+00
+                                                                        w : -0.15450E-01  0.52654E+00
+                                                                        =================================
+                 Calculate pressures    988.6735     10.8741      7.9254
+                                                                        raw    pressures : -0.10193E+01  0.81613E+00
+                 Smoothing pressures    988.7511     10.9517      0.0776
+                                                                        smooth pressures : -0.12929E+00  0.21467E-02
+                do leaf measurements    988.7633     10.9640      0.0123
+       compute convergence criterion    988.8035     11.0041      0.0402
+                                                                        velocity_diff_norm = 0.0055482 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube    988.8065     11.0071      0.0029
+Compute isostasy and adjust vertical    988.8067     11.0073      0.0002
+                  Compute divergence    988.8067     11.0073      0.0000
+                        wsmp_cleanup    988.8365     11.0371      0.0298
+              Reset surface geometry    988.8415     11.0422      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations     988.8426     11.0432      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces    989.0049     11.2056      0.1623
+                   Advect surface  1    989.0051     11.2057      0.0002
+                                                                        S. 1:    5 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1    989.0495     11.2501      0.0444
+                 Update cloud fields   1007.6192     29.8198     18.5697
+                    Advect the cloud   1008.4019     30.6025      0.7827
+                        Write output   1012.6695     34.8701      4.2676
+                    End of time step   1014.4482     36.6488      1.7787
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      28   1014.4484
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1014.4485      0.0001      0.0001
+                          surface 01   1014.4485      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface   1014.4486      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1014.4555      0.0070      0.0069
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1014.4585
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1014.4587      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1014.4745      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1014.4929      0.0344      0.0184
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.080282E+01
+                                                                        osolve%nleaves=   55049
+    Smoothen octree after refinement   1014.5144      0.0559      0.0215
+                                                                        osolve%nleaves=   55049
+                   Build osolve icon   1014.6499      0.1914      0.1355
+             Imbed surface in osolve   1014.8430      0.3845      0.1931
+                embedding surface  1   1014.8432      0.3847      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1017.4750      3.0165      2.6318
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    54951 leaves
+                                                                        nleaves_new=    55049 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1017.4773      3.0188      0.0022
+        Interpolate velo onto osolve   1019.7671      5.3086      2.2898
+                                                                        ov%pressure:     -0.1379E+00 0.1184E-01
+                                                                        osolve%pressure: -0.1293E+00 0.2147E-02
+                      Find bad faces   1019.9973      5.5389      0.2302
+                           Find void   1020.2146      5.7562      0.2173
+                                                                        whole leaf in fluid    24954
+          Define boundary conditions   1020.2276      5.7691      0.0129
+                         wsmp setup1   1020.2328      5.7743      0.0052
+                                                                        n =   103488
+                                                                        nz=  3831960
+                         wsmp setup2   1020.9817      6.5232      0.7490
+ -----------------------------------
+ start of non-linear iteratio      1   1021.0144
+ -----------------------------------
+                        build system   1021.0146      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3717        4188
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29259
+                          wsmp solve   1023.9782      2.9638      2.9636
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0940750 s
+                                                                        wsmp: ordering time:               2.5208681 s
+                                                                        wsmp: symbolic fact. time:         0.5125270 s
+                                                                        wsmp: Cholesky fact. time:         3.8894970 s
+                                                                        wsmp: Factorisation            15915.4675091 Mflops
+                                                                        wsmp: back substitution time:      0.0751019 s
+                                                                        wsmp: from b to rhs vector:        0.0056100 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.0980289 s
+                                                                        =================================
+                                                                        u : -0.33589E-03  0.10014E+01
+                                                                        v : -0.13688E+00  0.18349E+00
+                                                                        w : -0.15933E-01  0.52122E+00
+                                                                        =================================
+                 Calculate pressures   1031.0902     10.0758      7.1120
+                                                                        raw    pressures : -0.10395E+01  0.83535E+00
+                 Smoothing pressures   1031.1685     10.1541      0.0783
+                                                                        smooth pressures : -0.12934E+00  0.21519E-02
+                do leaf measurements   1031.1831     10.1687      0.0146
+       compute convergence criterion   1031.2219     10.2075      0.0387
+                                                                        velocity_diff_norm = 0.0087143 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1031.2247     10.2103      0.0029
+Compute isostasy and adjust vertical   1031.2249     10.2105      0.0002
+                  Compute divergence   1031.2250     10.2106      0.0000
+                        wsmp_cleanup   1031.2543     10.2399      0.0293
+              Reset surface geometry   1031.2598     10.2454      0.0055
+ -----------------------------------------------------------------------
+           Temperature calculations    1031.2608     10.2464      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1031.4246     10.4102      0.1638
+                   Advect surface  1   1031.4247     10.4103      0.0002
+                                                                        S. 1:    4 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1031.4689     10.4545      0.0442
+                 Update cloud fields   1049.1957     28.1813     17.7268
+                    Advect the cloud   1049.9792     28.9647      0.7835
+                        Write output   1054.2468     33.2324      4.2676
+                    End of time step   1056.5334     35.5190      2.2866
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      29   1056.5336
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1056.5337      0.0001      0.0001
+                          surface 01   1056.5337      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface   1056.5337      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1056.5407      0.0071      0.0070
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1056.5438
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1056.5439      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1056.5597      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1056.5777      0.0339      0.0179
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.086066E+01
+                                                                        osolve%nleaves=   55154
+    Smoothen octree after refinement   1056.5993      0.0555      0.0216
+                                                                        osolve%nleaves=   55154
+                   Build osolve icon   1056.7351      0.1913      0.1358
+             Imbed surface in osolve   1056.9264      0.3826      0.1913
+                embedding surface  1   1056.9266      0.3828      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1059.5585      3.0148      2.6320
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55049 leaves
+                                                                        nleaves_new=    55154 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1059.5608      3.0170      0.0023
+        Interpolate velo onto osolve   1061.8677      5.3239      2.3069
+                                                                        ov%pressure:     -0.1379E+00 0.1241E-01
+                                                                        osolve%pressure: -0.1293E+00 0.2152E-02
+                      Find bad faces   1062.0513      5.5076      0.1837
+                           Find void   1062.2686      5.7249      0.2173
+                                                                        whole leaf in fluid    24979
+          Define boundary conditions   1062.2818      5.7380      0.0132
+                         wsmp setup1   1062.2875      5.7437      0.0057
+                                                                        n =   103572
+                                                                        nz=  3835701
+                         wsmp setup2   1063.0364      6.4926      0.7489
+ -----------------------------------
+ start of non-linear iteratio      1   1063.0688
+ -----------------------------------
+                        build system   1063.0690      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3720        4193
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29286
+                          wsmp solve   1066.0381      2.9692      2.9691
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0939021 s
+                                                                        wsmp: ordering time:               2.5703552 s
+                                                                        wsmp: symbolic fact. time:         0.5157051 s
+                                                                        wsmp: Cholesky fact. time:         3.9013648 s
+                                                                        wsmp: Factorisation            15714.3598020 Mflops
+                                                                        wsmp: back substitution time:      0.0748510 s
+                                                                        wsmp: from b to rhs vector:        0.0057700 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.1623049 s
+                                                                        =================================
+                                                                        u : -0.30830E-03  0.10012E+01
+                                                                        v : -0.15537E+00  0.17451E+00
+                                                                        w : -0.16958E-01  0.51797E+00
+                                                                        =================================
+                 Calculate pressures   1073.2141     10.1453      7.1761
+                                                                        raw    pressures : -0.10479E+01  0.84440E+00
+                 Smoothing pressures   1073.2930     10.2241      0.0788
+                                                                        smooth pressures : -0.12946E+00  0.22953E-02
+                do leaf measurements   1073.3059     10.2370      0.0129
+       compute convergence criterion   1073.3464     10.2775      0.0405
+                                                                        velocity_diff_norm = 0.0041218 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1073.3494     10.2805      0.0030
+Compute isostasy and adjust vertical   1073.3496     10.2807      0.0002
+                  Compute divergence   1073.3496     10.2808      0.0000
+                        wsmp_cleanup   1073.3788     10.3100      0.0292
+              Reset surface geometry   1073.3840     10.3152      0.0052
+ -----------------------------------------------------------------------
+           Temperature calculations    1073.3851     10.3162      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1073.5492     10.4804      0.1641
+                   Advect surface  1   1073.5494     10.4805      0.0002
+                                                                        S. 1:    5 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1073.5937     10.5248      0.0443
+                 Update cloud fields   1092.2181     29.1492     18.6244
+                    Advect the cloud   1092.9977     29.9289      0.7796
+                        Write output   1097.3438     34.2750      4.3461
+                    End of time step   1099.1295     36.0606      1.7857
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      30   1099.1297
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1099.1297      0.0001      0.0001
+                          surface 01   1099.1298      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface   1099.1298      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1099.1363      0.0066      0.0065
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4356points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1099.1403
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1099.1405      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1099.1563      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1099.1741      0.0338      0.0178
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.091828E+01
+                                                                        osolve%nleaves=   55182
+    Smoothen octree after refinement   1099.1956      0.0553      0.0215
+                                                                        osolve%nleaves=   55182
+                   Build osolve icon   1099.3315      0.1912      0.1359
+             Imbed surface in osolve   1099.5464      0.4061      0.2149
+                embedding surface  1   1099.5466      0.4063      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1102.1573      3.0170      2.6108
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55154 leaves
+                                                                        nleaves_new=    55182 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1102.1595      3.0192      0.0022
+        Interpolate velo onto osolve   1104.4527      5.3124      2.2932
+                                                                        ov%pressure:     -0.1380E+00 0.1276E-01
+                                                                        osolve%pressure: -0.1295E+00 0.2295E-02
+                      Find bad faces   1104.6621      5.5218      0.2094
+                           Find void   1104.8794      5.7391      0.2173
+                                                                        whole leaf in fluid    24995
+          Define boundary conditions   1104.8924      5.7521      0.0129
+                         wsmp setup1   1104.8978      5.7575      0.0054
+                                                                        n =   103659
+                                                                        nz=  3840060
+                         wsmp setup2   1105.6434      6.5031      0.7457
+ -----------------------------------
+ start of non-linear iteratio      1   1105.6785
+ -----------------------------------
+                        build system   1105.6787      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3723        4198
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29317
+                          wsmp solve   1108.6557      2.9772      2.9770
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0869372 s
+                                                                        wsmp: ordering time:               2.6043441 s
+                                                                        wsmp: symbolic fact. time:         0.5165439 s
+                                                                        wsmp: Cholesky fact. time:         4.1940742 s
+                                                                        wsmp: Factorisation            14856.4113720 Mflops
+                                                                        wsmp: back substitution time:      0.0772722 s
+                                                                        wsmp: from b to rhs vector:        0.0057631 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.4852622 s
+                                                                        =================================
+                                                                        u : -0.27764E-03  0.10013E+01
+                                                                        v : -0.12748E+00  0.17243E+00
+                                                                        w : -0.17717E-01  0.51363E+00
+                                                                        =================================
+                 Calculate pressures   1116.1551     10.4766      7.4994
+                                                                        raw    pressures : -0.10743E+01  0.86780E+00
+                 Smoothing pressures   1116.2331     10.5546      0.0780
+                                                                        smooth pressures : -0.12959E+00  0.24653E-02
+                do leaf measurements   1116.2476     10.5691      0.0145
+       compute convergence criterion   1116.2870     10.6085      0.0394
+                                                                        velocity_diff_norm = 0.0064208 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1116.2899     10.6114      0.0029
+Compute isostasy and adjust vertical   1116.2901     10.6116      0.0002
+                  Compute divergence   1116.2901     10.6116      0.0000
+                        wsmp_cleanup   1116.3200     10.6415      0.0299
+              Reset surface geometry   1116.3254     10.6469      0.0054
+ -----------------------------------------------------------------------
+           Temperature calculations    1116.3264     10.6479      0.0010
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1116.4922     10.8137      0.1658
+                   Advect surface  1   1116.4924     10.8139      0.0002
+                                                                        S. 1:    5 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1116.5371     10.8586      0.0448
+                 Update cloud fields   1135.1439     29.4654     18.6067
+                    Advect the cloud   1135.9483     30.2698      0.8045
+                        Write output   1140.2598     34.5812      4.3114
+                    End of time step   1142.0499     36.3714      1.7902
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      31   1142.0501
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1142.0502      0.0001      0.0001
+                          surface 01   1142.0502      0.0001      0.0000
+                                                                        S. 1:     4356points
+                      refine surface   1142.0502      0.0001      0.0000
+                                                                        S. 1:   1 added ptcls in refine_surface
+                      check delaunay   1142.0567      0.0066      0.0065
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                      refine surface   1142.0605      0.0104      0.0038
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1142.0640      0.0139      0.0036
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4357points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1142.0672
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1142.0673      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1142.0831      0.0159      0.0157
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1142.1010      0.0338      0.0179
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.097005E+01
+                                                                        osolve%nleaves=   55196
+    Smoothen octree after refinement   1142.1225      0.0554      0.0216
+                                                                        osolve%nleaves=   55196
+                   Build osolve icon   1142.2584      0.1913      0.1359
+             Imbed surface in osolve   1142.4718      0.4046      0.2133
+                embedding surface  1   1142.4719      0.4048      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1145.1077      3.0405      2.6357
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55182 leaves
+                                                                        nleaves_new=    55196 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1145.1099      3.0427      0.0022
+        Interpolate velo onto osolve   1147.4139      5.3468      2.3041
+                                                                        ov%pressure:     -0.1380E+00 0.1262E-01
+                                                                        osolve%pressure: -0.1296E+00 0.2465E-02
+                      Find bad faces   1147.6473      5.5802      0.2334
+                           Find void   1147.8647      5.7975      0.2174
+                                                                        whole leaf in fluid    25022
+          Define boundary conditions   1147.8776      5.8104      0.0129
+                         wsmp setup1   1147.8827      5.8155      0.0051
+                                                                        n =   103725
+                                                                        nz=  3843063
+                         wsmp setup2   1148.6293      6.5621      0.7466
+ -----------------------------------
+ start of non-linear iteratio      1   1148.6647
+ -----------------------------------
+                        build system   1148.6649      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3725        4201
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29338
+                          wsmp solve   1151.6155      2.9508      2.9506
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0940261 s
+                                                                        wsmp: ordering time:               2.5650909 s
+                                                                        wsmp: symbolic fact. time:         0.5179472 s
+                                                                        wsmp: Cholesky fact. time:         4.0105300 s
+                                                                        wsmp: Factorisation            15580.5382724 Mflops
+                                                                        wsmp: back substitution time:      0.0726080 s
+                                                                        wsmp: from b to rhs vector:        0.0057750 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.2663400 s
+                                                                        =================================
+                                                                        u : -0.17872E-03  0.10013E+01
+                                                                        v : -0.21225E+00  0.16580E+00
+                                                                        w : -0.19314E-01  0.51231E+00
+                                                                        =================================
+                 Calculate pressures   1158.8958     10.2311      7.2803
+                                                                        raw    pressures : -0.10907E+01  0.88641E+00
+                 Smoothing pressures   1158.9745     10.3098      0.0787
+                                                                        smooth pressures : -0.12968E+00  0.24838E-02
+                do leaf measurements   1158.9887     10.3240      0.0142
+       compute convergence criterion   1159.0285     10.3638      0.0399
+                                                                        velocity_diff_norm = 0.0171027 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1159.0314     10.3667      0.0029
+Compute isostasy and adjust vertical   1159.0316     10.3669      0.0002
+                  Compute divergence   1159.0316     10.3670      0.0000
+                        wsmp_cleanup   1159.0620     10.3974      0.0304
+              Reset surface geometry   1159.0671     10.4025      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations    1159.0682     10.4035      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1159.2332     10.5685      0.1650
+                   Advect surface  1   1159.2333     10.5686      0.0002
+                                                                        S. 1:    4 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1159.2779     10.6132      0.0446
+                 Update cloud fields   1178.5115     29.8469     19.2337
+                    Advect the cloud   1179.2930     30.6283      0.7815
+                        Write output   1183.6262     34.9615      4.3332
+                    End of time step   1185.4227     36.7580      1.7965
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      32   1185.4229
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1185.4230      0.0001      0.0001
+                          surface 01   1185.4230      0.0001      0.0000
+                                                                        S. 1:     4357points
+                      refine surface   1185.4231      0.0002      0.0001
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1185.4298      0.0069      0.0068
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4357points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1185.4329
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1185.4331      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1185.4489      0.0159      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1185.4664      0.0335      0.0176
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.102160E+01
+                                                                        osolve%nleaves=   55217
+    Smoothen octree after refinement   1185.4878      0.0549      0.0214
+                                                                        osolve%nleaves=   55217
+                   Build osolve icon   1185.6238      0.1909      0.1360
+             Imbed surface in osolve   1185.8237      0.3907      0.1998
+                embedding surface  1   1185.8238      0.3909      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1188.4207      2.9878      2.5969
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55196 leaves
+                                                                        nleaves_new=    55217 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1188.4230      2.9900      0.0023
+        Interpolate velo onto osolve   1190.7496      5.3167      2.3267
+                                                                        ov%pressure:     -0.1384E+00 0.1338E-01
+                                                                        osolve%pressure: -0.1297E+00 0.2484E-02
+                      Find bad faces   1191.0022      5.5692      0.2525
+                           Find void   1191.2194      5.7865      0.2173
+                                                                        whole leaf in fluid    25044
+          Define boundary conditions   1191.2320      5.7991      0.0126
+                         wsmp setup1   1191.2366      5.8037      0.0046
+                                                                        n =   103812
+                                                                        nz=  3847134
+                         wsmp setup2   1191.9884      6.5555      0.7518
+ -----------------------------------
+ start of non-linear iteratio      1   1192.0210
+ -----------------------------------
+                        build system   1192.0212      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3729        4206
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29367
+                          wsmp solve   1195.0068      2.9858      2.9857
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0884559 s
+                                                                        wsmp: ordering time:               2.5337071 s
+                                                                        wsmp: symbolic fact. time:         0.5225141 s
+                                                                        wsmp: Cholesky fact. time:         4.2182539 s
+                                                                        wsmp: Factorisation            14895.3451915 Mflops
+                                                                        wsmp: back substitution time:      0.0772550 s
+                                                                        wsmp: from b to rhs vector:        0.0057070 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.4462550 s
+                                                                        =================================
+                                                                        u : -0.14701E-03  0.10014E+01
+                                                                        v : -0.20775E+00  0.14845E+00
+                                                                        w : -0.19857E-01  0.51074E+00
+                                                                        =================================
+                 Calculate pressures   1202.4669     10.4459      7.4601
+                                                                        raw    pressures : -0.10969E+01  0.89242E+00
+                 Smoothing pressures   1202.5451     10.5241      0.0782
+                                                                        smooth pressures : -0.12972E+00  0.24774E-02
+                do leaf measurements   1202.5617     10.5407      0.0166
+       compute convergence criterion   1202.6007     10.5797      0.0390
+                                                                        velocity_diff_norm = 0.0124324 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1202.6036     10.5826      0.0029
+Compute isostasy and adjust vertical   1202.6038     10.5828      0.0002
+                  Compute divergence   1202.6038     10.5828      0.0000
+                        wsmp_cleanup   1202.6338     10.6128      0.0300
+              Reset surface geometry   1202.6390     10.6180      0.0052
+ -----------------------------------------------------------------------
+           Temperature calculations    1202.6402     10.6192      0.0012
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1202.8059     10.7849      0.1658
+                   Advect surface  1   1202.8061     10.7851      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1202.8508     10.8298      0.0447
+                 Update cloud fields   1222.6934     30.6724     19.8427
+                    Advect the cloud   1223.4850     31.4640      0.7915
+                        Write output   1227.8251     35.8041      4.3401
+                    End of time step   1229.6465     37.6255      1.8215
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      33   1229.6467
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1229.6468      0.0001      0.0001
+                          surface 01   1229.6468      0.0001      0.0000
+                                                                        S. 1:     4357points
+                      refine surface   1229.6469      0.0002      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1229.6546      0.0079      0.0078
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4357points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1229.6587
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1229.6588      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1229.6747      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1229.6920      0.0334      0.0174
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.107040E+01
+                                                                        osolve%nleaves=   55245
+    Smoothen octree after refinement   1229.7136      0.0550      0.0216
+                                                                        osolve%nleaves=   55245
+                   Build osolve icon   1229.8497      0.1910      0.1361
+             Imbed surface in osolve   1230.0527      0.3940      0.2030
+                embedding surface  1   1230.0529      0.3943      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1232.6882      3.0295      2.6353
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55217 leaves
+                                                                        nleaves_new=    55245 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1232.6906      3.0320      0.0024
+        Interpolate velo onto osolve   1235.0299      5.3712      2.3393
+                                                                        ov%pressure:     -0.1386E+00 0.1370E-01
+                                                                        osolve%pressure: -0.1297E+00 0.2477E-02
+                      Find bad faces   1235.2244      5.5657      0.1945
+                           Find void   1235.4439      5.7852      0.2195
+                                                                        whole leaf in fluid    25064
+          Define boundary conditions   1235.4578      5.7991      0.0139
+                         wsmp setup1   1235.4641      5.8054      0.0063
+                                                                        n =   103899
+                                                                        nz=  3851637
+                         wsmp setup2   1236.2264      6.5677      0.7623
+ -----------------------------------
+ start of non-linear iteratio      1   1236.2599
+ -----------------------------------
+                        build system   1236.2602      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3731        4209
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29395
+                          wsmp solve   1239.2795      3.0196      3.0193
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.1002500 s
+                                                                        wsmp: ordering time:               2.5710790 s
+                                                                        wsmp: symbolic fact. time:         0.5226541 s
+                                                                        wsmp: Cholesky fact. time:         3.7609160 s
+                                                                        wsmp: Factorisation            16694.3238941 Mflops
+                                                                        wsmp: back substitution time:      0.0747890 s
+                                                                        wsmp: from b to rhs vector:        0.0056911 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.0357392 s
+                                                                        =================================
+                                                                        u : -0.13036E-03  0.10014E+01
+                                                                        v : -0.18251E+00  0.14635E+00
+                                                                        w : -0.20598E-01  0.50819E+00
+                                                                        =================================
+                 Calculate pressures   1246.3290     10.0690      7.0495
+                                                                        raw    pressures : -0.11109E+01  0.90569E+00
+                 Smoothing pressures   1246.4076     10.1477      0.0787
+                                                                        smooth pressures : -0.12983E+00  0.24966E-02
+                do leaf measurements   1246.4198     10.1598      0.0121
+       compute convergence criterion   1246.4605     10.2006      0.0407
+                                                                        velocity_diff_norm = 0.0102941 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1246.4633     10.2034      0.0028
+Compute isostasy and adjust vertical   1246.4634     10.2035      0.0002
+                  Compute divergence   1246.4635     10.2036      0.0000
+                        wsmp_cleanup   1246.4932     10.2333      0.0297
+              Reset surface geometry   1246.4983     10.2384      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations    1246.4994     10.2395      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1246.6638     10.4039      0.1645
+                   Advect surface  1   1246.6640     10.4041      0.0002
+                                                                        S. 1:    3 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1246.7082     10.4483      0.0442
+                 Update cloud fields   1266.0914     29.8315     19.3833
+                    Advect the cloud   1266.8686     30.6087      0.7772
+                        Write output   1271.1845     34.9246      4.3159
+                    End of time step   1273.0121     36.7522      1.8276
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      34   1273.0123
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1273.0124      0.0001      0.0001
+                          surface 01   1273.0124      0.0001      0.0000
+                                                                        S. 1:     4357points
+                      refine surface   1273.0124      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1273.0189      0.0066      0.0065
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4357points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1273.0221
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1273.0223      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1273.0381      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1273.0556      0.0335      0.0175
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.112065E+01
+                                                                        osolve%nleaves=   55287
+    Smoothen octree after refinement   1273.0772      0.0551      0.0216
+                                                                        osolve%nleaves=   55287
+                   Build osolve icon   1273.2134      0.1913      0.1362
+             Imbed surface in osolve   1273.4180      0.3959      0.2046
+                embedding surface  1   1273.4182      0.3961      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1276.0731      3.0510      2.6550
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55245 leaves
+                                                                        nleaves_new=    55287 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1276.0755      3.0534      0.0023
+        Interpolate velo onto osolve   1278.4024      5.3803      2.3269
+                                                                        ov%pressure:     -0.1389E+00 0.1418E-01
+                                                                        osolve%pressure: -0.1298E+00 0.2497E-02
+                      Find bad faces   1278.5659      5.5438      0.1635
+                           Find void   1278.7855      5.7634      0.2196
+                                                                        whole leaf in fluid    25086
+          Define boundary conditions   1278.7999      5.7778      0.0144
+                         wsmp setup1   1278.8063      5.7842      0.0064
+                                                                        n =   103974
+                                                                        nz=  3855270
+                         wsmp setup2   1279.5701      6.5480      0.7638
+ -----------------------------------
+ start of non-linear iteratio      1   1279.6065
+ -----------------------------------
+                        build system   1279.6066      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3735        4212
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29419
+                          wsmp solve   1282.6048      2.9983      2.9982
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0964370 s
+                                                                        wsmp: ordering time:               2.6314709 s
+                                                                        wsmp: symbolic fact. time:         0.5209210 s
+                                                                        wsmp: Cholesky fact. time:         4.1882761 s
+                                                                        wsmp: Factorisation            14983.5625827 Mflops
+                                                                        wsmp: back substitution time:      0.0752611 s
+                                                                        wsmp: from b to rhs vector:        0.0057662 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.5184700 s
+                                                                        =================================
+                                                                        u : -0.10575E-03  0.10015E+01
+                                                                        v : -0.17926E+00  0.14982E+00
+                                                                        w : -0.21407E-01  0.50651E+00
+                                                                        =================================
+                 Calculate pressures   1290.1375     10.5311      7.5327
+                                                                        raw    pressures : -0.11197E+01  0.91320E+00
+                 Smoothing pressures   1290.2170     10.6106      0.0795
+                                                                        smooth pressures : -0.12994E+00  0.25695E-02
+                do leaf measurements   1290.2290     10.6225      0.0119
+       compute convergence criterion   1290.2693     10.6628      0.0403
+                                                                        velocity_diff_norm = 0.0061375 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1290.2721     10.6657      0.0028
+Compute isostasy and adjust vertical   1290.2723     10.6658      0.0002
+                  Compute divergence   1290.2724     10.6659      0.0000
+                        wsmp_cleanup   1290.3020     10.6956      0.0297
+              Reset surface geometry   1290.3070     10.7005      0.0050
+ -----------------------------------------------------------------------
+           Temperature calculations    1290.3081     10.7016      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1290.4727     10.8663      0.1646
+                   Advect surface  1   1290.4729     10.8664      0.0002
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1290.5171     10.9107      0.0443
+                 Update cloud fields   1310.1013     30.4948     19.5841
+                    Advect the cloud   1310.8765     31.2701      0.7753
+                        Write output   1315.2007     35.5943      4.3242
+                    End of time step   1317.0111     37.4047      1.8104
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      35   1317.0113
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1317.0114      0.0001      0.0001
+                          surface 01   1317.0114      0.0001      0.0000
+                                                                        S. 1:     4357points
+                      refine surface   1317.0115      0.0001      0.0000
+                                                                        S. 1:   1 added ptcls in refine_surface
+                      check delaunay   1317.0185      0.0072      0.0071
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                      refine surface   1317.0222      0.0109      0.0036
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1317.0256      0.0143      0.0034
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4358points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1317.0290
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1317.0292      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1317.0451      0.0161      0.0159
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1317.0627      0.0337      0.0176
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.116854E+01
+                                                                        osolve%nleaves=   55329
+    Smoothen octree after refinement   1317.0845      0.0554      0.0218
+                                                                        osolve%nleaves=   55329
+                   Build osolve icon   1317.2221      0.1931      0.1377
+             Imbed surface in osolve   1317.4167      0.3876      0.1945
+                embedding surface  1   1317.4169      0.3878      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1320.0712      3.0422      2.6544
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55287 leaves
+                                                                        nleaves_new=    55329 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1320.0735      3.0445      0.0023
+        Interpolate velo onto osolve   1322.3806      5.3516      2.3071
+                                                                        ov%pressure:     -0.1394E+00 0.1459E-01
+                                                                        osolve%pressure: -0.1299E+00 0.2569E-02
+                      Find bad faces   1322.5912      5.5622      0.2106
+                           Find void   1322.8108      5.7817      0.2196
+                                                                        whole leaf in fluid    25105
+          Define boundary conditions   1322.8244      5.7954      0.0136
+                         wsmp setup1   1322.8304      5.8014      0.0060
+                                                                        n =   104025
+                                                                        nz=  3858189
+                         wsmp setup2   1323.5951      6.5661      0.7647
+ -----------------------------------
+ start of non-linear iteratio      1   1323.6313
+ -----------------------------------
+                        build system   1323.6315      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3736        4215
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29437
+                          wsmp solve   1326.5480      2.9167      2.9165
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.1014030 s
+                                                                        wsmp: ordering time:               2.7387211 s
+                                                                        wsmp: symbolic fact. time:         0.5198281 s
+                                                                        wsmp: Cholesky fact. time:         4.5010691 s
+                                                                        wsmp: Factorisation            14181.2814091 Mflops
+                                                                        wsmp: back substitution time:      0.0810130 s
+                                                                        wsmp: from b to rhs vector:        0.0057769 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.9481850 s
+                                                                        =================================
+                                                                        u : -0.72309E-04  0.10019E+01
+                                                                        v : -0.17383E+00  0.13958E+00
+                                                                        w : -0.21875E-01  0.50592E+00
+                                                                        =================================
+                 Calculate pressures   1334.5098     10.8785      7.9618
+                                                                        raw    pressures : -0.11315E+01  0.92421E+00
+                 Smoothing pressures   1334.5885     10.9572      0.0787
+                                                                        smooth pressures : -0.13008E+00  0.26737E-02
+                do leaf measurements   1334.6051     10.9738      0.0166
+       compute convergence criterion   1334.6467     11.0154      0.0416
+                                                                        velocity_diff_norm = 0.0086979 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1334.6495     11.0182      0.0028
+Compute isostasy and adjust vertical   1334.6497     11.0184      0.0002
+                  Compute divergence   1334.6497     11.0184      0.0000
+                        wsmp_cleanup   1334.6794     11.0481      0.0297
+              Reset surface geometry   1334.6848     11.0535      0.0054
+ -----------------------------------------------------------------------
+           Temperature calculations    1334.6859     11.0546      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1334.8508     11.2195      0.1649
+                   Advect surface  1   1334.8509     11.2196      0.0002
+                                                                        S. 1:    3 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1334.8953     11.2640      0.0444
+                 Update cloud fields   1355.5055     31.8742     20.6102
+                    Advect the cloud   1356.2951     32.6638      0.7896
+                        Write output   1360.6214     36.9900      4.3263
+                    End of time step   1362.4332     38.8019      1.8118
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      36   1362.4333
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1362.4334      0.0001      0.0001
+                          surface 01   1362.4334      0.0001      0.0000
+                                                                        S. 1:     4358points
+                      refine surface   1362.4335      0.0001      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1362.4405      0.0071      0.0070
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4358points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1362.4438
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1362.4439      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1362.4599      0.0161      0.0160
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1362.4774      0.0336      0.0175
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.121731E+01
+                                                                        osolve%nleaves=   55357
+    Smoothen octree after refinement   1362.4992      0.0554      0.0218
+                                                                        osolve%nleaves=   55357
+                   Build osolve icon   1362.6370      0.1932      0.1378
+             Imbed surface in osolve   1362.8304      0.3867      0.1935
+                embedding surface  1   1362.8306      0.3868      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1365.4658      3.0221      2.6352
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55329 leaves
+                                                                        nleaves_new=    55357 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1365.4682      3.0244      0.0023
+        Interpolate velo onto osolve   1367.7783      5.3345      2.3101
+                                                                        ov%pressure:     -0.1397E+00 0.1491E-01
+                                                                        osolve%pressure: -0.1301E+00 0.2674E-02
+                      Find bad faces   1368.0076      5.5638      0.2293
+                           Find void   1368.2250      5.7812      0.2174
+                                                                        whole leaf in fluid    25127
+          Define boundary conditions   1368.2375      5.7938      0.0125
+                         wsmp setup1   1368.2421      5.7983      0.0046
+                                                                        n =   104076
+                                                                        nz=  3861216
+                         wsmp setup2   1368.9929      6.5492      0.7508
+ -----------------------------------
+ start of non-linear iteratio      1   1369.0267
+ -----------------------------------
+                        build system   1369.0269      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3738        4216
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29453
+                          wsmp solve   1371.9886      2.9619      2.9617
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.1038649 s
+                                                                        wsmp: ordering time:               2.5754771 s
+                                                                        wsmp: symbolic fact. time:         0.5202422 s
+                                                                        wsmp: Cholesky fact. time:         3.8781161 s
+                                                                        wsmp: Factorisation            16423.8071603 Mflops
+                                                                        wsmp: back substitution time:      0.0743330 s
+                                                                        wsmp: from b to rhs vector:        0.0056958 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.1580641 s
+                                                                        =================================
+                                                                        u : -0.61750E-04  0.10019E+01
+                                                                        v : -0.15571E+00  0.14076E+00
+                                                                        w : -0.22057E-01  0.50932E+00
+                                                                        =================================
+                 Calculate pressures   1379.1608     10.1340      7.1721
+                                                                        raw    pressures : -0.11429E+01  0.93215E+00
+                 Smoothing pressures   1379.2424     10.2157      0.0817
+                                                                        smooth pressures : -0.13016E+00  0.27814E-02
+                do leaf measurements   1379.2570     10.2302      0.0145
+       compute convergence criterion   1379.2965     10.2697      0.0395
+                                                                        velocity_diff_norm = 0.0141149 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1379.2993     10.2726      0.0029
+Compute isostasy and adjust vertical   1379.2995     10.2728      0.0002
+                  Compute divergence   1379.2995     10.2728      0.0000
+                        wsmp_cleanup   1379.3316     10.3049      0.0321
+              Reset surface geometry   1379.3366     10.3099      0.0050
+ -----------------------------------------------------------------------
+           Temperature calculations    1379.3377     10.3110      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1379.5015     10.4747      0.1638
+                   Advect surface  1   1379.5016     10.4749      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1379.5459     10.5192      0.0443
+                 Update cloud fields   1399.7158     30.6891     20.1699
+                    Advect the cloud   1400.5090     31.4823      0.7932
+                        Write output   1404.8559     35.8292      4.3469
+                    End of time step   1406.6876     37.6609      1.8317
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      37   1406.6878
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1406.6878      0.0001      0.0001
+                          surface 01   1406.6879      0.0001      0.0000
+                                                                        S. 1:     4358points
+                      refine surface   1406.6879      0.0001      0.0000
+                                                                        S. 1:   1 added ptcls in refine_surface
+                      check delaunay   1406.6950      0.0072      0.0071
+                                                                        S. 1:    2 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                      refine surface   1406.6988      0.0110      0.0038
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1406.7021      0.0143      0.0033
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4359points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1406.7055
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1406.7056      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1406.7214      0.0159      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1406.7392      0.0338      0.0178
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.126029E+01
+                                                                        osolve%nleaves=   55399
+    Smoothen octree after refinement   1406.7609      0.0554      0.0216
+                                                                        osolve%nleaves=   55399
+                   Build osolve icon   1406.8974      0.1919      0.1365
+             Imbed surface in osolve   1407.0881      0.3826      0.1907
+                embedding surface  1   1407.0883      0.3829      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1409.7043      2.9988      2.6159
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55357 leaves
+                                                                        nleaves_new=    55399 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1409.7065      3.0010      0.0022
+        Interpolate velo onto osolve   1412.0179      5.3124      2.3114
+                                                                        ov%pressure:     -0.1400E+00 0.1526E-01
+                                                                        osolve%pressure: -0.1302E+00 0.2781E-02
+                      Find bad faces   1412.2969      5.5914      0.2790
+                           Find void   1412.5142      5.8087      0.2173
+                                                                        whole leaf in fluid    25139
+          Define boundary conditions   1412.5267      5.8212      0.0125
+                         wsmp setup1   1412.5315      5.8260      0.0048
+                                                                        n =   104127
+                                                                        nz=  3863667
+                         wsmp setup2   1413.2831      6.5776      0.7516
+ -----------------------------------
+ start of non-linear iteratio      1   1413.3176
+ -----------------------------------
+                        build system   1413.3178      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3741        4219
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29468
+                          wsmp solve   1416.2667      2.9492      2.9490
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0964942 s
+                                                                        wsmp: ordering time:               2.5908401 s
+                                                                        wsmp: symbolic fact. time:         0.5198591 s
+                                                                        wsmp: Cholesky fact. time:         4.0952530 s
+                                                                        wsmp: Factorisation            15495.3344037 Mflops
+                                                                        wsmp: back substitution time:      0.0747411 s
+                                                                        wsmp: from b to rhs vector:        0.0056782 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.3832228 s
+                                                                        =================================
+                                                                        u : -0.52428E-04  0.10019E+01
+                                                                        v : -0.21296E+00  0.22621E+00
+                                                                        w : -0.23112E-01  0.50698E+00
+                                                                        =================================
+                 Calculate pressures   1423.6638     10.3462      7.3971
+                                                                        raw    pressures : -0.11624E+01  0.95484E+00
+                 Smoothing pressures   1423.7444     10.4269      0.0806
+                                                                        smooth pressures : -0.13055E+00  0.34832E-02
+                do leaf measurements   1423.7626     10.4450      0.0181
+       compute convergence criterion   1423.8004     10.4828      0.0378
+                                                                        velocity_diff_norm = 0.0072876 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1423.8032     10.4856      0.0028
+Compute isostasy and adjust vertical   1423.8033     10.4858      0.0002
+                  Compute divergence   1423.8034     10.4858      0.0000
+                        wsmp_cleanup   1423.8341     10.5166      0.0308
+              Reset surface geometry   1423.8392     10.5216      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations    1423.8403     10.5227      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1424.0053     10.6877      0.1650
+                   Advect surface  1   1424.0055     10.6879      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1424.0498     10.7323      0.0444
+                 Update cloud fields   1445.0763     31.7588     21.0265
+                    Advect the cloud   1445.8995     32.5820      0.8232
+                        Write output   1450.2527     36.9351      4.3531
+                    End of time step   1452.0798     38.7622      1.8271
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      38   1452.0800
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1452.0801      0.0001      0.0001
+                          surface 01   1452.0801      0.0001      0.0000
+                                                                        S. 1:     4359points
+                      refine surface   1452.0802      0.0002      0.0000
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1452.0853      0.0053      0.0052
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4359points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1452.0896
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1452.0898      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1452.1056      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1452.1232      0.0337      0.0177
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.130573E+01
+                                                                        osolve%nleaves=   55427
+    Smoothen octree after refinement   1452.1448      0.0552      0.0215
+                                                                        osolve%nleaves=   55427
+                   Build osolve icon   1452.2813      0.1918      0.1366
+             Imbed surface in osolve   1452.4687      0.3791      0.1873
+                embedding surface  1   1452.4689      0.3793      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1455.0621      2.9725      2.5932
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55399 leaves
+                                                                        nleaves_new=    55427 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1455.0643      2.9747      0.0022
+        Interpolate velo onto osolve   1457.3774      5.2879      2.3131
+                                                                        ov%pressure:     -0.1406E+00 0.1523E-01
+                                                                        osolve%pressure: -0.1306E+00 0.3483E-02
+                      Find bad faces   1457.6693      5.5797      0.2919
+                           Find void   1457.8867      5.7971      0.2174
+                                                                        whole leaf in fluid    25147
+          Define boundary conditions   1457.8996      5.8100      0.0129
+                         wsmp setup1   1457.9048      5.8153      0.0052
+                                                                        n =   104154
+                                                                        nz=  3865422
+                         wsmp setup2   1458.6593      6.5697      0.7545
+ -----------------------------------
+ start of non-linear iteratio      1   1458.6953
+ -----------------------------------
+                        build system   1458.6956      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3740        4223
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29479
+                          wsmp solve   1461.6593      2.9640      2.9638
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0939209 s
+                                                                        wsmp: ordering time:               2.6639061 s
+                                                                        wsmp: symbolic fact. time:         0.5217431 s
+                                                                        wsmp: Cholesky fact. time:         4.4999781 s
+                                                                        wsmp: Factorisation            14167.2108120 Mflops
+                                                                        wsmp: back substitution time:      0.0802050 s
+                                                                        wsmp: from b to rhs vector:        0.0056829 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.8658011 s
+                                                                        =================================
+                                                                        u : -0.44528E-04  0.10019E+01
+                                                                        v : -0.20827E+00  0.22524E+00
+                                                                        w : -0.23423E-01  0.56272E+00
+                                                                        =================================
+                 Calculate pressures   1469.5385     10.8431      7.8791
+                                                                        raw    pressures : -0.12137E+01  0.10013E+01
+                 Smoothing pressures   1469.6199     10.9245      0.0814
+                                                                        smooth pressures : -0.13058E+00  0.36789E-02
+                do leaf measurements   1469.6380     10.9426      0.0181
+       compute convergence criterion   1469.6766     10.9813      0.0387
+                                                                        velocity_diff_norm = 0.0116445 (tol = 0.0100000)
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1469.6794     10.9841      0.0028
+Compute isostasy and adjust vertical   1469.6796     10.9843      0.0002
+                  Compute divergence   1469.6797     10.9843      0.0000
+                        wsmp_cleanup   1469.7102     11.0149      0.0306
+              Reset surface geometry   1469.7153     11.0200      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations    1469.7164     11.0211      0.0011
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1469.8808     11.1855      0.1644
+                   Advect surface  1   1469.8810     11.1856      0.0002
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1469.9244     11.2291      0.0435
+                 Update cloud fields   1491.0820     32.3866     21.1576
+                    Advect the cloud   1491.8729     33.1776      0.7909
+                        Write output   1496.2516     37.5562      4.3787
+                    End of time step   1498.1254     39.4301      1.8738
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+               Start of Step      39   1498.1256
+ ----------------------------------------------------------------------- 
+ ----------------------------------------------------------------------- 
+                  Loop over surfaces   1498.1256      0.0001      0.0001
+                          surface 01   1498.1257      0.0001      0.0000
+                                                                        S. 1:     4359points
+                      refine surface   1498.1257      0.0002      0.0001
+                                                                        S. 1:   0 added ptcls in refine_surface
+                      check delaunay   1498.1318      0.0062      0.0061
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                                                                        S. 1:     4359points
+ ----------------------------------------------------------------------- 
+          Start of iteration       1   1498.1361
+ ----------------------------------------------------------------------- 
+                 Create octree solve   1498.1363      0.0002      0.0002
+                                                                        (1) params%griditer < 0                F
+                                                                        (2) current_level==params%levelmax_oct T
+                                                                        (3) increase_current_level             T
+                                                                        converge = (1) & (2) & (3) ->  F
+        compute ref. criterion on ov   1498.1521      0.0160      0.0158
+                                                                        osolve%nleaves=   32768
+        improve osolve for criterion   1498.1702      0.0341      0.0181
+                                                                        max(crit)        =   3.200000E+01
+                                                                        max(crit,w_l_i_f)=   3.134700E+01
+                                                                        osolve%nleaves=   55448
+    Smoothen octree after refinement   1498.1917      0.0556      0.0215
+                                                                        osolve%nleaves=   55448
+                   Build osolve icon   1498.3284      0.1923      0.1367
+             Imbed surface in osolve   1498.5181      0.3820      0.1897
+                embedding surface  1   1498.5183      0.3822      0.0002
+                                                                        osolve%nleaves=   61440
+          Assessing octree stability   1501.1170      2.9809      2.5987
+                                                                        current refine level     6
+                                                                        After criterion based refinement:
+                                                                        nleaves_old=    55427 leaves
+                                                                        nleaves_new=    55448 leaves
+                                                                        After surfaces embedding:
+                                                                        nleaves_old=    61440 leaves
+                                                                        nleaves_new=    61440 leaves
+                                                                        C2: authorise increase of refine level =  T
+                     Change 3D cloud   1501.1192      2.9831      0.0022
+        Interpolate velo onto osolve   1503.4172      5.2811      2.2980
+                                                                        ov%pressure:     -0.1414E+00 0.1570E-01
+                                                                        osolve%pressure: -0.1306E+00 0.3679E-02
+                      Find bad faces   1503.7087      5.5726      0.2915
+                           Find void   1503.9260      5.7899      0.2173
+                                                                        whole leaf in fluid    25161
+          Define boundary conditions   1503.9387      5.8026      0.0127
+                         wsmp setup1   1503.9435      5.8074      0.0047
+                                                                        n =   104193
+                                                                        nz=  3867345
+                         wsmp setup2   1504.6986      6.5625      0.7551
+ -----------------------------------
+ start of non-linear iteratio      1   1504.7317
+ -----------------------------------
+                        build system   1504.7319      0.0002      0.0002
+                                                                        nelem per proc (min/max)        3742        4224
+                                                                        viscosity range   1.00000E-05  1.30880E+03
+                                                                        viscosity capped in     29493
+                          wsmp solve   1507.6977      2.9660      2.9658
+                                                                        =================================
+                                                                        wsmp: initialisation time:         0.0949180 s
+                                                                        wsmp: ordering time:               2.5654900 s
+                                                                        wsmp: symbolic fact. time:         0.5213399 s
+                                                                        wsmp: Cholesky fact. time:         4.0147929 s
+                                                                        wsmp: Factorisation            15953.0428377 Mflops
+                                                                        wsmp: back substitution time:      0.0798790 s
+                                                                        wsmp: from b to rhs vector:        0.0058281 s
+                                                                        ---------------------------------
+                                                                        wsmp: total solve time:            7.2825990 s
+                                                                        =================================
+                                                                        u : -0.34036E-04  0.10019E+01
+                                                                        v : -0.22691E+00  0.21529E+00
+                                                                        w : -0.23665E-01  0.56372E+00
+                                                                        =================================
+                 Calculate pressures   1514.9943     10.2626      7.2966
+                                                                        raw    pressures : -0.12096E+01  0.99530E+00
+                 Smoothing pressures   1515.0755     10.3438      0.0812
+                                                                        smooth pressures : -0.13066E+00  0.39648E-02
+                do leaf measurements   1515.0923     10.3606      0.0168
+       compute convergence criterion   1515.1319     10.4002      0.0396
+                                                                        velocity_diff_norm = 0.0070541 (tol = 0.0100000)
+                                                                         
+                                                                        convergence on this octree reached
+                                                                        increase_current_level -> T
+                                                                        current_level -> 6
+                    slicing the cube   1515.1347     10.4030      0.0028
+Compute isostasy and adjust vertical   1515.1349     10.4032      0.0002
+                  Compute divergence   1515.1349     10.4032      0.0000
+                        wsmp_cleanup   1515.1663     10.4346      0.0314
+              Reset surface geometry   1515.1714     10.4397      0.0051
+ -----------------------------------------------------------------------
+           Temperature calculations    1515.1726     10.4409      0.0012
+ -----------------------------------------------------------------------
+                                                                        =================================
+                                                                        skip temperature calculation
+                                                                        =================================
+                 Advect the surfaces   1515.3386     10.6069      0.1660
+                   Advect surface  1   1515.3387     10.6070      0.0002
+                                                                        S. 1:    1 flips in delaunay2
+                                                                        S. 1:    0 flips in delaunay2
+                    Erode surface  1   1515.3836     10.6520      0.0449
+Killed by signal 2.
+Killed by signal 2.
+Killed by signal 2.
+Killed by signal 2.
+Killed by signal 2.
-- 
GitLab