Differences

This shows you the differences between two versions of the page.

--- atmos:citation:research:soda_proecessing_1d_and_2d_hail_2d_hvps1 [2024/05/20 19:14] – klinman
+++ atmos:citation:research:soda_proecessing_1d_and_2d_hail_2d_hvps1 [2024/05/21 14:23] (current) – klinman
@@ Line 48: / Line 48: @@
 {{:atmos:citation:research:soda_conc_comparison.png?900|}}
+The SODA concentrations agree fairly well with the DISP concentrations.
+Notes to self on information to include, mention, explore, etc,
+  * Counts rounded. Ceiling or round up/down? Which is better representation of data?
+  * Time counts forward from time denoted for interval
+  * Idk why the counts are 5x higher. Idk why the correction factor. Idk why the divide by active time. Idk
+  * Sometimes, although I believe this is fixed in my math now, my calculated 2D concentrations are 1.1x higher. But again, I think this isn't an issue anymore.
+  * Sample volumes check out. Airspeed checks out. SODA uses one airspeed for whole time interval. Include a comparison of the numbers
+  * Include a comparison of the different probe specs used in 1D and 2D calculation
+  *
+to get rid of sec from airspeed, i like to think of it as counts per time interval. so 100 counts per 5 second interval would be 100/5s which gives the seconds on bottom to cancel out the seconds in the airspeed. 10 counts per 1 second likewise gives 10/1s. and all of this normalizes the concentration to one second, which is why the 5 second soda intervals and 1 second disp intervals produce the same looking numbers
-Notes for self
-  * np.ceil or np.round
-        * for cases where 1 count, but soda does correction of * 1.1 / 5
-              * so the 1 count would be 1 *1.1 /5 = ~.2     does the .2 count as a count since the /5 is a result of the 5 seconds?
-  * hail counts from disp count forward. so 241841 is 41 42 43 44 45
-  * cant say i understand it but in aarons code, the counts are multiplied by the 1.1 correction factor. then the conc is /5 which is the active time.
-        * what doesn't make sense is that the counts are ~5x higher than the disp file, which is wrong because when the concentration is calculated, it takes into account that 5x factor and becomes decently close to the disp calculated conc
-        * what also doesnt make sense is that the concentration from calculating it myself, using that /5 and *1.1 is that the concentration i calculate is consistently ~1.1 times higher than the concentration aaron gets.
-              * could i be looking at the wrong code
-  * ok update to previous bullet point. just not multiplying by 1.1 seems to be the move???
-        * i think im missing where in aarons code he actually makes the concentration file because something is off here
-              * and im still confused about why the counts should be /5 but thats not the step that I see taken in the code
-  * aaron takes conc1d and prints it to the ascii file and does nothing with it except format it
-Ok so, notes for making this page.
-  * all the counts check out. would be good to show this and include where the time starts
-  * the sample volume checks out. show this too, but not all the soda details. but a comparison of numbers would be good
-  * explain the differene between sodas counts and the disp counts and  how if you divide by 5 theyre the saame but im not sure why. however, it works out in the concentration with the divide by 5
-    * the counts are multipliued by a correction factor 1.1
-  * with all that considered the concs match well and the math was reproducable using the counts from soda and the numbers, and the conc ultimately lines up decently with the 1D
-  * the corr fac is the poisson corr factor
-    * not sure why this is needed tho
-  * the active time correction is for how ever long the chunk is
-    * not sure why it should be divided by 5? figure it should be a sum so no need for division but oh well it works so whatever
-  * why the 126/64 diode change?
-  *