I am using OpenIFS CY43R1 to simulate the AMORE-France region. I set up the model using the EartH2Observe dataset and run it with ERAI forcing for the period 1995 to 2004. Some months into the simulation, at the beginning of the first April, evaporation starts to explode over the course of a few days at a few grid cells. This seems to be caused by the variable PCFQTI having values larger than 4. in surfexcdriver_ctl_mod.F90 after subroutine VEXCS is called in line 711. Right before that subroutine is called the number of iterations for the solver in VEXCS is set to one (IITT = 1 in line 707). Increasing this number to 10 solves this problem. If I understand correctly, this means that the first solution that is calculated in VEXCS is not appropriate and iteration of the solver is required.

My question is why a value of 1 for IITT has been chosen here and if there is any argument against increasing this number.

Thanks

Stephan

p.s.: If you wish for a minimal working example, please indicate to me the files that you would like me to share.