diff --git a/JET/loadJETdata.m b/JET/loadJETdata.m index 0b02b05e639b0c48ad07af9fb8a1becae1b71bc7..2eae47999370765aada099ac45a90fbcf0742ea4 100644 --- a/JET/loadJETdata.m +++ b/JET/loadJETdata.m @@ -297,20 +297,32 @@ switch JETkeywrdcase{index} tracename=tracename(ij); [a,x,t,d,e]=rda_eff(shot,ppftype,tracename); switch tracename - case 'efit/btpd?uid=jetppf+seq=0' + case {'efit/btpd','efit/btpd?uid=jetppf+seq=0'} if isstr(t); disp('t is a string, assumes means btpd not defined'); - [wdia,x1,t1,d,e]=rda_eff(shot,'ppf','efit/wdia'); - [rgeo,x,t,d,e]=rda_eff(shot,'ppf','efit/rgeo'); - [xip,x3,t3,d,e]=rda_eff(shot,'ppf','efit/xip'); + [xip,x,t,d,e]=rda_eff(shot,'ppf','efit/xip'); + shot_mg3_list=[47274 47275 47276 47280 47281 47282 47283 47284 47285 47286 47287 47290 47295 47296 47301]; + if isempty(find(shot_mg3_list==shot)) + [wdia,x1,t1,d,e]=rda_eff(shot,'ppf','efit/wdia'); + else + [wdia,x1,t1,d,e]=rda_eff(shot,'ppf','mg3/wpd'); + wdia=interp1(t1,wdia,t); + end + [rgeo,x3,t3,d,e]=rda_eff(shot,'ppf','efit/rgeo'); a=2.122e6 .* wdia ./xip.^2 ./ rgeo; end - case 'efit/btnd?uid=jetppf+seq=0' + case {'efit/btnd','efit/btnd?uid=jetppf+seq=0'} if isstr(t); disp('t is a string, assumes means btnd not defined'); - [wdia,x1,t1,d,e]=rda_eff(shot,'ppf','efit/wdia'); - [rgeo,x,t,d,e]=rda_eff(shot,'ppf','efit/rgeo'); - [xip,x3,t3,d,e]=rda_eff(shot,'ppf','efit/xip'); + [xip,x,t,d,e]=rda_eff(shot,'ppf','efit/xip'); + shot_mg3_list=[47274 47275 47276 47280 47281 47282 47283 47284 47285 47286 47287 47290 47295 47296 47301]; + if isempty(find(shot_mg3_list==shot)) + [wdia,x1,t1,d,e]=rda_eff(shot,'ppf','efit/wdia'); + else + [wdia,x1,t1,d,e]=rda_eff(shot,'ppf','mg3/wpd'); + wdia=interp1(t1,wdia,t); + end + [rgeo,x3,t3,d,e]=rda_eff(shot,'ppf','efit/rgeo'); [cr0,x3,t3,d,e]=rda_eff(shot,'ppf','efit/cr0'); [bvac,x3,t3,d,e]=rda_eff(shot,'ppf','efit/bvac'); [volm,x3,t3,d,e]=rda_eff(shot,'ppf','efit/volm');