From GoogleEarth to loading effect (gravity) at altitude

in six steps
=================================================================

(1) Take a snapshot of GoogleEarth using PaintShopPro,

(2) Turn the sea cells: 
white
      the origo point:     
red

(3) Draw a 1-km precisely horizontal
     straight-line:         
yellow.
     Recolor the endpoint so that the origo remains the only red point.
     Turn non-sea cells into: neither white nor red nor yellow.

(4) Save in uncompressed tiff
     and move to linux (
Oload/),

(5) Use 
convert

       
convert Image3-uncompressed.tif Image3-uncompressed.txt

     Find the red point

       
fgrep red Image3-uncompressed.txt    =>  (io, jo)

     Measure out the 1.0 km in units of pixels:

       
fgrep yellow Image3-uncompressed.txt

    or
    awk -F, '/yellow/{sub(/:/,"",$2);print $1,$2}' Image3-uncompressed.txt |\
      minmax |\
        awk '{gsub(/\//,"",$0);print $7-$6}'  
=> 
r  [pixels/km]

(6) Calculation
    If the station is z meters above sea level, multiply z·r/1000   => dimensionless height z'

    Write a main program: Catch all white points,
    use subroutine
dtile_grav  (Oload/afor/p/tilegrav.f)  to sum up the gravity effect.
    Box size = (1,1), load location = (i -
io, j jo)
    Multiply with  G·ρ  and a tide coefficient [m]  =>  Loading effect at height.

    Program exists: 
~/OTEQ/PROG/OLOAD/hiresgsam.f
      program hi_resolution ! Newtionan gravity for site at altitude
c

c USAGE: hiresgam -a<alt> [-c] -f <file>
c
c -a<alt>   - station altitude in m
c -c        - correction for spherical earth
c -f <file> - The converted tiff file, see
c             http://holt.oso.chalmers.se/hgs/OTEQ/hires-loading-HOW.TO.html
c
c setenv ABSOFT_RT_FLAGS '-defaultcarriage'
c hiresgsam -a7.0 -f Image3-uncompressed.txt
c unsetenv ABSOFT_RT_FLAGS
c
c
DONE