SUBROUTINE ld_geom(success)
c============================================================================
c PROCEDURE NAME: 	ld_geom
c PROCEDURE TYPE: 	Fortran subroutine
c VERSION DATE:		23.VII.98
c
c PURPOSE:
c   Loads data into common block /geom/, which contains the geometric
c   data used by the program.  These are mostly copied or calculated
c   from the data in common block /instdata/.
c
c CALLING PARAMETERS:
c   Name	Type	i/o	units	Description
c   --------	----	---	-----	--------------
c   success      l       o              Successful completion flag
c
c EXTERNAL ROUTINES REQUIRED:
c
c NOTES & CAVEATS:
c
c REVISION HISTORY:
c
c============================================================================
      IMPLICIT NONE
c----------------------------------------------------------------------------
c Compile parameters
c

c----------------------------------------------------------------------------
c Calling parameters
c
      LOGICAL success

c----------------------------------------------------------------------------
c Common blocks
c
      INCLUDE "geom.cmn"
      INCLUDE "instdata.cmn"

c----------------------------------------------------------------------------
c Local variables
c
      INTEGER*4 i, j
      REAL*8 radeg      /57.295779513082323/    !deg-rad conversion

c============================================================================


c--Parameters for field-of-view rotations:  
c--         Rotation= asin(cos(180-i)/cos(lat)) ascending side, and 180 minus
c--                          this quantity on the descending side
c--   where lat is the latitude of the satellite (not of the center of FOV).
c--   and i  is the inclination.  
c--
c--   Latitude must be limited to between the maximum values.

      cosi= COS((180.-i_inclin)/radeg)
      latlimit= i_inclin
      IF(i_inclin .GT. 90.) latlimit= 180. - i_inclin

      chari= sign(1.,i_inclin -90.)


c-- Load the table of weights to use in the cost function.

      DO i=0, 69
        szweight(i)= 1.d0
      END DO

      DO i=70, 89
        szweight(i)= (1./3./COS(float(i)/radeg))**2
      END DO


c-- Copy the FOV corner data from the /instdata/ to /geom/ variables

      DO i= 1, i_nscans
        DO j= 1, 4
          fov_lat(j,i)= i_fovcorn(1,j,i)
          fov_lon(j,i)= i_fovcorn(2,j,i)
        END DO
      END DO

      success= .TRUE.
      RETURN
      END