informatics-lab / precip_rediagnosis

Project to use ML to re-diagnose precipitation fields from ensemble model fields
0 stars 0 forks source link

Add additional model precip variables to extraction script. #35

Closed stevehadd closed 2 years ago

stevehadd commented 2 years ago

We have found some documentation, which links stage output to stash codes, and it looks like the rainfall values we're currently getting from the model do not include convective rainfall. They also do not include snowfall. So we should get all of the following variables as described below from MASS and combine them in the following ways.

So for rainfall rates we should get of these

The sum is the total rainfall.

SImilarly for 3 hour accumulations, we should add these together to get total rainfall

Rainfall is not total precipitation, so we need to also get snowfall, both dynamic and convective

So our main 2 targets should be precip rate and precip accumulation, which has 4 components summed together as our target: precip_rate = rainfall_rate + rainfall_rate_from_convection + snowfall_rate + snowfall_from_convection precip_acc_3hr = rainfall_accumulation-PT03H + rainfall_accumulation_from_convection-PT03H + snowfall_accumulation-PT03H + snowfall_accumulation_from_convection-PT03H

We are also interested in max rainfall rates in a 3 hour period, as indicative of inputs. It might be interesting to get these for snow and rain and doing comparative processing for radar to use these as a target.