/[escript]/trunk/finley/test/python/runcoalgas.py
ViewVC logotype

Diff of /trunk/finley/test/python/runcoalgas.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3496 by gross, Tue Apr 5 04:12:08 2011 UTC revision 3497 by gross, Wed Apr 6 06:11:56 2011 UTC
# Line 202  print "fracture permeability in y direct Line 202  print "fracture permeability in y direct
202  print "fracture permeability in z direction= %f mD"%(PERM_F_Z/(U.mDarcy))  print "fracture permeability in z direction= %f mD"%(PERM_F_Z/(U.mDarcy))
203  print "initial porosity in fractured rock= %f"%EQUIL["DATUM_PRESS"]  print "initial porosity in fractured rock= %f"%EQUIL["DATUM_PRESS"]
204    
205  print PERM_F_X  
 1/0  
206  mkDir(OUTPUT_DIR)  mkDir(OUTPUT_DIR)
207    
208  domain=Rectangle(N_X, N_Y, l0=L_X, l1=L_Y)  domain=Rectangle(N_X, N_Y, l0=L_X, l1=L_Y)
# Line 260  for dt in DT: Line 259  for dt in DT:
259        
260    FN=os.path.join(OUTPUT_DIR, "state.%d.vtu"%(n_t+1))    FN=os.path.join(OUTPUT_DIR, "state.%d.vtu"%(n_t+1))
261    saveVTK(FN,p=p, S_fg=S_fg, C_mg=C_mg)    saveVTK(FN,p=p, S_fg=S_fg, C_mg=C_mg)
262    print "DDD", (t+dt)/U.day, well_P1.getBHP()/U.psi    print "DDD", (t+dt)/U.day, well_P1.getBHP()/U.psi, well_P1.getGasProduction(), well_P1.getWaterProduction()/U.Barrel*U.day
263    print "<%s>State %s saved to file %s."%(time.asctime(),n_t+1,FN )    print "<%s>State %s saved to file %s."%(time.asctime(),n_t+1,FN )
264    
265    n_t+=1    n_t+=1

Legend:
Removed from v.3496  
changed lines
  Added in v.3497

  ViewVC Help
Powered by ViewVC 1.1.26