Define variable(s) for Funcfit

Function Fitting_Engine(p,x)
    wave p
    variable x
   
    nvar ci
    nvar uc_occ_Y, uc_occ_Zr,  uc_occ_O, ci, u_Y, u_Zr, u_O
    wave ASF_Y, ASF_Zr, ASF_O, uc_position_Y, uc_position_Zr, uc_position_O

// ATOMIC POSITION //

// PERFECT UNIT CELL //
    variable/c FF_uc, FF_uc_Y, FF_uc_Zr, FF_uc_O, FF_ctr
       
    FF_uc_Y = atomic_frac_Y * ASF_Y * uc_occ_Zr * exp(ci*x*(temp_first_Y_position))*exp(-0.5*(x*u_Y)^2)
    FF_uc_Zr = atomic_frac_Zr * ASF_Zr * uc_occ_Zr * exp(ci*x*(temp_first_Zr_position))*exp(-0.5*(x*u_Zr)^2)
    FF_uc_O = 2*atomic_frac_O * ASF_O * uc_occ_O * exp(ci*x*(temp_first_O_position))*exp(-0.5*(x*u_O)^2)
   
    FF_uc = FF_uc_Y + FF_uc_Zr + FF_uc_O   
    FF_ctr = 1 / (1-exp(-ci*x*c))

// DISTORTED UNIT CELLS //
    variable i_F_duc
    variable/c FF_duc_Y, FF_duc_Zr, FF_duc_O, FF_duc   
    wave duc_occ_O
   
    FF_duc_Y = 0
    FF_duc_Zr = 0
    FF_duc_O = 0
   
    for (i_F_duc=0;i_F_duc<numpnts(temp_duc_position_O_reset);i_F_duc+=1)
        FF_duc_Y += atomic_frac_Y * ASF_Y* p[i_F_duc] * exp(ci*x*temp_duc_position_Y_reset[i_F_duc])*exp(-0.5*(x*u_Y)^2)
        FF_duc_Zr += atomic_frac_Zr *ASF_Zr * p[i_F_duc] * exp(ci*x*temp_duc_position_Zr_reset[i_F_duc])*exp(-0.5*(x*u_Zr)^2)
        FF_duc_O += 2*atomic_frac_O *ASF_O * duc_occ_O[i_F_duc] * exp(ci*x*temp_duc_position_O_reset[i_F_duc])*exp(-0.5*(x*u_O)^2)
    endfor

    FF_duc = FF_duc_Y + FF_duc_Zr + FF_duc_O
   

// INTERFACIAL MOLECULES //    
    variable/c FF_ads1, FF_ads2, FF_ads3, FF_ads
    nvar d_spacing_water

    FF_ads1 = ASF_Y*p[8]*exp(ci*x*(p[9]))*exp(-0.5*(x*p[10])^2)
    FF_ads2 = ASF_Y*p[11]*exp(ci*x*(p[12]))*exp(-0.5*(x*p[13])^2)
    FF_ads3 = ASF_Y*p[14]*exp(ci*x*(p[15]))*exp(-0.5*(x*p[16])^2)

    FF_ads = FF_ads1 + FF_ads2 + FF_ads3
   
    variable/c FF_adw1, FF_adw2, FF_adw, FF_lw
    wave ASF_H2O

    FF_adw1 = ASF_H2O * p[17] * exp(ci*x*(p[18]))*exp(-0.5*(x*p[19])^2)
    FF_adw2 = ASF_H2O * p[20] * exp(ci*x*(p[21]))*exp(-0.5*(x*p[22])^2)
    FF_adw = FF_adw1 + FF_adw2
    FF_lw = (ASF_H2O * exp(ci*x*p[23]) * exp(-0.5*(x*p[24])^2)) / (1-exp(-0.5*(x*p[25])^2)*exp(ci*x*d_spacing_water))

// MAGNIFICATION FACTORS //
    variable FF_mag, FT_Cell, FSpecular, FB
    nvar cell_water_thickness, photon_wavelength, attenuation_length, r_e, suc_area
   
    FT_cell = exp((-8*pi*cell_water_thickness)/(x*photon_wavelength*attenuation_length))
    FB = (1-p[26])^2/(1+p[26]^2-2*p[26]*cos(x*c))
    FSpecular = (4*pi*r_e/(x*suc_area))^2

    FF_mag = p[27] * FB * FT_Cell * FSpecular
   
    return FF_mag * magsqr(FF_uc*FF_ctr + FF_duc + FF_ads + FF_adw + FF_lw)
End