滑动平均法

,REFIRST)Р?END SUBROUTINEР!读取XYZ文件Р?SUBROUTINE READXYZ(INFILE,MLINE,XY,FIRST)Р?CHARACTER *(*)INFILEР?REAL FIRST(MLINE),XY(2,MLINE)Р?OPEN(10,FILE=INFILE,STATUS='OLD',ACCESS='SEQUENTIAL',Р *?FORM='FORMATTED')Р?DO I=1,MLINE,1Р READ(10,*) XY(1,I),XY(2,I),FIRST(I)Р?ENDDOР?CLOSE(10)Р?ENDSUBROUTINEР!对XYZ文件排序Р?SUBROUTINE ORDER_XYZ(XY,FIRST,MLINE)Р?REAL XY(2,MLINE),FIRST(MLINE)Р?INTEGER FLAGР?FLAG=1Р?DO WHILE(FLAG==1)Р?FLAG=0Р DO I=1,MLINE-1,1Р IF(XY(1,I+1)<XY(1,I)) THENР FLAG=1Р TEMP=XY(1,I)Р XY(1,I)=XY(1,I+1)Р XY(1,I+1)=TEMPР TEMP=XY(2,I)Р XY(2,I)=XY(2,I+1)Р XY(2,I+1)=TEMPР TEMP=FIRST(I)Р FIRST(I)=FIRST(I+1)Р FIRST(I+1)=TEMPР ELSEIF(XY(1,I+1)==XY(1,I) .AND. XY(2,I+1)<XY(2,I)) THENР FLAG=1Р TEMP=XY(2,I)Р XY(2,I)=XY(2,I+1)Р XY(2,I+1)=TEMPР TEMP=FIRST(I)Р FIRST(I)=FIRST(I+1)Р FIRST(I+1)=TEMP Р ENDIF