Fail to find perfect foresight solution

Dear all,

I am trying to run a deterministic simulation with a slightly extended version of Faria-e-Castro (2021) COVID-19 model. I believe I am working with correct equiilibrium conditions and plausible parameters. However, Dynare still fails to find the perfect foresight solution even with quite small (0.01) shocks to any exogenous varibale. I would really appreciate it if anyone could help me figure out what is going wrong.

The files I am working with are attached below. Please run main.m to conduct the simulation:
baseline_alterfec.mod (9.6 KB)
calibration_alterfec_alterw_alterck.m (6.3 KB)
main.m (483 Bytes)

I cannot run your code due to different file namings.

I am sorry. I have deteled redundant information and renamed the files. It should work now.
baseline_alterfec.mod (9.6 KB)
calibration.m (6.3 KB)
main.m (74 Bytes)

Your model does not satisfy the Blanchard-Kahn conditions:

There are 11 eigenvalue(s) larger than 1 in modulus
for 10 forward-looking variable(s)

The rank condition ISN’T verified!
and therefore shows explosive behavior.

Really appreciate your help! Now I have got the problem solved.