Home > utilities > grv.m

grv

PURPOSE

SYNOPSIS

function g=grv(lat)

DESCRIPTION

CROSS-REFERENCE INFORMATION

This function calls:

This function is called by:

cor25 x=[6.3 0 28.7 26 24.2 18.5 141 419 0 600 1010 15 15 15 0 1];
cor26_5a Version with shortened iteration; modified Rt and Rq

SOURCE CODE

0001 function g=grv(lat)
0002 gamma=9.7803267715;
0003 c1=0.0052790414;
0004 c2=0.0000232718;
0005 c3=0.0000001262;
0006 c4=0.0000000007;
0007 
0008 phi=lat*pi/180;
0009 x=sin(phi);
0010 g=gamma*(1+c1*x.^2+c2*x.^4+c3*x.^6+c4*x.^8);
0011