From 3a6dd0f86f97efe8a1d246b8f4634377ca9dd632 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Thu, 6 Jun 2019 19:24:07 +0000
Subject: [PATCH] add LD_LIBRARY_PATH for running torbeam from matlab
git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@12004 d63d8f72-b253-0410-a779-e742ad2e26cf
---
crpptbx/AUG/torbeam_plot_output_singlerun.m | 8 ++++++--
crpptbx/AUG/torbeam_prepare_inputs_and_run.m | 17 +++++++++++------
2 files changed, 17 insertions(+), 8 deletions(-)
diff --git a/crpptbx/AUG/torbeam_plot_output_singlerun.m b/crpptbx/AUG/torbeam_plot_output_singlerun.m
index a7a34ec8..b4424caf 100644
--- a/crpptbx/AUG/torbeam_plot_output_singlerun.m
+++ b/crpptbx/AUG/torbeam_plot_output_singlerun.m
@@ -66,8 +66,12 @@ for i=1:length(zzz.tb_out_files_to_cp_var)
catch
end
end
-zzz.tbeg_abs=load('cntpr.dat');
-zzz.tbeg_abs(2) = zzz.tbeg_abs(2) + 1; % counted from 0?
+try
+ zzz.tbeg_abs=load('cntpr.dat');
+ zzz.tbeg_abs(2) = zzz.tbeg_abs(2) + 1; % counted from 0?
+catch
+ warning('no cntpr.dat?')
+end
% plots
diff --git a/crpptbx/AUG/torbeam_prepare_inputs_and_run.m b/crpptbx/AUG/torbeam_prepare_inputs_and_run.m
index 50294aa3..e4b55e65 100644
--- a/crpptbx/AUG/torbeam_prepare_inputs_and_run.m
+++ b/crpptbx/AUG/torbeam_prepare_inputs_and_run.m
@@ -117,10 +117,10 @@ end
r0exp = torbeam_in.eqd_all.eqdsk{1}.r0;
if doload; disp(['time for reading data: ' num2str(toc)]); end
-
+
if doload,
nb_points_prof = 46;
- rhofit = linspace(0,1,nb_points_prof);
+ rhofit = linspace(0,1,nb_points_prof)';
tb_par_fields=fieldnames(torbeam_in.inbeam.tb_par);
@@ -141,7 +141,7 @@ if doload,
fprintf(fid,'%d\n',nb_points_prof);
% tefit = interpos(torbeam_in.nete.fit.rhopolnorm(:,it_te) , torbeam_in.nete.fit.te.data(:,it_te),rhofit,-1,[1 2],[0 0]);
tefit = interpos(torbeam_in.nete.fit.rhopolnorm(:,it_te) , torbeam_in.nete.fit.te.data(:,it_te),rhofit,-0.1,[1 0],[0 0]);
- aa=[rhofit(end:-1:1)' , abs(tefit(end:-1:1))/1e3]; % in keV
+ aa=[rhofit(end:-1:1) , abs(tefit(end:-1:1))/1e3]; % in keV
fprintf(fid,'%.7f %.7f\n',aa');
fclose(fid);
% ne
@@ -150,7 +150,7 @@ if doload,
fid=fopen(torbeam_in.fname{it}.ne,'w'); % ne.dat file etc per time
fprintf(fid,'%d\n',nb_points_prof);
nefit = interpos(torbeam_in.nete.fit.rhopolnorm(:,it_ne) , torbeam_in.nete.fit.ne.data(:,it_ne),rhofit,-0.1,[1 0],[0 0]);
- aa=[rhofit(end:-1:1)' , abs(nefit(end:-1:1))/1e19]; % in 1e19
+ aa=[rhofit(end:-1:1) , abs(nefit(end:-1:1))/1e19]; % in 1e19
fprintf(fid,'%.6f %.6f\n',aa');
fclose(fid);
% topfile
@@ -210,7 +210,7 @@ if doload,
end
disp(['time for creating all input files: ' num2str(toc)])
end
-
+
% run TORBEAM
current_dir = pwd;
eval(['cd ' torbeam_in.dir_torbeam])
@@ -228,6 +228,11 @@ torbeam_out.shotfile.shot = shot;
torbeam_out.shotfile.t = torbeam_out.t;
tic
icount = 0;
+ld_library_path=getenv('LD_LIBRARY_PATH'); % need /afs/@cell/common/soft/intel/ics2013/14.0/compiler/lib/intel64 (ok on toks but not toki now)
+ij_intel = findstr(ld_library_path,'intel');
+if isempty(ij_intel)
+ warning('no intel in LD_LIBRARY_PATH, torbeam may not run properly');
+end
for it_eff=1:length(times_for_torbeam_eff)
it = it_inb_beam_on(it_eff);
unix(['cp ' torbeam_in.fname{it}.te ' Te.dat']);
@@ -236,7 +241,7 @@ for it_eff=1:length(times_for_torbeam_eff)
disp([' t=' num2str(torbeam_in.inbeam.t(it_inb(it_eff))) ' beam=' num2str(find(torbeam_in.inbeam.beam_on(:,it_inb(it_eff)))','%d')]);
for ibeam=find(torbeam_in.inbeam.beam_on(:,it_inb(it_eff)))',
unix(['cp ' torbeam_in.fname{it}.inbeam{ibeam} ' inbeam.dat']);
- [a,b]=unix('/afs/ipp-garching.mpg.de/home/o/osauter/TORBEAM/lib-OUT/a.out');
+ [a,b]=unix(['setenv LD_LIBRARY_PATH ' ld_library_path '; /afs/ipp-garching.mpg.de/home/o/osauter/TORBEAM/lib-OUT/a.out']);
icount = icount + 1;
imw = findstr('(MW)',b);
if ~isempty(imw)
--
GitLab