diff --git a/crpptbx/AUG/gdat_aug.m b/crpptbx/AUG/gdat_aug.m
index 9eccee661f086ab68724e2115480ad9433ff5117..5b772e88377874e20435ae660e2696f442f87993 100644
--- a/crpptbx/AUG/gdat_aug.m
+++ b/crpptbx/AUG/gdat_aug.m
@@ -918,7 +918,7 @@ elseif strcmp(mapping_for_aug.method,'switchcase')
       [a,error_status]=rdaAUG_eff(shot,'VTA','Z_edge',exp_name_eff);
       gdat_data.(lower(node_child_nameeff_e)).z = repmat(a.data',1,size(gdat_data.(lower(node_child_nameeff_e)).data,2));
       nb_edge = size(gdat_data.(lower(node_child_nameeff_e)).z,1);
-      iaaa=iroundos(gdat_data.(lower(node_child_nameeff_e)).t,gdat_data.(lower(node_child_nameeff)).t);
+      iaaa=iround_os(gdat_data.(lower(node_child_nameeff_e)).t,gdat_data.(lower(node_child_nameeff)).t);
       gdat_data.data(nb_core+1:nb_core+nb_edge,:) = gdat_data.(lower(node_child_nameeff_e)).data(1:nb_edge,iaaa);
       gdat_data.dim{1}(nb_core+1:nb_core+nb_edge,:)=gdat_data.(lower(node_child_nameeff_e)).z(1:nb_edge,iaaa);
       [alow,e]=rdaAUG_eff(shot,'VTA',[upper(data_request_eff(1)) 'elow_e'],exp_name_eff);
@@ -1019,7 +1019,7 @@ elseif strcmp(mapping_for_aug.method,'switchcase')
     gdat_data.(lower(node_child_nameeff_e)).rhovolnorm = rhovolnorm_out_edge;
     % put values of rhopolnorm for dim{1} by default, all radial mesh for combined core, edge in grids_1d
     gdat_data.x = gdat_data.(lower(node_child_nameeff)).rhopolnorm;
-    iaaa=iroundos(gdat_data.(lower(node_child_nameeff_e)).t,gdat_data.(lower(node_child_nameeff)).t);
+    iaaa=iround_os(gdat_data.(lower(node_child_nameeff_e)).t,gdat_data.(lower(node_child_nameeff)).t);
     gdat_data.x(inb_chord_core+1:inb_chord_core+inb_chord_edge,:) = gdat_data.(lower(node_child_nameeff_e)).rhopolnorm(:,iaaa);
     gdat_data.dim{1} = gdat_data.x;
     gdat_data.dimunits{1} = 'rhopolnorm';