C C [HEMISHERICAL VIEW] (GRAEBER) C PROGRAM TO PRODUCE A PERSPECTIVE DRAWING OF TWO SINGLE VALUED C FUNCTIONS DEFINED OVER A HEMISPHERICAL SURFACE, IN SUCH A WAY C AS TO EXHIBIT THE ARCS OF LTITUDE AND LONGITUDE. C F1(NP,NT) ARRAY CONTAINING FIRST FUNCTION C F2(NP,NT) ARRAY CONTAINING SECOND FUNCTION C NP NUMBER (=4*N+1) OF POINTS ON ONE LATITUDE C NT NUMBER OF POINTS ON ONE LONGITUDE C O(3,3) ORTHOGONAL ROTATION MATRIX C PL PEN MOVEMENT SUBROUTINE C [14-MARCH-78] C EXTERNAL PL LOGICAL*1 B DIMENSION F1(1),F2(1),O(3,3) COMMON/VIS/N0 C CALL VISNP(PH,TH,JP,IT,NP,NT,O) NH=(NP-1)/2 IF(TH.GT.(1.57079)) GOTO 10 I1=MIN0(IT+1,NT) I2=NT GOTO 12 10 I1=1 I2=MAX0(IT-1,1) 12 J1=JP J2=JP+NH J3=JP-NH-1 J4=JP CALL VISSP(R,P,1.0,TH,PH+0.05,O) B=((-0.25).LT.P).AND.(P.LE.(0.25)) N1=0 N2=0 CALL VISHS(F1,F2,J1,J2,NP,I1,I2,NT,N1,N2,1,1,B,O,PL) CALL VISHS(F1,F2,J1,J2,NP,I1,I2,NT,N1,N2,1,-1,B,O,PL) N1=0 N2=0 CALL VISHS(F1,F2,J3,J4,NP,I1,I2,NT,N1,N2,-1,1,(.NOT.B),O,PL) CALL VISHS(F1,F2,J3,J4,NP,I1,I2,NT,N1,N2,-1,-1,(.NOT.B),O,PL) C RETURN END