Graflib:

NCAR:

Open Graphics System

call gnglxx (wrkarray, length)

wrkarray(real array) - work array

length (int) - size of work array

call gogks (ier, idum)

ierr (int) - fortran unit number for error messages

idum(int) - dummy, not used.

Close Graphics System:

call geglxx (dummy)

dummy(int) - dummy variable

call gclks

Open device:

call gndevice (0)

where device can be:

fr80 - creates a cgm file.

tekx - writes to a Tektronixs device.

etc.

call gopwk (iwkid, iconid, iwktype)

iwkid (int) - workstation identifier (similar to vs_num in GRAFLIB gnvsxx)

iconid (int) - Depends on the device specified by iwktype. In most cases it is ignored. If used (e.g., CGM), it is the fortran logical unit that is used internally for writing to a file. Be careful not to use this same unit number elsewhere in your program for reads or writes.

iwktype (int) - workstation type

Examples:

Create a cgm binary file name gmeta:

call gopwk (1,2, 1)

Create a X11 window:

call gopwk (2, idum, 8)

Create a postscript file named gmeta1.ps:

call gopwk (1, idum, 20)

call gnvsxx (vs_num, ....)

vs_num (int) - view surface number
(see gopwk for how this identifier is defined in NCAR.)

... rest of arguments dependent on device opened.

call gacwk (iwkid)

iwkid - workstation identifier specified in gopwk.

Deactivate device:

call gdusxx (vs_num)

vs_num (int) - view surface number

call gdawk (iwkid)

iwkid - workstation identifier specified in gopwk.

Reactivate device:

call gcvsxx (vs_num)

vs_num (int) - view surface number

See "gacwk" for reconnecting the workstation.

call gevsxx (vs_num)

vs_num (int) - view surface number

call gclwk (iwkid)

iwkid - workstation identifier specified in gopwk

Misc. Routines:

call gsvsfn - FR80 file name.

call gsvskp - keep /give file status.

There is no correspondence in NCAR for these routines. Comment them out.

Graflib:

Define coordinate systems:

call gswd2d(type,xmin,xmin,ymin,ymax)
call gsvp2d(vxmin,vxmax,vymin,vymax)

xmin,... defines the user coordinate system.

vxmin,...defines where on the frame the data is plotted.

type defines the type of transformation

linlin$ or 0

loglin$ or 1

linlog$ or 2

loglog$ or 3

call set(vxmn,vxmx,vymn,vymx,xmn,xmx,ymn,ymx,itype)

vxmn,...defines where on the frame the data is plotted..

xmn,.... defines the user coordinate system.

itype(int) defines the transformation.
1 = linlin
2 = linlog
3 = loglin
4 = loglog

call gscpvs(x,y) set current position
call gscp2d (x,y) set current position

The safest way to specify coordinate positions in NCAR, is to use the coordinate system defined in the call to set (…,xmn,xmx,ymn,ymx,...). Thus, if you want to duplicate gscp2d in GRAFLIB, you need to define both the limits of the data and the area to plot in to be from 0.0 to 1.0. The best way to do this in NCAR is to save the current coordinate system, change to the unity system, and then when done, reset to the original coordinates.

call getset ( x1,x2,x3,x4,x5,x6,x7,x8,ix)

call set (0.,1.,0.,1.,0.,1.,0.,1.,1)

< do your calls >

call set(x1,x2,x3,x4,x5,x6,x7,x8,ix)

Misc. Routines:

call gswdnx(nice)
call gswdny(nice)

nice (int)
0=nothing,
1=modify the limits of the user coords. to create "nice" values.

The nice feature in GRAFLIB modified the limits of your data so that when you called the routine to plot a grid, nice tick labels were produced starting at the endpoints of the grid. If you use thc method described in "Grids and Labels", a somewhat similar capability is available. See "Grids and Labels" for how to do this.

Graflib:

Paste:

Advance Frame:

call gxglfr(dummy)
call gxvsfr (ilvs)

ilvs (int) - advances frame for only the specified view surface.

call frame

advances to the next frame for all activated workstations.

Clipping:

call gsvtcl(iclip)

iclip (int)
"on$" or 1 to clip
"off$" or 0 to not clip

