Matching multiple parameters to target moments

I know some steady state values (for example, Lambda_S = 1, rgov_S = r_S or rk_S = (Gamma/beta + 1 + deltaK), and others), but it is not enough. For capital, for example, I used kk_S = Kkmean; and kb_S = Kbmean as initial values (drawn from data).

  1. In the mod file, should I create a block with steady_state_model for the variables I know and initval block for the variables I do not know but have an approximation?
  2. I still do not understand the loop according to steady state guesses. I was looking at Parameter loop in Dynare 5.x - #3 by jpfeifer, Loop over parameters - #67 by cmarch and Looping over initial values - guide/example code - #6 by jpfeifer (with code matlab/evaluate_steady_state.m · master · Dynare / dynare · GitLab), to test the ideas presented. I tried to simplify it to just see the changes in the steady state value of consumption c, but c_storage remains empty. That loop starts in line 414 of Loop_ss_eval.m. In line 291, I have the sketch of steady states and, in 411, I invoke dynare.
    SteadyStateLoop.zip (45.7 KB)