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://barre.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