Modis Chlorophyll product command line creation steps

Modis Chlorophyll product command line creation steps

Command line is usually quicker to perform, and avoid any GUI issues I am having in SeaDAS 7.0. I corrupt window default setting with ease. Once I have repeatable command line steps, automation scripts are simple to create.

1. identify clear days via http://oceancurrent.imos.org.au/ (Australia)
2. Download from Ocean Color Level 1 and 2 Browser
1. if area  available in a single swarth download the L1A_LAC file
2. otherwise download L0_LAC file.
3. For L0_LAC files
• Concatenate
• cat A2011098040500.L0_LAC A2011098041000.L0_LAC > A2011098.L0_LAC
• result is A2011098.L0_LAC
• Create L1A_LAC
• modis_L1A.py A2011098.L0_LAC
• result is A2011098040500.L1A_LAC
4. Create GEO file for L1A file
• modis_GEO.py A2011098040500.L1A_LAC
• If input product result of concatenation consider if -r option is required.
• result
• A2011098040500.GEO
5. Extract area slightly large than area of interest. Further subsection performed later.
• modis_L1A_extract.py A2011098040500.L1A_LAC  -n -40 -s -45 -w 142 -e 147 --verbose --log
• results include
• A2011098045000.L1A_LAC.sub
• A2011098045000.GEO.sub
6. Create L1B file for area of interest
• modis_L1B.py A2011098sub.L1A_LAC 
• I renamed my files 
• A2011098sub.L1A_LAC
• A2011098sub.GEO
• Results include
• A2011098sub.L1B_HKM
• A2011098sub.L1B_LAC
• A2011098sub.L1B_QKM
7. Adjust the names in my l2gen par files

for example my gsm.par

ifile=A2011098sub.L1B_LAC
geofile=A2011098sub.GEO
ofile=gsm.hdf
l2prod="adg_443_gsm,bbp_443_gsm,chl_gsm"

      8. Run l2gen par=gsm.par

Only need 3 dates at the moment so using SeaDAS 7.0 GUI and ERMapper.