wadewei Just to confirm that md.smb.precipitation should be filled with each months' precipitation? i.e. sum(md.smb.precipitation(i,1:12))=one year's precipitation at vertex i
mathieumorlighem If you type md.smb you will see that the unit of precipitation is m/year so you need to use the function mean and not sum
wadewei mathieumorlighem That makes sense. However I am getting calculated smb that is greater than precipitation by as much as 20%. Here's a plot showing a large domain has smb./precipitation>1 shorturl.at/aGKXY smb1=md1.results.TransientSolution(1).SmbMassBalance; [value pos]=max(smb1./mean(md1.smb.precipitation(1:end-1,1:12),2)) smb1(pos) md1.smb.precipitation(pos,1:12) value = 1.1909 pos = 6207 ans = 1.9118 ans = 1.9332 0.9310 1.8275 1.1984 0.9085 2.1312 2.4081 2.2663 1.6579 1.6181 1.7517 0.6334
mathieumorlighem Hi Wade. I am not an expert in PDD but the precipitations are in m/yr water equivalent while the SMB is in m/yr ice equivalent. The difference in densities may explain what you are observing.
wadewei mathieumorlighem Hi Mathieu. That's a good point. However with default density difference of 10%, it's not enough to explain smb/precipitation of approx. 120%
