C PROGRAM P O L A R P L O T C *************************** C C PROGRAM POLARPLOT IS DESIGNED FOR PLOTTING OF PARTICLE MOTION C DIAGRAMS FROM FILES LU3 GENERATED IN THE PROGRAM SYNTPL FOR DIF- C FERENT COMPONENTS OF DISPLACEMENT VECTOR. IT CAN BE ALSO USED C TO PLOT OTHER PARTICLE MOTION DIAGRAMS GENERATED ELSEWHERE BUT C STORED IN THE SAME FORM AS THOSE GENERATED IN PROGRAM SYNTPL. C C ***************************************************************** C C A SHORT DESCRIPTION OF THE PROGRAM C ********************************** C C PARTICLE MOTION DIAGRAMS ARE PLOTTED IN A FRAME, IN WHICH THE C X- AND Y-AXES CORRESPOND TO TWO SELECTED COMPONENTS OF DISPLA- C CEMENT VECTOR. THE DIAGRAMS MAY BE PLOTTED IN A NORMALIZED FORM, C WITH MAXIMUM AMPLITUDE REACHING THE BOUNDARIES OF THE FRAME FOR C EACH RECEIVER, OR WITH REAL CHANGES OF AMPLITUDES FROM RECEIVER C TO RECEIVER. A SELECTION OF RECEIVERS AT WHICH THE PARTICLE MO- C TION DIAGRAMS ARE TO BE COMPUTED, CAN BE PERFORMED. C CALCOMP PLOTTING ROUTINES ARE REQUIRED. C C ****************************************************************** C C INPUT DATA C ********** C C INPUT DATA CONSIST PARTIALLY OF THE DATA GENERATED BY A PROGRAM C SYNTPL AND STORED IN A FORMATTED FORM IN TWO FILES LU3 (EACH WITH C ONE COMPONENT), AND PARTIALLY OF THE ADDITIONAL INPUT DATA CONTROL- C LING PLOTTING, INTRODUCED DIRECTLY (CARD READER, TERMINAL). C C DATA STORED IN FILES LU3 C ************************ C IN A FILE LU3, THE COMPUTED RAY SYNTHETIC SEISMOGRAMS FOR INDI- C VIDUAL RECEIVER POSITIONS ARE STORED, TOGETHER WITH SOME OTHER C RELEVANT IMPORTANT INFORMATION. C THE DATA IN LU3 ARE STORED IN THE FOLLOWING ORDER: C 1) MTEXT FORMAT(17A4) C ARBITRARY ALPHANUMERIC TEXT DESCRIBING THE COMPUTATIONS. THIS C TEXT MAY APPEAR UNDER THE PLOTS. IT WAS READ IN IN THE PRO- C GRAM SEIS88. C 2) MDIST,MRED,MCOMP,ITPR,VRED,RSTEP,XSOUR,DT FORMAT(4I5,4F10.5) C MDIST... THE NUMBER OF RECEIVER POSITIONS. C MRED... MRED=0:NON REDUCED TRAVEL TIMES, C MRED=1: REDUCED TRAVEL TIMES. C MCOMP... MCOMP=0: VERTICAL COMPONENT. C MCOMP=1: RADIAL COMPONENT. C MCOMP=2: TRANSVERSE COMPONENT. C ITPR... ITPR=0: RECEIVERS ALONG THE EARTH'S SURFACE. C ITPR=1: RECEIVERS ALONG A VERTICAL PROFILE. C ITPR.GT.100: RECEIVERS ALONG AN INTERFACE. C VRED... REDUCTION VELOCITY C RSTEP... AVERAGE DIFFERENCE BETWEEN COORDINATES OF C NEIGHBOURING RECEIVER POSITIONS. C XSOUR... X-COORDINATE OF THE SOURCE. C DT... TIME STEP IN SYNTHETIC SEISMOGRAM. C 3) XMX,SMAXIM FORMAT(22X,F10.5,9X,E15.9) C XMX... COORDINATE OF OF THE RECEIVER AT WHICH THE C MAXIMUM AMPLITUDE 'SMAXIM' WAS RECORDED C (X-COORDINATE FOR RECEIVERS SITUATED ALONG THE C EARTH'S SURFACE OR ALONG AN INTERFACE, Z-COOR- C DINATE FOR RECEIVERS SITUATED ALONG A VERTICAL C PROFILE). C SMAXIM... THE MAXIMUM AMPLITUDE OF ALL TRACES. C THE FOLLOWING CARDS 4,5 ARE SUCCESSIVELY REPEATED FOR ALL C RECEIVER POSITIONS: C 4) XX,SMAX,TMIN,NPTS FORMAT(F10.5,E15.9,F10.5,I5) C XX... COORDINATE OF THE RECEIVER (X-COORDINATE WHEN C RECEIVERS ARE SITUATED ALONG THE EARTH'S SURFACE C OR ALONG AN INTERFACE AND Z-COORDINATE WHEN RE- C CEIVERS ARE SITUATED ALONG A VERTICAL PROFILE. C SMAX... MAXIMUM AMPLITUDE IN THE SYNTHETIC SEISMOGRAM C AT XX. C TMIN... THE TIME CORRESPONDING TO THE FIRST NON-ZERO C POINT IN THE SYNTHETIC SEISMOGRAM. C NPTS... NUMBER OF POINTS IN SYNTHETIC SEISMOGRAM. C 5) SYNTHETIC SEISMOGRAM FORMAT(20I4) C THE MAXIMUM LENGTH OF ONE SYNTHETIC SEISMOGRAM IS 3001 POINTS. C ONLY THE PARTS OF THE SEISMOGRAMS SITUATED BETWEEN THE FIRST C AND THE LAST NON-ZERO ARRIVAL ARE STORED. EACH SEISMOGRAM IS C STORED IN A NORMALIZED FORM, SEE THE FORMAT, WITH MAXIMUM C AMPLITUDE EQUAL 999. C C ADDITIONAL INPUT DATA C ********************* C C THE FOLLOWING DATA SHOULD BE INTRODUCED DIRECTLY (CARDS, TERMI- C NAL). THEY CONTROL THE PLOTTING OF SYNTHETIC SEISMOGRAMS. THE C DATA FROM LU3A AND LU3B (EACH FOR ONE COMPONENT) ARE READ IN C AFTER ALL THE ADDITIONAL INPUT DATA ARE READ IN. C C 1) ONE CARD C LU3A,LU3B,ISHIFT,IPRINT FORMAT(16I5) C LU3A,LU3B... NUMBERS OF FILES IN WHICH THE SYNTHETIC C SEISMOGRAMS FOR THE TWO CONSIDERED COMPONENTS C ARE STORED. C ISHIFT...THE SHIFT BETWEEN INDIVIDUAL PICTURES, IN CM. C DEFAULT VALUE, ISHIFT=10. C IPRINT...CONTROLS THE PRINTOUT OF SOME DATA ON THE LINE C PRINTER. SEE DETAILS IN THE SECTION ON OUTPUT C TABLES. C C 2) ONE CARD, VARIOUS SWITCHES C MCONT,MEPIC,NTICX,NTEXT FORMAT(16I5) C MCONT... MCONT=0: THE COMPUTATION TERMINATES. C MCONT=1: A NEW PLOT OF SYNTHETIC SEISMOGRAMS C WILL BE MADE FROM THE SAME FILE AS THE PRECE- C DING PLOT. C MCONT=-1: A NEW PLOT OF SYNTHETIC SEISMOGRAMS C WILL BE MADE FROM A NEW FILE. GO TO THE INPUT C DATA CARD NO.1. C MEPIC... CONTROLS THE SELECTION OF RECEIVER POSITIONS C MEPIC=0: NO SELECTION C MEPIC=1: SELECTION APPLIED, SEE CARD NO.3 C NTICX... THE NUMBER OF MARKED INTERVALS BETWEEN TWO AD- C JOINING TICS DENOTED BY COORDINATE VALUES ALONG C X- AND Y-AXES. C NTEXT... SELECTION OF THE TEXT WHICH SHOULD APPEAR UNDER C THE PLOTS: C NEXT=0: THE TEXT 'TEXT' STORED IN FILE LU3, IS C TO BE PLOTTED. C NTEXT=1: THE TEXT 'TXT' READ FROM INPUT DATA C CARD NO.5, SEE BELOW, IS TO BE PLOTTED. C C 3) ONE CARD, SELECTION OF RECEIVER POSITIONS. C INCLUDED ONLY WHEN MEPIC.NE.0. C NEPIC,(IEP(I),I=1,NEPIC) FORMAT(16I5) C NEPIC... NUMBER OF SELECTED RECEIVER POSITIONS AT WHICH C SYNTHETIC SEISMOGRAMS ARE TO BE CONSTRUCTED. C IEP(1),IEP(2),...,IEP(NEPIC)... SEQUENTIAL NUMBERS OF C SELECTED RECEIVER POSITIONS (IN INCREASING C ORDER). C NOTE: SIMILAR SELECTION CAN BE DONE ALREADY IN THE PROGRAM C SYNTPL. THE SELECTION HERE IS MADE AMONG THOSE RECEIVERS WHICH C REMAINED AFTER THE SELECTION IN SYNTPL. C C 4) ONE CARD, DESCRIPTION OF AXES IN THE PLOT C XLEN,DTICX,SC,TSTART,TFIN,AMP,B1... FORMAT(8F10.5) C XLEN... LENGTH OF X- AND Y-AXES (IN CM). THE MINIMUM C AND MAXIMUM VALUES ON BOTH AXES ARE -1.0 AND C 1.0 RESPECTIVELY. C DTICX... THE DISTANCE BETWEEN TWO ADJOINING TICS WHICH C ARE DENOTED BY CORRESPONDING COORDINATE VALUES C ON BOTH AXES. C SC... CONTROLS THE SCALES OF TICS AND ALPHANUMERICAL C TEXTS. FOR SC=1.0, THE TICS ARE 0.15 CM LONG C AND COORDINATES AND TEXTS DESCRIBING THE PLOTS C ARE 0.4 AND 0.45 CM HIGH, RESPECTIVELY. DEFAULT C VALUE, SC=1. C TSTART,TFIN... MINIMUM AND MAXIMUM VALUES OF THE TIME IN- C TERVAL FOR WHICH PARTICLE MOTION DIAGRAMS ARE C TO BE PLOTTED. THE VALUES APPEAR IN THE PLOT. C AMP,B1... PARAMETERS OF THE AMPLITUDE SCALING OF PARTICLE C MOTION DIAGRAMS. THEY CONTROL THE SCALING OF C PARTICLE MOTION DIAGRAMS. DEFAULT VALUE OF B1, C B1=1. C DENOTE DDX=0.5*XLEN. THEN THE FOLLOWING SYSTEMS C OF SCALING CAN BE USED: C C 1) FACTOR=B1*DDX/SMAXI, C WHEN AMP=0. C NORMALIZATION OF MAXIMUM AMPLITUDES IN ALL DIAGRAMS: C FOR B1=1, IN EACH FRAME THE PARTICLE MOTION DIAGRAM TOUCHES C THE FRAME OF THE DIAGRAM WITH ITS MAXIMUM AMPLITUDE. C 2) FACTOR=B1*DDX/SMAXIM, C WHEN AMP=-1, C REAL CHANGES OF AMPLITUDES BETWEEN DIAGRAMS: C FOR B1=1, THE PARTICLE MOTION DIAGRAM FOR THE RECEIVER WITH C MAXIMUM AMPLITUDE TOUCHES THE FRAME OF THE DIAGRAM. ALL THE C OTHER DIAGRAMS ARE WITHIN THEIR FRAMES. C 3) FACTOR=B1, C WHEN AMP=1, C REAL AMPLITUDES: MANUAL DETERMINATION OF THE SCALING FACTOR C BY SPECIFYING B1. C C 5) ONE CARD, ALPHANUMERICAL TEXT DESCIBING THE PLOTS. IT APPEARS C UNDER THE PLOTS WHEN NTEXT=1 IN INPUT DATA CARD NO.2. C TXT FORMAT(17A4) C C C **************************************************************** C C TERMINATION OF COMPUTATIONS C *************************** C C IF MCONT=1 IN CARD NO.2, THE CARDS 2-5 CAN BE REPEATED AN ARBI- C TRARY NUMBER OF TIMES TO PLOT PARTICLE MOTION DIAGRAMS WITH DIF- C FERENT AMPLITUDE SCALING, DIFFERENT TIME INTERVALS, ETC. C IF MCONT=-1 IN CARD NO.2, THE FILES LU3A AND LU3B ARE CLOSED, AND C NUMBERS OF NEW FILES ARE READ IN FROM INPUT DATA CARD NO.1. THE C FILES ARE OPENED AND THE COMPUTATIONS PROCEED AS DESCRIBED ABOVE. C THE COMPUTATION TERMINATES WHEN MCONT=0 IN THE INPUT DATA CARD C NO.2. C C ****************************************************************** C C OUTPUT ON A LINE PRINTER C ************************ C C ALL THE ADDITIONAL INPUT DATA ARE REPRODUCED ON THE LINE PRINTER. C THE PRINTOUT OF OTHER DATA IS CONTROLLED BY THE PARAMETER IPRINT, C SEE INPUT DATA CARD NO.1. C FOR IPRINT=0: ALSO THE DATA LU3/1,LU3/2 AND LU3/3 FROM BOTH FILES C LU3A AND LU3B ARE PRINTED. C FOR IPRINT=1: THE SAME AS FOR IPRINT=0. IN ADDITION, ONE LINE IS C PRINTED FOR EACH RECEIVER POSITION: C XX,SMAX,FACTOR,SFMAX C WHERE XX IS THE COORDINATE OF THE RECEIVER (X-COORDINATE FOR RECE- C IVERS SITUATED ALONG THE EARTH'S SURFACE OR ALONG AN INTERFACE, C Z-COORDINATE FOR RECEIVERS SITUATED ALONG A VERTICAL PROFILE), C SMAX IS THE MAXIMUM AMPLITUDE IN THE CORRESPONDING PARTICLE MOTION C DIAGRAM, FACTOR IS THE AMPLITUDE SCALING FACTOR CORRESPONDING TO C THE SPECIFIED RECEIVER, AND SFMAX IS THE MAXIMUM AMPLITUDE IN THE C PARTICLE MOTION DIAGRAM AFTER SCALING. C C ***************************************************************** C C