%  DoPrintonMaps.m  - create a series of plots with values for a given
%                     density or depth written on a map
%                   Saves figures to postscript files. 
%    needs variables set:  GEOMINFO CRUISEINFO DEFAULTS
%    sets variables: CruiseName Line FixGrid DepthResolution
%                    xrange yrange UseDepth Depth2Plot Dens2Plot cl
%    calls: SetupPlots printTonMap printSonMap printO2onMap printSigonMap

MakePlot=1;
SkipPlot=0;

%SetDefaults;

DepthResolution='OneMeter';
CRUISEINFO=0;
Line='grid';
FixGrid=0;
%CruiseName='NBP0103';
CruiseName='NBP0104';
SetupPlots;
%                    Whole area
xrange=[0 550];
yrange=[-50 300];
%                 Marguerite Bay and Adelaide Island
%xrange=[150 500];
%yrange=[-75 175];

if SkipPlot
%if MakePlot
  RefPress=300;
  printDynTopoMap;
%  print('-dpsc2',['printDynTopoMap.' CruiseName '.' num2str(RefPress) '.ps'])
end

%if SkipPlot
if MakePlot
  UseDepth=1;
%  Depth2Plot=300;
  Depth2Plot=50;
  printTonMap;
  print('-dps2',['printTonMap.' CruiseName '.' num2str(Depth2Plot) '.ps'])
end

%if SkipPlot
if MakePlot
  Depth2Plot=50;
  printSonMap;
  print('-dps2',['printSonMap.' CruiseName '.' num2str(Depth2Plot) '.ps'])
end

if SkipPlot
%if MakePlot
  Depth2Plot=300;
  printO2onMap;
%  print('-dpsc2',['printO2onMap.' CruiseName '.' num2str(Depth2Plot) '.ps'])
end

if SkipPlot
%if MakePlot
  Depth2Plot=300;
  printSigonMap;
%  print('-dpsc2',['printSigonMap.' CruiseName '.' num2str(Depth2Plot) '.ps'])
end




