%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
close all
clc
warning off all 
%
  TC_FILE		=	'/data/olvac/travail_en_cours/graffin/Climate_Index/data_tc/wptracks.nc';


% Period:
% ------
F_year		=	2000; % first year minimum 1949 for data West Pac
L_year		=	2017; % last year  maximum 2012 for data West Pac (2015 for EPac)
%
% Box area for calculation:


		lon_min = 120; %BASSIN PACIFIQUE OUEST
		lon_max = 140;
		lat_min = 15;
		lat_max = 25;
        
%       lon_min = 109; %BASSIN MER DE CHINE
% 		lon_max = 120;
% 		lat_min = 12;
% 		lat_max = 23;
        
% -------------------------
% lon_min		=	120; % East Degree        %
% lon_max		=	150; % East Degree        % ==> WEST
% lat_min		=	12;   % North Degree       % ==>      PACIFIC
% lat_max		=	30;  % North Degree       %
%
  % lon_min		=	140; % East Degree    %
  % lon_max		=	220; % East Degree    % ==> South
  % lat_min		=	-40;   % North Degree   % ==>      PACIFIC
  % lat_max		=	5;  % North Degree   %
%
% lon_min		=	260; % East Degree    %
% lon_max		=	360; % East Degree    % ==> Atlantic
% lat_min		=	5;   % North Degree   % ==>
% lat_max		=	40;  % North Degree   %

%  lon_min		=	180; % East Degree    %
%  lon_max		=	260; % East Degree    % ==> Atlantic
%  lat_min		=	5;   % North Degree   % ==>
%  lat_max		=	40;  % North Degree   %


  WS_EP    = ncread(TC_FILE,'vsmc');
  LAT_EP   = ncread(TC_FILE,'latmc');
  LON_EP   = ncread(TC_FILE,'longmc');
  HOUR_EP  = ncread(TC_FILE,'hourmc');
  DAY_EP   = ncread(TC_FILE,'daymc');
  MONTH_EP = ncread(TC_FILE,'monthmc');
  YEAR_EP  = ncread(TC_FILE,'yearic');
  %
  YEAR_EP=YEAR_EP';
  WS_EP    = WS_EP';
  LAT_EP   = LAT_EP';
  LON_EP   = LON_EP';
  HOUR_EP  = HOUR_EP';
  DAY_EP   = DAY_EP';
  MONTH_EP = MONTH_EP';
%


WS_EP=WS_EP*1.944012;       % convert unit : m/s ==> knot

 [ACEd2,INTENSITYd2,DENSITYalld2,ACEm2,INTENSITYm2,DENSITYallm2,DENSITYmajorm2,ACEm_anom2,INTENSITYm_anom2,DENSITYallm_anom2,DENSITYmajorm_anom2,Seasonal_ACEm,Seasonal_INTENSITYm,Seasonal_DENSITYallm,Seasonal_DENSITYmajorm] = calc_ACE_box_w_daily(F_year,L_year,lon_min,lon_max,lat_min,lat_max,WS_EP,LAT_EP,LON_EP,HOUR_EP,DAY_EP,MONTH_EP,YEAR_EP);

save('/data/olvac/travail_en_cours/graffin/Coastline/GEE/Scripts/test_ROIg1-2_2000_2017/TC_anomalies_wp.mat','ACEm_anom2','DENSITYmajorm_anom2');
save('/data/olvac/travail_en_cours/graffin/Coastline/GEE/Scripts/test_ROIg1-2_2000_2017/coordo_TC_box.mat','lon_max','lon_min','lat_max','lat_min');