call gsclip (clip)

clip (int )

0 = clip

1 = no clip

call gxglup(dummy)
call gxvsup(ilvs)

ilvs (int) - updates only the specified view surface.

call sflush

Updates all open workstations.

gsuskp(<get/keep status of file>)

This routine has no correspondence in NCAR. Comment it out.

Graflib:

NCAR:

Draw a grid:

call gpgr80 (type)

type - type of transformation
linlin$ or 0
loglin$ or 1
linlog$ or 2
loglog$ or 3

See the examples at the end of this section for how to duplicate drawing grids, labels, etc.

Paramter setting routines:

call setgr("xx=value, ...$")

where xx is the parameter to change, and value is the new value.

or

Individual parameter setting routines

call agsetc ("name", "char_value")
call agseti ("name", int_value)
call agsetr ("name", real_value)

Depending on the value of the parameter (char, int, or real), one of the above routines is called. NCAR has numerous parameters to be changed. For a complete list of the parameters to change, do a "man autograph_params".

xm, xn (int) - type of major and minor lines for x.

or call gsgrxm(value) and call gsgrnx(value)

where value: Ticks = 0, Grid = 1

call agsetr ('BACKGROUND.' , val)

val (real) where:

1.0 = perimeter background (default)

2.0 = grid background

3.0 = half-axis background

4.0 = no background.

ym, yn (int) - type of major and minor lines for y.

or call gsgrym(value) and call gsgrny(value)

where value is an integer: Ticks = 0, Grid = 1

Defined by the agsetr call immediately above.

nx, ny (int) - number of minor lines between majors for x and y.

 

or call gsgrny (value) and gsgrny (value)

 

 call agsetr ('AXIS/TICKS/MINOR/SPACING.).

value)

where value (real) means:

-1. Suppresses minor ticks

>=1. Specifies the number of minor ticks

0. specifies the NCAR will choose the number of minor ticks.

fx, fy - format of labels on x and y axes (default is scientific notation). The format is of the form "x$", where x can be S, F, E, or I optionally followed by the precision. The '$' is mandatory. Default is 'S$'

or

gsgrfx (format) or gsgrfy(format)

XLF, YLF(char) - format of labels on x and y. The format is of the form '(xprec)' where x can be E, F, G, and I followed by the precision. The character string must be surrounded by parenthesis, ala FORTRAN format. Default format is '(E10.3)'

call gswdnx(nice)
call gswdny(nice)

nice (int)
0=nothing,
1=modify user coords. to create "nice" values.

 

call agsetr ('X/NICE.', value)

call agsetr ('Y/NICE.', value)

where value (real)

