### Ex ##

################################################################################################
######          The follwoing commented lines will extract              ########################
######          subgrids from the global grids of EGM08.                ########################
######          The subgrids are provided along with the exercises      ########################
################################################################################################
############################# extract EGM08 subGRIDs #############################
# GRAVSOFT/geoip <<!
# /home/oa/PAVLIS/egm08_2190nx.bin
# egm08_nx.grd
# 4 0 0 f
# f 14.00833333333 20.00833333333 188.00833333333 195.00833333333 0.0166666666667 0.0166666666667
# !

# GRAVSOFT/geoip <<!
# /home/oa/PAVLIS/egm08_2190gx.bin
# egm08_gx.grd
# 4 0 0 f
# f 14.00833333333 20.00833333333 188.00833333333 195.00833333333 0.0166666666667 0.0166666666667
# !

# GRAVSOFT/geoip <<!
# /home/oa/PAVLIS/dot_min1x1_N180.bin
# egm08_MDT.grd
# 4 0 0 f
# f 14.00833333333 20.00833333333 188.00833333333 195.00833333333 0.0166666666667 0.0166666666667
# !

############################  Reduce altimetry data ################################
GRAVSOFT/geoip <<!
egm08_nx.grd
DATA/rads_SLA.alt
13 4 0 f
1
DATA/rads_SSH.dat
!

GRAVSOFT/geoip <<!
egm08_MDT.grd
DATA/ReducedAlt_C2_Rads.dat
13 4 0 f
1
DATA/rads_SLA.alt
!


#####################  Remove trend,Crossover adjustment  ########################
GRAVSOFT/aline_rads <<!
DATA/ReducedAlt_C2_Rads.dat
DATA/dddd.g
DATA/dddd
   15.800000000000001          18.20000000000001          189.80000000000001          193.20000000000001     
!
rm crsdif.dat crspar.dat 
GRAVSOFT/crsadj_ad <<! >DATA/res.dat                  
DATA/dddd
DATA/ddd1
-90.0  90.0 -200.0 362.0
55 0.00001
f
!

mv DATA/ddd1 DATA/dddd 

GRAVSOFT/aline_rads <<!
DATA/dddd
DATA/dddd.g
DATA/ddd1
   15.800000000000001          18.20000000000001          189.80000000000001          193.20000000000001     
!
cd DATA
mv ddd1 dddd
cd ..
rm crsdif.dat crspar.dat 
GRAVSOFT/crsadj_ad <<! >>DATA/res.dat                 
DATA/dddd
DATA/ddd1
-90.0  90.0 -200.0 362.0
55 0.00001
f
!

### final output ###
cp DATA/ddd1 DATA/geo.016.190.alt              


rm fort.*
rm crsdif.dat crspar.dat