anchor.diff:              ! Calculate location of station relative to fault center
               if (coord_type.eq.'geographic') then
                   call tri_center(center,tri(:,1),tri(:,2),tri(:,3))
  +                call lola2distaz(center(1),center(2),sta(1),sta(2),dist,az)
  +                dist = dist*radius_earth_m
  +                sta_new(1) = dist*dsin(az*d2r)
  +                sta_new(2) = dist*dcos(az*d2r)
                   do iTri = 1,3
                       call lola2distaz(center(1),center(2),tri(1,iTri),tri(2,iTri),dist,az)
                       dist = dist*radius_earth_m
  
positive.code: ! Calculate stress Green's function for each fault-fault pair
  do i = 1,fault%nrows
  
      ! Station coordinates
      if (gf_model.eq.'okada_rect') then
          sta(1) = fault%array(i,1)
          sta(2) = fault%array(i,2)
  
negative.code:         if (coord_type.eq.'geographic') then
              call tri_geo2cart(pt1,pt2,pt3,fault%array(i,1:3),fault%array(i,4:6), &
                                fault%array(i,7:9),'m')
              call tri_geometry(nvec,tvec1,tvec2,pt1,pt2,pt3)
          elseif (coord_type.eq.'cartesian') then
              call tri_geometry(nvec,tvec1,tvec2,fault%array(i,1:3),fault%array(i,4:6), &
                                fault%array(i,7:9))
          endif
      endif
  
      do j = 1,fault%nrows
  