-1. x/y forces major-tick positions at the endpoints of the axis (default - similar to GRAFLIB's nice).

+1. Forces major tick positions at the endpoints of the top (x) or right (y) axis.

0.0 specifies that the major ticks are not forced at the ends of either axis (same as GRAFLIB's no nice).

Graflib:

NCAR:

Plot Text:

call gptx2d ("string$")

where string contains the characters to plot. The string must end with a "$".

before the call to gptx2d, calls could be made to:
call gstxan(angle) angle to plot text
call gstxjf(just,just) justification of string
call gstxno(size) # of characters per line.

call plchhq(xpos,ypos,chrs,size,ang,cntr)

xpos,ypos (real) defines the position to plot the text

chrs (string) is the character string to plot

size (real) is the approximate size of the character string. To get the approximate width as in GRAFLIB, pass 1/(characters per line). For example, GRAFLIB's default was 72 character per line, which would be 1/72. This will give you about the same height. The width will be different, because NCAR uses proportional fonts.

ang (real) is the angle to plot the string. Same as in GRAFLlB.

cntr (real) defines the position of the string, where:
-1. left justified; 1. right justified; 0. centered

Paramter setting routines:

Parameters can be set by calling settx, or the individual grid parameter routines.

call settx("xx=value, ...$")

where xx is the parameter to change, and value is the new value.

or

Individual text parameter setting routines

 Similar NCAR Parameters

Font:

'ift' - font number or

call gstxft(ifont)

where ifont is the number of the font to use.

 

call gstxfp(ifont,prec)

 

where ifont is the font to use.

prec (int) is the precision of the text

0 = good, 1= better, 2 = best

The higher the quality of prec, the larger the CGM or PostScript file.

Text line width:

See gslwsc in Lines and Markers:

'fan' (real) - angle to plot text
or call gstxan(angle)

 Passed as an argument in the call to plchhq.

'ixj', 'iyj' (int) - justification of string
or call gstxjf(just,just)

where ixj can be 'left', 'right', or 'center' and iyj can be 'top', 'bottom', or 'center'.

Passed as an argument in the call to plchhq.

'fno' (int) - number of characters per line, i.e., character width.

or call gstxno(size)

Passed as an argument in the call to plchhq.

Plot a line between two points:

 

call gpln2d (x1, x2, y2, y2)

x1,x2,y1,y2 defines the two points to draw a line between.

 

call dpline (x1, x2, y2, y2)

x1,x2,y1,y2 defines the two points to draw a line between.

Plot call gpcv2d(xarray,yarray,length)

xarray,yarray contains the array of x,y coordinates to plot of length.

 

call dpcurv(xarray,yarray,length)

xarray,yarray contains the array of x,y coordinates to plot of length.

Paramter setting routines:

Parameters can be set by calling setcv, or the individual grid parameter routines.

call setcv("xx=value, xx=value....$")

where xx is the two character string which specifies the parameter to change and value is the new value.

 

call dpsetc(pnam, cval)

call dpsetr(pnam, rval)

call dpseti(pnam, ival)

Depending on the type of the parameter to change (character, real, or integer) one of the above routines is called. pnam is a quoted two character string which specifies the parameter to change. xval is the new value.

 

Some Parameters & alternate calls

Similar NCAR Parameters

md (int) - 0=lines, 1=marker, 2=both

One way to draw markers is to call 'dpsetc ('DPT', 'value'), where value is the marker you want.

cn (int) - 0=disjointed lines, 1= connected

Disjointed lines are not supported. You will need to break up the calls.

lb (int) - 0= no labels, l=labels

tx (int) - 0= no text , address of text

(You still need to call gncvlb to initialize labels as well as setting the above parameters.)

or

call gncvlb (0)

call gscvlb (1)

call gscvtx (loc(text))

 

call dpsetc ('DPT', 'string')

String defines the solid lines by a '$', spaces by '_', and alpanumeric characters specify the label.

sl (int) - 0= no smoothing , l= smoothing

tn (?) tension factor.

(You still need to call gncvsl to initialize smoothing as well as setting the above parameters.)

or

call gncvsl (0)

call gscvsl (1)

 

call dpsetr ("TCS", tension)

tension is the tension on cubic splines.

<0. Smoothing off

>=0. Desired tension to be used on the cubic splines used to do the smoothing.

nc(int) - 0= off , nice orientation of labels

NCAR alreadys does a really nice orientation of labels. The parameter 'SAF' gives you quite a bit of control over orientation, but is different then the GRAFLIB nice factor.

call gslnst(value)

sets line style (int), which varies from 1-20, where 1=solid line.

call gsln (itype)

where itype can be:

l = solid

2 = dashed

3 = dotted

4 = dash-dotted

If you need more line styles, you can create almost any style you want by setting various parameters. Check out "man dashpack" or the example at the end of this section.

call gslnsz (size)

set line size (width), where 0. is the default.

call gslwsc (size)

Sets the line width scale factor. A value of 1.0 is the default line width. A value of 0.0 produces no lines, where 2.0 will double the line width on most devices. That is, the displayed line width is dependent on the hardware capabilities of the graphics device being used.

A simple way of going from graflib to ncar size is:

ncar_size = graf_size*10. + 1.

Contour on a Rectangular Grid (No Labels):

call gpctv1 (…)

 

 

! initialize contour

call cprect (…)

! turn off the labeling of contours

call cpseti ('LLP - Line Label Position', 0)

call cpback(…) ! draw background

call cpcldr (…) ! draw contour lines

call gpcts1(…)

 

 

call cprect (…)

! set tension using cubic splines under tension

call cpsetr ('T2D - Tension', .0001)

call cpback (…)

call cpcldr (…)

Rectangular Mesh with X & Y spacing supplied (No Labels):

call gpctv2(…)

 

! X,Y must be strictly increasing

! 'cpsps2' interpolates from an array of data on

! rectangular mesh with X & Y supplied

! to an array of data on a "dense" rectangular

! grid and contour "dense" array

call cpsps2 (…)

call cpseti('LLP - Line Label Position',0)

call cpback (…)

call cpcldr (…)

or

call cpsetr ('XC1 - X COORDINATE AT I = 1',1.)

call cpsetr ('XCM - X COORDINATE AT I =M',<float(n)>)

call cpsetr ('YC1 - Y COORDINATE AT J = 1',1.) call cpsetr ('YCN - Y COORDINATE AT J = N',<float(n)>)

! NOTE: include the routine CPMPXY at the

!end of this section

call cpseti('MAP -- Mapping Function',3)

call cpseti('LLP - Line Label Position',0)

call cprect (…)

call cpback (…)

call cpcldr (…)

Smooth Labeled Contours for Rectangular Mesh

call gpcts2 (…)

.

 

call cpsps2 (…)

call cpsetr('T2D - Tension 2-dim spline',.0001)

call cpback (…)

call cpcldr (…)

or

call cpsetr ('XC1 - X COORDINATE AT I = 1',1.)

call cpsetr ('XCM - X COORDINATE AT I = M',<float(n)>)

call cpsetr ('YC1 - Y COORDINATE AT J = 1',1.)

call cpsetr ('YCN - Y COORDINATE AT J = N',<float(n)>)

! NOTE: include the routine CPMPXY at the

!end of this section

call cpseti('MAP -- Mapping Function',3)

call cpsetr('T2D - Tension 2-dim spline',.0001)

call cprect (…)

call cpback (…)

call cpcldr (…)

Contour from random (irregular) mesh:

gpctv3 (…)

 

 

 

call cpsetr ('XC1 - X COORDINATE AT I = 1',1.)

call cpsetr ('XCM - X COORDINATE AT I = M',<float(n)>)

call cpsetr ('YC1 - Y COORDINATE AT J = 1',1.)

call cpsetr ('YCN - Y COORDINATE AT J = N',<float(n)>)

! NOTE: include the routine CPMPXY at the

!end of this section

call cpseti('MAP -- Mapping Flag',4)

call cpseti('LLP - Line Label Position',0)

call cprect (…)

call cpback (…)

call cpcldr (…)

Smooth Contour from random (irregular) mesh:

call gpcts3 (…)

call cpsetr ('XC1 - X COORDINATE AT I = 1',1.)

call cpsetr ('XCM - X COORDINATE AT I = M',<float(n)>)

call cpsetr ('YC1 - Y COORDINATE AT J = 1',1.)

call cpsetr ('YCN - Y COORDINATE AT J = N',<float(n)>)

! NOTE: include the routine CPMPXY at the

!end of this section

call cpseti('MAP -- Mapping Flag',4)

call cpsetr('T2D - Tension 2-dim spline',.0001)

call cprect (…)

call cpback (…)

call cpcldr (…)

call gpctvx (…)

! Call IDSFFT to obtain a regular grid of values

! on the fitted surface.

call idsfft (…)

call cpsetr ('XC1 - X COORDINATE AT I = 1',1.)

call cpsetr ('XCM - X COORDINATE AT I = M',<float(n)>)

call cpsetr ('YC1 - Y COORDINATE AT J = 1',1.)

call cpsetr ('YCN - Y COORDINATE AT J = N',<float(n)>)

call cpseti('LLP - Line Label Position',0)

call cprect (…)

call cpback (…)

call cpcldr (…)

call gpctsx ()

call idsfft (…)

call cpsetr ('XC1 - X COORDINATE AT I = 1',1.)

call cpsetr ('XCM - X COORDINATE AT I = M',<float(n)>)

call cpsetr ('YC1 - Y COORDINATE AT J = 1',1.)

call cpsetr ('YCN - Y COORDINATE AT J = N',<float(n)>)

call cprect (…)

call cpsetr('T2D - Tension 2-dim spline',.0001) Default 0.

call cpback (…)

call cpcldr (…)

Font:

'ilb' - contour label switch (0=off, 1=on)

 

call cpseti ('LLP ', iflag)

if iflag = 0, no labels. A value of 1 gives reasonable contour labels, a value of 3 is better.

 

fmn,fmx - min & max label for label mode 4. (See the description of imd immediately following.)

See the code segment in label mode 4 below.

imd - format of labels if ilb is turned on. The possible values of imd and their meanings are:

1 = label the contour lines alphabetically (A-Z)

2 = use E format

3 = use F format (default)

4 = use G format, where fmn & fmx define the range. Labels above and below fmn & fmx are plotting in E format. Otherwise, labels are plotted in F format.

5 = use integer format

0 = if 0, then use the text labels defined in itx

= the format to use, e.g., "xf.w"

To simulate the GRAFLIB behavior:

Alphabetical Labels:

     data clab/ 'A', 'B', 'C', 'D'/
     call cpseti ('LLP', 2)
     do I = 1, ncl 
        call cpseti ('PAI', I)
        call cpsetc ('LLT', clab(I))
        call cpseti ('CLU', 3)
     endo

E format labels:

     call cpseti ('NET', 0)

F format labels

default

G format

     call cpseti('LLP', 2)
     do i = 1,ncl
        call cpseti('PAI',I)
	call cpgetr('CLV', clv)
	if ( clv .gt. fmx or
           clv .lt. fmn) then
           call cpseti('NET',0)	
        endif
     enddo

Integer Format

     call cpseti ('SFS', -4)

Specified format

     character*10 clab
     call cpseti('LLP',2)
     do i = 1,ncl
        call cpseti('PAI',I)
        call cpgetr('CLV',clv)	
	write(clab,imd)clv
	call cpsetc('LLT',clab)
	call CPSETI('CLU',3)
     enddo

itx - the text for the labels if imd = 0. For example,

itx = "lab1:lab2:lab3$"

To simulate "lab1:lab2:lab3$" do:
     data clab/"lab1","lab2","lab3"/
     call cpseti('LLP',2)
     do i = 1,ncl
        call cpseti('PAI',I)
        call cpsetc('LLT', clab(i))
        call CPSETI('CLU',3)
     enddo

isl, ftn - smoothing switch (0=off, 1=on) and smoothing tension factor

call cpsetr ('T2D', ftn)

inc - nice orientation of labels

0 = use user's current justification

1= on

To use the current justification, do nothing.

To set justification do:

     call cpseti('LLP', 2)
     call cpseti ('LLO', 1)

fft - adjustment to the first solid pattern defined by fdl.

     call cpseti ('LLP', 2)
     call cpsetr ('RC1', fft*fdl)

fdl - length of lines between labels

     cpsetr ('RC2', fdl)
     cpsetr ('RC3',0.)

irt - rotate labels, if 1, the labels with follow the lines

To simulate the value of 1, do:

     call cpseti('LLP', 2)
     call cpseti ('LLO',1)

ipz, fpv - "phony" value switch and the value to ignore

0 = off

1 = on (when on, fpv is the value to ignore)

To turn off this feature:

call cpsetr ('SPV - Special value',0.)

To turn on this feature and supply the value to ignore:

call cpsetr ('SPV',fpv)

Note: In NCAR, you cannot ignore the value of 0.0.

izc, isv - zone centering switch and zone centering array mode.

izc = 0, switch is off, ignore isv

izc = 1, switch is on, and isv informs the main contour routines to adjust their indexing mode.

To simulate izc = 1,

! Adjust X,Y coordinates 
! Compute data values at the center of the grid box
     MM1=M-1
     JM1=N-1
     DO 20 J=1,N
     DO 10 I=1,MM1
	Z(I,J)=0.5*(Z(I,J)+Z(I+1,J))
 10   CONTINUE
     Z(M,J)=Z(MM1,J)
 20   CONTINUE

fsc - scale applied to contour levels before labeled

     call cpsetr ('SFS', fsc)

isk - step for labeled contour lines

     call cpseti ('LIU', isk)

fst - stagger labels from level to level

     cpseti ('LLP',2)
     do i = 1,ncl
       call cpseti('PAI',I)
       if (mod(i,3).eq.0) then
          call cpsetr 
              ('RC1',fst+fft*fdl)
       end if
     enddo