The W3FB00 routine converts the coordinates of a location on earth
from the natural coordinate system of latitude/longitude
to the grid (i,j) coordinate system overlaid on the polar
stereographic map projection true at 60 N.
USAGE: CALL W3FB00 (ALAT, ALONG, XMESHL, XI, XJ)
ALAT LATITUDE IN DEG. (-20.0(S. HEMIS) ) ALAT ) 90.0)
ALONG WEST LONGITUDE IN DEGREES
XMESHL MESH LENGTH OF GRID IN KILOMETERS AT 60N
XI I OF THE POINT RELATIVE TO NORTH POLE
XJ J OF THE POINT RELATIVE TO NORTH POLE
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. The subroutine grid is oriented so that longitude
80W is parallel to the gridlines of i=constant. the radius of
the earth is taken to be 6371.2 km.
All parameters in the call statement must be real.
This code will not vectorize on a cray. You will have put
it line to vectorize it.
||Library contains Fortran 90 decoder/encoder
routines for GRIB edition 1. (Fortran90)
Date posted: 2/22/2007