Review index · Download raw source

reference/legacy/run_full_air_grid.jl

1include("photometer_model.jl")2ENV["CLOUDSCAT_DEBUG_REDUCTION"]="1"3out=joinpath(MCROOT,"full_air_grid"); mkpath(out)4open(joinpath(out,"summary.csv"),"w") do io5    println(io,"case_id,tau_over,batch,N,direct_fluence,total_fluence,seconds")6    for (i,τ) in enumerate((0.,.1,.3,1.,3.,10.,30.,100.))7        id="full_tau_"*replace(string(τ),'.'=>'p')8        for b in 1:49            N=τ>=100 ? 25000 : 5000010            r,secs,direct,_,_=makecase(N=N,beta=τ/2000,seed=20263000+10i+b,11                   domainbottom=0.,domainradius=400000.,domaintop=120000.)12            F=sum(r.photo)*DT13            writecurve(joinpath(out,id*"_batch"*string(b)*".csv"),r)14            println(io,join((id,τ,b,N,direct,F,secs),','));flush(io)15            println("FULL AIR ",id," batch",b," seconds=",secs," F=",F);flush(stdout)16        end17    end18end19