move private flux points psi_norm to max to be outside
2 unresolved threads
2 unresolved threads
Merge request reports
Activity
1997 1997 warning('Found negative values of psi_norm at TS positions, replacing them with 0'); 1998 1998 psi_norm(mask) = 0; 1999 1999 end 2000 params_eff = gdat_data.gdat_params; 2001 params_eff.doplot=0; 2002 params_eff.data_request='z_axis'; 2003 z_axis = gdat(gdat_data.shot,params_eff); 2004 ii=iround_os(psiscatvol.dim{1},z_axis.data); 2005 is_psi_plasma = psi_norm < 1; 2006 is_in_lcfs = zeros(size(psi_norm)); 2007 for it=1:size(is_psi_plasma,2) 2008 is_in_lcfs(ii(it):-1:1,it)=cumprod(is_psi_plasma(ii(it):-1:1,it)); 2009 is_in_lcfs(ii(it):end,it)=cumprod(is_psi_plasma(ii(it):end,it)); 1998 1998 psi_norm(mask) = 0; 1999 1999 end 2000 params_eff = gdat_data.gdat_params; 2001 params_eff.doplot=0; 2002 params_eff.data_request='z_axis'; 2003 z_axis = gdat(gdat_data.shot,params_eff); 2004 ii=iround_os(psiscatvol.dim{1},z_axis.data); 2005 is_psi_plasma = psi_norm < 1; 2006 is_in_lcfs = zeros(size(psi_norm)); 2007 for it=1:size(is_psi_plasma,2) 2008 is_in_lcfs(ii(it):-1:1,it)=cumprod(is_psi_plasma(ii(it):-1:1,it)); 2009 is_in_lcfs(ii(it):end,it)=cumprod(is_psi_plasma(ii(it):end,it)); 2010 end 2011 is_private_flux = is_psi_plasma & ~is_in_lcfs; 2012 % at this stage set private flux "rho" to max, could use 1.3^2 2013 psi_norm(is_private_flux) = max(max(psi_norm));
Please register or sign in to reply