code: def getvec(gh, lat, lon):
      
      sv = []
      pad = 120 - len(gh)
      for x in range(pad):
          gh.append(0.)
      for x in range(len(gh)):
          sv.append(0.)
  #! convert to colatitude for MB routine
      itype = 1
      colat = 90. - lat
      date, alt = 2000., 0.  # use a dummy date and altitude
      x, y, z, f = magsyn(gh, sv, date, date, itype, alt, colat, lon)
      vec = cart2dir([x, y, z])
      vec[2] = f
      return vec
docs: Evaluates the vector at a given latitude and longitude for a specified
      set of coefficients
  
      Parameters
      ----------
      gh : a list of gauss coefficients
      lat : latitude of location
      long : longitude of location
  
      Returns
      -------
      vec : direction in [dec, inc, intensity]