clc; close all; clear all; format long; axis square; tic;

fid = fopen('RadSurvey.2017_2_2_8_44_23_IN.dat');
data = textscan(fid, '%f %f %f');
fclose(fid);
r = data{:,1}; theta = data{:,2}; rad = data{:,3};
r = r./12;
%theta = theta - pi./6;
theta = theta.*pi./180;
time = 0.5*60./rad; %time till dose exceeds daily guideline
%rad = log(rad);
%rad = ln(rad);
[x,y] = pol2cart(theta,r);
[X,Y] = meshgrid(linspace(min(r),max(r),length(r)));
N = 100;
[x1,y1] = meshgrid(linspace(min(x),max(x),N),linspace(min(y),max(y),N));

vq = griddata(x,y,rad,x1,y1);
vt = griddata(x,y,time,x1,y1);

ax = gca;
ax.FontSize = 14;

subplot(2,2,1)
[Crad,hrad]=contour(x1,y1,vq);%,N2)
c = colorbar;
c.Label.String = 'Radiation Field (mSv/h)';
c.FontSize = 14;
colormap jet
hrad.LevelStep = 0.025;
axis square;
hrad.Fill = 'on';
hrad.LevelList = [min(rad):0.1:4];
ax = gca;
ax.FontSize = 14;
ax.XTick = -25:5:25;
ax.YTick = -25:5:25;
ylabel('Feet from Center (ft.)');
xlabel('Feet from Center (ft.)');
title('Dataset Cut >= 4mSv/h to Show low-field Contours')

subplot(2,2,2)
axis square
hold on
[Ctime,htime]=contour(x1,y1,vt);%,N2)
c = colorbar;
c.Label.String = 'Time to Exceed Daily Guideline (min.)';
c.FontSize = 14;
colormap jet
htime.LevelStep = 0.5;
set(gca,'ydir','normal');axis square;
htime.Fill = 'on';
htime.LevelList = [0.5, 1, 1.5, 2, 5, 10, 15, 20,30, 40, 60];
ax = gca;
ax.FontSize = 14;
ax.XTick = -25:5:25;
ax.YTick = -25:5:25;
ylabel('Feet from Center (ft.)');
xlabel('Feet from Center (ft.)');
title('Dataset Cut >= 4mSv/h to Show low-field Contours')

subplot(2,2,3)
[Crad,hrad]=contour(x1,y1,vq);%,N2)
c = colorbar;
c.Label.String = 'Radiation Field (mSv/h)';
c.FontSize = 14;
colormap jet
hrad.LevelStep = 0.5;
axis square;
hrad.Fill = 'on';
hrad.LevelList = [4:0.1:max(rad)];
ax = gca;
ax.FontSize = 14;
ax.XTick = -25:5:25;
ax.YTick = -25:5:25;
ylabel('Feet from Center (ft.)');
xlabel('Feet from Center (ft.)');
title('Dataset Cut <= 4mSv/h to Show high-field Contours')

subplot(2,2,4)
axis square
hold on
[Ctime,htime]=contour(x1,y1,vt);%,N2)
c = colorbar;
c.Label.String = 'Time to Exceed Daily Guideline (min.)';
c.FontSize = 14;
colormap jet
htime.LevelStep = 0.5;
set(gca,'ydir','normal');axis square;
htime.Fill = 'on';
htime.LevelList = 1:10;
ax = gca;
ax.FontSize = 14;
ax.XTick = -25:5:25;
ax.YTick = -25:5:25;
ylabel('Feet from Center (ft.)');
xlabel('Feet from Center (ft.)');
title('Dataset Cut <= 4mSv/h to Show high-field Contours')

toc;