Sample in sceaua

[XiaoXiaoLeisure]

Hi, I am using SPOTY's SCEUA algorithm to optimize the VIC model. I want to keep the parameter combinations of each calculation, but I find that the number of parameter combinations in the output csv file is less than iterations, for example, iterations=500, len(params)=208.

sampler = spotpy.algorithms.sceua(cal_setup, dbname=name, dbformat='csv')
sampler.sample(iterations)

results_df = pd.DataFrame(np.array(sampler.getdata()))
columns_to_save = results_df.columns[:7] # 获取前7列的列名
filtered_df = results_df[columns_to_save]
filtered_df.columns = ['like1', 'binfil', 'Ds','Ds_max', 'Ws', 'Soil_D2', 'Soil_D3']
filtered_df.to_excel(outCal + paramfile, sheet_name='param', index=False)

How can I get the VIC parameters for each calculation?

[krishnakafle]

I also faced the same issue and am trying to figure out why this is happening. But for the moment, I am writing a script to save parameters for each iteration. Let us know if you have any updates or a better method for this.

    def simulation(self,vector):
        simulations= self.vicmodel.run_vic(
            binfilt=vector[0],
            Ws=vector[1],
            DsMax=vector[2],
            Ds=vector[3],
            soil_d2=vector[4],
            soil_d3=vector[5], 

            c=vector[6],

            exp1=vector[7],
            exp2=vector[8],
            exp3=vector[9],
            
            )
        global iter_count
        rep = iter_count  # This will now reflect the current iteration

        iter_folder = f'output/SAVE_{rep}'
        if not os.path.exists(iter_folder):
            os.makedirs(iter_folder)
        content_to_write = 'Paramters Used for the simulation ...\n \
            {0}, {1}, {2}, {3}, {4}, {5}, {6}, {7}, {8}, {9}'.format(vector[0], vector[1], vector[2], vector[3], vector[4], vector[5], vector[6], vector[7], vector[8], vector[9])

        with open(f"{iter_folder}/param.txt", "w") as file:
            file.write(content_to_write)

        return simulations

Edit:
I am keeping track of the iterations with a global variable (iter_count) and updating its value by 1 inside the objective function.