Project

General

Profile

Biologically effective degree days calculation

Added by Antonio Beltran almost 2 years ago

Hello everyone,

I have a file .nc, with daily mean temperature data (tas) for 10 years and I want to calculate for each year this equation (Biologically effective degree days - Β°C):
𝐡𝐸𝐷𝐷 = βˆ‘π‘šπ‘–π‘›[π‘šπ‘Žπ‘₯[𝑇𝐺𝑖𝑗 βˆ’ π‘‡π‘™π‘œπ‘€ , 0], π‘‡β„Žπ‘–π‘”β„Ž βˆ’ π‘‡π‘™π‘œπ‘€] where TGii is the daily mean temperature at day i of period j. The π‘‡β„Žπ‘–π‘”β„Ž and π‘‡π‘™π‘œπ‘€ are constants: Thigh=25Β°C and Tlow=0Β°C. Is that possible with CDO? because it is difficult for me to calculate spatial data (200 grid points) with R or Python.

Thanks,
Antonio B.


Replies (1)

RE: Biologically effective degree days calculation - Added by Estanislao Gavilan almost 2 years ago

Hi Antonio,

I think it is possible. In fact, there are many examples on this forum about data processing. You might need to divide the equation in small steps though. I am really surprise you cannot use python, 200 points is a small domain. If it is a ncfile, you can loop the time axis of the ncfile, and then you just need to do the sum up. Therefore, you only store one dimesional grid (BEDD). You can easily open 720*1440 grid points in an ordinary computer if you only call one time step.

Kind regards,

Estanislao

    (1-1/1)