"merra_smb_harmonics.py", line 141,
if np.any(MASKS): ---> if bool(MASKS)
2."merra_smb_harmonics.py", line 407
merra_Ylms.month[t] =... ---> merra_Ylms.month =...
3.after running "merra_smb_harmonics.py", it'll gives out 2 files (MERRA2_401_tavgM_2d_SMB_CLM_L60_225.nc & MERRA2_400_tavgM_2d_SMB_CLM_L60_225) for 2020/09 (the reprocessed month). Need to manually delete the 400 one. (and the path in index file)
Hey Tyler, could you add 1 more variable in --product (i.e. MERRA2 hybrid product to calculate) in /model-harmonics/SMB/merra_hybrid_regrid.py?
I need "density". So either add firn height anomaly (h_a) in --product options and I can calculate density myself, or add density option is good. Many thanks
*Now I think about it, maybe adding h_a is better so I have more flexibility