Loading computation.

   This is an auxiliary dataset that provides the contribution of 
   the degree-1 spherical harmonics to loading displacement: 
   the difference in displacements between the motion with respect 
   to the center of frame, i.e. the center of the solid Earth mass 
   excluding the contribution of the ocean. When this displacement 
   is added to the displacement computed with respect to the center 
   of mass of the total Earth including the mass of the ocean, 
   the sum becomes the displacement with respect to the center 
   of mass of the solid Earth, excluding the ocean.


   Displacements are computed in frequency domain. The original data are
   re-gridded to the uniform equiangular grid D2047 (4097x8192) with
   refined land-sea mask. The points of the ocean tides model that fall
   in land according to GTOPO30 are considered as land. The points of 
   the ocean tide model that fall in ocean according to GTOPO30 but 
   have the fill value (i.e. considered as land in odenan tide model) 
   are computed from adjacent points by interpolation with with 
   exponentially decaying weights.

   Then for every tide, cosine and sine components of ocean tide heights 
   is expanded into spherical harmonic of degree/order 2047. The harmonics
   are scaled by factors that depends on loading Love numbers. These scaled 
   spherical harmonics are used for computing 3D tidal displacement: 
   up, east, and north components. The up component is defined
   as a displacement towards direction orthogonal to the reference ellipsoid.

   Adopted density of sea water: 1025 kg/meter^3.

   The displacement is computed with respect to the center of the total 
   Earth's mass, including the mass of the ocean.
