00001 /* zzphph.f -- translated by f2c (version 19961017).
00002    You must link the resulting object file with the libraries:
00003         -lf2c -lm   (in that order)
00004 */
00005 
00006 #include "f2c.h"
00007 
00008 /* Common Block Declarations */
00009 
00010 extern struct {
00011     real xpgmin, ypgmin, xpgmax, ypgmax, xclbot, yclbot, xcltop, ycltop, xbot,
00012              ybot, xtop, ytop, xmin, ymin, xmax, ymax;
00013     integer ixcoor, iycoor;
00014     real alphxx, betaxx, alphyy, betayy, tmajx, tminx, tmajy, tminy;
00015     integer majrx, minrx, majry, minry, isizx, isizy;
00016     real xphold, yphold;
00017 } zzzplt_;
00018 
00019 #define zzzplt_1 zzzplt_
00020 
00021 extern struct {
00022     real xphmax, yphmax;
00023     integer ixpmax, iypmax;
00024     real xpscal, ypscal;
00025     integer iflip, nplotr;
00026     char cfile[64];
00027 } zzpltr_;
00028 
00029 #define zzpltr_1 zzpltr_
00030 
00031 
00032 
00033 
00034 /* Subroutine */ int zzphph_(real *x1, real *y1, real *x2, real *y2)
00035 {
00036     extern /* Subroutine */ int zzmpli_(real *, real *, real *, real *);
00037     static real px1, px2, py1, py2;
00038 
00039 
00040 /*  Plot a physical coordinate line. */
00041 
00042 /*  Everything in this package eventually funnels into this routine, */
00043 /*  which depends on the hardware being drawn upon.  The device */
00044 /*  is chosen by the NPLOTR variable (in 'plotpak.inc' COMMON /ZZPLTR/): 
00045 */
00046 
00047 /*    1 = Los Alamos CGS SunView window */
00048 /*    2 = Los Alamos CGS Metafile */
00049 /*    3 = Write all line coordinates to an ASCII file for later work */
00050 /*    4 = Write all line coordinates to a Unix plot format file */
00051 /*    5 = Use Microsoft GRAPHICS.LIB */
00052 /*    6 = Write PostScript code to an output file */
00053 /*    7 = Call the C "memplot" routine */
00054 /* .......................................................................
00055  */
00056 /* +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
00057  */
00058 
00059 /*  Internal Data for PLOTPAK */
00060 
00061     if (zzpltr_1.iflip == 1) {
00062 /*  Flipped coordinates */
00063         px1 = zzpltr_1.xpscal * (zzpltr_1.yphmax - *y1);
00064         px2 = zzpltr_1.xpscal * (zzpltr_1.yphmax - *y2);
00065         py1 = zzpltr_1.ypscal * (zzpltr_1.xphmax - *x1);
00066         py2 = zzpltr_1.ypscal * (zzpltr_1.xphmax - *x2);
00067     } else if (zzpltr_1.iflip == 2) {
00068 /*  Flipped coordinates */
00069         px1 = zzpltr_1.xpscal * (zzpltr_1.yphmax - *y1);
00070         px2 = zzpltr_1.xpscal * (zzpltr_1.yphmax - *y2);
00071         py1 = zzpltr_1.xpscal * *x1;
00072         py2 = zzpltr_1.xpscal * *x2;
00073     } else {
00074 /*  Normal coordinates */
00075         px1 = zzpltr_1.xpscal * *x1;
00076         px2 = zzpltr_1.xpscal * *x2;
00077         py1 = zzpltr_1.ypscal * *y1;
00078         py2 = zzpltr_1.ypscal * *y2;
00079     }
00080 
00081 /*  Plot it! */
00082 
00083 /* cc      IF( NPLOTR.EQ.1 .OR. NPLOTR.EQ.2 )THEN */
00084 /* cc         CALL GMOVA2( PX1 , PY1 ) */
00085 /* cc         CALL GLINA2( PX2 , PY2 ) */
00086 /* cc      ENDIF */
00087 
00088 /* cc      IF( NPLOTR .EQ. 3 )THEN */
00089 /* cc         WRITE(99,101) PX1,PY1,PX2,PY2 */
00090 /* cc101      FORMAT(4(F6.4,',')) */
00091 /* cc      ENDIF */
00092 
00093 /* cc      IF( NPLOTR .EQ. 4 )THEN */
00094 /* cc         CALL ZZUPLI( NINT(PX1),NINT(PY1) , NINT(PX2),NINT(PY2) ) */
00095 /* cc      ENDIF */
00096 
00097 /* cc      IF( NPLOTR .EQ. 5 )THEN */
00098 /* cc         CALL ZZPCLI( PX1,PY1 , PX2,PY2 ) */
00099 /* cc      ENDIF */
00100 
00101 /* cc      IF( NPLOTR .EQ. 6 )THEN */
00102 /* cc         CALL ZZPSLI( NINT(PX1),NINT(PY1) , NINT(PX2),NINT(PY2) ) */
00103 /* cc      ENDIF */
00104 
00105     if (zzpltr_1.nplotr == 7) {
00106         zzmpli_(&px1, &py1, &px2, &py2);
00107     }
00108 
00109 /* cc      WRITE(*,999) PX1,PY1,PX2,PY2 */
00110 /* cc999   FORMAT('ZZPHPH:',4(1X,1PG10.3)) */
00111 
00112     return 0;
00113 } /* zzphph_ */
00114