W3FB08

The W3FB08 routine converts a location on earth given in
           the coordinate system of latitude/longitudei
           to an (i,j) coordinate system overlaid on a
           mercator map projection.

USAGE:  CALL W3FB08 (ALAT,ALON,ALAT1,ALON1,ALATIN,DX,XI,XJ)

   Input argument list:
      ALAT     - LATITUDE IN DEGREES (NEGATIVE IN SOUTHERN HEMIS)
      ALON     - EAST LONGITUDE IN DEGREES, REAL*4
      ALAT1    - LATITUDE  OF LOWER LEFT CORNER OF GRID (POINT (1,1))
      ALON1    - LONGITUDE OF LOWER LEFT CORNER OF GRID (POINT (1,1))
                 ALL REAL*4
      ALATIN   - THE LATITUDE AT WHICH THE MERCATOR CYLINDER
                 INTERSECTS THE EARTH
      DX       - MESH LENGTH OF GRID IN METERS AT ALATIN

   Output argument list:
      XI       - I COORDINATE OF THE POINT SPECIFIED BY ALAT, ALON
      XJ       - J COORDINATE OF THE POINT; BOTH REAL*4

REMARKS:
        Formulae and notation loosely based on Hoke, Hayes,
        and Renninger's " map projections and grid systems...",
        March 1981 AFGWC/TN-79/003

 
W3lib.tar Library contains Fortran 90 decoder/encoder routines for GRIB edition 1. (Fortran90)
Date posted: 2/22/2007