Procedure for calculating three components of vector field from spherical harmonic using either geocentric or geodetic coordinates as input and output is described. Three subroutines of computer program are explained. Program is written in FORTRAN for IBM 360 computer