Name: GRAFMP 24 The routine GRAFMP plots a geographical axis system. The call is: CALL GRAFMP (XA, XE, XOR, XSTP, YA, YE, YOR, YSTP) or: void graf3 ( float xa, float xe, float xor, float xstp, float ya, float ye, float yor, float ystp); XA, XE are the lower and upper limits of the X-axis. XOR, XSTP are the first X-axis label and the step between labels. YA, YE are the lower and upper limits of the Y-axis. YOR, YSTP are the first Y-axis label and the step between labels. Notes: - GRAFMP must be called from level 1 and sets the level to 2. - The axes must be linear and scaled in ascending order. In general, X-axes must be scaled between -180 and 180 degrees and Y-axes between -90 and 90 degrees. - For elliptical projections, the plotting of an axis system will be suppressed. This will also be done for azimuthal projections with YE - YA > 90. - The statement CALL GRIDMP (I, J) overlays an axis system with a grid of longitude and latitude where I and J are the number of grid lines between la- bels in the X- and Y-direction. Name: GRIDMP 24 The routine GRIDMP overlays an axis system with a longitude and latitude grid. The call is: CALL GRIDMP (I, J) level 1, 2, 3 or: void gridmp (int i, int j); I, J are the number of grid lines between labels in the X- and Y-direction. Name: XAXMAP 24 The routine XAXMAP plots a secondary X-axis. The call is: CALL XAXMAP (A, B, OR, STEP, CSTR, NT, NY) level 2 or: void xaxmap (float a, float b, float or, float step, char *cstr, int nt, int ny); A, B are the lower and upper limits of the X-axis. OR, STEP are the first label and the step between labels. CSTR is a character string containing the axis name. NT indicates how ticks, labels and the axis name are plotted. If NT = 0, they are plotted in a clock- wise direction. If NT = 1, they are plotted in a counter-clockwise direction. NY defines the horizontal position of the X-axis. A secondary axis must be located inside the axis system. Name: YAXMAP 24 The routine YAXMAP plots a secondary Y-axis. The call is: CALL YAXMAP (A, B, OR, STEP, CSTR, NT, NX) level 2 or: void yaxmap (float a, float b, float or, float step, char *cstr, int nt, int nx); A, B are the lower and upper limits of the Y-axis. OR, STEP are the first label and the step between labels. CSTR is a character string containing the axis name. NT indicates how ticks, labels and the axis name are plotted. If NT = 0, they are plotted in a clock- wise direction. If NT = 1, they are plotted in a counter-clockwise direction. NX defines the vertical position of the Y-axis. A secondary axis must be located inside the axis system. Name: PROJCT 24 The routine PROJCT selects the geographical projection. The call is: CALL PROJCT (CTYPE) level 1 or: void projct (char *ctype); CTYPE is a character string defining the projection. = 'CYLI' defines a cylindrical equidistant projection. = 'MERC' selects the Mercator projection. = 'EQUA' defines a cylindrical equal-area projection. = 'HAMM' defines the elliptical projection of Hammer. = 'AITO' defines the elliptical projection of Aitoff. = 'WINK' defines the elliptical projection of Winkel. = 'SANS' defines the elliptical Mercator-Sanson projection. = 'CONI' defines a conical equidistant projection. = 'ALBE' defines a conical equal-area projection (Albert). = 'CONF' defines a conical conformal projection. = 'AZIM' defines an azimuthal equidistant projection. = 'LAMB' defines an azimuthal equal-area projection. = 'STER' defines an azimuthal stereographic projection. = 'ORTH' defines an azimuthal orthographic projection. Default: CTYPE = 'CYLI'. Notes: - For cylindrical equidistant projections, the sca- ling parameters in GRAFMP must be in the range: -540 <= XA <= XE <= 540 -180 <= YA <= YE <= 180 For Mercator procektions: -540 <= XA <= XE <= 540 - 85 <= YA <= YE <= 85 For cylindrical equal-area projections: -540 <= XA <= XE <= 540 - 90 <= YA <= YE <= 90 For elliptical projections: -180 <= XA <= XE <= 180 - 90 <= YA <= YE <= 90 For conical projections: -180 <= XA <= XE <= 180 0 <= YA <= YE <= 90 or - 90 <= YA <= YE <= 0 For azimuthal projections with YE - YA > 90, the hemisphere around the map pole is projected onto a circle. Therefore, the hemisphere can be selected with the map pole. The plotting of the axis system will be suppressed. If YE - YA <= 90, the part of the globe defined by the axis scaling is projected onto a rectangle. The map pole will be set by GRAFMP to ((XA+XE)/2, (YE+YA)/2). The scaling parameters must be in the range: -180 <= XA <= XE <= 180 and XE - XA <= 180 - 90 <= YA <= YE <= 90 - For all projections except the default projection, longitude and latitude coordinates will be projec- ted with the same scaling factor for the X- and Y-axis. The scaling factor is determined by the scaling of the Y-axis while the scaling of the X-axis is used to centre the map. The longitude (XA+XE)/2 always lies in the centre of the axis system. Name: WORLD 24 The routine WORLD plots coastlines and political borders in low resolution. The call is: CALL WORLD level 2 or: void world (); Name: SHDMAP 24 The routine SHDMAP plots shaded continents. The call is: CALL SHDMAP (CMAP) level 2 or: void shdmap (char *cmap); CMAP is a character string defining the continent. = 'AFRI' means Africa. = 'ANTA' means the Antarctic. = 'AUST' means Australia. = 'EURA' means Europe and Asia. = 'NORT' means North America. = 'SOUT' means South America. = 'ALL' means all continents. Note: Shading patterns can be selected with SHDPAT and MYPAT. Colours can be defined with COLOR and SET- CLR. Name: SHDEUR 24 The routine SHDEUR plots shaded European countries. The call is: CALL SHDEUR (IRAY, IPRAY, ICRAY, N) level 2 or: void shdeur (int *iray, long *ipray, int *icray, int n); IRAY is an integer array containing the countries to be shaded. IRAY can have the following values: 1: Albania 13: Iceland 24: Portugal 2: Andorra 14: Italy 25: Romania 3: Belgium 15: Yugoslavia 26: Sweden 4: Bulgaria 16: Liechtenstein 27: Switzerland 5: Germany 17: Luxembourg 28: Spain 6: Denmark 18: Malta 29: CSFR 8: England/GB 19: Netherlands 30: Turkey 9: Finland 20: North Ireland 31: USSR 10: France 21: Norway 32: Hungary 11: Greece 22: Austria 12: Ireland 23: Poland IPRAY is an integer array containing shading patterns. ICRAY is an integer array containing colour numbers. N is the number of countries to be shaded. Name: CURVMP 24 The routine CURVMP plots curves through data points or marks them with symbols. The call is: CALL CURVMP (XRAY, YRAY, N) level 2 or: void curvmp (float *xray, float *yray, int n); XRAY, YRAY are real arrays containing the data points. N is the number of data points. Notes: - CURVMP is similar to CURVE except that only a linear interpolation can be used. - In general, a line between two points on the globe will not be projected as a straight line. There- fore, CURVMP interpolates lines linearly by small steps. An alternate plotting mode can be set with MAPMOD. Name: MAPBAS 24 M A P B A S The routine MAPBAS defines the map data file used in the rou- tine WORLD. A DISLIN map file and some external map files compi- led by Paul Wessel can be used. The external map files can be copied via FTP anonymous from the servers ftp://ftp.ngdc.noaa.gov/MGG/shorelines/ ftp://kiawe.soest.hawaii.edu/pub/wessel/gshhs/ The external map files 'gshhs_l.b', 'gshhs_i.b', 'gshhs_h.b' and 'gshhs_f.b' must be copied to the map subdirectory of the DISLIN directory. The call is: CALL MAPBAS (CBAS) level 1, 2 CBAS is a character string defining the map data file: = 'DISLIN' defines the DISLIN base map. = 'GSHL' defines 'gshhs_l.b' as base map. = 'GSHI' defines 'gshhs_i.b' as base map. = 'GSHH' defines 'gshhs_h.b' as base map. = 'GSHF' defines 'gshhs_f.b' as base map. Default: CBAS = 'DISLIN'. Name: MAPLEV 24 The routine MAPLEV defines land or lake coordinates for WORLD if the external map files from Paul Wessel are used. The call is: CALL MAPLEV (COPT) level 1, 2 COPT is a character string that can have the values 'BOTH', 'LAND' and 'LAKE'. Default: COPT = 'BOTH'. Name: MAPPOL 24 MAPPOL defines the map pole used for azimuthal projections. The call is: CALL MAPPOL (XPOL, YPOL) level 1 or: void mappol (float xpol, float ypol); XPOL, YPOL are the longitude and latitude coordinates in de- grees where: -180 <= XPOL <= 180 and - 90 <= YPOL <= 90 Default: (0., 0.) Note: For an azimuthal projection with YE - YA <= 90, the map pole will be set by GRAFMP to ((XA+XE)/2, (YA+YE)/2). Name: MAPSPH 24 For an azimuthal projection with YE - YA > 90, DISLIN automati- cally projects a hemisphere around the map pole on a circle. The hemisphere can be reduced using MAPSPH. The call is: CALL MAPSPH (XRAD) level 1 or: void mapshd (float xrad); XRAD defines the region around the map pole that will be projected onto a circle (0 < XRAD <= 90). Default: XRAD = 90. Name: MAPREF 24 The routine MAPREF defines, for conical projections, two lati- tudes where the cone intersects or touches the globe. The call is: CALL MAPREF (YLW, YUP) level 1 or: void mapref (float ylw, float yup); YLW, YUP are the lower and upper latitudes where: 0 <= YLW <= YUP <= 90 or - 90 <= YLW <= YUP <= 0 Default: YLW = YA + 0.25 * (YE - YA) YUP = YA + 0.75 * (YE - YA) Note: YLW and YUP can have identical values and lie out- side of the axis scaling. Name: MAPMOD 24 The routine MAPMOD determines how data points will be connected by CURVMP. The call is: CALL MAPMOD (CMODE) level 1, 2, 3 or: void mapmod (char *cmode); CMODE is a character string defining the connection mode. = 'STRAIGHT' defines straight lines. = 'INTER' means that lines will be interpolated linearly. Default: CMODE = 'INTER'. Name: POS2PT 24 The routine POS2PT converts map coordinates to plot coordinates. The call is: CALL POS2PT (XLONG, YLAT, XP, YP) level 2, 3 or: void pos2pt (float xlong, float ylat, float *xp, float *yp); XLONG, YLAT are the map coordinates in degrees. XP, YP are the plot coordinates calculated by POS2PT. Name: X2DPOS 24 The funtion X2DPOS converts map coordinates to plot coordinates. The call is: XP = X2DPOS (XLONG, YLAT) level 2, 3 or: float x2dpos (float xlong, float ylat); XP is the plot coordinate calculated by X2DPOS. XLONG, YLAT are the map coordinates in degrees. Name: Y2DPOS 24 The funtion Y2DPOS converts map coordinates to plot coordinates. The call is: YP = Y2DPOS (XLONG, YLAT) level 2, 3 or: float y2dpos (float xlong, float ylat); YP is the plot coordinate calculated by Y2DPOS. XLONG, YLAT are the map coordinates in degrees.