diff --git a/matlab/AUG/gdat_aug.m b/matlab/AUG/gdat_aug.m
index 3ae36d18548e769f6137de5bd55790c1474094ff..3cc7d98244b1a1f4de6e549ebef36c7c41c55f8c 100644
--- a/matlab/AUG/gdat_aug.m
+++ b/matlab/AUG/gdat_aug.m
@@ -2849,9 +2849,17 @@ elseif ischar(NTIME) || ischar(NTIME_Lpf) || (NTIME > NTIME_Lpf)
   end
   return
 end
-Lpf_tot = Lpf_par.value(1:NTIME); % nb of points: 100000*nb_SOL points + nb_core
-Lpf_SOL = fix(Lpf_tot/100000);
-Lpf1_t = mod(Lpf_tot,100000)+1; % nb of Lpf points
+i_ref = 100000;
+Lpf_tot = Lpf_par.value(1:NTIME); % nb of points: i_ref*nb_SOL points + nb_core
+Lpf_SOL = fix(Lpf_tot/i_ref);
+Lpf1_t = mod(Lpf_tot,i_ref)+1; % nb of Lpf points
+if max(Lpf1_t) > 3000
+  % most probably i_ref wrong
+  i_ref = 10^fix(log10(max(Lpf1_t)));
+  Lpf_tot = Lpf_par.value(1:NTIME); % nb of points: i_ref*nb_SOL points + nb_core
+  Lpf_SOL = fix(Lpf_tot/i_ref);
+  Lpf1_t = mod(Lpf_tot,i_ref)+1; % nb of Lpf points
+end
 
 [equil_time,e]=rdaAUG_eff(shot,DIAG,'time',exp_name_eff,[],gdat_data.gdat_params.extra_arg_sf2sig);
 gdat_data.t = equil_time.value(1:NTIME);