% program to make line plots of new particle formation events
clear all;format long g;

% select station
stationname='Hyytiala';

% Neccessary input files for the line maps
% 1. MNC files (start/end time of formation at lowest diameter of dmps spectra)
%    Put the MNC files in subdirectory  called 'mncdata'
% 2. EOG files (End of growth rate files)
%    Put the EOG files in subdirectory  called 'eogdata'
% 3. dat (trajectory files)
%   Put the dat files in subdirectory called 'trajdata'

% Select period
period=[datenum(2005,5,11,0,0,0):datenum(2005,5,14,0,0,0)];

% setting the time-shift parameter (time for growth between 1.5 nm diameter
% and the lowest size bin of the instrument)
timeshift=0.75; % h

% setting the position of the station (decimal fraction, not actual
% minutes)
poslat=61.850;%Latitude of measurement station
poslon=24.283; %longitude of measurement station

% setting the latitude and longitude limits of the map projection area
latlim=[54 72]; % latitude
lonlim=[0 40]; % longitude

% setting the resolution of the gridded map
resolution_lat=0.4; % degrees in latitude direction
resolution_lon=0.8; % degrees in longitude direction

% select the time shift of the GMT time compared to your station (if your
% station is to the east of the Greenwich meridian, the timeshift should be
% negative).
GMT_shift=-2; % select the GMT shift in hours

% loading data from this period
number_map_plotting_loaddata;

% plotting line maps
number_map_plotting_lineplot

% calculating the NPF counts in grid boxes 0.1 grades in size
number_map_plotting_gridcount
