%date : 17/02/2005 
%auteur : Gael Forget 
%version : beta

%domaine.m : charge les parametres principaux associes 
%	au domaine (spatio_temporel)

clear tmp*;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%0.1) declaration des variables en global
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
domaine_global_def;

rep_domaine='/net/ross/raid0/gforget/1x1_50levels/GRID/';
%rep_domaine='/net/altix3700/raid5/gforget/exp_CTRL/exp_OFFLINE/exp1/';

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%0.2) valeurs de reference
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r0=6.37e6; %rayon terrestre
deg2rad = 2*pi/360; %conversion angles
gravity=9.81; %acceleration gravitationnelle 
%rho0=999.8; %densite de ref pour l'eau
%rho0=1035;
rho0=1029;

%%%%%%%%%%%%%%%%%%%%%%%%%
%1) domaine géographique : 
%%%%%%%%%%%%%%%%%%%%%%%%%
%dimensions du probleme :
jpi=360; jpj=160; jpk=50;

%masks:
fid=fopen([rep_domaine 'maskCtrlC.data'],'r','b');
tmask3D=fread(fid,[jpi*jpj jpk],'float32'); fclose(fid);
tmask3D=reshape(tmask3D,[jpi jpj jpk]);tmask3D(find(tmask3D==0))=NaN;
fid=fopen([rep_domaine 'maskCtrlW.data'],'r','b');
umask3D=fread(fid,[jpi*jpj jpk],'float32'); fclose(fid);
umask3D=reshape(umask3D,[jpi jpj jpk]);umask3D(find(umask3D==0))=NaN;
fid=fopen([rep_domaine 'maskCtrlS.data'],'r','b');
vmask3D=fread(fid,[jpi*jpj jpk],'float32'); fclose(fid);
vmask3D=reshape(vmask3D,[jpi jpj jpk]);vmask3D(find(vmask3D==0))=NaN;

%positions:
fid=fopen([rep_domaine 'XC.data'],'r','b');
lon2D_t=fread(fid,[jpi jpj],'float32'); fclose(fid);
fid=fopen([rep_domaine 'YC.data'],'r','b');
lat2D_t=fread(fid,[jpi jpj],'float32'); fclose(fid);
fid=fopen([rep_domaine 'XG.data'],'r','b');
lon2D_u=fread(fid,[jpi jpj],'float32'); fclose(fid);
fid=fopen([rep_domaine 'YG.data'],'r','b');
lat2D_v=fread(fid,[jpi jpj],'float32'); fclose(fid);

%depths:
fid=fopen([rep_domaine 'RC.data'],'r','b');
gdept=-fread(fid,jpk,'float32'); fclose(fid);
fid=fopen([rep_domaine 'RF.data'],'r','b');
gdepw=-fread(fid,jpk+1,'float32'); fclose(fid);

%restriction du mask a 1 bassin :
%load MASK_pac.mat; 
%load /data4/gforget/FIRSTdiags/various_mat/MASK_atl.mat;
%load /data4/gforget/FIRSTdiags/MASK_indpac.mat;
%load /data4/gforget/FIRSTdiags/MASK_south.mat;
%for kcur=1:jpk; tmask3D(:,:,kcur)=squeeze(tmask3D(:,:,kcur)).*mask_boites; end; umask3D=umask3D.*tmask3D; vmask3D=vmask3D.*tmask3D;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%2) facteurs d'echelle :
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%fid=fopen([rep_domaine 'DXG.data'],'r','b');
fid=fopen([rep_domaine 'DXC.data'],'r','b');
e1t=fread(fid,[jpi jpj],'float32'); fclose(fid);
fid=fopen([rep_domaine 'DYG.data'],'r','b');
e2t=fread(fid,[jpi jpj],'float32'); fclose(fid);
%fid=fopen([rep_domaine 'DXC.data'],'r','b');
fid=fopen([rep_domaine 'DXG.data'],'r','b');
e1v=fread(fid,[jpi jpj],'float32'); fclose(fid);
%e2v=e2t;
fid=fopen([rep_domaine 'DYC.data'],'r','b');
e2v=fread(fid,[jpi jpj],'float32'); fclose(fid);
%fid=fopen([rep_domaine 'DYC.data'],'r','b');
%e2u=fread(fid,[jpi jpj],'float32'); fclose(fid);
e2u=e2t;
e1u=e1t;


fid=fopen([rep_domaine 'DRF.data'],'r','b');
e3t=fread(fid,jpk,'float32'); fclose(fid);
fid=fopen([rep_domaine 'DRC.data'],'r','b');
e3w=fread(fid,jpk+1,'float32'); fclose(fid);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%3) parametres temporels :
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%1) temps 0 du modele en jours julens :
model_time0=jul_0h([2002 01 01 06 00 00]);  %en jours juliens
%model_time0=jul_0h([1992 01 01 06 00 00]);  %en jours juliens
%TEMPO%
%model_time0=0;
%2) pas de temps du modele :
dtime=3600;
%explication : ensuite pour passer d'un temps modele a une date :
%date_tmp=greg_0h(model_time0+model_time/86400);	


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%4) diffusivites/viscosites 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


% Coefficient de viscosite (U,V) du modele
Vh=1e4; % horizontale
Vv=1e-3; % vertical

% Coefficient de diffusivite (T,S) du modele
Kh=1e2; %horizontale
Kv=1e-5; %verticale