C     CHUONG TRINH DU BAO K/N XH SM BANG FF XAC SUAT
$LARGE
$DEBUG
      PROGRAM Plop
      DIMENSION T(350,8),R(350,8),TD(350,8),P(200),V(350),XT(8)
      DIMENSION TMIN(350),MA(80),TNX(350,24),Hg(350)
      REAL X(350,30),MAY(350),X1(5500)
      INTEGER MSM(350,4),NG(350,4),IDX1(2,30),IDX(30),TSMTT(350,8)
      CHARACTER*12 FN,FIN,FN1,Fts
  905 CONTINUE
      WRITE(*,'(A\)')' FILE SO LIEU S.MU='
      READ(*,'(A)')FN
C
      OPEN(1,FILE=FN)

      I=1
    5 READ(1,*,END=100)(NG(I,J),J=1,4)
      WRITE(*,*)(NG(I,J),J=1,4)
      READ(1,*)(T(I,j),j=1,8)
      READ(1,*)(TD(I,j),j=1,8)
      READ(1,*)(R(I,j),j=1,8)
      READ(1,*)P(I)
      READ(1,*)MAY(I)
      READ(1,*)H
      READ(1,*)H
      READ(1,*)H
      READ(1,*)H
      READ(1,'(a)')H
      READ(1,*)V(I)
      READ(1,*)(tnx(i,j),j=1,24)
C     READ(1,*)H
C     READ(1,*)H
      READ(1,*)H
      READ(1,*)(MSM(I,J),J=1,4)
C     READ(1,*)H
      READ(1,*)(TSMTT(I,J),J=1,8)
      I=I+1
      IF(I-350) 5,5,100
  100 N=I-1
      write(*,*)' Vao xong so lieu !'
      DO 105 I=1,N
      DO 105 J=1,8
      IF((T(I,J).LE.0.OR.R(I,J).LE.0.)) GOTO 110
      A=ALOG(100./R(I,J))
      B=T(I,J)/10.+273.
      TD(I,J)=5443.2*B/(5443.2+A*B)-273.
      GOTO 105
  110 TD(I,J)=-1
  105 CONTINUE
C      amsm=4
  275 CONTINUE
      WRITE(*,'(A\)')' TINH CHO T.H M+SM (4) HAY SM (1)?'
      READ(*,*) AMSM
      WRITE(*,'(A\)')' FILE IN RA='
      READ(*,'(A)')FIN
      OPEN(3,FILE=FIN,STATUS='NEW')
C      do 789 kk=1,9
C      read(4,*)a1,a2,v1,v2
C
C     WRITE(*,'(A\)')' TINH VOI CAP GIO (V1,V2):'
C     READ(*,*)V1,V2
C     WRITE(*,'(A\)')' TINH VOI CAP L.MAY TRONG KHOANG (A1,A2):'
C     READ(*,*) A1,A2
      N1=0
      DO 20 I=1,N-1
      IF(NG(I,3)-NG(I+1,3))20,15,20
   15 IF(NG(I,1)-NG(I+1,1)+1)20,30,20
   30 DO 40 J=1,8
      if(j.le.2)xt(j)=t(i,j)
      if(j.gt.2)xt(j)=t(i+1,j)
   40 continue
   50 SM=999.
      DO 200 JJ=1,24
      if(jj.gt.18)goto 201
      IF(TNX(I+1,JJ).LE.0) GOTO 200
      IF(SM.LE.TNX(I+1,JJ)) GOTO 200
      SM=TNX(I+1,JJ)
      goto 200
  201 if(tnx(i,jj).le.0)goto 200
      if(sm.le.tnx(i,jj))goto 200
      sm=tnx(i,jj)
  200 CONTINUE
      IF(SM.LE.0.) GOTO 20
C      IF(SM.GT.AMSM)GOTO 20
      IF(TD(I+1,1).LE.5.) GOTO 20
      IF(T(I+1,1).LE.50.) GOTO 20
      IF(R(I+1,1).LE.30.) GOTO 20
      IF(TNX(I+1,19).LE.0.) GOTO 20
      IF(TD(I,7).LE.5.) GOTO 20
      IF(T(I,7).LE.50.) GOTO 20
      IF(may(I+1).LT.0.) GOTO 20
      IF(V(I+1).LT.0.) GOTO 20
      IF(P(I+1).LE.0.) GOTO 20
C      DO 20 I=2,N
      IF(T(I,7).LT.0.OR.TD(I,7).LT.0)GOTO 20
      IF(T(I+1,1).LT.0.OR.TD(I+1,1).LT.0)GOTO 20
      N1=N1+1
      if(Sm.LE.amsm) goto 185
      X(N1,1)=0.
      GOTO 190
  185 X(N1,1)=1.
  190 X(N1,2)=T(I+1,1)/10
      X(N1,3)=TD(I+1,1)
      X(N1,4)=T(I,7)/10-T(I+1,1)/10
      X(N1,5)=TD(I,7)-TD(I+1,1)
      X(N1,6)=R(I,7)-R(I+1,1)
      X(N1,7)=P(I+1)-P(I)
      X(N1,8)=MAY(I+1)
C      X(N1,9)=Hg(I+1)
      X(N1,9)=V(I+1)
      write(3,65)(x(n1,jj),jj=1,9)
   65 FORMAT(1X,10F7.1)
   20 CONTINUE
      if(n1.le.3)goto 2110
C      IF(AMSM.EQ.4)WRITE(3,392)
C      IF(AMSM.EQ.1)WRITE(3,393)
  392 FORMAT(//30X,'Tinh cho ca t.hop M+SM')
  393 FORMAT(30X,'Chi tinh cho t.hop SM')
      goto 789
 2110 write(*,*)' Dung luong mau qua nho (<=3 )'
 789  continue
      close(3)
      WRITE(*,'(A\)')' CO TINH TIEP VOI T.HOP KHAC ? (0/1)'
      READ(*,*)IL
      IF(IL.EQ.1) GOTO 275
      close(1)
  790 WRITE(*,'(A\)')' TINH VOI FILE SL KHAC ? (0/1)'
      READ(*,*)IFILE
      IF(IFILE.EQ.1)GOTO 905
      STOP
      END