clear all
close all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  Extrapole and interpole biology data and write in climatology file
%
%  Data source : IRI/LDEO Climate Data Library (World Ocean Atlas 1998)
%    http://ingrid.ldgo.columbia.edu/
%    http://iridl.ldeo.columbia.edu/SOURCES/.NOAA/.NODC/.WOA2001/
%
%  Pierrick Penven, IRD, 2005.
%  Olivier Aumont the master, IRD, 2006.
%  Patricio Marchesiello, chief, IRD, 2007.
%  Christophe Eugene Raoul Menkes, the slave, IRD, 2007.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
close all
%%%%%%%%%%%%%%%%%%%%% USERS DEFINED VARIABLES %%%%%%%%%%%%%%%%%%%%%%%%
%
%  Switches for selecting what to process (1=ON)
%
crocotools_param
%
% Climatology files
%
%
no3_seas_data  = [woapisces_dir,'no3_seas.cdf'];
no3_ann_data   = [woapisces_dir,'no3_ann.cdf'];
po4_seas_data  = [woapisces_dir,'po4_seas.cdf'];
po4_ann_data   = [woapisces_dir,'po4_ann.cdf'];
o2_seas_data   = [woapisces_dir,'o2_seas.cdf'];
o2_ann_data    = [woapisces_dir,'o2_ann.cdf'];
sio3_seas_data = [woapisces_dir,'sio3_seas.cdf'];
sio3_ann_data  = [woapisces_dir,'sio3_ann.cdf'];
dic_seas_data  = [woapisces_dir,'dic_seas.cdf'];
dic_ann_data   = [woapisces_dir,'dic_ann.cdf'];
talk_seas_data = [woapisces_dir,'talk_seas.cdf'];
talk_ann_data  = [woapisces_dir,'talk_ann.cdf'];
doc_seas_data  = [woapisces_dir,'doc_seas.cdf'];
doc_ann_data   = [woapisces_dir,'doc_ann.cdf'];
fer_seas_data  = [woapisces_dir,'fer_seas.cdf'];
fer_ann_data   = [woapisces_dir,'fer_ann.cdf'];
dust_seas_data = [woapisces_dir,'dust_seas.cdf'];
dust_ann_data  = [woapisces_dir,'dust_ann.cdf'];
%
cycle=woa_cycle;
NO3min=1;
if strcmp(climato_dir,cars2009_dir);
    err_msg=sprintf(['Error : you need to use woadir when creating the croco_oa.nc (Z) \n'...
                     'files to be compatible with PISCES'])
    error(err_msg)
end
    
%
%%%%%%%%%%%%%%%%%%% END USERS DEFINED VARIABLES %%%%%%%%%%%%%%%%%%%%%%%
%
% Add variables in the files
%
disp('')
disp('====================================================== ')
disp('=> You need the croco_oa.nc file created by oct_make_clim.m ')
disp('=> with makeoa=1 from crocotools_param.m                ')
disp('====================================================== ')
oct_add_no3(oaname,clmname,ininame,grdname,no3_seas_data,...
        no3_ann_data,cycle,makeoa,makeclim,makeini)

oct_add_dic(oaname,clmname,ininame,grdname,dic_seas_data,...
       dic_ann_data,cycle,makeoa,makeclim,makeini)

oct_add_po4(oaname,clmname,ininame,grdname,po4_seas_data,...
        po4_ann_data,cycle,makeoa,makeclim,makeini)

oct_add_sio3(oaname,clmname,ininame,grdname,sio3_seas_data,...
        sio3_ann_data,cycle,makeoa,makeclim,makeini)

oct_add_o2(oaname,clmname,ininame,grdname,o2_seas_data,...
        o2_ann_data,cycle,makeoa,makeclim,makeini)

oct_add_talk(oaname,clmname,ininame,grdname,talk_seas_data,...
       talk_ann_data,cycle,makeoa,makeclim,makeini)

oct_add_doc(oaname,clmname,ininame,grdname,doc_seas_data,...
       doc_ann_data,cycle,makeoa,makeclim,makeini)

if makequota
   oct_add_don(oaname,clmname,ininame,grdname,doc_seas_data,...
       doc_ann_data,cycle,makeoa,makeclim,makequota,makeini)
   oct_add_dop(oaname,clmname,ininame,grdname,doc_seas_data,...
       doc_ann_data,cycle,makeoa,makeclim,makequota,makeini)
end
            
oct_add_fer(oaname,clmname,ininame,grdname,fer_seas_data,...
       fer_ann_data,cycle,makeoa,makeclim,makeini)

