search

mrghg / py12box_invert

MIT License
3 stars 0 forks source link

Add firn data #14

Open mrghg opened 3 years ago

mrghg commented 3 years ago

Need to add a non-linear solver to cope with this. Will add some example input files.

The relevant snippets of IDL code:

# For one firn measurement, interpolating the concentration Greens functions to our coordinates:
y_greens=interpol(concGreens[fi, *], year_firn, parameters['STARTY'] + $
          indgen(parameters['NYEARS']*12)/12.)

# Modelled firn concentration:
total(y['MODEL', why]*y_greens)/total(y_greens) 

# Emissions sensitivity:
total(H[wh[whi], why]*y_greens)/total(y_greens)

# Approximate weighted timestep of firn measurement:
total(findgen(n_elements(y_greens))*y_greens)/total(y_greens)
mrghg commented 3 years ago

Have added an example input to py12box_agage repo under CFC-115.