W3FB01

The W3FB01 routine converts the coordinates of a location from the
           grid(i,j) coordinate system overlaid on the polar
           stereographic map projection true at 60N to the natural
           coordinate system of latitude/longitude on the earth. 

USAGE:  CALL W3FB01 (XI, XJ, XMESHL, ALAT, ALONG)

   Input variables:
     XI            I OF THE POINT RELATIVE TO NORTH POLE
     XJ            J OF THE POINT RELATIVE TO NORTH POLE
     XMESHL        MESH LENGTH OF GRID IN KILOMETERS AT 60N

   Output variables:
     ALAT          LATITUDE IN DEG. (-20.0(S. HEMIS) < ALAT < 90.0)
     ALONG         WEST LONGITUDE IN DEGREES

REMARKS: 
        The grid used in this subroutine has its origin (i=0,j=0)
        at the north pole, so if the user's grid has its origin at a
        point other than the north pole, a translation is required to
        get i and j for input into w3fb01. The subroutine grid is
        oriented so that longitude 80w is parallel to gridlines of
        i=constant. The earth's radius is taken to be 6371.2 km.
        All parameters in the call statement must be reaL.

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