anchor.diff:  Write(3,*) "! --- Number of tide gauges"
   Write(3,*) "INTEGER, PARAMETER :: NTIDEGAUGES=", NTIDEGAUGES
   !
  -If(option_3d=='y')then 
  +!If(option_3d=='y')then 
   Write(3,*) "!"
   Write(3,*) "!"
   Write(3,*) "!------------------- 8- Predictions at geodetic sites ----------------"
  
positive.code: !
  !
  ! Set default parameters
  RSL_FILE = "'none'"
  RSL_DATABASE = "'none'"
  RSL_DATABASE_FORMAT = '-1'
  NRSL = 0
  TGAUGES_DATABASE="'none'"
  NTIDEGAUGES=0
  degree_st_min = "0"
  degree_st_max = "0"
  
  ! ************************************************************************
  ! ************************************************************************
  ! ************************************************************************
  
negative.code: 						   GEODETIC_DATABASE
  Write(3,*) "! --- Number of tide gauges"
  Write(3,*) "INTEGER, PARAMETER :: NGEOD=", NGEOD
  !Endif 
  !
  If(option_3d_regions=='y')then 
  Write(3,*) "!"
  