equil_out_ids%time_slice(it)%boundary%psi_norm=equil_in_itm(it)%global_param%psi_bound! since value at separatrix not known, may need to use namelist choice for shaving
equil_out_ids%time_slice(it)%boundary%psi_norm=equil_in_itm(it)%global_param%psi_bound! since value at separatrix not known, may need to use namelist choice for shaving