Subroutine MFGELO
c    ===================
c
c Loads the handy geometry common block.
c
c /MFGEOX/  spatial coordinates and geometrical quantities of DC/LT/LW etc.
c
      Implicit None
c
#include "mfgeox.inc"
c
      Integer p,o,is,q
      Real    R
c
      Real rd /57.29577951/
c
c--- geometry
c
C ZLT
      ZLT(1,1)=465.1600
      ZLT(1,2)=465.6400
      ZLT(1,3)=465.1600
      ZLT(1,4)=465.6400
C
      ZLT(2,1)=788.1600
      ZLT(2,2)=788.6400
      ZLT(2,3)=788.1600
      ZLT(2,4)=788.6400
C
      ZLT(3,1)=887.2600
      ZLT(3,2)=887.7400
      ZLT(3,3)=887.2600
      ZLT(3,4)=887.7400
C
      ZLT(4,1)=992.5100
      ZLT(4,2)=992.9900
      ZLT(4,3)=992.5100
      ZLT(4,4)=992.9900
C
C ZLW
      ZLW(1,1)=673.8600
      ZLW(1,2)=659.8600
      ZLW(1,3)=673.8600
      ZLW(1,4)=659.8600
      ZLW(1,5)=673.8600
      ZLW(1,6)=659.8600
      ZLW(1,7)=673.8600
      ZLW(1,8)=659.8600
C
      ZLW(2,1)=769.8600
      ZLW(2,2)=755.8600
      ZLW(2,3)=769.8600
      ZLW(2,4)=755.8600
      ZLW(2,5)=769.8600
      ZLW(2,6)=755.8600
      ZLW(2,7)=769.8600
      ZLW(2,8)=755.8600
C ZDC
      ZDC(1,1)= 476.7
      ZDC(1,2)= 492.3
      ZDC(1,3)= 492.3
      ZDC(1,4)= 476.7
      ZDC(1,5)= 476.7
      ZDC(1,6)= 492.3
      ZDC(1,7)= 492.3
      ZDC(1,8)= 476.7
      do o=1,8
       if(mod(o,2).eq.1) then
        ZDC(2,o)=815.8
        ZDC(3,o)=975.8
        ZDC(4,o)=1025.8
       else
        ZDC(2,o)=800.2
        ZDC(3,o)=960.2
        ZDC(4,o)=1010.2
       endif
      enddo
      Do o=1,8
       q=(o+1)/2
       ZFM( 1,o)=ZLT(1,q)
       ZFM( 2,o)=ZDC(1,o)
       ZFM( 3,o)=ZLW(1,o)
       ZFM( 4,o)=ZLW(2,o)
       ZFM( 5,o)=ZLT(2,q)
       ZFM( 6,o)=ZDC(2,o)
       ZFM( 7,o)=ZLT(3,q)
       ZFM( 8,o)=ZDC(3,o)
       ZFM( 9,o)=ZLT(4,q)
       ZFM(10,o)=ZDC(4,o)
      Enddo
C
C-- ltp
      fsel_ltp=0.7031250
      spsel_ltp=1.406250
C-- ltr
      fsel_ltr=75.25000
      spsel_ltr=1.9
C-- lwp
      fsel_lwp=0.3415625
      spsel_lwp=0.7031250
C-- lwr
      fsel_lwr=190.3554
      spsel_lwr=1.8
c--> fastor (16 strip)
       fs_fastor=fsel_lwr+7.5*spsel_lwr
       ss_fastor=16*spsel_lwr
c--> 4-plets (4 strip)
       fs_4plet=fsel_lwr+1.5*spsel_lwr
       ss_4plet=4*spsel_lwr
C-- dc
      fsel_dci=69.00736
      spsel_dci=7.295160
      fsel_dco=68.99223
      spsel_dco=7.295160
      inorout(1)=1
      inorout(2)=2
      inorout(3)=2
      inorout(4)=2
      fsdc(1)=fsel_dci
      fsdc(2)=fsel_dco
      spsdc(1)=spsel_dci
      spsdc(2)=spsel_dco
      psw(1)=69.007362366
      psw(2)=68.992225647
      End
c