REAL AB, BC,CD,DE,EF,FA,BE,BF
      REAL D1,D2,D3
      AB=-.047
      FA=.01
      BF=-.01              
      
     
      EF=-.048
      BE=.01
      
      BC=-.055
      CD=-.07
      DE=.1
      OPEN(UNIT=2, FILE='outwat1.txt',STATUS='UNKNOWN')  
      
      
      
      
    
      PRINT*,'=======================================================' 
      PRINT*,'|                  Problem #  1                       |' 
      PRINT*,'|                   FINAL EXAM                        |'
      PRINT*,'|  CIVE 6383 Adv.Planning and Des. of Water Resource  |'
      PRINT*,'|              By  Pracha Vorarittinapa               |'
      PRINT*,'======================================================='   
              
      WRITE(2,400)
  400 FORMAT('=======================================================')  
      WRITE(2,401)
  401 FORMAT('|                  Problem #  1                       |')
      WRITE(2,402)
  402 FORMAT('|                     RESULT                          |') 
     
      WRITE(2,403)
  403 FORMAT('|              By  Pracha Vorarittinapa               |')
      WRITE(2,400)    
  404 FORMAT('')    
      WRITE(2,404)
      WRITE(2,404)     
      WRITE(2,600)  
  600 FORMAT('-------------------------------------') 
       
      WRITE(2,700)
  700 FORMAT('       PIPE        FLOW(m^3/s)  ')
      WRITE(2,600)
      
   10 D1=(-1/1.85)*(4633*(ABS(FA))**1.85*FA/ABS(FA)+3470*(ABS(AB))**
     +1.85*AB/ABS(AB)+1286*(ABS(BF))**1.85*BF/ABS(BF))/(4633*(ABS(FA))
     +**.85+3470*(ABS(AB))**.85+1286*(ABS(BF))**.85) 
     
      D2=(-1/1.85)*(1286*(ABS(BF))**1.85*BF/ABS(BF)+23510*(ABS(EF))**
     +1.85*EF/ABS(EF)+47020*(ABS(BE))**1.85*BE/ABS(BE))/(1286*(ABS(BF))
     +**.85+23510*(ABS(EF))**.85+47020*(ABS(BE))**.85)       
     
      D3=(-1/1.85)*(47020*(ABS(BE))**1.85*BE/ABS(BE)+43.99*(ABS(DE))**
     +1.85*DE/ABS(DE)+482.3*(ABS(BC))**1.85*BC/ABS(BC)+7.633*(ABS(CD))**
     +1.85*CD/ABS(CD))/(47020*(ABS(BE))**.85+43.99*(ABS(DE))**.85+482.3*
     +(ABS(BC))**.85+7.633*(ABS(CD))**.85)  
      
      
      AB=AB+D1   
      FA=FA+D1
      BF=BF+D1+D2
      
      EF=EF+D2
      BE=BE+D2+D3
      
      BC=BC+D3
      CD=CD+D3
      DE=DE+D3
      IF((ABS(D1).GT.0.000001).OR.(ABS(D2).GT.0.000001).OR.
     +(ABS(D3).GT.0.0000001)) THEN    
      
      
      GOTO 10
      END IF
      
      PRINT*,'AB = ',AB
      PRINT*,'BC = ',BC
      PRINT*,'CD = ',CD
      PRINT*,'DE = ',DE
      PRINT*,'EF = ',EF
      PRINT*,'FA = ',FA
      PRINT*,'BE = ',BE
      PRINT*,'BF = ',BF
      
      WRITE(2,500) AB  
  500 FORMAT(5X,'AB       ',F10.5) 
      WRITE(2,501) BC  
  501 FORMAT(5X,'BC       ',F10.5) 
      WRITE(2,502) CD  
  502 FORMAT(5X,'CD       ',F10.5) 
      WRITE(2,503) DE 
  503 FORMAT(5X,'DE       ',F10.5) 
      WRITE(2,504) EF 
  504 FORMAT(5X,'EF       ',F10.5) 
      WRITE(2,505) FA 
  505 FORMAT(5X,'FA       ',F10.5) 
      WRITE(2,506) BE 
  506 FORMAT(5X,'BE       ',F10.5) 
      WRITE(2,507) BF 
  507 FORMAT(5X,'BF       ',F10.5)
      WRITE(2,600)
      PRINT*,'D1=', D1  
      PRINT*,'D2=', D2  
      PRINT*,'D3=', D3  
      END