## compute gravity anomalies from Height Anomalies ###


rm DATA/dddd.grd
rm DATA/dddd.tmp

### grid the altimetry data and derive residule gravity anomaly grid using geofour ###
GRAVSOFT/geogrid <<! 
DATA/geo.016.190.alt              
DATA/geo.016.190.geo                              
DATA/geo.016.190.err                              
2 -1 20 0 1
  14.000 0.01
1 8.0
   15.508333333333333        18.491666666666667        189.50833333333333        193.49166666666667     
  1.66666666666666664E-002  1.66666666666666664E-002
!
GRAVSOFT/geofour <<! 
DATA/geo.016.190.geo                              
DATA/ddd
DATA/dddd.grd
 3   14.000
  15.5083 189.5083 256 512 8
!

########### select grid, using geoip ############
rm DATA/dddd.tmp
GRAVSOFT/geoip <<!
DATA/geo.016.190.geo                              
DATA/dddd.tmp
4 4 0 f
f
   15.841666666666667        18.158333333333335        189.84166666666667        193.15833333333333       1.66666666666666664E-002  1.66666666666666664E-002
!
mv DATA/dddd.tmp DATA/geo.016.190.geo                              
rm DATA/dddd.tmp
GRAVSOFT/geoip <<!
DATA/geo.016.190.err                              
DATA/dddd.tmp
4 4 0 f
f
   15.841666666666667        18.158333333333335        189.84166666666667        193.15833333333333       1.66666666666666664E-002  1.66666666666666664E-002
!
mv DATA/dddd.tmp DATA/geo.016.190.err                              
GRAVSOFT/geoip<<!
DATA/dddd.grd
DATA/dddd.tmp
4 4 0 f
f
   15.841666666666667        18.158333333333335        189.84166666666667        193.15833333333333       1.66666666666666664E-002  1.66666666666666664E-002
!
cp DATA/dddd.tmp DATA/geo.016.190.est


############ restore the full wavelength  #############################    
GRAVSOFT/geoip <<!
egm08_gx.grd
DATA/geo.016.190.gra                              
16 4 0 f
DATA/dddd.tmp
!


########## Interpolate along the ship track ####################

GRAVSOFT/geoip <<!
DATA/geo.016.190.gra
DATA/dg_diff.dat
13 4 0 f
1
Pacific_Ship_track_dg.dat
!