
%Units
KM = 1000;
PW = 1.e15;
DAY = 24*3600;
MONTH = 30*DAY; 
YEAR = 12*MONTH;
MB = 100; %millibar  (Pa)
SV = 1.e09; %Sverdrup, mass transport unit (kg/s)

%Planetary & Physical constants
RADIUS = 6371*KM;       %Earth Radius (km)
GRAVI = 9.81;           %Earth gravity (ms-2)
SIGMA = 5.67*1.e-08;    %Stefan Boltzmann (Wm-2K-4)
So = 1367;              %Solar constant (Wm-2) 
ALPHAp = .3;            %Planetary albedo (non dim.)

%Ocean
RHO = 1000;             %density of water (kg m-3)
CPO = 4000;             %specific heat capacity of water (J kg-1 K-1)
HM  = 50;               %mixed layer thickness
HO = 500;               %Thermocline thickness (m)
HCM = RHO*CPO*HM;       %heat capacity of mixed layer (J K-1 m-2)
HCO = RHO*CPO*HO;       %heat capacity of thermocline (J K-1 m-2)

%Amosphere
CPA = 1000; %specific heat capacity of dry air (J kg-1 K-1)
CAheat = 2*CPA; %pseudo-heat capacity (to include moisture effects)
PATM = 1000*MB; %surface atm. pressure (Pa)
LV = 2.5*1.e06; %Latent heat of vaporization
HCA = CAheat*PATM/GRAVI; %atmospheric heat capacity (J K-1 m-2)

%Convective & circulation parameterizations
DTcrit_conv = 40; %moist adiabatic lapse rate (K)
PSIfrac = 0.1; % ratio PSIo/PSIa of circulation intensity
BIGONE = 100;
KEFF = 100*SV/15; %linear parameterization

%Radiation
TEMI = (So*(1-ALPHAp)/(4*SIGMA))^.25; %Global mean emission temp. (K)
TEMI1 = ((So*(1-ALPHAp)*(pi/3+sqrt(3)/2))/(2*pi*SIGMA))^.25; %Box1 Emission temp.(K)
TEMI2 = ((So*(1-ALPHAp)*(2*pi/3-sqrt(3)/2))/(2*pi*SIGMA))^.25; %Box2 Emission temp.(K)
GAMMA = 1.25; %non dim. factor for H2o emissivity calc. (See Emanuel, 2002)
ALPHA = .0012; %non dim. factor for CO2 emissivity calc. (see Emanuel, 2002)