Commit: 548001fc24e9a87ce5d7497b3576575a37595639 - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revision: 548001fc24e9a87ce5d7497b3576575a37595639 (tree)
• Zeit: 2007-06-13 17:48:31
• Autor: iselllo
• Commiter: iselllo

Log Message

I fixed some bug in lognormal_binning.R. Now I correctly generate a set
of diameters in [D_min,D_max] which are log spaced.

Ändern Zusammenfassung

• modified: R-codes/lognormal_binning.R

Diff

diff -r 001a27c28856 -r 548001fc24e9 R-codes/lognormal_binning.R
--- a/R-codes/lognormal_binning.R	Tue Jun 12 16:48:44 2007 +0000
+++ b/R-codes/lognormal_binning.R	Wed Jun 13 08:48:31 2007 +0000
@@ -2,11 +2,11 @@
 #first I fix a grid in the volume which will then become a diamter grid
 #the maximum diameter can be 1000nm and careful about the units we use!
 #1000nm=1e-6m
-D_max<-600 #in nm
+D_max<-1200 #in nm
 #Now I am going to create a sequence of 100 volumes
-N_vol<-100
+N_vol<-250
 # Now I introduce the value of the fractal dimension D_f
-D_f<-2.9
+D_f<-2.25
 print ("the chosen fractal dimension is")
 print (D_f)
 write.table(D_f,"fractal_dimension",col.names=FALSE,row.names=FALSE)
@@ -16,7 +16,9 @@
 print("the diameter of the basic monomer is [nm]")
 print(D_0)
 print ("and consequently the volume is [nm^3]")
-V_0<-pi/6*D_0^3
+V_0<-pi/6*D_0^3 # NB: the fundamental particle does not have a fractal structure
+# so its solid volume is simply the one of a sphere (or equivalently the effective density
+# is equal to the physical density of the material it is made up of).
 print(V_0)
 # I now save the volume of the fundamental monomer in m
 write.table(D_0/1e9/2,"monomer_radius",col.names=FALSE,row.names=FALSE)
@@ -26,7 +28,8 @@
 print("D_min is [nm]")
 print (D_min)
 
-
+#In the end, the lowest bin must correspond to a D_min which has to be at least as large as
+# the one of the primary particles
 
 linear_bin<-0
 
@@ -68,21 +71,23 @@
 }
 if (linear_bin !=1)
 {# then I want to introduce some log-binning
+# first I work on the diameter grid, then I will translate it into the 
+# grid on the volume of solid
 print ("chosen a log binning")
-D_seq<-seq(len=(N_vol+1))
+D_seq<-seq(0,len=(N_vol+1))
 # klog<-log(D_max)
 
 
 klog<-log(D_max)-log(D_min)
 delta_log<-klog/(length(D_seq)-1)
 
-# D_seq[1]<-1
-D_seq[1]<-D_min
-# delta_log_seq<-exp(seq(1:(length(D_seq)-1))*delta_log)
-delta_log_seq<-exp(seq(1:(length(D_seq)-1))*delta_log)*D_min
-D_seq[2:length(D_seq)]<-delta_log_seq # this way I created a log-spaced diameter binning
+
+# this way I create a log-spaced diameter binning
+D_seq<-exp(D_seq*delta_log)*D_min
 d_temp<-D_seq
 print("OK D_seq_log")
+# now that I have my diameter grid, I can create a corresponding grid of solid volumes
+# which turns out to be a function of the fractal dimension
 V_seq<-pi/6*D_seq^3*(D_seq/D_0)^(D_f-3) #New sequence of volumes
 }
 #now I have a list of diameters corresponding to a constant shape of the vectors
@@ -106,7 +111,7 @@
 N_meas<-1.061e8*1e6 #I multiplied by 1e6 in order to go from cm^-3 
 # to m^-3 
 
-correction<-0.67 # correction factor to account for the conservation of the total volume
+correction<-0.640 # correction factor to account for the conservation of the total volume
 
 bin_pop<-bin_pop*N_meas*correction 
 
@@ -394,16 +399,20 @@
 
 #Now I work out the diameter sequence
 # D_seq[1]<-2*D_seq[2]-D_seq[3] #or maybe I set the first value by hand
-D_seq[1]<-10 #now I have the problem I mentioned above. Since I cannot take the 
+# D_seq[1]<-10 #now I have the problem I mentioned above. Since I cannot take the 
 # log(0), I need to re-set the first value of the diameter sequence
 
 
 
+#probably I do not need this bit any more, now I am setting D_seq[1]=D_min
+
 #carefuk with this bit!!!!!!!!
-if (linear_bin != 1)
-{
-D_seq[1]<-1
-}
+# if (linear_bin != 1)
+# {
+# D_seq[1]<-1
+# }
+
+
 
 write.table(D_seq,"diameter_sequence",col.names=FALSE,row.names=FALSE)