% gfd_flightlevel_set_variables
% sets variables for core data
% called from: gfd_run_plot or gfd_profile_plot

switch i(k)     % Set variable to be plotted in panel j, line k
    case 1
        yd=lat; y_label='latitude (deg)'; y_label2=''; y_lim=[59 61]; y_tick=1;
    case 2
        yd=lon; y_label='longitude (deg)'; y_label2=''; y_lim=[-43 -40]; y_tick=1;
    case 3
        yd=heading_inu; y_label='heading (deg)'; y_label2=''; y_lim=[0 360]; y_tick=60;
    case 4
        yd=radar_alt; y_label='height (m)'; y_label2='radar alt.'; y_lim=[0 200]; y_tick=50;
    case 5
        yd=pressure_alt; y_label='height (m)'; y_label2='pressure alt.'; y_lim=[0 200]; y_tick=50;
    case 6
        yd=gps_alt; y_label='gps height (m)'; y_label2='gps alt.'; y_lim=[70 120]; y_tick=5;
    case 7 
        yd=static_pressure; y_label='static pressure (hPa)'; y_label2=''; y_lim=[980 1000]; y_tick=10;
    case 8
        yd=pressure_alt - gps_alt'; y_label='height (m)'; y_label2='diff alt.'; y_lim=[50 100]; y_tick=50;
        
    case 10
        yd=air_temp; y_label='Temperature (K)'; y_label2='T_{air} non-deiced'; y_lim=[273 278]; y_tick=1;
    case 11
        yd=deiced_air_temp; y_label='Temperature (K)'; y_label2='T_{air} deiced'; y_lim=[273 278]; y_tick=1;
    case 12
        yd=sfc_temp; y_label='Temperature (K)'; y_label2='T_{surface}'; y_lim=[273 278]; y_tick=1;
    case 13
        yd=sfc_temp_lwu; y_label='Temperature (K)'; y_label2='T_{sfc-lwu}'; y_lim=[273 279]; y_tick=1;
    case 14
        yd=dewpoint_ge; y_label='Temperature'; y_lim=[273 278]; y_label2='T_{dew}'; y_lim=[273 278]; y_tick=1;
    case 15
        yd=air_theta; y_label='\theta (K)'; y_label2='\theta'; y_lim=[275 280]; y_tick=1;
    case 16
        yd=air_thetae; y_label='\theta _e (K)'; y_label2='\theta _e'; y_lim=[285 291]; y_tick=1; 
        %y_lim=[284 306]; y_tick=4; % use for profiles
        
    case 20
        yd=ws; y_label='wind speed (m s^{-1})'; y_label2=''; y_lim=[15 35]; y_tick=5;
    case 21
        yd=wd; y_label='wind direction (deg.)'; y_label2=''; y_lim=[0 120]; y_tick=30;
    case 22
        yd=u; y_label='u (m s^{-1})'; y_label2=''; y_lim=[-30 0]; y_tick=5;
    case 23 
        yd=v; y_label='v (m s^{-1})'; y_label2=''; y_lim=[-30 0]; y_tick=5;
    case 24
        yd=w; y_label='w (m s^{-1})'; y_label2=''; y_lim=[-3 3]; y_tick=1;
        
    case 30
        yd=q_ge; y_label='specific humidity (g kg^{-1})'; y_label2='q'; y_lim=[3 6]; y_tick=1;
    case 31
        yd=qsat; y_label='specific humidity (g kg^{-1})'; y_label2='q_{sat}'; y_lim=[3 6]; y_tick=1;
    case 32
        yd=qsati; y_label='specific humidity (g kg^{-1})'; y_label2='q_{sati}'; y_lim=[3 6]; y_tick=1;
    case 33
        yd=rh_ge; y_label='relative humidity (%)'; y_label2='RH'; y_lim=[20 140]; y_tick=20;
    case 34
        yd=rhi_ge; y_label='relative humidity (g kg^{-1})'; y_label2='RH_{ice}'; y_lim=[20 140]; y_tick=20;
        
    case 40
        yd=lwc_jw.*rho; y_label='liquid water content (g m^{-3})'; y_label2='JW'; y_lim=[0 0.8]; y_tick=0.2;
    case 41
        yd=lwc_nev; y_label='water content (g m^{-3})'; y_label2='LWC-NEV'; y_lim=[0 2]; y_tick=0.2; % Not working!
    case 42
        yd=twc_nev; y_label='water content (g m^{-3})'; y_label2='TWC-NEV'; y_lim=[0 0.8]; y_tick=0.2; 
    case 43
        yd=twc_lyman; y_label='liquid water content - Lyman (g kg^{-1})'; y_label2='Lyman'; y_lim=[0 0.1]; y_tick=0.02;
    case 44
        yd=part_conc_2dc; y_label='particle conc. - 2DC (m^{-3})'; y_label2='2DC'; y_lim=[0 3e6]; y_tick=0.5e6;
    case 45
        yd=lwc_2dc; y_label='liquid water content (g m^{-3})'; y_label2='2DC'; y_lim=[0 0.8]; y_tick=0.2;
    case 46
        yd=part_conc_2dp; y_label='particle conc. - 2DP (m^{-3})'; y_label2='2DP'; y_lim=[0 0.5e5]; y_tick=0.1e5;
    case 47
        yd=lwc_2dp; y_label='liquid water content - 2DP (g m^{-3})'; y_label2='2DP'; y_lim=[0 3]; y_tick=0.5;
    case 48
        yd=part_conc_pcasp; y_label='particle conc. - PCASP (m^{-3})'; y_label2='PCASP'; y_lim=[0 2e3]; y_tick=0.5e3;
    case 49
        yd=mass_content_pcasp; y_label='mass content - PCASP (g m^{-3})'; y_label2='PCASP'; y_lim=[0 2.0]; y_tick=0.4;
        
    case 50
        yd=total_scat_blue_neph; y_label='Total Scattering - Neph'; y_label2='Blue'; y_lim=[0 3.5e-4]; y_tick=1e-4;
    case 51
        yd=total_scat_red_neph; y_label='Total Scattering - Neph'; y_label2='Red'; y_lim=[0 3.5e-4]; y_tick=1e-4;
    case 52
        yd=total_scat_green_neph; y_label='Total Scattering - Neph'; y_label2='Green'; y_lim=[0 3.5e-4]; y_tick=1e-4;
        
    case 60
        yd=swd_clear; y_label='SW Radiation (W m^{-2})'; y_label2='SW down'; y_lim=[0 80]; y_tick=20;
        y_lim=[0 450]; y_tick=50; % for profiles
    case 61
        yd=lwd; y_label='LW Radiation (W m^{-2})'; y_label2='LW down'; y_lim=[300 380]; y_tick=20;
        y_lim=[200 400]; y_tick=50; % for profiles
    case 62
        yd=swu_clear; y_label='SW Radiation (W m^{-2})'; y_label2='SW up'; y_lim=[0 80]; y_tick=20;
        y_lim=[0 450]; y_tick=50; % for profiles
    case 63
        yd=lwu; y_label='LW Radiation (W m^{-2})'; y_label2='LW up'; y_lim=[300 380]; y_tick=20;
        y_lim=[200 400]; y_tick=50; % for profiles
        
    case 70
        yd=shf_bulk; y_label='surface heat flux (W m^{-2})'; y_label2='shf bulk'; y_lim=[0 200]; y_tick=50;
    case 71
        yd=lhf_bulk; y_label='surface heat flux (W m^{-2})'; y_label2='lhf bulk'; y_lim=[0 200]; y_tick=50;
    case 72
        yd=tau_bulk; y_label='surface momentum flux (N m^{-2})'; y_label2='\tau bulk'; y_lim=[0 2]; y_tick=0.5;
        
    case 80
        yd=ozone; y_label='ozone (ppb)'; y_label2=''; y_lim=[0 100]; y_tick=5;
        
end;