search

uataq / X-STILT

Source codes for X-STILT(NOx) [Wu et al., 2018, 2023] based on STILT-R version 2 [Fasoli et al., 2018] and STILT [Lin et al., 2003].
16 stars 11 forks source link

M1 Boundary Conditions with TCCON footprint #10

Open JPMastrogiacomo opened 1 week ago

JPMastrogiacomo commented 1 week ago

Hi, I am relatively new to X-STILT but I was able to successfully run it with TCCON measurements and get footprint and trajectory files. I am attempting to calculate the background with the M1 method and having trouble with trajec.edpt.ctnrt.r specifically. Part 1 works great and it successfully gets the end points from carbon tracker. Part 2 is supposed to get co2 from carbon tracker for oco-2/tccon levels above the released particles. But it is based on output$wgt.prof which does not seem to exist in the trajectory files anymore. Perhaps the variables that we need for this part (press, ak.pwf, pwf, ak.norm, etc) should be attached to the output$receptor or output$qt_prof variables? If you point me to the right file I can try to make those changes.

One more question. I always seem to get fewer particles than set by numpar in run_xstilt_tccon.r. And each trajectory has a different number of particles. Is this intended behaviour? Thanks so much!

JPMastrogiacomo commented 1 week ago

I have figured out a solution. In wgt.trajec.foot.tccon.r I am saving the output of get.wgt.tccon.func as the variable output$combine.prof and that gives me the information I need to calculate the background.

For most of my my trajectories I am getting the correct number of particles, but for some I get less than half which is strange.