anchor.diff:          ! load parameters from namelist
           call surface_physics_par_load(surface%par,trim(nml_file))
           ! open file for CRMSD output
  -        write(output,"(A,I3.3,A)") trim(prefix), i,  ".out"
  +        write(output,"(A,I6.6,A)") trim(prefix), i,  ".out"
           open(2,file=trim(output),form='formatted')
  +        if (mod(i,100) == 0) write(*,*) "\x1B[32mrun_particle:\x1B[0m Load file ", trim(nml_file)
   
   
           do k=1,nloop ! re-iterate 'nloop' times
  
positive.code: 
      ! name list to drive model
      namelist /driver/ nloop, ntime, input_forcing, output, validation 
  
      write(*,*) "\x1B[32mrun_particles:\x1B[0m start particle calculation."
      prog_name = "run_particles"
  
negative.code:         cost_smb = calculate_crmsd(state%smb,vali%smb,ntime)
          cost_alb = calculate_crmsd(state%alb,vali%alb,ntime)
          cost = dsqrt(cost_alb**2+cost_stt**2+cost_smb**2)
          write(2,*) cost_smb
          close(2)
      end do
          
  