Nares Strait Modis Processing as of July-14, 2009
(MacPro Powerbook, OS X 10.5.7)
<ul>
<li>1. Browse oceancolor.gsfc.nasa using:
  <ul>
  <li>1.1 Region [80,83]N latitude, [72,54]W longitude require 100% coverage within swath/scenes
  <li>1.2 Download list of filenames
  <li>1.3 Edit this file for desired scenes (list.txt containing file names $name)
  </ul>

<li>2. Process with <a href="../snipets/step-0.txt">step-0.csh</a> to get ascii-file for band-1
  <ul>
  <li>2.1 Using list.txt to generate step-1.csh 
  <li>2.2 Execute step-1.csh that calls step-2.csh with $name as argument
  <li>2.3 Execute <a href="../snipets/step-2.txt">step-2.csh</a> consisting of
    <ul>
    <li>2.3.1 Download Level-0 (LO_LAC.bz2) file from oceancolor.gsfc.nasa
    <li>2.3.2 Unzip L0_LAC file
    <li>2.3.3 Process Level-0 through Level-2 with seabatch1.1 script (requires SeaDas5.4)
    <li>2.3.4 Extract 250-m resolution band-1 from Level-1B (L1B_QKM) and geolocation (.GEO) files
      <ul>
      <li>2.3.4.1 Create command script .cmd for embedded IDL processing
        <ul>
        <li>2.3.4.1.1 "load" .hdf file into memory
        <li>2.3.4.1.2 "out" binary bin1, bin2, and nav files to current directory
        </ul>
      <li>2.3.4.2 Execute command script by invoking "seadas -em test.cmd >log" (requires SeaDas5.4)
      </ul>
    <li>2.3.5 Create control file log.dat from hdfdump -h filename.L1B_QKM
    <li>2.3.6 Run fortran-77 <a href="../snipets/getbin.txt">getbin.f</a> to read binaries, apply scaling coefficients, sub-sample into an area of interest using log.dat as control, and produce 1-line summary
     <li>2.3.7 Clean up files
        <ul>
        <li>2.3.8.1 rename sub-sampled ascii file as .dat
        <li>2.3.8.2 move .LO_LAC, .L1B_QKM and .GEO files to archive directory
        <li>2.3.8.3 remove all other files
        </ul>
      </ul>
    </ul>
<li>3. Create regular grid for display and physical analyses
  <ul>
  <li>3.1 Using log.all to generate gpath1.csh
  <li>3.2 Execute gpath1.csh that calls GMT script file gpath.g (or gpathmed.g for larger area)
    <ul>
    <li>3.2.1 Filter .dat ascii data with "blockmean" (GMT command) with 0.01/0.002 long/lat width
    <li>3.2.2 Generate .grd gridded data with "surface" (GMT commant) with same resolution as 3.2.1
    <li>3.2.3 Color contour gridded field with "grdimage" (GMT command)
    <li>3.2.4 Create .grd, .ps, and .pdf files
    <li>3.2.5 Create .track ascii file of properties along a specified track line
    <li>3.2.6 Create and save .track ascii file of properties along a specified track line
    </ul>
  <li>3.3 Create .gif movie from saved .pdf files using <a href="../images/movie.txt">movie.csh</a>
  <li>3.4 Create <a href="../images/Modis2009/test.html">Image Link Table</a>
  </ul>
<li>4. Analyze data along track lines as data that varies along track and time
   <ul>
   <li>4.1 Read sequence of binary gridded 1801x1501 images at 2-byte/pixel (3.2.4)
   <li>4.2 Calculate <a href="../images/GridCor.jpg">correlation between 2 successive images</a> as a function of spatial lag
      </ul>
    </ul>
</ul>
<ul>
<li>Software needs:
  <ul>
  <li>SeaDas5.4 (from oceancolor.gsfc.nasa)
  <li>SeaBatch1.1 (from seabatch.com)
  <li>g77 (from fink)
  <li>GMT (from fink)
  <li>ImageMagick (from fink)
  <li>wget (from fink)
  <li>hdf (from fink)
  <li>X11 (from apple)
  </ul>
</ul>