% Final version of the code,
now tracking brittle-injected crust (BIC), 
ductile-injected crust (DIC) and mantle.
Olive (8/08)
Also sediment-basement surface (de Sagazan, 6/18)

The code that resets the phase numbers during remeshing 
is written in the rem_interpolate subroutine, which is part of remesh.f90.


Tracers:
- initiated in init_tracer.f90
- move in tracer_move.f90
- are interpolated during remeshing in remesh_tracers.f90

Sediment-basement surface:
- initiated in init_surf.f90
- moves in surf_move.f90
- is reinterpolated at every time step (to make sure elements do not leave the domain) in interp_surf.f90
- surface elements are relocated after interpolation and remeshings in remesh_surf.f90

M cutoff:
in subroutine cutM.f90

+ obvious modifications in read_params.f90



%%%%%%%%%%%%%%%% MODIFICATIONS IN THE DREZINA FILE, in the magma injection section :

TRACERS

;Initial position of tracer
0 1000 -15000

=tracer or not ? 0 or 1
 initial position of tracer (x and y, in m)


TRACKING PHASES

;T_solidus, T_liquidus,  rate_inject_brittle (m/s), rate_inject_ductile (m/s), inj_period, ext_period, inj_count, ynreset, ynmove, Mcutoff,  whencut (Myr),duration of transition (Myr),  rate_inject_brittle2 (m/s), rate_inject_ductile2 (m/s)
1000,         1300,        7.922e-10,   7.922e-10 ,          1,        0,        0,       1,     0,    1,  1.0, 0.3, 3.169e-10, 3.169e-10


ynreset: resets the phase numbers in the injection zone: (=1, tracks crust and mantle, =2, tracks BIC, DIC and mantle, = 0, doesn't track anything)
CAREFUL ! if you track 3 phases you must initiate them !
           


ynmove: the injection zone can move freely or not
M value in the brittle or ductile layers = rate_inject_ductile or rate_inject_brittle
Mcutoff: implements a sudden change in the M value or not.
IF yes, the cutoff happens after 'whencut' (in Myrs) and the transition between both M values takes 'duration of transition' (in Myrs, 0 gives a sudden change)  
The new M values are rate_inject_ductile2 and rate_inject_brittle2.

