Matlab basic Spherical Harmonic Sythesis tools

Example:

• load the coefficients:

[i,header]=load_icgem('ggm05g.gfc.txt');

• convert to geoid height:

g=mod_convert(i,'non-dim','geoid',header.earth_gravity_constant,header.radius);

• perform synthesis with 1 degree samplint at the equator:

[long,lat,grid_out]=mod_sh_synth(g,0,360);

• plot the geoid height above the Earth mean radius of the model:

plot(long*180/pi,grid_out-header.radius);