function [elev_z,az_N]=aircraft2earth_tail(elev_tail, az_a, pitch, roll, drift, trackp)


sinpitch = sin(deg2rad(pitch));
cospitch = cos(deg2rad(pitch));
sinrollaz  = sin(deg2rad(roll+az_a));
cosrollaz  = cos(deg2rad(roll+az_a));
sindrift = sin(deg2rad(drift));
cosdrift = cos(deg2rad(drift));
coselev = cos(deg2rad(elev_tail));
sinelev = sin(deg2rad(elev_tail));


for azi=1:length(az_a)
    rg=1;
    
    xt=cosrollaz(azi)*sindrift(azi)*coselev(azi)*sinpitch(azi) + cosdrift(azi)*sinrollaz(azi)*coselev(azi) - sindrift(azi)*cospitch(azi)*sinelev(azi);
    yt=-cosrollaz(azi)*cosdrift(azi)*coselev(azi)*sinpitch(azi) + sindrift(azi)*sinrollaz(azi)*coselev(azi)  +cospitch(azi)*cosdrift(azi)*sinelev(azi);
    zt=cospitch(azi)*cosrollaz(azi)*coselev(azi) + sinpitch(azi)*sinelev(azi);
    
    Xt_over_r = [xt;yt;zt];
            
            
    
    
    
    tilt_t(rg,azi)=asin((Xt_over_r(2)));
    tilt_t(rg,azi)=rad2deg(tilt_t(rg,azi));
    
    elev_t(rg,azi)=asin((Xt_over_r(3)));
    elev_z(rg,azi)=rad2deg(elev_t(rg,azi));
    
    az_t(rg,azi)=atan2(Xt_over_r(1),Xt_over_r(2));
    az_t(rg,azi)=rad2deg(az_t(rg,azi));
    az_N(rg,azi)=az_t(rg,azi) + trackp(rg,azi);
    
end

    u_neg=find(az_N <0);
    az_N(u_neg)=az_N(u_neg)+360;
    u_360=find(az_N>=360);
    az_N(u_360)=az_N(u_360)-360;