diff --git a/lib/misc/C2K.m b/lib/misc/C2K.m
new file mode 100644
index 00000000..79b91e4c
--- /dev/null
+++ b/lib/misc/C2K.m
@@ -0,0 +1,27 @@
+%% C2K
+% computes Kelvin from Celsius
+
+%%
+function K = C2K(C)
+  % created by Bas Kooijman 2015/01/30; modified 2015/07/28
+  
+  %% Syntax
+  % K = <../C2K.m *C2K*>(C)
+  
+  %% Description
+  % Obtains Kelvin from temperatures defined in Celsius
+  %
+  % Input
+  %
+  % * C: scalar or matrix in temperatures in degrees Celsius
+  %  
+  % Output
+  %
+  % * K: temperature(s) in Kelvin
+  
+  %% Example 
+  % C2K(20)
+
+  K = C + 273.15; % K, temperature in Kelvin
+end
+
diff --git a/lib/misc/K2C.m b/lib/misc/K2C.m
new file mode 100644
index 00000000..b6f7fc28
--- /dev/null
+++ b/lib/misc/K2C.m
@@ -0,0 +1,27 @@
+%% K2C
+% computes Celsius from Kelvin
+
+%%
+function C = K2C(K)
+  % created by Bas Kooijman 2015/01/30; modified 2015/07/28
+  
+  %% Syntax
+  % K = <../C2K.m *C2K*>(C)
+  
+  %% Description
+  % Obtains Kelvin from temperatures defined in degrees Celsius
+  %
+  % Input
+  %
+  % * K: scalar or matrix in temperatures in Kelvin
+  %  
+  % Output
+  %
+  % * C: temperature(s) in degrees Celsius
+  
+  %% Example 
+  % K2C(280)
+
+  C = K - 273.15; % K, temperature in Kelvin
+end
+
diff --git a/lib/misc/html/C2K.html b/lib/misc/html/C2K.html
new file mode 100644
index 00000000..092ce532
--- /dev/null
+++ b/lib/misc/html/C2K.html
@@ -0,0 +1,90 @@
+
+<!DOCTYPE html
+  PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
+<html><head>
+      <meta http-equiv="Content-Type" content="text/html; charset=utf-8">
+   <!--
+This HTML was auto-generated from MATLAB code.
+To make changes, update the MATLAB code and republish this document.
+      --><title>C2K</title><meta name="generator" content="MATLAB 7.14"><link rel="schema.DC" href="http://purl.org/dc/elements/1.1/"><meta name="DC.date" content="2015-07-28"><meta name="DC.source" content="C2K.m"><style type="text/css">
+html,body,div,span,applet,object,iframe,h1,h2,h3,h4,h5,h6,p,blockquote,pre,a,abbr,acronym,address,big,cite,code,del,dfn,em,font,img,ins,kbd,q,s,samp,small,strike,strong,sub,sup,tt,var,b,u,i,center,dl,dt,dd,ol,ul,li,fieldset,form,label,legend,table,caption,tbody,tfoot,thead,tr,th,td{margin:0;padding:0;border:0;outline:0;font-size:100%;vertical-align:baseline;background:transparent}body{line-height:1}ol,ul{list-style:none}blockquote,q{quotes:none}blockquote:before,blockquote:after,q:before,q:after{content:'';content:none}:focus{outine:0}ins{text-decoration:none}del{text-decoration:line-through}table{border-collapse:collapse;border-spacing:0}
+
+html { min-height:100%; margin-bottom:1px; }
+html body { height:100%; margin:0px; font-family:Arial, Helvetica, sans-serif; font-size:10px; color:#000; line-height:140%; background:#fff none; overflow-y:scroll; }
+html body td { vertical-align:top; text-align:left; }
+
+h1 { padding:0px; margin:0px 0px 25px; font-family:Arial, Helvetica, sans-serif; font-size:1.5em; color:#d55000; line-height:100%; font-weight:normal; }
+h2 { padding:0px; margin:0px 0px 8px; font-family:Arial, Helvetica, sans-serif; font-size:1.2em; color:#000; font-weight:bold; line-height:140%; border-bottom:1px solid #d6d4d4; display:block; }
+h3 { padding:0px; margin:0px 0px 5px; font-family:Arial, Helvetica, sans-serif; font-size:1.1em; color:#000; font-weight:bold; line-height:140%; }
+
+a { color:#005fce; text-decoration:none; }
+a:hover { color:#005fce; text-decoration:underline; }
+a:visited { color:#004aa0; text-decoration:none; }
+
+p { padding:0px; margin:0px 0px 20px; }
+img { padding:0px; margin:0px 0px 20px; border:none; }
+p img, pre img, tt img, li img { margin-bottom:0px; } 
+
+ul { padding:0px; margin:0px 0px 20px 23px; list-style:square; }
+ul li { padding:0px; margin:0px 0px 7px 0px; }
+ul li ul { padding:5px 0px 0px; margin:0px 0px 7px 23px; }
+ul li ol li { list-style:decimal; }
+ol { padding:0px; margin:0px 0px 20px 0px; list-style:decimal; }
+ol li { padding:0px; margin:0px 0px 7px 23px; list-style-type:decimal; }
+ol li ol { padding:5px 0px 0px; margin:0px 0px 7px 0px; }
+ol li ol li { list-style-type:lower-alpha; }
+ol li ul { padding-top:7px; }
+ol li ul li { list-style:square; }
+
+.content { font-size:1.2em; line-height:140%; padding: 20px; }
+
+pre, tt, code { font-size:12px; }
+pre { margin:0px 0px 20px; }
+pre.error { color:red; }
+pre.codeinput { padding:10px; border:1px solid #d3d3d3; background:#f7f7f7; }
+pre.codeoutput { padding:10px 11px; margin:0px 0px 20px; color:#4c4c4c; }
+
+@media print { pre.codeinput, pre.codeoutput { word-wrap:break-word; width:100%; } }
+
+span.keyword { color:#0000FF }
+span.comment { color:#228B22 }
+span.string { color:#A020F0 }
+span.untermstring { color:#B20000 }
+span.syscmd { color:#B28C00 }
+
+.footer { width:auto; padding:10px 0px; margin:25px 0px 0px; border-top:1px dotted #878787; font-size:0.8em; line-height:140%; font-style:italic; color:#878787; text-align:left; float:none; }
+.footer p { margin:0px; }
+
+  </style></head><body><div class="content"><h1>C2K</h1><!--introduction--><p>computes Kelvin from Celsius</p><!--/introduction--><h2>Contents</h2><div><ul><li><a href="#3">Syntax</a></li><li><a href="#4">Description</a></li><li><a href="#5">Example</a></li></ul></div><h2>Syntax<a name="3"></a></h2><p>K = <a href="../C2K.m"><b>C2K</b></a>(C)</p><h2>Description<a name="4"></a></h2><p>Obtains Kelvin from temperatures defined in Celsius</p><p>Input</p><div><ul><li>C: scalar or matrix in temperatures in degrees Celsius</li></ul></div><p>Output</p><div><ul><li>K: temperature(s) in Kelvin</li></ul></div><h2>Example<a name="5"></a></h2><p>C2K(20)</p><p class="footer"><br>
+      Published with MATLAB&reg; 7.14<br></p></div><!--
+##### SOURCE BEGIN #####
+%% C2K
+% computes Kelvin from Celsius
+
+%%
+function K = C2K(C)
+  % created by Bas Kooijman 2015/01/30; modified 2015/07/28
+  
+  %% Syntax
+  % K = <../C2K.m *C2K*>(C)
+  
+  %% Description
+  % Obtains Kelvin from temperatures defined in Celsius
+  %
+  % Input
+  %
+  % * C: scalar or matrix in temperatures in degrees Celsius
+  %  
+  % Output
+  %
+  % * K: temperature(s) in Kelvin
+  
+  %% Example 
+  % C2K(20)
+
+  K = C + 273.15; % K, temperature in Kelvin
+end
+
+
+##### SOURCE END #####
+--></body></html>
\ No newline at end of file
diff --git a/lib/pet/html/check_my_pet.html b/lib/pet/html/check_my_pet.html
index 18efe617..7000ee2f 100644
--- a/lib/pet/html/check_my_pet.html
+++ b/lib/pet/html/check_my_pet.html
@@ -6,7 +6,7 @@
    <!--
 This HTML was auto-generated from MATLAB code.
 To make changes, update the MATLAB code and republish this document.
-      --><title>check_my_pet</title><meta name="generator" content="MATLAB 7.14"><link rel="schema.DC" href="http://purl.org/dc/elements/1.1/"><meta name="DC.date" content="2015-07-18"><meta name="DC.source" content="check_my_pet.m"><style type="text/css">
+      --><title>check_my_pet</title><meta name="generator" content="MATLAB 7.14"><link rel="schema.DC" href="http://purl.org/dc/elements/1.1/"><meta name="DC.date" content="2015-07-28"><meta name="DC.source" content="check_my_pet.m"><style type="text/css">
 html,body,div,span,applet,object,iframe,h1,h2,h3,h4,h5,h6,p,blockquote,pre,a,abbr,acronym,address,big,cite,code,del,dfn,em,font,img,ins,kbd,q,s,samp,small,strike,strong,sub,sup,tt,var,b,u,i,center,dl,dt,dd,ol,ul,li,fieldset,form,label,legend,table,caption,tbody,tfoot,thead,tr,th,td{margin:0;padding:0;border:0;outline:0;font-size:100%;vertical-align:baseline;background:transparent}body{line-height:1}ol,ul{list-style:none}blockquote,q{quotes:none}blockquote:before,blockquote:after,q:before,q:after{content:'';content:none}:focus{outine:0}ins{text-decoration:none}del{text-decoration:line-through}table{border-collapse:collapse;border-spacing:0}
 
 html { min-height:100%; margin-bottom:1px; }
@@ -55,21 +55,7 @@
 .footer { width:auto; padding:10px 0px; margin:25px 0px 0px; border-top:1px dotted #878787; font-size:0.8em; line-height:140%; font-style:italic; color:#878787; text-align:left; float:none; }
 .footer p { margin:0px; }
 
-  </style></head><body><div class="content"><h1>check_my_pet</h1><!--introduction--><p>Checks my_data, pars_init and predict files for field inconsistencies</p><!--/introduction--><h2>Contents</h2><div><ul><li><a href="#3">Syntax</a></li><li><a href="#4">Description</a></li><li><a href="#5">Example of use</a></li><li><a href="#8">Checking the my_data file</a></li><li><a href="#9">Checking the pars_init file</a></li><li><a href="#10">Checking the predict file</a></li></ul></div><h2>Syntax<a name="3"></a></h2><p><a href="../check_my_pet.m"><b>check_my_pet</b></a> (speciesnm)</p><h2>Description<a name="4"></a></h2><p>Checks my_data, pars_init and predict files for field inconsistencies. Checking points for my_data:</p><pre class="language-matlab">* existence of <span class="string">standard</span> <span class="string">metadata</span> <span class="string">fields</span>
-* existence of <span class="string">temp</span>
-* existence of <span class="string">pseudodata</span>
-* existence of <span class="string">weights</span>
-* existence of <span class="string">units</span> <span class="string">(and number consistence with type of data)</span>
-* existence of <span class="string">labels</span> <span class="string">(and number consistence with type of data)</span>
-* existence of <span class="string">bibkeys</span> <span class="string">and</span> <span class="string">the</span> <span class="string">corresponding</span> <span class="string">entrances</span> <span class="string">in</span> <span class="string">biblist</span>
-</pre><p>Checking points for pars_init:</p><pre class="language-matlab">* existence of <span class="string">standard</span> <span class="string">metapar</span> <span class="string">fields</span>
-* model is <span class="string">one</span> <span class="string">of</span> <span class="string">the</span> <span class="string">predefined</span> <span class="string">models</span>
-* existence of <span class="string">free</span>
-* existence of <span class="string">units</span>
-* existence of <span class="string">labels</span>
-</pre><p>Checking points for predict:</p><pre class="language-matlab">* existence of <span class="string">the</span> <span class="string">same</span> <span class="string">fields</span> <span class="string">as</span> <span class="string">in</span> <span class="string">data</span>
-* length of <span class="string">prediction</span> <span class="string">results</span> <span class="string">matches</span> <span class="string">the</span> <span class="string">length</span> <span class="string">of</span> <span class="string">data</span>
-</pre><p>Input:</p><div><ul><li>string with species name or cell vector with multiple species names</li></ul></div><p>Output:</p><div><ul><li>warnings are printed to screen</li></ul></div><h2>Example of use<a name="5"></a></h2><p>See <a href="../run_my_pet.m"><b>run_my_pet</b></a></p><h2>Checking the my_data file<a name="8"></a></h2><h2>Checking the pars_init file<a name="9"></a></h2><h2>Checking the predict file<a name="10"></a></h2><p class="footer"><br>
+  </style></head><body><div class="content"><h1>check_my_pet</h1><!--introduction--><p>Checks my_data, pars_init and predict files for field inconsistencies</p><!--/introduction--><h2>Contents</h2><div><ul><li><a href="#3">Syntax</a></li><li><a href="#4">Description</a></li><li><a href="#5">Example of use</a></li><li><a href="#8">%%%%%%%%%%%%%%%%%%%%%%%%%%%</a></li><li><a href="#10">Checking the my_data file</a></li><li><a href="#12">%%%%%%%%%%%%%%%%%%%%%%%%%%%</a></li><li><a href="#13">checking the existence of metadata fields</a></li><li><a href="#14">checking if metadata.species matches speciesnm</a></li><li><a href="#15">checking the existence of temp in the data structure</a></li><li><a href="#16">checking the existence of psd in the data structure</a></li><li><a href="#17">checking the existence of weight in the data structure</a></li><li><a href="#18">checking the existence of units in the txt_data structure</a></li><li><a href="#19">checking the existence of labels in the txt_data structure</a></li><li><a href="#20">checking the existence of facts</a></li><li><a href="#21">checking the existence of bibkeys in the txt_data structure</a></li><li><a href="#22">checking the existence of bibkeys in the biblist structure</a></li><li><a href="#23">%%%%%%%%%%%%%%%%%%%%%%%%%%%</a></li><li><a href="#25">Checking the pars_init file</a></li><li><a href="#27">%%%%%%%%%%%%%%%%%%%%%%%%%%%</a></li><li><a href="#28">checking the existence of metapar fields</a></li><li><a href="#29">checking the existence of metapar fields</a></li><li><a href="#30">checking the existence of par fields</a></li><li><a href="#31">checking the existence of free in the par structure and if it filled with either 0 or 1</a></li><li><a href="#32">checking the existence of units in the txt_par structure</a></li><li><a href="#33">checking the existence of labels in the txt_par structure</a></li><li><a href="#34">checking the realism of the par set</a></li><li><a href="#35">%%%%%%%%%%%%%%%%%%%%%%%%%%%</a></li><li><a href="#37">Checking the predict file</a></li><li><a href="#39">%%%%%%%%%%%%%%%%%%%%%%%%%%%</a></li><li><a href="#40">checking the existence of psd in the Prd_data structure</a></li><li><a href="#41">checking for the bijection data-predictions</a></li><li><a href="#42">checking for the pseudodate predictions in data</a></li></ul></div><h2>Syntax<a name="3"></a></h2><p><a href="../check_my_pet.m"><b>check_my_pet</b></a> (speciesnm)</p><h2>Description<a name="4"></a></h2><p>Checks my_data, pars_init and predict files for field inconsistencies Checking points for my_data:   - existence of standard metadata fields   - existence of temp   - existence of pseudodata   - existence of weights   - existence of units (and number consistence with type of data)   - existence of labels (and number consistence with type of data)   - existence of bibkeys               and the corresponding entrances in biblist Checking points for pars_init:   - existence of standard metapar fields   - model is one of the predefined models   - existence of free   - existence of units   - existence of labels Checking points for predict:   - existence of the same fields as in data   - length of prediction results matches the length of data</p><p>Input</p><div><ul><li>speciesnms: string with species name or cell vector with multiple species names</li></ul></div><h2>Example of use<a name="5"></a></h2><p>check_files('my_pet')</p><h2>%%%%%%%%%%%%%%%%%%%%%%%%%%%<a name="8"></a></h2><h2>Checking the my_data file<a name="10"></a></h2><h2>%%%%%%%%%%%%%%%%%%%%%%%%%%%<a name="12"></a></h2><h2>checking the existence of metadata fields<a name="13"></a></h2><h2>checking if metadata.species matches speciesnm<a name="14"></a></h2><h2>checking the existence of temp in the data structure<a name="15"></a></h2><h2>checking the existence of psd in the data structure<a name="16"></a></h2><h2>checking the existence of weight in the data structure<a name="17"></a></h2><h2>checking the existence of units in the txt_data structure<a name="18"></a></h2><h2>checking the existence of labels in the txt_data structure<a name="19"></a></h2><h2>checking the existence of facts<a name="20"></a></h2><h2>checking the existence of bibkeys in the txt_data structure<a name="21"></a></h2><h2>checking the existence of bibkeys in the biblist structure<a name="22"></a></h2><h2>%%%%%%%%%%%%%%%%%%%%%%%%%%%<a name="23"></a></h2><h2>Checking the pars_init file<a name="25"></a></h2><h2>%%%%%%%%%%%%%%%%%%%%%%%%%%%<a name="27"></a></h2><h2>checking the existence of metapar fields<a name="28"></a></h2><h2>checking the existence of metapar fields<a name="29"></a></h2><h2>checking the existence of par fields<a name="30"></a></h2><h2>checking the existence of free in the par structure and if it filled with either 0 or 1<a name="31"></a></h2><h2>checking the existence of units in the txt_par structure<a name="32"></a></h2><h2>checking the existence of labels in the txt_par structure<a name="33"></a></h2><h2>checking the realism of the par set<a name="34"></a></h2><h2>%%%%%%%%%%%%%%%%%%%%%%%%%%%<a name="35"></a></h2><h2>Checking the predict file<a name="37"></a></h2><h2>%%%%%%%%%%%%%%%%%%%%%%%%%%%<a name="39"></a></h2><h2>checking the existence of psd in the Prd_data structure<a name="40"></a></h2><h2>checking for the bijection data-predictions<a name="41"></a></h2><h2>checking for the pseudodate predictions in data<a name="42"></a></h2><p class="footer"><br>
       Published with MATLAB&reg; 7.14<br></p></div><!--
 ##### SOURCE BEGIN #####
 %% check_my_pet
@@ -77,46 +63,39 @@
 
 %%
 function check_my_pet(speciesnms)
-  % created 2015/05/05 by Goncalo Marques; modified 2015/05/12, 2015/06/05, 2015/07/02 and 2015/07/08 by Goncalo Marques
+  % created 2015/05/05 by Goncalo Marques; modified 2015/05/12, 2015/06/05, 2015/07/02, 2015/07/08 and 2015/07/22 by Goncalo Marques
   
   %% Syntax 
   % <../check_my_pet.m *check_my_pet*> (speciesnm)
 
   %% Description
-  % Checks my_data, pars_init and predict files for field inconsistencies.
+  % Checks my_data, pars_init and predict files for field inconsistencies
   % Checking points for my_data:
-  %
-  %   * existence of standard metadata fields
-  %   * existence of temp
-  %   * existence of pseudodata
-  %   * existence of weights
-  %   * existence of units (and number consistence with type of data)
-  %   * existence of labels (and number consistence with type of data)
-  %   * existence of bibkeys and the corresponding entrances in biblist
-  %
+  %   - existence of standard metadata fields
+  %   - existence of temp
+  %   - existence of pseudodata
+  %   - existence of weights
+  %   - existence of units (and number consistence with type of data)
+  %   - existence of labels (and number consistence with type of data)
+  %   - existence of bibkeys
+  %               and the corresponding entrances in biblist
   % Checking points for pars_init:
-  %
-  %   * existence of standard metapar fields
-  %   * model is one of the predefined models
-  %   * existence of free 
-  %   * existence of units 
-  %   * existence of labels 
-  %
+  %   - existence of standard metapar fields
+  %   - model is one of the predefined models
+  %   - existence of free 
+  %   - existence of units 
+  %   - existence of labels 
   % Checking points for predict:
+  %   - existence of the same fields as in data
+  %   - length of prediction results matches the length of data
   %
-  %   * existence of the same fields as in data
-  %   * length of prediction results matches the length of data
-  %
-  % Input:
+  % Input
   %
-  % * string with species name or cell vector with multiple species names
+  % * speciesnms: string with species name or cell vector with multiple species names
   %  
-  % Output:
-  %
-  % * warnings are printed to screen
 
   %% Example of use
-  % See <../run_my_pet.m *run_my_pet*> 
+  % check_files('my_pet') 
 
 if iscell(speciesnms)
   k = length(speciesnms);
@@ -133,15 +112,21 @@
 info = 0;
 
 if ~isempty(strfind(speciesnm, ' '))
-  fprintf(['The species name in input should not have spaces.\n']);
-  fprintf(['The standard species name follow the form ''Genus_species''.\n']);
+  fprintf('The species name in input should not have spaces.\n');
+  fprintf('The standard species name follow the form ''Genus_species''.\n');
   return;
 end
   
 
+%% %%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%
 %% Checking the my_data file
- 
-if exist(['mydata_', speciesnm], 'var')~=2
+%%
+%% %%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+bibNotToCheck = {'Kooy2010'};    % list of bib references that does not need to have an explicit reference
+
+if exist(['mydata_', speciesnm], 'file')~=2
   fprintf(['There is no mydata_', speciesnm,' file.\n']);
   return;
 end
@@ -151,7 +136,7 @@
 datafields = fields(data);
 txtdatafields = fields(txt_data);
 
-% checking the existence of metadata fields
+%% checking the existence of metadata fields
 mtdtfields = {'phylum', 'class', 'order', 'family', 'species', 'species_en', 'T_typical', 'data_0', 'data_1'};
 
 for i = 1:length(mtdtfields)
@@ -160,12 +145,12 @@
   end
 end
 
-% checking if metadata.species matches speciesnm
+%% checking if metadata.species matches speciesnm
 if ~strcmp(speciesnm, metadata.species)
-  fprintf(['The species name in metadata.species does not match the species name in the mydata file name.\n']);
+  fprintf('The species name in metadata.species does not match the species name in the mydata file name.\n');
 end
 
-% checking the existence of temp in the data structure
+%% checking the existence of temp in the data structure
 if sum(strcmp(datafields, 'temp')) == 0
   fprintf('The data structure does not include temperature data. \n');
 else
@@ -179,7 +164,7 @@
   end
 end
 
-% checking the existence of psd in the data structure
+%% checking the existence of psd in the data structure
 if sum(strcmp(datafields, 'psd')) == 0
   fprintf('The data structure does not include the pseudodata for the regression. \n');
   psdexist = 0;
@@ -189,7 +174,7 @@
   psdexist = 1;
 end
 
-% checking the existence of weight in the data structure
+%% checking the existence of weight in the data structure
 if sum(strcmp(datafields, 'weight')) == 0
   fprintf('The data structure does not include the weights for the regression. \n');
   info = 1;
@@ -232,7 +217,7 @@
   end
 end
 
-% checking the existence of units in the txt_data structure
+%% checking the existence of units in the txt_data structure
 if sum(strcmp(txtdatafields, 'units')) == 0
   fprintf('The txt_data structure does not include the data units. \n');
 else
@@ -281,7 +266,7 @@
   end
 end
   
-% checking the existence of labels in the txt_data structure
+%% checking the existence of labels in the txt_data structure
 if sum(strcmp(txtdatafields, 'label')) == 0
   fprintf('The txt_data structure does not include the data labels. \n');
 else
@@ -330,40 +315,50 @@
   end
 end
 
-% checking the existence of bibkeys in the txt_data structure
+%% checking the existence of facts
+if isfield(metadata, 'facts')
+  factsfields = fields(metadata.facts);
+  factsfields = factsfields(~strcmp(factsfields, 'bibkey'));
+else
+  factsfields = {};
+end
+
+%% checking the existence of bibkeys in the txt_data structure
 if sum(strcmp(txtdatafields, 'bibkey')) == 0
   fprintf('The txt_data structure does not include the bibkeys. \n');
 else
   bibkeyfields = fields(txt_data.bibkey);
+  referencedfields = [datafields; factsfields];
 
-  if length(bibkeyfields) > length(datafields)
+  if length(bibkeyfields) > length(referencedfields)
     for i = 1:length(bibkeyfields)
-      if sum(strcmp(datafields, bibkeyfields(i))) == 0
-        fprintf(['There is a bibkey defined for ', bibkeyfields{i}, ' but there is no corresponding data. \n']);
+      if sum(strcmp(referencedfields, bibkeyfields(i))) == 0
+        fprintf(['There is a bibkey defined for ', bibkeyfields{i}, ' but there is no corresponding data or fact. \n']);
       end
     end
   else
-    for i = 1:length(datafields)
-      if sum(strcmp(datafields(i), bibkeyfields)) == 0
-        fprintf(['There is no bibkey defined for data point/set ', datafields{i}, '. \n']);
+    for i = 1:length(referencedfields)
+      if sum(strcmp(referencedfields(i), bibkeyfields)) == 0
+        fprintf(['There is no bibkey defined for data point/set or fact ', referencedfields{i}, '. \n']);
       end
     end
   end
 end
 
-% checking the existence of bibkeys in the biblist structure
+%% checking the existence of bibkeys in the biblist structure
 if sum(strcmp(fields(metadata), 'biblist')) == 0
   fprintf('The metadata structure does not include the biblist. \n');
 else
   biblistfields = fields(metadata.biblist);
+  biblistfields = biblistfields(~strcmp(biblistfields, bibNotToCheck));    % bibNotToCheck defined at the beginning of mydata section
   for i = 1:length(bibkeyfields)
-    eval(['bibkeynm = txt_data.bibkey.', bibkeyfields{i}, ';']);
+    bibkeynm = eval(['txt_data.bibkey.', bibkeyfields{i}]);
     if ~iscell(bibkeynm) 
       if sum(strcmp(biblistfields, cellstr(bibkeynm))) == 0
         fprintf(['The bibkey ', bibkeynm, ' defined for ', bibkeyfields{i}, ' has no corresponding reference in biblist. \n']);
       end
     else
-      [lin k] = size(bibkeynm);
+      [lin, k] = size(bibkeynm);
       if lin ~= 1 % check if for multiple references they are in a line vector
         fprintf(['The bibkey defined for ', bibkeyfields{i}, ' should be a line vector - use commas (,) to separate multiple references. \n']);
       else
@@ -398,9 +393,13 @@
 end
 
 
+%% %%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%
 %% Checking the pars_init file
+%%
+%% %%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-if exist(['pars_init_', speciesnm], 'var')~=2
+if exist(['pars_init_', speciesnm], 'file')~=2
   fprintf(['There is no pars_init_', speciesnm,' file.\n']);
   return;
 end
@@ -410,7 +409,7 @@
 parfields = fields(par);
 txtparfields = fields(txt_par);
 
-% checking the existence of metapar fields
+%% checking the existence of metapar fields
 mtparfields = {'model', 'T_ref'};
 
 for i = 1:length(mtparfields)
@@ -419,15 +418,26 @@
   end
 end
 
-% checking the existence of metapar fields
-possible_models = {'std', 'stf', 'stx', 'abj', 'hex'};
+%% checking the existence of metapar fields
+possible_models = {'std', 'stf', 'stx', 'ssj', 'abj', 'asj', 'hex'};
 
 if sum(strcmp(metapar.model, possible_models)) == 0
   fprintf(['The model ', metapar.model, ' is not one of the predefined models. \n']);
 end
 
-% checking the existence of par fields
-Eparfields = {'z', 'F_m', 'kap_X', 'kap_P', 'v', 'kap', 'kap_R', 'p_M', 'p_T', 'k_J', 'E_G', 'E_Hb', 'E_Hp', 'h_a', 's_G'};
+%% checking the existence of par fields
+switch metapar.model
+  case {'std', 'stf'}
+    Eparfields = {'z', 'F_m', 'kap_X', 'kap_P', 'v', 'kap', 'kap_R', 'p_M', 'p_T', 'k_J', 'E_G', 'E_Hb', 'E_Hp', 'h_a', 's_G'};
+  case 'stx'
+    Eparfields = {'z', 'F_m', 'kap_X', 'kap_P', 'v', 'kap', 'kap_R', 'p_M', 'p_T', 'k_J', 'E_G', 'E_Hb', 'E_Hx', 'E_Hp', 'h_a', 's_G', 'a_0'};
+  case 'abj'
+    Eparfields = {'z', 'F_m', 'kap_X', 'kap_P', 'v', 'kap', 'kap_R', 'p_M', 'p_T', 'k_J', 'E_G', 'E_Hb', 'E_Hj', 'E_Hp', 'h_a', 's_G'};
+  case 'asj'
+    Eparfields = {'z', 'F_m', 'kap_X', 'kap_P', 'v', 'kap', 'kap_R', 'p_M', 'p_T', 'k_J', 'E_G', 'E_Hb', 'E_Hs', 'E_Hj', 'E_Hp', 'h_a', 's_G'};
+  case 'hex'
+    Eparfields = {'z', 'F_m', 'kap_X', 'kap_P', 'v', 'kap', 'kap_R', 'p_M', 'p_T', 'k_J', 'E_G', 'E_Hb', 's_s', 'E_He', 'h_a', 's_G'};
+end
 
 for i = 1:length(Eparfields)
   if ~isfield(par, Eparfields{i})
@@ -435,19 +445,7 @@
   end
 end
 
-if strcmp(metapar.model, 'abj')
-  if ~isfield(par, 'E_Hj')
-    fprintf('The parameter E_Hj is missing in the par structure. \n');
-  end
-end
-
-if strcmp(metapar.model, 'stx')
-  if ~isfield(par, 'E_Hx')
-    fprintf('The parameter E_Hx is missing in the par structure. \n');
-  end
-end
-
-% checking the existence of free in the par structure
+%% checking the existence of free in the par structure and if it filled with either 0 or 1
 if sum(strcmp(parfields, 'free')) == 0
   fprintf('The par structure does not include the free substructure. \n');
 else
@@ -457,19 +455,29 @@
   if length(freefields) > length(parfields)
     for i = 1:length(freefields)
       if sum(strcmp(parfields, freefields(i))) == 0
-        fprintf(['There are units defined for ', freefields{i}, ' but there is no corresponding data. \n']);
+        fprintf(['There is free defined for ', freefields{i}, ' but there is no corresponding data. \n']);
+      else
+        freeval = eval(['par.free.', freefields{i}]);
+        if freeval ~= 0 && freeval ~= 1
+          fprintf(['The value in free for ', freefields{i}, ' should be either 0 or 1. \n']);
+        end
       end
     end
   else
     for i = 1:length(parfields)
       if sum(strcmp(parfields(i), freefields)) == 0
-        fprintf(['There are no units defined for data point/set ', parfields{i}, '. \n']);
+        fprintf(['There is no free defined for data point/set ', parfields{i}, '. \n']);
+      else
+        freeval = eval(['par.free.', parfields{i}]);
+        if freeval ~= 0 && freeval ~= 1
+          fprintf(['The value in free for ', parfields{i}, ' should be either 0 or 1. \n']);
+        end
       end
     end
   end
 end
 
-% checking the existence of units in the txt_par structure
+%% checking the existence of units in the txt_par structure
 if sum(strcmp(txtparfields, 'units')) == 0
   fprintf('The txt_par structure does not include the parameter units. \n');
 else
@@ -490,7 +498,7 @@
   end
 end
 
-% checking the existence of labels in the txt_par structure
+%% checking the existence of labels in the txt_par structure
 if sum(strcmp(txtparfields, 'label')) == 0
   fprintf('The txt_par structure does not include the parameter units. \n');
 else
@@ -512,7 +520,7 @@
 end
 
 
-% checking the realism of the par set
+%% checking the realism of the par set
 filternm = ['filter_', metapar.model];
 [pass, flag]  = feval(filternm, par, chem);
 if ~pass 
@@ -527,9 +535,13 @@
 end
 
 
+%% %%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%
 %% Checking the predict file
+%%
+%% %%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-if exist(['predict_', speciesnm], 'var')~=2
+if exist(['predict_', speciesnm], 'file')~=2
   fprintf(['There is no predict_', speciesnm,' file.\n']);
   return;
 end
@@ -543,6 +555,16 @@
 
 Prd_datafields = fields(Prd_data);
 
+%% checking the existence of psd in the Prd_data structure
+if sum(strcmp(Prd_datafields, 'psd')) == 0
+  Prd_psdexist = 0;
+else
+  Prd_datafields = Prd_datafields(~strcmp(Prd_datafields, 'psd'));
+  Prd_psdfields = fields(Prd_data.psd);
+  Prd_psdexist = 1;
+end
+
+%% checking for the bijection data-predictions
 if length(Prd_datafields) > length(datafields)
   for i = 1:length(Prd_datafields)
     if sum(strcmp(datafields, Prd_datafields(i))) == 0
@@ -569,6 +591,20 @@
   end
 end
 
+%% checking for the pseudodate predictions in data
+if Prd_psdexist 
+  for i = 1:length(Prd_psdfields)
+    if sum(strcmp(psdfields, Prd_psdfields(i))) == 0
+      fprintf(['There is a prediction defined for psd.', Prd_psdfields{i}, ' but there is no corresponding data. \n']);
+    else 
+      eval(['[ldt, k] = size(data.psd.', Prd_psdfields{i}, ');']);  % number of data points per set
+      eval(['[lprdt, k] = size(data.psd.', Prd_psdfields{i}, ');']);  % number of data points per set
+      if ldt ~= lprdt
+        fprintf(['There is a prediction defined for psd.', Prd_psdfields{i}, ' has a length of ', num2str(lprdt), ' but the corresponding data has a length of ', num2str(ldt), '. \n']);
+      end
+    end
+  end
+end
 
 
 ##### SOURCE END #####
diff --git a/lib/pet/mydata_my_pet.asv b/lib/pet/mydata_my_pet.asv
deleted file mode 100644
index 13ca7a9f..00000000
--- a/lib/pet/mydata_my_pet.asv
+++ /dev/null
@@ -1,212 +0,0 @@
-%% mydata_my_pet
-% Sets referenced data
-
-%%
-function [data, txt_data, metadata] = mydata_my_pet 
-  % created by Starrlight Augustine, Bas Kooijman, Dina Lika, Goncalo Marques and Laure Pecquerie 2015/03/31
-  % last modified: 07/07/2015 by starrlight
-  
-  %% Syntax
-  % [data, txt_data, metadata] = <../mydata_my_pet.m *mydata_my_pet*>
-  
-  %% Description
-  % Sets data, pseudodata, metadata, explanatory text, weight coefficients.
-  % Meant to be a template in add-my-pet
-  %
-  % Output
-  %
-  % * data: structure with data
-  % * txt_data: text vector for the presentation of results
-  % * metadata: structure with info about this entry
-  
-  %% To do (remove these remarks after editing this file)
-  % * copy this template; replace 'my_pet' by the name of your species
-  % * fill in metadata fields with the proper information
-  % * insert references for the data (an example is given), for multiple references, please use commas to separate references
-  % * edit fact-list for your species, where you can add species and/or data properties
-  % * edit real data; remove all data that to not belong to your pet
-  % * complete reference list
-  % * OPTIONAL : add discussion points / comments before the reference list
-
-%% set metadata
-
-T_C = 273.15; % K, temperature at 0 degrees C (used in T_typical)
-
-metadata.phylum     = 'phylum'; 
-metadata.class      = 'class'; 
-metadata.order      = 'order'; 
-metadata.family     = 'family';
-metadata.species    = 'my_pet'; 
-metadata.species_en = 'my_pet_english_name'; 
-metadata.T_typical  = T_C + 20; % K
-metadata.data_0     = {'ab'; 'ap'; 'am'; 'Lb'; 'Lp'; 'Li'; 'Wdb'; 'Wdp'; 'Wdi'; 'Ri'};  % tags for different types of zero-variate data
-metadata.data_1     = {'t-L'; 'L-W'}; % tags for different types of uni-variate data
-
-metadata.COMPLETE = 2.5; % using criteria of LikaKear2011
-
-metadata.author   = {'FirstName1 LastName1'};            % put names as authors as separate strings:  {'FirstName1 LastName2','FirstName2 LastName2'} , with corresponding author in first place 
-metadata.date_subm = [2015 04 20];                       % [year month day], date at which the entry is submitted
-metadata.email    = {'myname@myuniv.univ'};              % e-mail of corresponding author
-metadata.address  = {'affiliation, zipcode, country'};   % affiliation, postcode, country of the corresponding author
-
-% uncomment and fill in the following fields when the entry is updated:
-% metadata.author_mod_1  = {'FirstName3 LastName3'};          % put names as authors as separate strings:  {'author1','author2'} , with corresponding author in first place 
-% metadata.date_mod_1    = [2017 09 18];                      % [year month day], date modified entry is accepted into the collection
-% metadata.email_mod_1   = {'myname@myuniv.univ'};            % e-mail of corresponding author
-% metadata.address_mod_1 = {'affiliation, zipcode, country'}; % affiliation, postcode, country of the corresponding author
-
-% for curators only ------------------------------
-% metadata.curator     = {'FirstName LastName'};
-% metadata.email_cur   = {'myname@myuniv.univ'}; 
-% metadata.date_acc    = [2015 04 22]; 
-%-------------------------------------------------
-
-%% set data
-% zero-variate data;
-% typically depend on scaled functional response f.
-% here assumed to be equal for all real data; the value of f is specified in pars_init_my_pet.
-
-% age 0 is at onset of embryo development
-data.ab = 15;      units.ab = 'd';    label.ab = 'age at birth';  bibkey.ab = 'MollCano2010';   comment.ab  = 'mean value taken from several measurements'; 
-  temp.ab = T_C + 20;  % K, temperature 
-  % observed age at birth is frequently larger than ab, because of diapauzes during incubation
-data.ap = 261;     units.ap = 'd';    label.ap = 'age at puberty'; bibkey.ap = 'Anon2015';
-  temp.ap = T_C + 20;  % K, temperature 
-  % observed age at puberty is frequently larger than ap, 
-  %   because allocation to reproduction starts before first eggs appear
-data.am = 591;     units.am = 'd';    label.am = 'life span';     bibkey.am = 'Wiki';   
-  temp.am = T_C + 20;  % K, temperature 
-% (accounting for aging only) 
-
-% Please specify what type of length measurement is used for your species.
-% We put here snout-to-vent length, but you should change this depending on your species:
-% add an optional comment structure to give any additional explanations on
-% how the value was chosen, see the last column of the Lb data set for an
-% example
-data.Lb  = 0.45;   units.Lb  = 'cm';   label.Lb  = 'snout to vent length at birth';    bibkey.Lb  = 'Anon2015'; 
-data.Lp  = 2.36;   units.Lp  = 'cm';   label.Lp  = 'snout to vent length at puberty';  bibkey.Lp  = {'Anon2015','Wiki'}; % for multiple references, please use commas to separate references
-data.Li  = 6.25;   units.Li  = 'cm';   label.Li  = 'ultimate snout to vent length';    bibkey.Li  = 'Wiki';
-data.Wdb = 5.8e-5; units.Wdb = 'g';    label.Wdb = 'dry weight at birth';              bibkey.Wdb = 'Anon2015';
-data.Wdp = 8e-3;   units.Wdp = 'g';    label.Wdp = 'dry weight at puberty';            bibkey.Wdp = 'Anon2015';
-data.Wdi = 0.15;   units.Wdi = 'g';    label.Wdi = 'ultimate dry weight';              bibkey.Wdi = 'Wiki';
-data.Ri  = 2.3;    units.Ri  = '#/d';  label.Ri  = 'maximum reprod rate';              bibkey.Ri  = 'Wiki';   
-  % for an individual of ultimate length Li 
-  temp.Ri = T_C + 20;  % K, temperature 
- 
-% uni-variate data
-
-% uni-variate data at f = 0.8 (this value should be added in pars_init_my_pet as a parameter f_tL) 
-% snout-to-vent length and wet weight were measured at the same time
-data.tL = [0     50  100 200 300 400 500 600;    % d, time since birth
-           0.45  1.1 1.7 2.7 3.4 4.0 4.5 4.9]';  % cm, snout-to-vent length at f and T
-units.tL = {'d', 'cm'};     label.tL = {'time since birth', 'snout to vent length'};  bibkey.tL = 'Anon2015';
-  temp.tL = T_C + 25;  % K, temperature
-comment.tL = 'Put here any remarks about the experimental protocol'; % optional field
-
-  
-data.LW = [0.9 1.8 3.2 4.3 5.0;      % cm, snout-to-vent length at f
-           0.004 0.03 0.20 0.55 1.1]';   % g, wet weight at f and T
-units.LW = {'cm', 'g'};     label.LW = {'snout to vent length', 'wet weight'};  bibkey.LW = 'Anon2015';
-comment.LW = 'Put here any remarks about the experimental protocol'; % optional field
-
-%% set weights for all real data
-weight = setweights(data, []);
-
-%% overwriting weights (remove these remarks after editing the file)
-% the weights were set automatically with the function setweigths,
-% if one wants to ovewrite one of the weights it should always present an explanation example:
-%
-% zero-variate data:
-% weight.Wdi = 100 * weight.Wdi; % Much more confidence in the ultimate dry
-%                                % weight than the other data points
-% uni-variate data: 
-% weight.tL = 2 * weight.tL;
-
-%% set pseudodata and respective weights
-% (pseudo data are in data.psd and weights are in weight.psd)
-[data, units, label, weight] = addpseudodata(data, units, label, weight);
-
-%% overwriting pseudodata and respective weights (remove these remarks after editing the file)
-% the pseudodata and respective weights were set automatically with the function setpseudodata
-% if one wants to overwrite one of the values it should always present an explanation
-% example:
-% data.psd.p_M = 1000;                    % my_pet belongs to a group with high somatic maint 
-% weight.psd.kap = 10 * weight.psd.kap;   % I need to give this pseudo data a higher weight
-
-%% Facts
-% list facts: F1, F2, etc.
-% make sure each fact has a corresponding bib key
-% do not put any DEB modelling assumptions here, only relevant information on
-% biology and life-cycles etc.
-F1 = 'The larval stage last 202 days and no feeding occurs';
-metadata.bibkey.F1 = 'Wiki'; % optional bibkey
-metadata.facts = struct('F1',F1);
-
-%% Discussion points
-D1 = 'Author_mod_1: I found information on the number of eggs per female as a function of length in Anon2013 that was much higher than in Anon2015 but chose to not include it as the temperature was not provided';
-% optional bibkey: bibkey.D1 = 'Anon2013';
-D2 = 'Author_mod_1: I was surprised to observe that the weight coefficient for ab changed so much the parameter values';     
-% optional bibkey: bibkey.D2 = 'Kooy2010';
-metadata.discussion = struct('D1', D1, 'D2', D2);
-
-
-%% pack data and txt_data for output
-data.weight = weight;
-data.temp = temp;
-txt_data.units = units;
-txt_data.label = label;
-txt_data.bibkey = bibkey;
-if exist('comment','var')
-  txt_data.comment = comment;
-end
-
-%% References
-  % the following two references should be kept-----------------------------------------------------------
-  bibkey = 'Kooy2010'; type = 'Book'; bib = [ ...  % used in setting of chemical parameters and pseudodata
-  'author = {Kooijman, S.A.L.M.}, ' ...
-  'year = {2010}, ' ...
-  'title  = {Dynamic Energy Budget theory for metabolic organisation}, ' ...
-  'publisher = {Cambridge Univ. Press, Cambridge}, ' ...
-  'pages = {Table 4.2 (page 150), 8.1 (page 300)}, ' ...
-  'howpublished = {\url{http://www.bio.vu.nl/thb/research/bib/Kooy2010.html}}'];
-  eval(['metadata.biblist.' bibkey, '= ''@', type, '{', bibkey, ', ' bib, '}'';']);
-  %
-  bibkey = 'LikaKear2011'; type = 'Article'; bib = [ ...  % used for the estimation method
-   'author = {Lika, K. and Kearney, M.R. and Freitas, V. and van der Veer, H.W. and van der Meer, J. and Wijsman, J.W.M. and Pecquerie, L. and Kooijman, S.A.L.M.},'...
-   'year = {2011},'...
-   'title = {The ''''covariation method'''' for estimating the parameters of the standard Dynamic Energy Budget model \textrm{I}: Philosophy and approach},'...
-   'journal = {Journal of Sea Research},'...
-   'volume = {66},'...
-   'number = {4},'...
-   'pages = {270-277},'...
-   'DOI = {10.1016/j.seares.2011.07.010},'...
-   'howpublished = {\url{http://www.sciencedirect.com/science/article/pii/S1385110111001055}}'];
-  eval(['metadata.biblist.' bibkey, '= ''@', type, '{', bibkey, ', ' bib, '}'';']);
-  %------------------------------------------------------------------------------------------------------
-
-  % References for the data, following BibTex rules
- % author names : author = {Last Name, F. and Last Name2, F2. and Last Name 3, F3. and Last Name 4, F4.}
- % latin names in title e.g. \emph{Pleurobrachia pileus}
-
-  bibkey = 'Wiki'; type = 'Misc'; bib = [...
-  'howpublished = {\url{http://en.wikipedia.org/wiki/my_pet}},'...% replace my_pet by latin species name
-  'note = {Accessed : 2015-04-30}'];
-  eval(['metadata.biblist.' bibkey, '= ''@', type, '{', bibkey, ', ' bib, '}'';']);
-  %
-  bibkey = 'MollCano2010'; type = 'Article'; bib = [ ... % meant as example; replace this and further bib entries
-  'author = {M{\o}ller, L. F. and Canon, J. M. and Tiselius, P.}, ' ... 
-  'year = {2010}, ' ...
-  'title = {Bioenergetics and growth in the ctenophore \emph{Pleurobrachia pileus}}, ' ...
-  'journal = {Hydrobiologia}, ' ...
-  'volume = {645}, ' ...
-  'number = {4}, '...
-  'pages = {167-178}'];
-  eval(['metadata.biblist.' bibkey, '= ''@', type, '{', bibkey, ', ' bib, '}'';']);
-  %
-  bibkey = 'Anon2015'; type = 'Misc'; bib = [ ...
-  'author = {Anonymous}, ' ...
-  'year = {2015}, ' ...
-  'howpublished = {\url{http://www.fishbase.org/summary/Rhincodon-typus.html}}'];
-  eval(['metadata.biblist.' bibkey, '= ''@', type, '{', bibkey, ', ' bib, '}'';']);
-
-
diff --git a/lib/pet/mydata_my_pet.m b/lib/pet/mydata_my_pet.m
index e755ae0e..80a6b008 100644
--- a/lib/pet/mydata_my_pet.m
+++ b/lib/pet/mydata_my_pet.m
@@ -11,7 +11,7 @@
   
   %% Description
   % Sets data, pseudodata, metadata, explanatory text, weight coefficients.
-  % Meant to be a template in add-my-pet
+  % Meant to be a template in add_my_pet
   %
   % Output
   %
@@ -29,7 +29,6 @@
   % * OPTIONAL : add discussion points / comments before the reference list
 
 %% set metadata
-%% blah bal albe
 
 T_C = 273.15; % K, temperature at 0 degrees C (used in T_typical)
 
@@ -68,14 +67,14 @@
 % here assumed to be equal for all real data; the value of f is specified in pars_init_my_pet.
 
 % age 0 is at onset of embryo development
-data.ab = 15;      units.ab = 'd';    label.ab = 'age at birth';  bibkey.ab = 'MollCano2010';   comment.ab  = 'mean value taken from several measurements'; 
+data.ab = 15;      units.ab = 'd';    label.ab = 'age at birth';                bibkey.ab = 'MollCano2010';    
   temp.ab = T_C + 20;  % K, temperature 
   % observed age at birth is frequently larger than ab, because of diapauzes during incubation
-data.ap = 261;     units.ap = 'd';    label.ap = 'age at puberty'; bibkey.ap = 'Anon2015';
+data.ap = 261;     units.ap = 'd';    label.ap = 'age at puberty';              bibkey.ap = 'Anon2015';
   temp.ap = T_C + 20;  % K, temperature 
   % observed age at puberty is frequently larger than ap, 
   %   because allocation to reproduction starts before first eggs appear
-data.am = 591;     units.am = 'd';    label.am = 'life span';     bibkey.am = 'Wiki';   
+data.am = 591;     units.am = 'd';    label.am = 'life span';                   bibkey.am = 'Wiki';   
   temp.am = T_C + 20;  % K, temperature 
 % (accounting for aging only) 
 
@@ -84,7 +83,7 @@
 % add an optional comment structure to give any additional explanations on
 % how the value was chosen, see the last column of the Lb data set for an
 % example
-data.Lb  = 0.45;   units.Lb  = 'cm';   label.Lb  = 'snout to vent length at birth';    bibkey.Lb  = 'Anon2015'; 
+data.Lb  = 0.45;   units.Lb  = 'cm';   label.Lb  = 'snout to vent length at birth';    bibkey.Lb  = 'Anon2015'; comment.Lb  = 'mean value taken from several measurements';
 data.Lp  = 2.36;   units.Lp  = 'cm';   label.Lp  = 'snout to vent length at puberty';  bibkey.Lp  = {'Anon2015','Wiki'}; % for multiple references, please use commas to separate references
 data.Li  = 6.25;   units.Li  = 'cm';   label.Li  = 'ultimate snout to vent length';    bibkey.Li  = 'Wiki';
 data.Wdb = 5.8e-5; units.Wdb = 'g';    label.Wdb = 'dry weight at birth';              bibkey.Wdb = 'Anon2015';
@@ -134,11 +133,22 @@
 % data.psd.p_M = 1000;                    % my_pet belongs to a group with high somatic maint 
 % weight.psd.kap = 10 * weight.psd.kap;   % I need to give this pseudo data a higher weight
 
+%% pack data and txt_data for output
+data.weight = weight;
+data.temp = temp;
+txt_data.units = units;
+txt_data.label = label;
+txt_data.bibkey = bibkey;
+if exist('comment','var')
+  txt_data.comment = comment;
+end
+
 %% Facts
 % list facts: F1, F2, etc.
 % make sure each fact has a corresponding bib key
 % do not put any DEB modelling assumptions here, only relevant information on
 % biology and life-cycles etc.
+
 F1 = 'The larval stage last 202 days and no feeding occurs';
 metadata.bibkey.F1 = 'Wiki'; % optional bibkey
 metadata.facts = struct('F1',F1);
@@ -150,17 +160,6 @@
 % optional bibkey: bibkey.D2 = 'Kooy2010';
 metadata.discussion = struct('D1', D1, 'D2', D2);
 
-
-%% pack data and txt_data for output
-data.weight = weight;
-data.temp = temp;
-txt_data.units = units;
-txt_data.label = label;
-txt_data.bibkey = bibkey;
-if exist('comment','var')
-  txt_data.comment = comment;
-end
-
 %% References
   % the following two references should be kept-----------------------------------------------------------
   bibkey = 'Kooy2010'; type = 'Book'; bib = [ ...  % used in setting of chemical parameters and pseudodata
diff --git a/lib/pet/petregr.m b/lib/pet/petregr.m
index 1aea96d8..96bdddc2 100644
--- a/lib/pet/petregr.m
+++ b/lib/pet/petregr.m
@@ -4,7 +4,7 @@
 %%
 function [q, info] = petregr(func, par, chem, T_ref, data)
 % created 2001/09/07 by Bas Kooijman; 
-% modified 2015/03/30 by Goncalo Marques
+% modified 2015/03/30, 2015/07/28 by Goncalo Marques
 
 %% Syntax
 % [q, info] = <../petregr.m *petregr*> (func, par, chem, T_ref, data)
@@ -36,7 +36,7 @@
 % The number of fields in data is variable
    
    
-  global report max_step_number max_fun_evals tol_simplex tol_fun
+  global report max_step_number max_fun_evals tol_simplex tol_fun simplex_size
 
   % option settings
   info = 1; % initiate info setting
@@ -98,6 +98,9 @@
   if ~exist('tol_fun','var') || isempty(tol_fun)
     nmregr_options('tol_fun', 1e-4);
   end
+  if ~exist('simplex_size','var') || isempty(simplex_size)
+    nmregr_options('simplex_size', 0.05);
+  end
   if ~exist('report','var') || isempty(report)
     nmregr_options('report', 1);
   end
@@ -115,7 +118,7 @@
   eval(['f = ', func, '(q, chem, T_ref, data4pred);']);
   eval(['fv(:,1) = W'' * ([', listf, '] - Y).^2;']);
   % Following improvement suggested by L.Pfeffer at Stanford
-  usual_delta = 0.05;             % 5 percent deltas for non-zero terms
+  usual_delta = simplex_size;     % 5 percent deltas is the default for non-zero terms
   zero_term_delta = 0.00025;      % Even smaller delta for zero elements of q
   for j = 1:n_par
     y = xin;
diff --git a/lib/pet/petregr_f.m b/lib/pet/petregr_f.m
index 96ea483c..d214442b 100644
--- a/lib/pet/petregr_f.m
+++ b/lib/pet/petregr_f.m
@@ -4,7 +4,7 @@
 %%
 function [q, info] = petregr_f(func, par, chem, T_ref, data, filter)
 % created 2001/09/07 by Bas Kooijman; 
-% modified 2015/01/29 by Goncalo Marques, 2015/03/21 by Bas Kooijman, 2015/03/30 by Goncalo Marques, 2015/04/27 by Goncalo Marques
+% modified 2015/01/29 by Goncalo Marques, 2015/03/21 by Bas Kooijman, 2015/03/30, 2015/04/27, 2015/07/28 by Goncalo Marques
 
 %% Syntax
 % [q, info] = <../petregr_f.m *petregr_f*> (func, par, chem, T_ref, data, filter)
@@ -40,7 +40,7 @@
 
    
    
-  global report max_step_number max_fun_evals tol_simplex tol_fun
+  global report max_step_number max_fun_evals tol_simplex tol_fun simplex_size
 
   % option settings
   info = 1; % initiate info setting
@@ -102,6 +102,9 @@
   if ~exist('tol_fun','var') || isempty(tol_fun)
     nmregr_options('tol_fun', 1e-4);
   end
+  if ~exist('simplex_size','var') || isempty(simplex_size)
+    nmregr_options('simplex_size', 0.05);
+  end
   if ~exist('report','var') || isempty(report)
     nmregr_options('report', 1);
   end
@@ -119,7 +122,7 @@
   eval(['f = ', func, '(q, chem, T_ref, data4pred);']);
   eval(['fv(:,1) = W'' * ([', listf, '] - Y).^2;']);
   % Following improvement suggested by L.Pfeffer at Stanford
-  usual_delta = 0.05;             % 5 percent deltas for non-zero terms
+  usual_delta = simplex_size;     % 5 percent deltas is the default for non-zero terms
   zero_term_delta = 0.00025;      % Even smaller delta for zero elements of q
   for j = 1:n_par
     y = xin;
diff --git a/lib/regr/nmregr_options.m b/lib/regr/nmregr_options.m
index 861aaeb4..b0dcb73f 100644
--- a/lib/regr/nmregr_options.m
+++ b/lib/regr/nmregr_options.m
@@ -3,7 +3,7 @@
 
 %%
 function nmregr_options (key, val)
-  %  created at 2002/02/10 by Bas Kooijman; modified 2015/01/16, 2015/02/27 Goncalo Marques
+  %  created at 2002/02/10 by Bas Kooijman; modified 2015/01/16, 2015/02/27, 2015/07/28 Goncalo Marques
   
   %% Syntax
   % <../nmregr_options .m *nmregr_options*> (key, val)
@@ -28,6 +28,8 @@ function nmregr_options (key, val)
   %       together to call them the same
   %    'tol_tun': tolerance for how close the loss-function values must be
   %       together to call them the same
+  %    'simplex_size': fraction added (subtracted if negative) to the 
+  %       free parameters when building the simplex
   %
   % Output
   %
@@ -36,7 +38,7 @@ function nmregr_options (key, val)
   %% Example of use
   % nmregr_options('default'); nmregr_options('report', 0)
 
-    global report max_step_number max_fun_evals tol_simplex tol_fun
+    global report max_step_number max_fun_evals tol_simplex tol_fun simplex_size
  
     if ~exist('key','var')
       key = 'inexistent';
@@ -50,6 +52,7 @@ function nmregr_options (key, val)
 	    max_fun_evals = 2000;
 	    tol_simplex = 1e-4;
 	    tol_fun = 1e-4;
+        simplex_size = 0.05;
 
       case 'report'
 	    if exist('val','var') == 0
@@ -105,6 +108,17 @@ function nmregr_options (key, val)
 	    else
 	      tol_fun = val;
         end
+
+      case 'simplex_size'
+	    if ~exist('val','var') 
+	      if numel(simplex_size) ~= 0
+	        fprintf(['simplex_size = ', num2str(simplex_size),' \n']);
+	      else
+	        fprintf('simplex_size = unknown \n');
+	      end
+	    else
+	      simplex_size = val;
+        end
         
       otherwise
         if ~strcmp(key, 'inexistent')   
@@ -135,4 +149,9 @@ function nmregr_options (key, val)
 	    else
 	      fprintf('tol_fun = unknown \n');
 	    end
+	    if numel(simplex_size) ~= 0
+	      fprintf(['simplex_size = ', num2str(simplex_size),' \n']);
+	    else
+	      fprintf('simplex_size = unknown \n');
+	    end
      end
\ No newline at end of file
