Added dilatation rate, bulk viscosity, initial e2d

src/module_definitions.f90

@@ -95,7 +95,7 @@ module definitions
       double precision,dimension(:),pointer::lode
       double precision,dimension(:),pointer::u,v,w,wiso,temp
       double precision,dimension(:,:),pointer::lsf
-      double precision,dimension(:),pointer::eviscosity,q
+      double precision,dimension(:),pointer::eviscosity,q,dilatr
       integer,dimension(:,:),pointer::icon,iface
       logical,dimension(:),pointer::is_plastic
       end type octreesolve
@@ -216,7 +216,10 @@ module definitions
       logical isostasy
       logical flexure
       logical isobc
+      logical bulkvisc
+      logical init_e2d
       integer nmpi
+      integer matrule
       integer levelcut
       integer levelapprox
       integer nstep
@@ -256,6 +259,7 @@ module definitions
       double precision elastic_plate_thickness
       double precision density_difference
       double precision vex
+      double precision e2d0
       integer baselevelx0
       integer baselevelx1
       integer baselevely0