%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  Add radswbio/swradbio (solar short-wave radiation without diurnal cycle) 
%  in frc/bulk climatological atmopheric forcing files. 
%  Variable needed by biogeochemical model PISCES
%
%
%  Further Information:
%  http://www.crocoagrif.org
%
%  This file is part of CROCOTOOLS
%
%  CROCOTOOLS is free software; you can redistribute it and/or modify
%  it under the terms of the GNU General Public License as published
%  by the Free Software Foundation; either version 2 of the License,
%  or (at your option) any later version.
%
%  CROCOTOOLS is distributed in the hope that it will be useful, but
%  WITHOUT ANY WARRANTY; without even the implied warranty of
%  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%  GNU General Public License for more details.
%
%  You should have received a copy of the GNU General Public License
%  along with this program; if not, write to the Free Software
%  Foundation, Inc., 59 Temple Place, Suite 330, Boston,
%  MA  02111-1307  USA
%
%
%  June 2014: G.Cambon (IRD/LEGOS)
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
close all

disp(['Add swradbio: daily averaged solar short wave radiation     '])
disp(['Creating variables and attributes for forcing and bulk ' ...
      'files'])

crocotools_param
%
% Open the frc file
% 
if (makefrc)
disp(['Add swradbio: daily averaged solar short wave radiation'])
disp(['No change needed in this climato case                  '])
disp(['Creating variables and attributes for frc files        '])
end
oct_nc_add_swradbio_frc(frcname,'r');
%
nw_id = netcdf.open(frcname, 'NC_WRITE');
swrad1=netcdf.getVar(nw_id, netcdf.inqVarID(nw_id, 'swrad'));
netcdf.putVar(nw_id, netcdf.inqVarID(nw_id, 'swradbio'), swrad1);
netcdf.close(nw_id);
%
% Open the bulk time
%
if (makeblk)
disp(['Add radswbio: daily averaged solar short wave radiation'])
disp(['No change needed in this climato case                  '])
disp(['Creating variables and attributes for bulk files       '])
end
oct_nc_add_swradbio_blk(blkname);
%
nw_id = netcdf.open(blkname, 'NC_WRITE');
radsw1=netcdf.getVar(nw_id, netcdf.inqVarID(nw_id, 'radsw'));
netcdf.putVar(nw_id, netcdf.inqVarID(nw_id, 'radswbio'), radsw1);
netcdf.close(nw_id);
%