W3FB07

The W3FB07 routine converts the coordinates of a location on earth
           given in a grid coordinate system overlaid on a polar
           stereographic map projection true at 60 degrees N or S
           latitude to the natural coordinate system of
           latitude/longitude.

USAGE:  CALL W3FB07(XI,XJ,ALAT1,ALON1,DX,ALONV,ALAT,ALON)

   Input argument list:
     XI       - I COORDINATE OF THE POINT  REAL*4
     XJ       - J COORDINATE OF THE POINT  REAL*4
     ALAT1    - LATITUDE  OF LOWER LEFT POINT OF GRID (POINT 1,1)
                LATITUDE <0 FOR SOUTHERN HEMISPHERE; REAL*4
     ALON1    - LONGITUDE OF LOWER LEFT POINT OF GRID (POINT 1,1)
                EAST LONGITUDE USED THROUGHOUT; REAL*4
     DX       - MESH LENGTH OF GRID IN METERS AT 60 DEG LAT
                MUST BE SET NEGATIVE IF USING
                SOUTHERN HEMISPHERE PROJECTION; REAL*4
                   190500.0 LFM GRID,
                   381000.0 NH PE GRID, -381000.0 SH PE GRID, ETC.
     ALONV    - THE ORIENTATION OF THE GRID.  I.E.,
                THE EAST LONGITUDE VALUE OF THE VERTICAL MERIDIAN
                WHICH IS PARALLEL TO THE Y-AXIS (OR COLUMNS OF
                THE GRID) ALONG WHICH LATITUDE INCREASES AS
                THE Y-COORDINATE INCREASES.  REAL*4
                   FOR EXAMPLE:
                   255.0 FOR LFM GRID,
                   280.0 NH PE GRID, 100.0 SH PE GRID, ETC.

   Output argument list:
     ALAT     - LATITUDE IN DEGREES (NEGATIVE IN SOUTHERN HEMI.)
     ALON     - EAST LONGITUDE IN DEGREES, 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