%Parameters of the model

%% Parameters
epsilon=3.8571;               % Elasticity of substitution between differentiated goods (NAWM)
sig=2;                        % risk aversion (NAWM)
ies=0.5;                      % inverse of Frisch elasticity (NAWM)
delta=0.025;                  % depreciation rate (NAWM)
chi=0.95;                     % survivial rate of bankers (NAWM)
csi=2;                        % elasticity of substitution btw green and brown capital (CHM 21) 
ps=0.4;                       % elasticity of emissions with respect to output (CHM 21)
deltax=0.0035;                % pollution depreciation (GH 21)
D0=-0.0076;                   % constant in damage function (GH 21)
D1tilde=8.10e-6;              % 1st order coefficient term in damage function, before adjustment  (GH 21)
D2tilde=1.05e-8;              % 2nd order coefficient term in damage function, before adjustment  (GH 21)
alfa=0.36;                    % elasticity of production wrt capital (NAWM)
pi=1.005;                     % steady state inflation (ECB target)
omega=0.32;                   % Labor Disutility
pl=1;
kappaFG=3;
theta1 = 0.0334;
theta2 = 2.6;

% Targets
rr=1.005;                     % real interest rate: 2% in the annualized steady state (NAWM)
RoW=15.31;                    % rr Row RoW emissions/EA emissions (our_world_in_data) 
KG=0.2;                       % green capital out of total capital (Giovanardi et. al)
levss=6;                        % steady-state leverage (NAWM)
spss=0.0217/4;                  % steady state corporate spread (NAWM)
bP_Y=0.0254;                       % GDP share of corporate bonds held by the public sector (ECB); 
G=0.215;                      % public spending/GDP ratio (GK11)
XGtC=851;                     % DICE 2016R2

% Steady State for bank values
betta=1/rr;                                 % discount factor
iota=(1-chi*(spss*levss+1/betta))/levss;          % transfer for new bankers


% Parameters that do not affect the ss
kappaI=10.78;             % investment adjustment cost (NAWM)
phipi=2.74;               % mp response to inflation (NAWM) 
phiy=0.03;                   % mp response to output (NAWM) 
rhor=0.93;                % monetary policy inertia (NAWM)
rhoa=0.92;                % tfp persistence (NAWM)
rhoc=0.95;                % climate shock persistence
rho_xi=0.66;
calvo=0.82;               % price rigidity in calvo framework (NAWM)
kappaP=(epsilon-1)*calvo/(pi^2*(1-calvo)*(1-betta*calvo)); % 
psi=10;

spGss=spss;                                         % green spread            
spBss=spss;                                         % brown spread
rG=rr+spGss;                                      % loan rate for green bonds
rB=rr+spBss;                                      % loan rate for brown bonds
rkG=rG-(1-delta);                               % rental rate of green capital
rkB=rB-(1-delta);                               % rental rate of brown capital
pI=(epsilon-1)/epsilon;                         % real marginal costs
r=pi/betta;                                     % nominal interest rate
D=D0+D1tilde*XGtC+D2tilde*XGtC^2;               % damage function
A=1-D;  
piss=1.005;
nuss=(1-chi)/(1-chi*betta*(spss*levss+1/betta));  % bank discount factor
theta=nuss/levss+betta*nuss*spss;                   % absconding share
muss=0.0135;
levGss=.2*levss;
levBss=levss-levGss;
muGss=0.2;
muBss=0.8;
bstar=levGss/levss;
erow=13.203483300117536;
sgss=muss*muGss;
D0=-0.0076;
D1=1.715202229340095e-06;
D2=4.708146047721220e-10;
g=0.226201639026692;
zeta=0.861750832965662;
alpha0=2.505570065159853e-04;
G=0.215;   
Abar=0.9931;
nuss=(1-chi)/(1-chi*betta*(spss*levss+1/betta));  % bank discount factor
theta=nuss/levss+betta*nuss*spss;                   % absconding share
rss=1.0100;
KG=.2;
erow=13.203483300117536;

save par_ce betta alfa epsilon delta sig ies chi csi iota theta ps deltax D0 D1 D2 erow omega pl  ...
       zeta muGss kappaI phipi phiy rhor  rhoa rho_xi kappaP A sgss kappaFG G  nuss spss spGss spBss...
      levss levGss levBss bstar bP_Y rhoc rss  alpha0 psi muss muBss theta theta1 theta2 piss rr 
save par_sp zeta g erow  