anchor.diff:          real*8,dimension(Dij_d0,Dij_d1),intent(out) :: Dij
           integer,dimension(dh_mask_d0,dh_mask_d1),intent(in) :: dh_mask
           integer,intent(in) :: n_dh
  -        integer,intent(in),optional :: sampling
  +        integer,intent(in) :: sampling
           integer,intent(in) :: lmax
           integer,intent(in) :: dh_mask_d0
           integer,intent(in) :: dh_mask_d1
  
positive.code:         integer,intent(out) :: exitstatus
          real*8,dimension(cindex_d0,cindex_d1),intent(in) :: cindex
          real*8,dimension(cilm_d0,cilm_d1,cilm_d2),intent(out) :: cilm
          integer,optional,intent(in) :: degmax
          integer,intent(in) :: cindex_d0
          integer,intent(in) :: cindex_d1
          integer,intent(in) :: cilm_d0
  
negative.code:         real*8,dimension(admit_d0),intent(out) :: admit
          real*8,dimension(corr_d0),intent(out) :: corr
          integer,intent(in) :: k
          real*8,dimension(admit_error_d0),intent(out) :: admit_error
          real*8,dimension(corr_error_d0),intent(out) :: corr_error
          real*8,dimension(taper_wt_d0),intent(in) :: taper_wt
          integer,intent(in) :: mtdef
          integer,intent(in) :: k1linsig
          integer,intent(in) :: taper_order_d0
          integer,intent(in) :: g_d0
          integer,intent(in) :: g_d1
  