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)
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
ALAT LATITUDE IN DEG. (-20.0(S. HEMIS) < ALAT < 90.0)
ALONG WEST LONGITUDE IN DEGREES
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.
||Library contains Fortran 90 decoder/encoder
routines for GRIB edition 1. (Fortran90)
Date posted: 2/22/2007