Merge pull request #26 from GeoscienceAustralia/single-file-writes

Share output files between chains

Project.toml

@@ -32,7 +32,6 @@ StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
-julia = "1.5.3"
 DataInterpolations = "3.6.1"
 Distances = "0.10.7"
 DistributedArrays = "0.6.6"
@@ -50,4 +49,4 @@ PyPlot = "2.10.0"
 Roots = "2.0.0"
 SpecialFunctions = "1.6 - 2.1.4"
 StatsBase = "0.33.16"
-
+julia = "1.5.3"

examples/1D/stat/02_run_lineregression_highλ.jl

@@ -39,7 +39,7 @@ opt = transD_GP.OptionsStat(nmin = nmin,
                         peskycholesky = true
                         )
 ## set up McMC
-nsamples, nchains, nchainsatone = 50001, 4, 1
+nsamples, nchains, nchainsatone = 5001, 4, 1
 Tmax = 2.50
 addprocs(nchains)
 @info "workers are $(workers())"

examples/AEM_MT_joint/AEM_only/02_make_aem_inversion_opts_s.jl

@@ -1,4 +1,3 @@
-using transD_GP, Distributed
 ## make options for the stationary GP
 fileprefix = "OrstUhn_"
 K = transD_GP.GP.OrstUhn()

examples/AEM_MT_joint/AEM_only/03_run_aem_inversion_s.jl

@@ -1,3 +1,4 @@
+using Distributed
 ## set up McMC
 nsamples, nchains, nchainsatone = 100001, 5, 1
 Tmax = 2.50

examples/AEM_MT_joint/MT_only/allfreqs/01_make_synthetic.jl

@@ -23,7 +23,10 @@ z, ρ, nfixed = transD_GP.makezρ(zboundaries; zfixed=zfixed, ρfixed=ρfixed)
 ρ[(z.>=120) .&(z.<140)] .= 2
 ρ[(z.>=140) .&(z.<250)] .= 100.
 ρ[(z.>=250) .&(z.<270)] .= 100.
-ρ[(z.>=270)]            .= 150
+ρ[(z.>=270) .&(z.<600)] .= 150
+ρ[(z.>=600) .&(z.<800)] .= 2
+ρ[(z.>=800) .&(z.<5000)] .= 500
+ρ[z.>=5000] .= 700
 # add jitter to model in log10 domain
 Random.seed!(11)
 ρ[1:50] = 10 .^(0.1*randn(length(ρ[1:50])) + log10.(ρ[1:50]))

examples/AEM_MT_joint/MT_only/allfreqs/03_run_MT_inversion_s.jl

@@ -1,6 +1,6 @@
 ## set up McMC
 using Distributed
-nsamples, nchains, nchainsatone = 300001, 4, 1
+nsamples, nchains, nchainsatone = 100001, 4, 1
 Tmax = 2.50
 addprocs(nchains)
 @info "workers are $(workers())"

examples/AEM_MT_joint/MT_only/limitedfreqs2/01_make_synthetic.jl → examples/AEM_MT_joint/MT_only/highfreqs/01_make_synthetic.jl

@@ -23,7 +23,10 @@ z, ρ, nfixed = transD_GP.makezρ(zboundaries; zfixed=zfixed, ρfixed=ρfixed)
 ρ[(z.>=120) .&(z.<140)] .= 2
 ρ[(z.>=140) .&(z.<250)] .= 100.
 ρ[(z.>=250) .&(z.<270)] .= 100.
-ρ[(z.>=270)]            .= 150
+ρ[(z.>=270) .&(z.<600)] .= 150
+ρ[(z.>=600) .&(z.<800)] .= 2
+ρ[(z.>=800) .&(z.<5000)] .= 500
+ρ[z.>=5000] .= 700
 # add jitter to model in log10 domain
 Random.seed!(11)
 ρ[1:50] = 10 .^(0.1*randn(length(ρ[1:50])) + log10.(ρ[1:50]))

examples/AEM_MT_joint/MT_only/limitedfreqs2/03_run_MT_inversion_s.jl → examples/AEM_MT_joint/MT_only/highfreqs/03_run_MT_inversion_s.jl

@@ -1,6 +1,6 @@
 ## set up McMC
 using Distributed
-nsamples, nchains, nchainsatone = 300001, 4, 1
+nsamples, nchains, nchainsatone = 100001, 4, 1
 Tmax = 2.50
 addprocs(nchains)
 @info "workers are $(workers())"

examples/AEM_MT_joint/MT_only/useful/03_run_MT_inversion_s.jl → examples/AEM_MT_joint/MT_only/lowfreqs/03_run_MT_inversion_s.jl

@@ -1,6 +1,6 @@
 ## set up McMC
 using Distributed
-nsamples, nchains, nchainsatone = 300001, 4, 1
+nsamples, nchains, nchainsatone = 100001, 4, 1
 Tmax = 2.50
 addprocs(nchains)
 @info "workers are $(workers())"