%
% Horizontal extrapolation
%
if (makeoa)
  oct_ext_tracers(oaname,no3_seas_data,no3_ann_data,...
              'nitrate','NO3','no3_time','Zno3',Roa);

  oct_ext_tracers(oaname,po4_seas_data,po4_ann_data,...
              'phosphate','PO4','po4_time','Zpo4',Roa);

  oct_ext_tracers(oaname,sio3_seas_data,sio3_ann_data,...
              'silicate','Si','si_time','Zsi',Roa);

  oct_ext_tracers(oaname,o2_seas_data,o2_ann_data,...
              'oxygen','O2','o2_time','Zo2',Roa);

  oct_ext_tracers(oaname,dic_seas_data,dic_ann_data,...
             'dic','DIC','dic_time','Zdic',Roa);

  oct_ext_tracers(oaname,talk_seas_data,talk_ann_data,...
              'talk','TALK','talk_time','Ztalk',Roa);

  oct_ext_tracers(oaname,doc_seas_data,doc_ann_data,...
             'doc','DOC','doc_time','Zdoc',Roa);

     if makequota
     oct_ext_tracers(oaname,doc_seas_data,doc_ann_data,...
             'doc','DON','don_time','Zdon',Roa);

     oct_ext_tracers(oaname,doc_seas_data,doc_ann_data,...
             'doc','DOP','dop_time','Zdop',Roa);
     end

  oct_ext_tracers(oaname,fer_seas_data,fer_ann_data,...
             'fer','FER','fer_time','Zfer',Roa);
end
%
% Vertical interpolations 
%
if (makeclim)
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' NO3...')
  oct_vinterp_clm(clmname,grdname,oaname,'NO3','no3_time','Zno3',0,'r');
%
%  PO4
%
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' PO4...')
  oct_vinterp_clm(clmname,grdname,oaname,'PO4','po4_time','Zpo4',0,'r');
 %
 %  Si
 %
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' Si...')
  oct_vinterp_clm(clmname,grdname,oaname,'Si','si_time','Zsi',0,'r');
%
%  O2
%
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' O2...')
  oct_vinterp_clm(clmname,grdname,oaname,'O2','o2_time','Zo2',0,'r');

%
% DIC 
%
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' DIC...')
  oct_vinterp_clm(clmname,grdname,oaname,'DIC','dic_time','Zdic',0,'r');
%
% TALK
%
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' TALK...')
  oct_vinterp_clm(clmname,grdname,oaname,'TALK','talk_time','Ztalk',0,'r');

%
% DOC
%
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' DOC...')
  oct_vinterp_clm(clmname,grdname,oaname,'DOC','doc_time','Zdoc',0,'r');

if makequota
%
% DON
%
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' DON...')
  oct_vinterp_clm(clmname,grdname,oaname,'DON','don_time','Zdon',0,'r');
%
% DOP
%
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' DOP...')
  oct_vinterp_clm(clmname,grdname,oaname,'DOP','dop_time','Zdop',0,'r');
end

%
% FER
%
  disp(' ')
  disp(' Vertical interpolations')
  disp(' ')
  disp(' FER...')
  oct_vinterp_clm(clmname,grdname,oaname,'FER','fer_time','Zfer',0,'r');
end

% %----------------------------------------------------------------------------
% % Initial file
% %----------------------------------------------------------------------------
% if (makeini)
% %
% %  NO3
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' NO3...')
%   oct_vinterp_clm(ininame,grdname,oaname,'NO3','no3_time','Zno3',tini,'r',1);
% %
% % DIC
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' DIC...')
%   oct_vinterp_clm(ininame,grdname,oaname,'DIC','dic_time','Zdic',tini,'r',1);
% %
% %  PO4
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' PO4...')
%   oct_vinterp_clm(ininame,grdname,oaname,'PO4','po4_time','Zpo4',tini,'r',1);
%  %
%  %  Si
%  %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' Si...')
%   oct_vinterp_clm(ininame,grdname,oaname,'Si','si_time','Zsi',tini,'r',1);
% %
% %  O2
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' O2...')
%   oct_vinterp_clm(ininame,grdname,oaname,'O2','o2_time','Zo2',tini,'r',1);
% %
% % TALK
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' TALK...')
%   oct_vinterp_clm(ininame,grdname,oaname,'TALK','talk_time','Ztalk',tini,'r',1);
% %
% % DOC
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' DOC...')
%   oct_vinterp_clm(ininame,grdname,oaname,'DOC','doc_time','Zdoc',tini,'r',1);

% if makequota
% %
% % DON
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' DON...')
%   oct_vinterp_clm(ininame,grdname,oaname,'DON','don_time','Zdon',tini,'r',1);
% %
% % DOP
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' DOP...')
%   oct_vinterp_clm(ininame,grdname,oaname,'DOP','dop_time','Zdop',tini,'r',1);
% end

% %
% % FER
% %
%   disp(' ')
%   disp(' Vertical interpolations')
%   disp(' ')
%   disp(' FER...')
%   oct_vinterp_clm(ininame,grdname,oaname,'FER','fer_time','Zfer',tini,'r',1);
% end

if (makeplot)
disp(' ')
disp(' Make a few plots...')
oct_test_clim(clmname,grdname,'NO3',1,coastfileplot)
figure
oct_test_clim(clmname,grdname,'PO4',1,coastfileplot)
figure
oct_test_clim(clmname,grdname,'Si',1,coastfileplot)
figure
oct_test_clim(clmname,grdname,'O2',1,coastfileplot)
figure
oct_test_clim(clmname,grdname,'DIC',1,coastfileplot)
figure
oct_test_clim(clmname,grdname,'TALK',1,coastfileplot)
figure
oct_test_clim(clmname,grdname,'DOC',1,coastfileplot)
figure
if makequota
oct_test_clim(clmname,grdname,'DON',1,coastfileplot)
figure
end
oct_test_clim(clmname,grdname,'FER',1,coastfileplot)
figure
end
%
% End
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%