關於干濕熱浪

Question

如果有人可以協助編寫代碼來計算“相對濕度 > 66% 和 < 33% 的特定熱浪天數”，我將不勝感激。 （然而，熱浪事件被定義為溫度分別至少連續三天超過每日平均溫度的 90% 的事件）。

Answer 1

好的，這是一個解決方案

# temperature percentile
cdo timpctl,90 infile -timmin infile -timmax t2m.nc t2m_pcen90.nc

# mask the temperature
cdo ge t2m.nc t2m_pcen90.nc mask.nc

# Need to make sure we have three consecutive days
cdo --timestat_date last runmean,3 mask.nc mask3.nc
cdo gec,1 mask3.nc heatwave_T.nc 

# Now filter for dry heatwaves, assuming RH is %, change X if fraction
cdo lec,33 rh.nc rhdry.nc
cdo mul heatwave_T.nc rhdry.nc heatwave_dry.nc 

# and same for wet
cdo gec,66 rh.nc rhwet.nc
cdo mul heatwave_T.nc rhwet.nc heatwave_wet.nc 

根據您的定義，當您處於熱浪中時，每個文件中的每個位置/時間都應該有一個 1。 當然，元數據適用於 T2m 而不是索引，如果需要，請使用 NCO 進行更改。 我有幾個視頻指南可以幫助解決這個問題，關鍵的一個是關於掩蔽的（雖然它不包括運行平均值部分）。 另請注意，僅在當天應用 RH 標准（沒有運行平均值），但這就是您在問題中編寫定義的方式。 如果需要，復制運行平均值部分。

ps：一般來說，最好表明您自己嘗試過解決方案，在詢問之前，SO 准則是問題具有調試性質，或者可以是單行請求，但不是像“寫給我”這樣的編碼請求執行 X 或 Y 的代碼”-我認為這就是您被否決的原因。