1 MODULE trace 2 3 4 ! Trace of D_s 5 ! Note same quantity in visc_s_mod: trd_s 6 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: trD_s_C 7 8 ! Trace of the square of D_s 9 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: trD_s2 10 11 ! Trace of D_s at previous timestep 12 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: trD_s_Co 13 14 ! Trace of D_s at previous time 15 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: trD_s_Co2 16 17 ! Trace of the dot of D_s (Mth solids phase) and 18 ! D_sl (Lth solids phase) 19 DOUBLE PRECISION, DIMENSION (:,:,:), ALLOCATABLE :: trD_s2_ip 20 21 END MODULE trace 22