Project

General

Profile

All events above 95th percentile

Added by Francis Nkrumah over 6 years ago

I would like to know how to generate a file containing all daily rainfall events above the 95th percentile, not only the 95th percentile value but all the days above that value


Replies (6)

RE: All events above 95th percentile - Added by Ralf Mueller over 6 years ago

Hi Francis!

If you have the 95th percentile field, you can substract this field from the input data. All points with values larger than zero are above the 95th percentile. If you want, you can create a 0-1-mask out of that, then each location has its own timeseries, being 1 for a day above the 95th percentile value at that locations.

hth
ralf

RE: All events above 95th percentile - Added by Francis Nkrumah over 6 years ago

Please Ralf, can you add a sample code on how to get this done?

RE: All events above 95th percentile - Added by Ralf Mueller over 6 years ago

cdo -gtc,0 -sub <dailyData> <95th percentile field> <outputMaskFile>

Selecting corresponding events - Added by Francis Nkrumah over 6 years ago

Hello

I'm having two netcdf files: one contains days with precipitation events above the 95th percentile and the other is daily zonal wind. I would like to select the corresponding days of the zonal wind when the precipitation events occurred. Please is there a code that can help me generate the corresponding events?

RE: All events above 95th percentile - Added by Ralf Mueller over 6 years ago

I suspect, that both files have the same horizontal grid.

if you have hte 0-1-mask files, you can multiply it with your zonal data. so that everything non-zero in the results file do belong to the desired events/days

cdo mul  <zonalWinds> <above-95th-perc-events-mask> <zonamWindsOnDaysWithEvents>

BUT: it can only be done, if your 0-1-mask file is on a the same horizontal grid as the zonal data. it this makes sense is up to you

hth
ralf

    (1-6/6)