Commit: 30739cf3253656d1aa3ac958117f7f083cc3d660 - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revisión: 30739cf3253656d1aa3ac958117f7f083cc3d660 (tree)
• Tiempo: 2007-08-06 08:53:30
• Autor: iselllo
• Commiter: iselllo

Log Message

I fixed a bug in 3D_R_plotting.R: since V_list is already the volume of
solid in each bin, it has to be multiplied by the primary (NOT the
effective) density to provide the mass in each bin.

Cambiar Resumen

• modified: R-codes/3D_R_plotting.R

Diferencia incremental

diff -r c22ec535b2ff -r 30739cf32536 R-codes/3D_R_plotting.R
--- a/R-codes/3D_R_plotting.R	Sun Aug 05 17:25:06 2007 +0000
+++ b/R-codes/3D_R_plotting.R	Sun Aug 05 23:53:30 2007 +0000
@@ -198,11 +198,14 @@
 V_150<-read.table("V_list", header=FALSE)
 V_150<-as.matrix(V_150)
 
-eff_dens<-read.table("./binning-plots/effective_density_bin",header=FALSE)
-eff_dens<-as.matrix(eff_dens)
+# eff_dens<-read.table("./binning-plots/effective_density_bin",header=FALSE)
+# eff_dens<-as.matrix(eff_dens)
 
-mass_fitting<-eff_dens*V_150*final_fit*log_factor*log(10) #since final fit is dN/dlogd
-mass_simu<-eff_dens*V_150*data_150[number_distr, ]*log_factor*log(10)
+
+rho_0<-1.3
+
+mass_fitting<-rho_0*V_150*final_fit*log_factor*log(10) #since final fit is dN/dlogd
+mass_simu<-rho_0*V_150*data_150[number_distr, ]*log_factor*log(10)
 
 my_error_mass<-sum((mass_fitting-mass_simu)^2)
 
@@ -214,12 +217,13 @@
 
 #initial mass distribution
 
-mass_ini<-mass_simu<-eff_dens*V_150*data_150[1, ]*log_factor*log(10)
+mass_ini<-rho_0*V_150*data_150[1, ]*log_factor*log(10)
 
 pdf("2D-mass_fitting.pdf")
 plot(D_mean,mass_ini,type="l",lty=2,lwd=2,col="black", log="x",
 ylab = list(expression(paste(mass," ["*cm^-3*"]"))),
-xlab = expression(paste(D[agg]," [nm]"))
+xlab = expression(paste(D[agg]," [nm]")),
+ylim=range(0,max(mass_fitting,mass_ini,mass_simu))
 )
 lines(D_mean,mass_fitting,type="b",lty=1,pch=4,lwd=2,col="blue")
  lines(D_mean,mass_simu,lwd=2,col="brown")