anchor.diff:      allocate( VMTR_RGAM       (ADM_gall,   ADM_kall,ADM_lall   ) )
       allocate( VMTR_RGAM_pl    (ADM_gall_pl,ADM_kall,ADM_lall_pl) )
   
  +    allocate( VMTR_RGAMH      (ADM_gall,   ADM_kall,ADM_lall   ) )
  +    allocate( VMTR_RGAMH_pl   (ADM_gall_pl,ADM_kall,ADM_lall_pl) )
  +
       allocate( VMTR_GAM2       (ADM_gall,   ADM_kall,ADM_lall   ) )
       allocate( VMTR_GAM2_pl    (ADM_gall_pl,ADM_kall,ADM_lall_pl) )
   
  
positive.code: 
               VMTR_RGSH_pl    (n,k,l) = 1.D0 / GSQRTH_pl(n,k)
               VMTR_RGAM_pl    (n,k,l) = 1.D0 / GAM
               VMTR_RGAM2_pl   (n,k,l) = 1.D0 / VMTR_GAM2_pl (n,k,l)
               VMTR_RGAM2H_pl  (n,k,l) = 1.D0 / VMTR_GAM2H_pl(n,k,l)
               VMTR_RGSGAM2_pl (n,k,l) = 1.D0 / VMTR_GSGAM2_pl (n,k,l)
  
negative.code:   !--- 1/gamma at the integer level
    real(8), public, allocatable, save :: VMTR_RGAM   (:,:,:)
    real(8), public, allocatable, save :: VMTR_RGAM_pl(:,:,:)
  
    !--- 1/gamma at the half-integer level
    real(8), public, allocatable, save :: VMTR_RGAMH   (:,:,:)
    real(8), public, allocatable, save :: VMTR_RGAMH_pl(:,:,:)
  
    !--- 1/gamma^2 at the integer level
    real(8), public, allocatable, save :: VMTR_RGAM2   (:,:,:)
  