A question about CWatM output variable: totalPotET

[ray-yicong]

I want to ask you a question about CWatM output, recently, I wanted to output varible called totalPotET, but failed, it came to a erros said that :
File "C:\ProgramData\Anaconda3\envs\ycr372\lib\site-packages\numpy\ma\core.py", line 3385, in setitem
_data[indx] = dval
TypeError: NumPy boolean array indexing assignment requires a 0 or 1-dimensional input, input has 2 dimensions

However, I tried other variables related to evapotranspiration are fine, for example：totalET, EWRef.

Can you help me find the problem?

[mikhailsmilovic]

Guten Tag @ray-yicong,

totalPotET is a land class-specific variable. That means it can be output as totalPotET[0] for the m/day over forests, totalPotET[1] for the m/day over grasslands, etc. Note then that these are not exactly in units of m, as they are m over the specific land class.

Outputting variables with square brackets do not work with all output types.
OUT_MAP_Daily and OUT_Map_MonthAvg will work fine.
OUT_Map_AnnualAvg and "total" maps (ex., MonthTot) will not work, and we must do something else, but perhaps these two output types are sufficient. This "issue" has been included here: #174

[ray-yicong]

Sorry i still can't output the totalPotET, i tried both OUT_MAP_Daily and OUT_Map_MonthAvg, the model can run, but it stopped at the second month after the spinup date with the same error: these are the related simulation settings

and this is the error

What is the problem? Can you please help me see again?

[mikhailsmilovic]

Guten Morgen, land class-specific variables are to be output indicating the land class, for example
OUT_MAP_Daily = totalPotET[0], totalPotET[1], totalPotET[2], totalPotET[3]
OUT_Map_MonthAvg = totalPotET[0], totalPotET[1], totalPotET[2], totalPotET[3]

Note then that these are not exactly in units of m, as they are m over the specific land class.
Specifically, m3/(m2 of land class)

0 = forest
1 = grassland
2 = paddy
3 = non-paddy irrigated land

[ray-yicong]

Problem solved, it worked！I didn't specify the type of landcover I wanted to output before. That's why the error said that can't has two dimensions.

Thank you again for your detailed answer!

[ray-yicong]

Thank you for your detailed answer! 易聪睿 ***@***.***  

…

------------------ 原始邮件 ------------------ 发件人: "iiasa/CWatM" ***@***.***>; 发送时间: 2024年12月5日(星期四) 晚上11:38 ***@***.***>; ***@***.******@***.***>; 主题: Re: [iiasa/CWatM] A question about CWatM output variable: totalPotET (Discussion #173) Guten Tag @ray-yicong, totalPotET is a land class-specific variable. That means it can be output as totalPotET[0] for the m/day over forests, totalPotET[1] for the m/day over grasslands, etc. Note then that these are not exactly in units of m, as they are m over the specific land class. Outputting variables with square brackets do not work with all output types. OUT_MAP_Daily and OUT_Map_MonthAvg will work fine. OUT_Map_AnnualAvg and "total" maps (ex., MonthTot) will not work, and we must do something else, but perhaps these two output types are sufficient. — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you were mentioned.Message ID: ***@***.***>