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

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