Hi
I am trying to set up a inverse calibration of a transient problem with Pest(standalone). I run Feflow with the following command
[font=courier]feflow61c -work ..\ -steps steps.pow -dar base_ascii.dar base_ascii.fem[/font]
and I use the -steps option to control the output times in the .dar-file. This works fine 99.9 % of the time, but eventually Feflow will miss an output time with a very small amount e.g. use 79.99999 days instead of 80.00000 days. This will cause Pest to stop, since it is looking for an observation at 80.00000 days.
My Pest instruction file looks like this:
[font=courier]pif %
%RESULTS% %TIME = 8.000000e+01%
%----------%
%----------%
l6 [obswell1_1]4:18[/font]
And Feflow produces a .dar-file like this:
[font=courier]***************************************************
RESULTS AT STEP = 64 AND TIME = 7.999999e+01 [d]:
---------------------------------------------------
Obs H [m] Vx [m2/d] Vy [m2/d]
---------------------------------------------
1 -1.452183e-01 -1.798366e-02 1.786912e-02
2 -3.418914e-01 -6.729764e-02 5.468340e-02
3 -5.286261e-01 -1.310629e-01 -8.419374e-02
4 -5.296815e-01 -1.313571e-01 -8.475472e-02[/font]
Time steps is supplied like this
[font=courier]# 1
! Timesteps
.
.
66.36 0.01
80 13.64
80.01 0.01
.
END
END[/font]
I cannot use step number to identify output at a specific time since I use automatic time step control and the number of steps can change from model run to model run.
Is there a way to force Feflow to use the timesteps specified with -steps? Alternatively can the output times in the .dar file be written with fewer decimals?
Best regards
Magnus