LNDMAP plots the earth land edge map using a user-supplied projection subroutine. This routine will plot the entire map surface. PLOT routine clipping option only parts of the map surface can be plotted if desired. Uses Fortran file unit 2.
CALL LNDMAP(proj,s) proj (EXTERNAL): user-supplied projection subroutine name s (R): longitude of left-mode edge of map (-180 to +180) proj is a subroutine with the call format CALL MAPPLT(long,lat,ip) long (R): shifted longitude in deg (0 to 360) where "zero" corresponds to the longitude "s" lat (R): lattitude in deg (-90 to 90) ip (I): "pen" control flag = 3: move to (long,lat) with pen up = 2: draw to (long,lat) with pen down
An example implementation of proj which uses a linear projection is:
SUBROUTINE MAPPLT(ALONG,ALAT,IP) DATA XLONG/0.02/ ! LONG. SCALE FACTOR (INCHES/DEG LONG) DATA YLAT/0.02/ ! LAT. SCALE FACTOR (INCHES/DEG LAT) X1=ALONG*XLONG Y1=(ALAT+90.)*YLAT CALL PLOT(X1,Y1,IP) RETURN END