function to run line command torbeam via tbdemo and to read output file

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@4420 d63d8f72-b253-0410-a779-e742ad2e26cf

crpptbx/AUG/run_torbeam_fullshot.m

+function [torbeam_out_struct, file_out_prefix]=run_torbeam_fullshot(shot,varargin);
+%
+% Run tbdemo, using line command execution of TORBEAM, to run torbeam acroos full shot on reduced time axis
+%
+% shot: shot number for which to run Torbeam. 
+%       If empty or not given, ask for it
+%       if negative, assume output file exists, do not re-run TORBEAM and just reads data into structure torbeam_out_struct
+%
+% varargin{1:2}: 'time',time_array: time points for which to run Torbeam (by default run every 50ms for each ECRH2 gyrotron having power>0)
+% varargin{3:4}: 'file_prefix', filename_prefix: results are saved on [filename '_gyroXX.dat'] (default: tbdemo_shot_gyroXX.dat with XX being the gyro nb)
+%
+% Examples:
+%      torbeam_out_struct = run_torbeam_fullshot(shot);   % standard call to run Torbeam
+%      torbeam_out_struct = run_torbeam_fullshot(-shot);  % standard call to just read the data from the local file from a previous run
+%
+%      torbeam_out_struct = run_torbeam_fullshot(shot,'time',linspace(tstart,tend,nbpoints)); % to run on given time array
+%      torbeam_out_struct = run_torbeam_fullshot(shot,'time',[t1 t2 t3...],'file_prefix','tbdemo_shot_fewtimes_gyro'); % to run on specific times
+
+if ~exist('shot')
+  shot=input('shot: ');
+end
+
+time_array=linspace(0,10,10/0.05);
+file_out_prefix = ['tbdemo_' num2str(abs(shot)) '_gyro'];
+
+if nargin>=3 
+  for i=1:2:length(varargin)-1
+    if ~isempty(varargin{i}) && ~isempty(varargin{i+1})
+      switch lower(varargin{i})
+       case {'time'}
+	time_array = varargin{i+1};
+       case {'file_prefix'}
+	file_out_prefix = varargin{i+1};
+       otherwise
+	disp(['case ' lower(varargin{i}) ' not yet implemented, ask Olivier.Sauter@epfl.ch'])
+	return
+      end
+    end
+  end
+end
+
+if shot>0
+  pgyro=gdat(shot,'pgyro',1);
+
+  for igyro=5:8
+    filename=[file_out_prefix num2str(igyro) '.dat'];
+    unix(['rm ' filename ' >& /dev/null']); % since cannot write over an existing file by default
+    unix(['rm ' '.' filename ' >& /dev/null']); % since cannot write over an existing file by default
+    unix(['touch ' '.' filename]);
+    itt=find(pgyro.data(:,igyro)>1e5);
+    if ~isempty(itt)
+      itt_tok=find(time_array>=pgyro.t(itt(1)) & time_array<=pgyro.t(itt(end)));
+      if ~isempty(itt_tok)
+	pgyro_tok=interp1(pgyro.data(itt,igyro),time_array(itt_tok));
+	disp(['running system ' num2str(igyro) '...']);
+	for j=1:length(itt_tok)
+	  if pgyro_tok(j)>1e5
+	    [a,b]=unix(['ssh sxaug21.aug.ipp.mpg.de ~rem/public/tbm_demo/tbmdemo -batch -silent -shot ' num2str(shot) ' -time ' num2str(time_array(itt_tok(j))) ' -system ' num2str(igyro) ' -usedcr -eqdiag EQI -betapol -usemirror -cdrive 1 >> ' ['.' filename]]);
+	  end
+	end
+	[a,b]=unix(['grep -v ''could not''' ' .' filename '| grep -v ''invalid'' > ' filename]);
+      end
+    end
+  end
+end
+
+for igyro=5:8
+  try
+    eval(['[gyro_tbnam,tbout{igyro}.gyro,shot_tbnam,tbout{igyro}.shot,time_tbnam,tbout{igyro}.time,pol_tbnam,tbout{igyro}.pol,tor_tbnam,tbout{igyro}.tor,' ...
+	  'pow_tbnam,tbout{igyro}.pow,r_tbnam,tbout{igyro}.r,z_tbnam,tbout{igyro}.z,pwr_dep_tbnam,tbout{igyro}.rhopol_pwr_dep,rhopol_cd_dep_tbnam,tbout{igyro}.rhopol_cd_dep,jcd_peak_tbnam,' ...
+	  'tbout{igyro}.jcd_peak,rhopol_c_wid_tbnam,tbout{igyro}.rhopol_c_wid,pni_tbnam,tbout{igyro}.pni,picr_tbnam,tbout{igyro}.picr,prad_tbnam,tbout{igyro}.prad,' ...
+	  'betan_tbnam,tbout{igyro}.betan,ipl_tbnam,tbout{igyro}.ipl,bt_tbnam,tbout{igyro}.bt,a1_tbnam,tbout{igyro}.a1,a2_tbnam,tbout{igyro}.a2,rhopolntm_tbnam,tbout{igyro}.rhopolntm]=' ...
+	  'textread(''' file_out_prefix  num2str(igyro) '.dat'',''%s%d%s%d%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f%s%f'');']);
+  catch
+    disp(['problems reading file ' file_out_prefix  num2str(igyro) '.dat: may be some errors with tbdemo, check file, skipped'])
+  end
+end
+
+torbeam_out_struct = tbout;