-
Notifications
You must be signed in to change notification settings - Fork 18
/
Copy pathstd2abj.m
44 lines (35 loc) · 1.08 KB
/
std2abj.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
%% std2abj
% writes table with abj and std values for tau_p, tau_b, l_p, l_b
%%
function std2abj(spec)
% created 2024/11/19 by Bas Kooijman
%% Syntax
% <../std2abj.m *std2abj*> (spec)
%% Description
% writes table with abj and std values for tau_p, tau_b, l_p, l_b and shows it in the bowser
%
% Input:
%
% * spec: entry name
%
% Output:
%
% * text-file with the name title.html is written and shown in browser
%% Example of use
% std2abj('Daphnia_pulex');
par = allStat2par(spec);
if ~isfield(par,'E_Hj')
par.E_Hj = par.E_Hb + 1e-3;
else
fprintf('Warning from std2abj: E_Hj is already in pars\n');
end
cPar = parscomp_st(par); vars_pull(par); vars_pull(cPar);
res = NaN(4,2); vars = {'tau_p';'tau_b';'l_p';'l_b'};
pars_tj = [g k l_T v_Hb v_Hj v_Hp]; % compose pars for get_tj
[tau_j, tau_p, tau_b, l_j, l_p, l_b] = get_tj(pars_tj, 1);
res(:,1) = [tau_p; tau_b; l_p; l_b];
pars_tp = [g k l_T v_Hb v_Hp]; % compose pars for get_tp
[tau_p, tau_b, l_p, l_b] = get_tp(pars_tp, 1);
res(:,2) = [tau_p; tau_b; l_p; l_b];
prt_tab({vars,res},{spec,'abj','std'})
end