%     PLOT 1-D FINITE DIFFERENCE RESULT AND COMPARE WITH MOC SOLUTION
%     MOC FROM ASG12.M
      clear all
      clf reset
      LX=1000
      OOIP=53433.67
      SWR=0.0
      SOR=0.0
      KRWEP=0.1667
      KROEP=1.0
      VISW=1.0
      VISO=12.0
      M=(KRWEP*VISO)/(KROEP*VISW)
      MPV=OOIP*(1-SWR-SOR)
      EXP=1.5
      Q=20
      tcar=MPV/Q
      load prod.dat
      tmoc1=prod(:,1);
      cutmoc=prod(:,2);
      tmoc2=prod(:,3);
      ermoc=prod(:,4);
      load prod1.dat
      time=prod1(:,1)./tcar;
      Q1=prod1(:,5);
      Q2=prod1(:,6);
      cut=Q1./(Q1+Q2);
      er=-prod1(:,7)./MPV;
      figure(1)
      plot(tmoc1,cutmoc,':',time,cut,'-')
      axis([0 2 0 1])
      title('M=2.0, n=1.5, NX=10, ')
      xlabel('DIMENSIONLESS TIME')
      ylabel('WATER CUT')
      in=input('hit ENTER to continue');
      figure(2)
      clf reset
      plot(tmoc2,ermoc,':',time,er,'-')
      axis([0 2 0 1])
      title('M=2.0, n=1.5, NX=10 ')
      xlabel('DIMENSIONLESS TIME')
      ylabel('RECOVERY EFFICIENCY')
