search

ciceroOslo / ciceroscm

Python version of the CICERO-SCM simple climate model/emulator
Apache License 2.0
19 stars 4 forks source link

Add read methods for LUC, volcanic and solar forcing #63

Closed benmsanderson closed 2 years ago

benmsanderson commented 2 years ago

Current build allows for passing of emissions, concentrations and gas parameters dataframes as specified in cfg using "concentrations_data", "emissions_data", "gaspam_data" - which is very useful.

Idealised runs are limited by hard coding of Land use change, Solar and volcanic forcing text files which have finite length, and for some idealised runs it is preferable to have these forcings as constants.

A similar interactive treatment would address this issue and allow - ie:

cscm=CICEROSCM({ "luc_data": df_luc, "volc_data": df_volc, "solar_data" df_solar, })

where df_luc, df_volc and df_solar are pandas dataframes (or whatever)

ragnhibs commented 2 years ago

And I have made AR6 timeseries input files that can be used for the three sets here: /div/qbo/utrics/RadiativeForcing/RFforSCM

ragnhibs commented 2 years ago

Still not able to run with own solar data:

"rf_sun_data":pd.read_csv("/div/qbo/utrics/RadiativeForcing/RFforSCM/solar_erf_ar6.txt", header=None, skiprows=1, index_col=0)

return series._values[index] IndexError: index 270 is out of bounds for axis 0 with size 270

For LUC and volc, I do not get any error message, but I do not think the data is used...

ragnhibs commented 2 years ago

Do not think I got the others to work either...

Could be useful with an error message if your input to the code is not an input that are used... Misspelling can happen...