I attempted to use the [rainfall_onset2] with my data but encountered some difficulties.
Please, I would appreciate it if you could take a look at it for me.

Hi,

Great work with the python code for climatological analysis of rainfall using satellite data; very helpful and relevant.

I am trying to use the code for a similar analysis, but first, I need to even repeat the code for the Ghana study area to understand exactly what each line of code is doing or trying to achieve.

In line 1 dailyClim function is defined but never used in the code. Can one ignore this function and still get results? Besides, line 4 in that function yields an error since the input data copied in line 3 is 3-d which is been reshaped to a 2-d

In line 9 of getOnsetCessation.py code, in the getDate function, you define a wInd parameter; could you please provide some info on this parameter? Later on, I saw you used wInd = 45 (in line 18) in the getOnsetCessation function. Is wInd an index value, a threshold value relating to date, or should wInd be calibrated per study location?

More questions might follow :)

I will like to commend the details and robustness of the Python code. It is quite useful for my current project.

The current implementation outputs onset and cessation dates as indexes. Is it possible to have both onset and cessation dates as Julian days (number of days from January 1)?