OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added R script to extract micromodel subset

Browse Source
This commit is contained in:
James E McClure 2018-02-19 10:44:49 -05:00
parent 1117059e86
commit 579689746a
2 changed files with 5487 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5466
example/Juanes/FullMicromodel.discs Normal file
View File

File diff suppressed because it is too large Load Diff

21
example/Juanes/ReadDiscPack.R Normal file
View File

@ -0,0 +1,21 @@
require("ggplot2")
Discs<-read.csv("FullMicromodel.discs",head=FALSE,sep=" ")
colnames(Discs)<-c("cx","cy","radius")
L=0.45
# Extract some subset from the interior of the discs
SubDiscs<-subset(Discs,Discs$cy>0.9-L)
SubDiscs<-subset(SubDiscs,SubDiscs$cy<0.9+L)
SubDiscs<-subset(SubDiscs,SubDiscs$cx>0.9-L)
SubDiscs<-subset(SubDiscs,SubDiscs$cx<0.9+L)
SubDiscs$cx<-SubDiscs$cx-0.9+L
SubDiscs$cy<-SubDiscs$cy-0.9+L
write.table(SubDiscs,file="DiscPack.in",quote=FALSE,row.names=FALSE,col.names=FALSE,sep=" ")
ShowPlot<-ggplot(SubDiscs)+geom_circle(aes(x0=cx,y0=cy,r=radius))