Set max_b1 in .mod header to peak b1, not sysGE.maxRF

On branch develop modified: writemod.m
curtcorum · Dec 20, 2023 · ffc37e3 · ffc37e3
1 parent 3878dcd
commit ffc37e3
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/+toppe/writemod.m b/+toppe/writemod.m
@@ -144,7 +144,8 @@ function writemod(sysGE, varargin)
 dtycyc        = length(find(abs(b1)>0.2236*max(abs(b1)))) / length(b1);
 maxpw         = dtycyc;
 num           = 1;
-max_b1        = sysGE.maxRF;                       	% Gauss. Full instruction amplitude (32766) should produce max_b1 RF amplitude,
+%max_b1        = sysGE.maxRF;        % Gauss. Full instruction amplitude (32766) should produce max_b1 RF amplitude,
+max_b1        = max(abs(b1(:)));     % Gauss
 																		% as long as other RF .mod files (if any) use the same max_b1.
 max_int_b1_sq = max( cumsum(abs(b1).^2)*dt*1e3 );   	% Gauss^2 - ms
 max_rms_b1    = sqrt(mean(abs(b1).^2));              	% Gauss