The magnetic measurements of declination (D), horizontal (H) and vertical (Z) components of earth’s magnetic field, collected from ground surveys between 1962 and 1966, are used to develop an analytical model of geomagnetic field variations over Indian region for the epoch 1965. In order to reflect spatial features with wavelengths of approximately 1000 km, sixth degree polynomial as a function of differential latitude and longitude is calculated by the method of least squares. The root mean square fit of the model to the input data is better than that accounted by the International Geomagnetic Reference Field for 1965.0. Isomagnetic charts drawn forD, H, Z and total force (F) reflect more details than that shown on world magnetic charts. Further, the values of the field at common repeat stations recorded between 1962 and 1974, after eliminating the field values for the epoch 1965.0, are used to get the secular variation as well as its spatial dependence again by means of polynomial which now includes coefficients which are functions of time and of geographical locations. The accuracy of coefficients is tested against the behaviour of secular variation at permanent magnetic observatories. The merits and limitations of the model are discussed.