Add python script for generating a TIPMIP CO2 emission file. #14
Conversation
| # unit correction sec month-1 | ||
| dt = (time_bnds[1, idx] - time_bnds[0, idx]) * 86400.0 | ||
| # unit correction - Gt C month-1 to kg CO2 s-1 | ||
| co2flx[idx, :, :] = Edata[idx] * 1e12 * 3.664 / dt / totarea / 12 |
There was a problem hiding this comment.
Why do you leave this block here? I think the protocol specifies a constant emission rate, doesn't it? Then case constant_em == False is not needed and it is confusing that it is there as an option.
(I don't think that there is a discernible difference between the two cases)
There was a problem hiding this comment.
Indeed, the block can be removed for the TIPMIP protocol. I thought it could be of interest for other applications as a reference. If that is not useful and too confusing, I will remove it.
There was a problem hiding this comment.
We have 30 years simulation with the constant_em == False option, I would like to see if this makes a difference at all for the model output (hopefully it doesn't). The script will in any case have to be changed to some extent if someone wants to use it to make a different emission input file.
I would like to add this example on how to create a CO2 emission file.
Some paths are hard coded, so a user will have to modify the script for their own use. However, I think this can be useful as a blue print for similar idealized emission experiments.
For scripts generating input data, I suggest to copy the directory structure used for
inputdataon the shared noresm directory on Betzy.