μ’μ, λλ μ΄κ²μ λν΄ Gitterμμ λ§μ μλ€λ₯Ό λ¨μμ§λ§ λ¬Έμ λ₯Ό μ΄μ΄μΌ ν νμκ° μλ κ² κ°μ΅λλ€. λλ μ¬κΈ°μ λ¬Έμ λ₯Ό λͺ νν νλ €κ³ λ Έλ ₯ν κ²μ΄μ§λ§ νμμ λ°λΌ λͺ νμ±μ μμ²νμμμ€.
λ¬Έμ λ λ€μκ³Ό κ°μ΅λλ€. κ°μκΈ° μ΅μ λ²μ μ κ°λ° λΆκΈ° W/μ΅μ dvmdostem μ½λλ₯Ό μ€νν λ μν¬νλ‘ λ‘κ·Έμ μΆκ°ν logger.debug() λ¬Έμ λ°λΌ μ¬λ°λ₯Έ μ¬ν( SLAμ κ²½μ°) λ° μ΄μ (μ¬κΈ°μμ λ³Ό μ μλ dvmdostem CMT04-Salix PFTμ λ€λ₯Έ λͺ¨λ νΉμ±: https://www.betydb.org/pfts/136). μλ₯Ό λ€μ΄:
2017-10-10 07:56:52 INFO [PEcAn.logger::logger.info] :
Selected PFT(s): CMT04-Salix
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'prior.distns' not found
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'post.distns' not found
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'trait.mcmc' not found
2017-10-10 07:56:52 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Salix has MCMC samples for: SLA
2017-10-10 07:56:52 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Salix will use prior distributions for: mort2
growth_resp_factor leaf_turnover_rate leaf_width nonlocal_dispersal
fineroot2leaf root_turnover_rate seedling_mortality stomatal_slope
r_fract root_respiration_rate Vm_low_temp Vcmax quantum_efficiency
2017-10-10 07:56:52 INFO [PEcAn.logger::logger.info] :
using 5004 samples per trait
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-10 07:56:52 INFO [PEcAn.logger::logger.info] :
No sampling method supplied, defaulting to uniform random sampling
μ΄μ dist λΆλΆμ λμ΄λ μλͺ»λ νΉμ±μ 무μμ λκΉ?
λ°λΌμ dbfilesμ μ¬ν IDμ PEcAn μ€νμ PFT ν΄λλ₯Ό 보면 λ λ€ μ¬λ°λ₯Έ μ΄μ /νλ°© νμΌ/μ 보μ μΌμΉλ₯Ό 보μ¬μ€λλ€.
μλ₯Ό λ€μ΄:
> load('prior.distns.Rdata')
> ls()
[1] "post.distns" "prior.distns" "q"
> prior.distns
distn parama paramb n
SW_albedo unif 0.010 0.300 NA
gcmax unif 0.001 0.005 NA
cuticular_cond weibull 1.120 54531.000 29
SLA norm 10.000 5.000 NA
frprod_perc_20 unif 0.100 99.000 NA
frprod_perc_10 unif 0.100 99.000 NA
extinction_coefficient_diffuse weibull 4.000 0.400 NA
>
그리κ³
> load('prior.distns.Rdata')
> ls()
[1] "prior.distns" "q"
> prior.distns
distn parama paramb n
SW_albedo unif 0.010 0.300 NA
gcmax unif 0.001 0.005 NA
cuticular_cond weibull 1.120 54531.000 29
SLA norm 10.000 5.000 NA
frprod_perc_20 unif 0.100 99.000 NA
frprod_perc_10 unif 0.100 99.000 NA
extinction_coefficient_diffuse weibull 4.000 0.400 NA
> getwd()
[1] "/data/sserbin/Modeling/dvmdostem/pecan_runs/run.56/pft/CMT04-Salix"
> load('post.distns.MA.Rdata')
> ls90
Error: object 'ls90' not found
> ls()
[1] "post.distns" "prior.distns" "q"
> post.distns
distn parama paramb n
SW_albedo unif 0.01000 3.000000e-01 NA
gcmax unif 0.00100 5.000000e-03 NA
cuticular_cond weibull 1.12000 5.453100e+04 29
SLA norm 10.91521 1.807541e-01 NA
frprod_perc_20 unif 0.10000 9.900000e+01 NA
frprod_perc_10 unif 0.10000 9.900000e+01 NA
extinction_coefficient_diffuse weibull 4.00000 4.000000e-01 NA
> load('post.distns.Rdata')
> ls()
[1] "post.distns" "prior.distns" "q"
> post.distns
distn parama paramb n
SW_albedo unif 0.01000 3.000000e-01 NA
gcmax unif 0.00100 5.000000e-03 NA
cuticular_cond weibull 1.12000 5.453100e+04 29
SLA norm 10.91521 1.807541e-01 NA
frprod_perc_20 unif 0.10000 9.900000e+01 NA
frprod_perc_10 unif 0.10000 9.900000e+01 NA
extinction_coefficient_diffuse weibull 4.00000 4.000000e-01 NA
> prior.distns
distn parama paramb n
SW_albedo unif 0.010 0.300 NA
gcmax unif 0.001 0.005 NA
cuticular_cond weibull 1.120 54531.000 29
SLA norm 10.000 5.000 NA
frprod_perc_20 unif 0.100 99.000 NA
frprod_perc_10 unif 0.100 99.000 NA
extinction_coefficient_diffuse weibull 4.000 0.400 NA
κ·Έλ° λ€μ μ°κΈ° κ΅¬μ± μ½λλ MCMC μνμ΄ μκΈ° λλ¬Έμ μ μΌν λͺ¨λΈ νΉμ±μ΄ SLAμΈ μ€νμ μ΄λνλ μλͺ»λ 맀κ°λ³μλ₯Ό μμ±νλ €κ³ μλν©λλ€. λλ ν¬ν
μ¬μ μ΄ μ΄λμμ μ€λμ§ μ μ μμ΅λλ€. UAFμ Tobey Carmanμ λμΌν μμ€μ κ°λ°/dvmdostem 컀λ°μ μλ VMμμ λμΌν λ¬Έμ λ₯Ό κ²ͺκ³ μμ΅λλ€.
https://github.com/PecanProject/pecan/issues/new
μ§κΈμ λ§νκΈ°κ° λ§€μ° μ΄λ ΅μ§λ§ μ΄ νΉμ PFTμλ§ κ΄λ ¨μ΄ μλ κ² κ°μ΅λλ€. κ·Έλ¬λ νμ€νμ§ μμ΅λλ€. PFT μ½€λ³΄λ‘ λ€λ₯Έ λͺ¨λλ₯Ό μ΄μ§λ½νμ μμ§ λ³Έ μ μ΄ μμ΅λλ€. λν pEcAn/betyμ νλͺ© λͺ¨λΈ λ±λ‘μλ§ ν΄λΉλ μ μμ΅λλ€.
λ€μμ λ΄ pecan.CONFIGS.xml μ€ νλμ μ€λν«μ λλ€.
?xml version="1.0"?>
<pecan>
<outdir>/data/sserbin/Modeling/dvmdostem/pecan_runs/run.56/</outdir>
<database>
<bety>
<user>bety</user>
<password>bety</password>
<host>localhost</host>
<dbname>bety</dbname>
<driver>PostgreSQL</driver>
<write>FALSE</write>
</bety>
<dbfiles>/data/home/sserbin/.pecan/dbfiles</dbfiles>
</database>
<pfts>
<pft>
<name>CMT04-Salix</name>
<constants>
<num>1</num>
</constants>
<posteriorid>2000001342</posteriorid>
<outdir>/data/sserbin/Modeling/dvmdostem/pecan_runs/run.56//pft/CMT04-Salix</outdir>
</pft>
</pfts>
κ·Έλμ get.parameter.samplesμ λͺ κ°μ§ μΈμ λ¬Έμ μΆκ°νμ΅λλ€.
### Identify PFTs in the input settings.xml file
get.parameter.samples <- function(settings,
posterior.files = rep(NA, length(settings$pfts)),
ens.sample.method = "uniform") {
## added by sps
PEcAn.logger::logger.debug(paste0("get.param.samples.R posterior file: ", posterior.files))
PEcAn.logger::logger.debug(paste0("priors: ", prior.distns))
# added by sps
pfts <- settings$pfts
num.pfts <- length(settings$pfts)
pft.names <- list()
outdirs <- list()
for (i.pft in seq_along(pfts)) {
pft.names[i.pft] <- settings$pfts[[i.pft]]$name
### If no PFT(s) are specified insert NULL to warn user
if (length(pft.names) == 0) {
pft.names[1] <- "NULL"
}
### Get output directory info
outdirs[i.pft] <- settings$pfts[[i.pft]]$outdir
PEcAn.logger::logger.info(paste0("PFT Output Directories: ", outdirs)) # added by sps
} ### End of for loop to extract pft names
PEcAn.logger::logger.info("Selected PFT(s): ", pft.names)
## Generate empty list arrays for output.
trait.samples <- sa.samples <- ensemble.samples <- env.samples <- runs.samples <- list()
## Load PFT priors and posteriors
for (i in seq_along(pft.names)) {
#rm(prior.distns, post.distns, trait.mcmc) # sps
## Load posteriors
if (!is.na(posterior.files[i])) {
# Load specified file
load(posterior.files[i])
if (!exists("prior.distns") & exists("post.distns")) {
prior.distns <- post.distns
}
} else {
# Default to most recent posterior in the workflow, or the prior if there is none
fname <- file.path(outdirs[i], "post.distns.Rdata")
PEcAn.logger::logger.debug(paste0("posterior fname: ",fname)) # added by sps
if (file.exists(fname)) {
load(fname)
prior.distns <- post.distns
} else {
load(file.path(outdirs[i], "prior.distns.Rdata"))
}
}
## added by sps
PEcAn.logger::logger.debug(paste0("get.param.samples.R posterior file: ", posterior.files))
PEcAn.logger::logger.debug(paste0("priors: ", prior.distns))
# added by sps
### Load trait mcmc data (if exists)
if ("trait.mcmc.Rdata" %in% dir(unlist(outdirs[i]))) {
ma.results <- TRUE
PEcAn.logger::logger.debug("*** Trait MCMC samples found ***") # added by sps
load(file.path(outdirs[i], "trait.mcmc.Rdata"))
} else {
ma.results <- FALSE
PEcAn.logger::logger.debug("*** Trait MCMC samples not found ***") # added by sps
}
pft.name <- unlist(pft.names[i])
λ‘κ·Έ μΆλ ₯
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-10 16:54:40 WARN [PEcAn.logger::logger.warn] :
Multiple Valid Files found on host machine. Returning last updated
record.
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-10 16:54:40 DEBUG [PEcAn.logger::logger.debug] :
get.param.samples.R posterior file:
/data/home/sserbin/.pecan/dbfiles/posterior/2000001348/post.distns.MA.Rdata
2017-10-10 16:54:40 DEBUG [PEcAn.logger::logger.debug] :
priors: c("gamma", "beta", "gamma", "gamma", "beta", "lnorm", "weibull",
"beta", "weibull", "beta", "weibull", "norm", "weibull", "weibull",
"gamma") priors: c(1.47, 2.63, 2.9, 6.53, 20.3, 0.811, 1.67, 3.61, 3.63,
2, 2.66, 10, 2.06, 3.55, 90.9) priors: c(0.0578, 6.52, 0.63, 1.49, 76.1,
0.843, 0.657, 0.433, 3.81, 4, 6.29, 1.02, 19, 24.7, 1580) priors: c(0,
0, 40, 17, 30, 0, 66, 0, 4, 0, 35, 0, 125, 97, 56)
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-10 16:54:40 INFO [PEcAn.logger::logger.info] :
PFT Output Directories:
/data/sserbin/Modeling/dvmdostem/pecan_runs/run.60//pft/CMT04-Salix
2017-10-10 16:54:40 INFO [PEcAn.logger::logger.info] :
Selected PFT(s): CMT04-Salix
2017-10-10 16:54:40 DEBUG [PEcAn.logger::logger.debug] :
get.param.samples.R posterior file:
/data/home/sserbin/.pecan/dbfiles/posterior/2000001348/post.distns.MA.Rdata
2017-10-10 16:54:40 DEBUG [PEcAn.logger::logger.debug] :
priors: c("gamma", "beta", "gamma", "gamma", "beta", "lnorm", "weibull",
"beta", "weibull", "beta", "weibull", "norm", "weibull", "weibull",
"gamma") priors: c(1.47, 2.63, 2.9, 6.53, 20.3, 0.811, 1.67, 3.61, 3.63,
2, 2.66, 10, 2.06, 3.55, 90.9) priors: c(0.0578, 6.52, 0.63, 1.49, 76.1,
0.843, 0.657, 0.433, 3.81, 4, 6.29, 1.02, 19, 24.7, 1580) priors: c(0,
0, 40, 17, 30, 0, 66, 0, 4, 0, 35, 0, 125, 97, 56)
2017-10-10 16:54:40 DEBUG [PEcAn.logger::logger.debug] :
*** Trait MCMC samples found ***
2017-10-10 16:54:40 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Salix has MCMC samples for: SLA
2017-10-10 16:54:40 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Salix will use prior distributions for: mort2
growth_resp_factor leaf_turnover_rate leaf_width nonlocal_dispersal
fineroot2leaf root_turnover_rate seedling_mortality stomatal_slope
r_fract root_respiration_rate Vm_low_temp Vcmax quantum_efficiency
2017-10-10 16:54:40 INFO [PEcAn.logger::logger.info] :
using 5004 samples per trait
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
μ΄μν μ μ 첫 λ²μ§Έ λ‘κ±° λ¬Έ μ§ν©μ΄ for 루νμμ μ¬μ /νμ ꡬ문 λΆμμ λλ¬νκΈ° μ μ μ€νμ μ΄λ―Έ posterior.files λ° prior.distnsκ° μ€μ λμ΄ μμμ 보μ¬μ€λλ€. κ·Έλμ λλ μ¬μ ν μ°μ μμκ° μ΄λμμ μμΌλ©° κ·Έ κΈ°λ₯μμ μ΄λ»κ² μ¬μ©λλμ§λ₯Ό 격리ν μ μμ΅λκΉ? μ½λλ₯Ό κΈ°λ°μΌλ‘ PFT ν΄λμ μ 보λ₯Ό μ¬μ©ν΄μΌ νλ κ²μ²λΌ 보μ΄μ§λ§ λ€λ₯Έ κ³³μμ λ€μ΄μ€λ κ²μ²λΌ 보μ λκΉ?
νμ€νλνμ§λ§μ΄ 'μλμ€'μ΄λ€ μ΄μ μ λν μ νν μ¬ν μ 보λ₯Ό μ»μμ΄ λΉνΈλ₯Ό runnignν΄μΌνλ€λ κ²
} else {
# Default to most recent posterior in the workflow, or the prior if there is none
fname <- file.path(outdirs[i], "post.distns.Rdata")
PEcAn.logger::logger.debug(paste0("posterior fname: ",fname)) # added by sps
if (file.exists(fname)) {
load(fname)
prior.distns <- post.distns
} else {
load(file.path(outdirs[i], "prior.distns.Rdata"))
}
}
λ‘κ±° μΆλ ₯μ΄ μκΈ° λλ¬Έμ μ€νλμ§ μλλ€λ κ²μ μκ³ μμ΅λλ€ ....
λ°λΌμ ν΄λΉ λΆλΆμ μ€ννλ λμ λ€λ₯Έ κ³³μμ λ€λ₯Έ prior.distnsλ₯Ό μ¬μ©νμ¬ λκ° μ΄λμ μ μ μμ΅λκΉ???
μ’μ, λ§€μ° μ΄μνκ³ λλ μ΄κ²μΌλ‘ κΈΈμ μμλ€. μ DVM-DOS-TEM PFTλ₯Ό μμ±νλλ° μ¬μ ν λμΌν λ¬Έμ κ° λ°μν©λλ€.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:28:17 INFO [PEcAn.logger::logger.info] :
Selected PFT(s): CMT04-Betula
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'prior.distns' not found
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'post.distns' not found
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'trait.mcmc' not found
2017-10-13 08:28:17 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Betula has MCMC samples for: SLA
2017-10-13 08:28:17 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Betula will use prior distributions for: mort2
growth_resp_factor leaf_turnover_rate leaf_width nonlocal_dispersal
fineroot2leaf root_turnover_rate seedling_mortality stomatal_slope
r_fract root_respiration_rate Vm_low_temp Vcmax quantum_efficiency
2017-10-13 08:28:17 INFO [PEcAn.logger::logger.info] :
using 5004 samples per trait
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
μ΄ μ€λ₯κ° dvm-dos-temκ³Ό κ΄λ ¨νμ¬ νΉμ λͺ¨λΈμΌ μ μλ μ΄μ κ° μμ΅λκΉ? μλ§λ 무μΈκ°κ° NAMESPACEμμ λΆμ μ νκ² λ§μ€νΉλκ³ μμ΅λκΉ? μ΄λ»κ² μ΄λ° μΌμ΄ μΌμ΄λλμ§ λ³Ό μ μμ΅λκΉ?
μκ² μ΅λλ€. λ¬Έμ κ° λ°μνμ΅λλ€. λ€μμ SIPNET μ€νμ λλ€.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:32:10 INFO [PEcAn.logger::logger.info] :
Selected PFT(s): tundra.deciduous.NGEE_Arctic
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'prior.distns' not found
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'post.distns' not found
Warning in rm(prior.distns, post.distns, trait.mcmc) :
object 'trait.mcmc' not found
2017-10-13 08:32:10 INFO [PEcAn.logger::logger.info] :
PFT tundra.deciduous.NGEE_Arctic has MCMC samples for: leafC Amax
Vm_low_temp leaf_respiration_rate_m2 SLA
2017-10-13 08:32:10 INFO [PEcAn.logger::logger.info] :
PFT tundra.deciduous.NGEE_Arctic will use prior distributions for: mort2
growth_resp_factor leaf_turnover_rate leaf_width nonlocal_dispersal
fineroot2leaf root_turnover_rate seedling_mortality stomatal_slope
r_fract root_respiration_rate Vcmax quantum_efficiency
2017-10-13 08:32:10 INFO [PEcAn.logger::logger.info] :
using 3004 samples per trait
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:32:10 INFO [PEcAn.logger::logger.info] :
No sampling method supplied, defaulting to uniform random sampling
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:32:10 INFO [PEcAn.logger::logger.info] :
Using uniform random sampling
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
λΆλͺ ν Vcmaxμ mort2λ μ€λ₯μ λλ€. μ΄ μ¬μ μ μ΄λμμ μ€λκ°!
λν λ€μμ meta.analysis μ μ λ΄ μν¬νλ‘.logμ μ΄μ μΆλ ₯μ λλ€.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:28:10 INFO [PEcAn.logger::logger.info] :
---------------------------------------------------------
2017-10-13 08:28:10 INFO [PEcAn.logger::logger.info] : SLA
Auto-disconnecting PostgreSQLConnection
2017-10-13 08:28:13 INFO [PEcAn.logger::logger.info] : Median SLA : 12.525
2017-10-13 08:28:13 INFO [PEcAn.logger::logger.info] :
---------------------------------------------------------
2017-10-13 08:28:13 INFO [PEcAn.logger::logger.info] :
Summary of Prior distributions for: CMT04-Betula
2017-10-13 08:28:13 INFO [PEcAn.logger::logger.info] :
distn parama paramb n
2017-10-13 08:28:13 INFO [FUN] :
cuticular_cond lnorm 8.400 0.900 0
2017-10-13 08:28:13 INFO [FUN] :
extinction_coefficient_diffuse gamma 5.000 10.000 NA
2017-10-13 08:28:13 INFO [FUN] :
gcmax unif 0.001 0.005 NA
2017-10-13 08:28:13 INFO [FUN] :
SLA norm 10.000 5.000 NA
2017-10-13 08:28:13 INFO [FUN] :
frprod_perc_20 unif 0.100 99.000 NA
2017-10-13 08:28:13 INFO [FUN] :
frprod_perc_10 unif 0.100 99.000 NA
2017-10-13 08:28:13 INFO [PEcAn.logger::logger.info] :
number of observations per trait for CMT04-Betula
2017-10-13 08:28:13 INFO [PEcAn.logger::logger.info] :
34 observations of SLA
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
μ¬λ°λ₯Έ νΉμ±μ΄ νμλλ κ²μ λ³Ό μ μμΌλ―λ‘ λ¬Έμ λ λ©ν λΆμ ν write.configsλ₯Ό μ€ννκΈ° μ μ λ°μν©λλ€.
μ’μ, κ·Έλ¦¬κ³ λλλ‘ μλμ κ°μ΄ μ μλν κ²μ΄λ€.
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:48:05 INFO [PEcAn.logger::logger.info] :
PFT Output Directories:
/data/sserbin/Modeling/dvmdostem/pecan_runs/run.61//pft/CMT04-Salix
2017-10-13 08:48:05 INFO [PEcAn.logger::logger.info] :
Selected PFT(s): CMT04-Salix
2017-10-13 08:48:05 DEBUG [PEcAn.logger::logger.debug] :
posterior fname:
/data/sserbin/Modeling/dvmdostem/pecan_runs/run.61//pft/CMT04-Salix/post.distns.Rdata
2017-10-13 08:48:05 DEBUG [PEcAn.logger::logger.debug] :
get.param.samples.R posterior file: NA
2017-10-13 08:48:05 DEBUG [PEcAn.logger::logger.debug] :
priors: c("weibull", "unif", "unif", "norm", "unif", "unif", "weibull")
priors: c(1.12, 0.01, 0.001, 10.9152064936782, 0.1, 0.1, 4) priors:
c(54531, 0.3, 0.005, 0.180754050808388, 99, 99, 0.4) priors: c(29, NA,
NA, NA, NA, NA, NA)
2017-10-13 08:48:05 DEBUG [PEcAn.logger::logger.debug] :
*** Trait MCMC samples found ***
2017-10-13 08:48:05 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Salix has MCMC samples for: SLA
2017-10-13 08:48:05 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Salix will use prior distributions for: cuticular_cond
SW_albedo gcmax frprod_perc_20 frprod_perc_10
extinction_coefficient_diffuse
2017-10-13 08:48:05 INFO [PEcAn.logger::logger.info] :
using 5004 samples per trait
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
μ¬κΈ° μ§μ§ λκ° μ§μκ° μλ κ² κ°μμ....
get.param.samples.R μ¬ν νμΌ μμ λ‘κ·Έ μΆλ ₯μμ ββμ°Έκ³ : NA λ° κ·Έκ²μ μ μλνμ§λ§ λ λμκ°μ κ·Έ print λ¬Έμ νμΌμ λν κ²½λ‘, μ¬λ°λ₯Έ μ λ³΄κ° μλ νμΌμ 보μ¬μ£Όμ§λ§ μλͺ»λ κ²μ μ¬μ©ν©λλ€. μ¬μ ?
μ, pecan.xmlμμ λμΌν μν¬νλ‘μ μ€μ μ μ€ννλ©΄, μ¦ μν¬νλ‘λ₯Ό λ€μ μμ±νλ©΄ λ€λ₯Έ κ²°κ³Όλ₯Ό μ»μ΅λλ€!!!
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:56:35 WARN [PEcAn.logger::logger.warn] :
Multiple Valid Files found on host machine. Returning last updated
record.
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:56:36 DEBUG [PEcAn.logger::logger.debug] :
get.param.samples.R posterior file:
/data/home/sserbin/.pecan/dbfiles/posterior/2000001353/post.distns.MA.Rdata
2017-10-13 08:56:36 DEBUG [PEcAn.logger::logger.debug] :
priors: c("gamma", "beta", "gamma", "gamma", "beta", "lnorm", "weibull",
"beta", "weibull", "beta", "weibull", "norm", "weibull", "weibull",
"gamma") priors: c(1.47, 2.63, 2.9, 6.53, 20.3, 0.811, 1.67, 3.61, 3.63,
2, 2.66, 10, 2.06, 3.55, 90.9) priors: c(0.0578, 6.52, 0.63, 1.49, 76.1,
0.843, 0.657, 0.433, 3.81, 4, 6.29, 1.02, 19, 24.7, 1580) priors: c(0,
0, 40, 17, 30, 0, 66, 0, 4, 0, 35, 0, 125, 97, 56)
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 08:56:36 INFO [PEcAn.logger::logger.info] :
PFT Output Directories:
/data/sserbin/Modeling/dvmdostem/pecan_runs/run.61//pft/CMT04-Salix
2017-10-13 08:56:36 INFO [PEcAn.logger::logger.info] :
Selected PFT(s): CMT04-Salix
2017-10-13 08:56:36 DEBUG [PEcAn.logger::logger.debug] :
μ΄λ²μλ print λ¬Έμ΄ NAκ° μλλ©° write.configsμμ μλͺ»λ νΉμ±μ μ»μμ μ μ μμ΅λλ€!!
λ°λΌμ μ¬λ°λ₯΄κ² μλνμ§ μμΌλ©΄ μ¬λ°λ₯Έ μ¬ν νμΌμ κ°λ¦¬ν€μ§λ§ κ·Έ νμΌμ΄λ PFTμ μ‘΄μ¬νμ§ μλ νΉμ±μ λν μ¬μ μ 보λ₯Ό μ΄λ»κ²λ μ»μ΅λλ€!?
κ·Έλμ μ΄κ²μ΄ λμμ΄ λλμ§ λͺ¨λ₯΄κ² μ§λ§ μ΄ λ¬Έμ λ₯Ό κ²ͺκ³ μλ @tobeycarman μ΄ μ€μλ‘ μλͺ»λ μμ λλ ν 리μμ μν¬νλ‘λ₯Ό μμνμ¬(μμ€ μΆμ²) PEcAn.xmlμ μ°Ύμ§ λͺ»νμ΅λλ€. κ·Έλ¬λ μ€λ₯ λ©μμ§λ ν₯λ―Έ λ‘μ΅λλ€.
>
>
> # ----------------------------------------------------------------------
> # PEcAn Workflow .... [TRUNCATED]
> if (is.na(args[1])){
+ settings <- PEcAn.settings::read.settings("pecan.xml")
+ } else {
+ settings.file = args[1]
+ settings <- PEcAn.setting .... [TRUNCATED]
2017-10-13 08:05:51 SEVERE [PEcAn.logger::logger.severe] :
Could not find a pecan.xml file
Warning in file(file, ifelse(append, "a", "w")) :
cannot open file '/home/dlebauer/out/ebifarm/post//STATUS': No such file or directory
Error in file(file, ifelse(append, "a", "w")) :
cannot open the connection
Error during wrapup: cannot open the connection
>
μ΄λ»κ²λ λ°μ΄ν°λ² μ΄μ€μ @dlebauer μ μ΄μ μ€νμ μ¬νκ° μκ³ ED2 μ€νμ²λΌ 보μ΄λ©° μ΄λ»κ²λ κΈ°λ³Έκ°μ΄ μ€μ λλ κ²μ΄ κ°λ₯ν©λκΉ? λλ R ν¨ν€μ§μ μ΄λκ°μ μμ μ μμ΅λλ€. λλ μ΄μ΄ μμ΄ μ΄λ¦¬μ 리 κ²μνλ€. λλ μ΄κ²μ΄ ν°λ¬΄λμλ λ‘± μ·μ΄λΌλ κ²μ μκ³ μμ§λ§ μμ΄λμ΄λ₯Ό μ΄λ°μν€λ κ²μ΄λ©΄ λ¬Όμ΄λ³Ό κ²μ΄λΌκ³ μκ°νμ΅λλ€.
VMκ³Ό κ΄λ ¨μ΄ μλ ν κ·Έμ μ€λ₯κ° μ΄λ»κ² κ·Έλ κ² λνλ μ§ μ μ μκΈ° λλ¬Έμ λλ€.
λ€μμ @tobeycarmanμ μΆλ ₯λ
λλ μΆκ°νλ€
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 09:57:01 WARN [PEcAn.logger::logger.warn] :
Multiple Valid Files found on host machine. Returning last updated
record.
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 09:57:01 DEBUG [PEcAn.logger::logger.debug] :
get.param.samples.R posterior file: NA
2017-10-13 09:57:01 DEBUG [PEcAn.logger::logger.debug] :
priors: c("gamma", "beta", "gamma", "gamma", "beta", "lnorm", "weibull",
"beta", "weibull", "beta", "weibull", "norm", "weibull", "weibull",
"gamma") priors: c(1.47, 2.63, 2.9, 6.53, 20.3, 0.811, 1.67, 3.61, 3.63,
2, 2.66, 10, 2.06, 3.55, 90.9) priors: c(0.0578, 6.52, 0.63, 1.49, 76.1,
0.843, 0.657, 0.433, 3.81, 4, 6.29, 1.02, 19, 24.7, 1580) priors: c(0,
0, 40, 17, 30, 0, 66, 0, 4, 0, 35, 0, 125, 97, 56)
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
See help("Deprecated")
2017-10-13 09:57:01 INFO [PEcAn.logger::logger.info] :
PFT Output Directories:
/data/sserbin/Modeling/dvmdostem/pecan_runs/run.64//pft/CMT04-Salix
2017-10-13 09:57:01 INFO [PEcAn.logger::logger.info] :
Selected PFT(s): CMT04-Salix
2017-10-13 09:57:01 DEBUG [PEcAn.logger::logger.debug] :
posterior fname:
/data/sserbin/Modeling/dvmdostem/pecan_runs/run.64//pft/CMT04-Salix/post.distns.Rdata
2017-10-13 09:57:01 DEBUG [PEcAn.logger::logger.debug] :
get.param.samples.R posterior file: NA
2017-10-13 09:57:01 DEBUG [PEcAn.logger::logger.debug] :
priors: c("weibull", "unif", "unif", "norm", "unif", "unif", "weibull")
priors: c(1.12, 0.01, 0.001, 10.9152064936782, 0.1, 0.1, 4) priors:
c(54531, 0.3, 0.005, 0.180754050808388, 99, 99, 0.4) priors: c(29, NA,
NA, NA, NA, NA, NA)
2017-10-13 09:57:01 DEBUG [PEcAn.logger::logger.debug] :
*** Trait MCMC samples found ***
2017-10-13 09:57:01 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Salix has MCMC samples for: SLA
2017-10-13 09:57:01 INFO [PEcAn.logger::logger.info] :
PFT CMT04-Salix will use prior distributions for: cuticular_cond
SW_albedo gcmax frprod_perc_20 frprod_perc_10
extinction_coefficient_diffuse
2017-10-13 09:57:01 INFO [PEcAn.logger::logger.info] :
using 5004 samples per trait
Warning: '$.SafeList' is deprecated.
Use 'PEcAn.settings::$.SafeList' instead.
λ°λΌμ 무μΈκ°λ‘ μΈν΄ νΌμΉΈμ΄ λ³κ²½ μ¬νμ ν½μ νλ©΄ μ λλ‘ μλνμ§λ§ κ·Έλ μ§ μμΌλ©΄ κ·Έλ μ§ μμ κ² κ°μ΅λλ€. κ·Έλ¦¬κ³ λ€μ get.param.samples.R μ¬ν νμΌ: μλν λ
μΉμμ κ·Έλ κ² νλ©΄ λμΌν κ²°κ³Όκ° λμ€μ§ μμΌλ―λ‘ μλΉν μΌκ΄μ±μ΄ μμ΅λλ€.
μΆκ°
@serbinsh μ¬κΈ° κΉμ§ μν¬νλ‘λ₯Ό ν μ€μ© μ€νν μ μμ΅λκΉ?
"νΈμ€νΈ μμ€ν
μμ μ¬λ¬ μ ν¨ν νμΌμ μ°Ύμμ΅λλ€. λ§μ§λ§μΌλ‘ μ
λ°μ΄νΈλ νμΌμ λ°νν©λλ€.
κΈ°λ‘. " κ²½κ³ κ° λ¨μκ° λ μ μμ΅λλ€. dbfile.check λ λλλ‘ μμΉ μλ νμΌ λͺ©λ‘μ λ°νν μ μμ΅λλ€.
κ·Έκ²μ΄ κ·νμ λ¬Έμ μ κ΄λ ¨μ΄ μλμ§ νμ€νμ§ μμ§λ§ μ΄κ²μ΄ PDAμμ λ¬Έμ λ₯Ό μΌμΌν€κ³ μμκΈ° λλ¬Έμ ν¨μμ return.all
μ΅μ
μ μΆκ°νμ΅λλ€.
DALECμ λμΌν λ¬Έμ κ° μμ΅λλ€.
κ°μ₯ μ΅κ·Όμ μ€νν 1000008435
λ‘ κ΅¬μ± νμΌμ μμ±νμ΅λλ€.
-mort2 19.5518192434762
-growth_resp_factor 0.269134877386544
-t5 0.00115746296987216
-leaf_width 4.21460770704145
-nonlocal_dispersal 0.208363236071724
-fineroot2leaf 2.28414274764454 -t7 0.00143050364071698
-seedling_mortality 0.95094009313231
-stomatal_slope 3.44801807715082
-r_fract 0.31488851365167
-root_respiration_rate 5.48891170784511
-Vm_low_temp 9.99793554087408
-SLA 0.0329788044168264
-Vcmax 22.4703903038444
-quantum_efficiency 0.0574557888428778
νμ€ν μλͺ»λ λ³μμ λλ€.
λ¦΄λ¦¬μ€ VMμμ Betsyμ μ€νμ 볡μ νλ €κ³ μλνμΌλ©° DALEC κ΅¬μ± νμΌμ
μ€νμ΄ μ±κ³΅νμΌλ―λ‘ μ΄λ»κ²λ κ°λ°μ μ΅μ μ½λκ° νΉμ νΉμ±μ κΈ°λ³Έκ°μΌλ‘ μ€μ νκ³ μμμ΄ λΆλͺ ν©λλ€.
@tonygardella λλ κ·Έκ²μ΄ κ°λ° λΆκΈ° λ¬Έμ λΌλ λ° λμν©λλ€. λ΄ μ§κ°μ κΈ°λ₯κ³Ό NAMESPACEλ₯Ό μ΄λ¦¬μ 리 μ΄λν μλνμ§ μμ κ²°κ³Όμ΄μ§λ§ μμ§ μ΄μ λ₯Ό μ μ μμ΅λλ€.
μ΄κ²μ write.configs λ¨κ³κΉμ§ μΆμ νμ΅λλ€. νΉν get.parameter.samples λ¨κ³μμ μ¬λ°λ₯Έ μ΄μ λ°/λλ μ¬ν λ°μ΄ν°λ₯Ό μ°Ύμ λ€μ ν΄λΉ νΉμ νΉμ±μ κΈ°λ³Έκ°μΌλ‘ μ€μ νλ κ²μ μλ§μΌλ‘ λ§λλλ€.
@tonygardella κ·Έ μ΄μ μ μ‘λ κ²½λ‘/μμΉλ₯Ό λ³Ό μ μμ΅λκΉ?
λλ μ΄κ²μ΄ λ¬Έμ κ° λλ νμΌμ΄λΌκ³ μκ°ν©λλ€: https://github.com/PecanProject/pecan/blob/develop/base/utils/data/post.distns.RData
get.parameter.samplesμμ μ΄ λΌμΈμ μ΄μ λλ¬ ν λ κ·Έλ° prior.distnsμ λμμ΄ μ κ±° μ΄ λΌμΈ μ prior.distnsμ posterior.distnsμ ν λΉ ν μ μμ§λ§,μ΄ κ²½μ°λ‘λ prior.distns.RData νμΌμ μ¬μ νλ‘λ PEcAn.utils ν¨ν€μ§κ° λ‘λλκΈ° λλ¬Έμ μ¬μ .distnsλ‘. ν΄λΉ κ°μ²΄κ° μμΌλ©΄ ν΄λΉ ν¨ν€μ§ νμΌμ κΈ°λ³Έκ°μΌλ‘ κ³μ μ¬μ©ν©λλ€.
μ΄κ²μ΄ VMμ΄ μλ λ¨Έμ μμ λ°μνλ μ΄μ λ VMμ μλ μ¬ν νμΌμ΄ μ‘΄μ¬νμ§ μλ ν μ½λμ ν΄λΉ λΆλΆμ΄ νΈλ¦¬κ±°λμ§ μκΈ° λλ¬Έμ λλ€.
ν μ€μ© μ€νν λ λ΄ νκ²½μμ prior.distnsλ₯Ό κ³μ μμ νλ €κ³ νμ§λ§ κ°μ²΄λ₯Ό ββμ°Ύμ μ μμμ§λ§ μ½λλ₯Ό μ€νν λ μ΄λ»κ²λ prior.distns κ°μ²΄κ° μμ΄μ ν΄λΉ νμΌλ‘ μ΄λνμ΅λλ€. utils ν¨ν€μ§μ μΌλΆμ λλ€.
νμΌμ΄ μ΄λκ°μμ ν μ€νΈμ μ¬μ©λ©λκΉ? https://github.com/PecanProject/pecan/blob/develop/base/utils/data/post.distns.RData
λΉ λ₯Έ κ²μμ μ΄λμλ μ¬μ©λ κ²μΌλ‘ νμλμ§ μμ΅λλ€.
Aaaand μ΄κ²μ μ κ·Έκ²μ΄ λνλ¬λμ§ μ€λͺ
ν©λλ€ -- ν¨ν€μ§ λ΄ post.distns
λ #1684μμ utilsλ₯Ό μ§μ° λ‘λ©μΌλ‘ μ νν λκΉμ§ μν¬νλ‘μμ μλμΌλ‘ νμλμ§ μμμ΅λλ€. λ©μ μΏ¨ν!
κ°μ₯ μ μ©ν λκΈ
Aaaand μ΄κ²μ μ κ·Έκ²μ΄ λνλ¬λμ§ μ€λͺ ν©λλ€ -- ν¨ν€μ§ λ΄
post.distns
λ #1684μμ utilsλ₯Ό μ§μ° λ‘λ©μΌλ‘ μ νν λκΉμ§ μν¬νλ‘μμ μλμΌλ‘ νμλμ§ μμμ΅λλ€. λ©μ μΏ¨ν!