calcrn_vect.m

function [albs,arns,brns,crns,arnv,brnv,crnv,cras,crav]=calcrn_vect(rg,ratm,emiss,emisv,albe,albv,sigmav,sigmava)


    albs=(albe-sigmav*albv)./((1-sigmav).^2+sigmav*albv.*albe-(sigmav*albv).^2);

	v1=1-albv*albs.*sigmav;
    v1a=1-albv*albs.*sigmava;
	v2=1-emisv;
	v3=1-emiss;
	v4=1-sigmav;
    v4a=1-sigmava;
	v5=1-sigmav*v2*v3;
    v5a=1-sigmava*v2*v3;
    
    arns=-(v4a*emiss+emisv*emiss*sigmava)./v5a;
    brns=emisv*emiss*sigmava./v5a;
    crns=(rg*(1-albs).*v4./v1)+(v4a*emiss*ratm./v5a);
    cras=(v4a*emiss*ratm./v5a);
    
    arnv=brns;
    brnv=-sigmava.*(emisv+(emisv*emiss+v4a.*v3*emisv)./v5a);
    crnv=rg*(1-albv)*sigmav.*(1+(albs.*v4)./v1)+sigmava*emisv*ratm.*(1+(v4a.*v3)./v5a);
    crav=sigmava*emisv*ratm.*(1+(v4a.*v3)./v5a);

end