Calibration of hydraulic heads

Hi everyone,

I am trying to model a contaminated site and simulate groundwater flow with feflow 6.1 for 60 years (1950-2010). It's a model with both unsaturated and saturated zone.
I have run the model for groundwater flow. Now I want to calibrate this.
I have 20 observation wells and hydraulic head for 1980 and 1990. I want to calibrate the calculated hydraulic head for this time.
I tried to read the values of calculated heads from Hydraulic-head chart plot. But as it is a very big number of simulation time, It's not very clear. Also, I want to see the scatter plot between actual and calculated heads.

Could anyone please suggest me how can I do this?


Thanks :)

You can derive calculated values stored as history charts by developing a plugin:

[b][font=courier][color=blue]  IfmResult  IfmGetHistoryDimensions (IfmDocument pDoc, IfmHistoryID nID, int nKind, int* nTimes, int* nHist);
  IfmResult  IfmGetHistoryValues (IfmDocument pDoc, IfmHistoryID nID, int nKind, int nTimes, int nHist, double* pTimes, double** pValues);[/color][/font][/b]

Here is a an overview about available histories. The ID is also indicated:

[b][font=courier][color=blue]typedef enum IfmHistoryID {
    IfmHIST_HEAD = 1, /* The 'head distribution' window */
    IfmHIST_AVER_M = 2, /* The 'average mass distribution' window */
    IfmHIST_LOCA_M = 3, /* The 'local mass distribution' window */
    IfmHIST_AVER_H = 4, /* The 'average heat distribution' window */
    IfmHIST_LOCA_H = 5, /* The 'local heat distribution' window */
    IfmHIST_FFLUX = 6, /* The 'fluid flux' window */
    IfmHIST_MFLUX = 7, /* The 'mass flux' window */
    IfmHIST_HFLUX = 8, /* The 'heat flux' window */
    IfmHIST_ITERA = 10, /* Iteration window index */
    IfmHIST_TIMES = 12, /* Time step history index */
    IfmHIST_RESIDUAL = 13, /* Residual window index */
    IfmHIST_VARIO = 15, /* Kriging variogram index */
    IfmHIST_ANA = 16, /* Analytical function index */
    IfmHIST_MULTW_FLUX = 17, /* Muli-layer well flux window index */
    IfmHIST_BHE = 18, /* Borehole heat exchanger window index */
    IfmHIST_IFM = 19, /* IFM based diagram window index */
    IfmHIST_FFLUX_I = 20, /* The 'integral fluid flux' window */
    IfmHIST_MFLUX_I = 21, /* The 'integral mass flux' window */
    IfmHIST_HFLUX_I = 22, /* The 'integral heat flux' window */
    IfmHIST_BHE_IO = 23 /* The Inlet/outlet bhe diagram */
} IfmHistoryID;[/color][/font][/b]

These API functions are shown in C, but can be easily migrated to Python. Please let me know if I may further assist you.