Versions of thermo-mechanical model code:

HTR HTRL HTRV20 HTRV45

      thermal grid    201 x 68
      Eulerian grid   201 x 41
      Lagrangian grid 501 x 41
      
    Code Versions: 
      
      jun15_201x68    old version, first port from SP2
                      [spmicx]
      nov21_201x68    cooling histories added to microfem_exhum_s
      jan11.01_201x68 saves peak metamorphic data for all
                      Lagrangian particles in microfem_peaks_s
      feb14.01_201x68 Extends Lagrangian grid to the left
                      [test_extend_spmicx]
      mar05_201x68    Fixed conservation of mass for advancing subduction
      may08_201x68    Fixed bug which assigned wrong initial temperatures
                      in extended Lagrangian grid.
                      [nonzero_vr]
                      

HTRR HTRRV
      thermal grid     401 x 68
      Eulerian grid    401 x 41
      Lagrangian grid 1001 x 41

    Code Versions:
       
      sep6_401x68    RR version of jun15_201x68 code
                     [rrspmicx]
      extend_401x68  shifts Lagrangian grid to the left by a specified
                     amount
      may08_401x68   The same as may08_201x68 code except for grid resolution
      jul18_401x68   When microfem_g1_s is written, the field "aa" (array nela0)
                     is now written out in I8 format instead of I5.  The corresponding
                     READ statement now uses free format instead of I5 format. 
                     [extend_rrspmicx]

HTRT
      thermal grid     201 x 88
      Eulerian grid    201 x 41
      Lagrangian grid  501 x 41
 
    Code Versions:
       
      mar05_201x88    RT version of mar05_201x68 
      may08_201x88    RT version of may08_201x68
            
HTR3V
      thermal grid     601 x 68
      Eulerian grid    601 x 41
      Lagrangian grid 1501 x 41

    Code Versions:

      jul18_601x68    When microfem_g1_s is written, the field "aa" (array nela0)
                      is now written out in I8 format instead of I5.  The corresponding
                      READ statement now uses free format instead of I5 format. 
                      (same as jul18_401x68 except for grid resolution).

      jul31_601x68    Added Lax stabilization to erosion code.
      aug29_601x68    Modified fead.f - negative temperatures are set to zero
                      within the thermal iteration loop.

????
      thermal grid     801 x 68
      Eulerian grid    801 x 41
      Lagrangian grid 2001 x 41

    Code Versions:

      jul18_801x68    Same as jul18_401x68 except for grid resolution.