call to gdat 'powers' not working if no ohmic power
The following call
powers = gdat(69866,'powers')
gives the following error
> In gdat (line 196)
One or more output arguments not assigned during call to "interpos".
Error in gdat_tcv (line 2441)
gdat_data.data(:,1) = interpos(-21,gdat_data.ohm.t(ij),gdat_data.ohm.data(ij),gdat_data.t);
Error in gdat (line 190)
[gdat_data,gdat_params,error_status,varargout] = feval(['gdat_' lower(machine_eff)],args{:});
It is given by the fact that
ptot_ohm = tdi('\results::conf:ptot_ohm');
is filled with NaNs
@osauter should we add a try/catch
on gdat
or should we solve the problem at a higher level (i.e. fill whichever node for the ohmic power gdat
is using)?