#!/usr/bin/perl # # makes time points for SAE expt # run in /vd directory ################################ DEFINE INPUT PARAMS BELOW # using exponential increase in timepoints below # t_start * scale**npts = t_max $outfile="stimecho"; $npts = 32; # number of t1 points $t_start=0.000010; # first time point $t_max=10.0000; # max time! $a = $t_max/$t_start; $b = log($a)/log(10); #need log(10) to convert to base-10 in perl $c = $b/($npts-1); $scale = (10**$c); ########################################### open(OUTPUT, "> $outfile"); for($i = 0; $i < $npts; $i++){ $t = ($scale**$i)*$t_start; printf OUTPUT " %3.8f \n", $t; } close (OUTPUT); ########################################### exit;