%% Generates example used hytgrid6ck.m
%%

close all
clear all

kp=3;

nop=50;                  %Number of points to break field line into.  
notp=48;                 %Number of magnetic local time (MLT) intervals to be used.
tint=1;                  %UT time iteration step size to be used (hours)

re=6371;                                               %Earth radi (km)
lshmx=9;                                               %maximum L-shell    
lshcp=6.8;                                             %maximum L-shell for calculation of profiles 
lshv=[(re+400)/re 1.1 1.15 1.2:0.1:6.8 7:0.2:lshmx];    %required L-shell vector  

noti=24/tint;                                          %number of time increments to be used
timmat=0:tint:(24-tint);                               %UT time vector
mltv=(0:(24/notp):(24-(24/notp)));                      %magnetic local time vector

load plsimeq1

mlti=0;
hytot=ones(length(lshv),length(mltv));

hytgrid6ck(mlti,hytot,lshv,mltv,kp,lsht,mltt,ratt,kpt)
