GCONTR draws contour lines of a 2-d array using a technique which produces long, connected contour lines. It assumes that the points of the input array are equally spaced in each dimension. Several options for contouring are provided. When cs, the label character size is positive, the origin point (1,1) is in the lower-left corner and the point (i,j) in the array is plotted at,
xplot = (i-1)*xl yplot = (j-1)*ylIf cs is negative, the x and y values plotted (xp,yp) for the point (i,j) in the array are computed using the PLT3D transformation common block PLT3B (see PLT3D),
xp = a1 * (yl * j) + a2 * (xl * i) + a3 yp = b1 * (yl * j) + b2 * (xl * i) + b4where the vertical height (z) is zero. In this case, xl and yl should be set to 1.0 and PLT3D should be called before GCONTR. GCONTR is used in the MASTER routine LCNTR and CVAX3DX.
CALL GCONTR(z,ndx,ndy,nx,ny,xl,yl,cv,nv,zm,iw,n,cs,m,i,ic,il) z (R): 2-d array of values dimensioned z(nx,ny) ndx,ndy (I): dimensions of data array nx,ny (I): number of points to use in array xl,yl (R): axis length scale factors (inches/array index) cv (R): array of contour levels dimensioned cl(nv) nv (I): number of contour levels (note: if nv < 0 then only one contour level is used. It will be labeled with the abs(nv)'th symbol) zm (R): maximum value of z for consideration. A z value which exceeds this value will be ignored. The cell edges which include this point will not be included in contouring. iw (I): workspace dimensioned at least (2*nx*ny*nv+1)/31 n (I): contour labeling option < 0 label with contour value (number with n digits to the right of the decimal point) = 0 no labelling of contours > 0 label with alphabet (nl should be less than 26) cs (R): size of labels < 0 : plot contours using PLT3B transformation xl and yl should then be set to 1.0 (see notes) > 0 : normal location specification m (I): minimum number of cells crossed by contour in order for contour to be labeled i (I): color and line type flag = 0 color and line type arrays not used = 1 color array used = 2 line type array used = 3 color and line type array used ic (I): color list for each contour (only requied for i>0) l (I): line type list for contours (accessed only for i>1)