Ocean Color Forum - Not logged in
Forum Ocean Color Home Help Search Login
Previous Next Up Topic SeaDAS / SeaDAS 6.x - General Questions / MODIS l1b to l2 processing problem (locked) (1760 hits)
By Sarat76 Date 2011-07-28 08:39 Edited 2011-07-28 08:44
Hi,
While processing MODIS level1b to level2 by using script in the command line, I could only create the following files!

A2006146044000.stdproc.L2
A2006146044000.stdproc.L2.par

However I am not able to create the following L2_LAC files!

A2006146044000.L2_LAC
A2006146044000.L2_LAC.par

I could create the
A2006146044000.L2_LAC
A2006146044000.L2_LAC.par files by interactive mode (seadas main menu) but not through using the script in command line.

Can anyone figure out the problem.Thank you

Sarat
By gwyn Date 2011-07-28 14:14
If this is in reference to the script posted in another thread, it's because you never call l2gen with ofile1='A2006146044000.L2_LAC'.
By gwyn Date 2011-07-29 14:03
Sarat -

You set the output filename as:
set name3=`echo $name|cut -c 1-14`.stdproc.L2``

Then use it in the call:
echo "l2gen,ifile='$dir/$name',geofile='$dir/$name2',ofile1='$dir/$name3',slat=32.25,elat=34.25,slon=130.75,elon=132.75,ctl_pt_incr=1,resolution=-1,/wait" >> $dir/load_seadas

You never set an output filename to be A2006146044000.L2_LAC.

Maybe that's what you wanted for this call, commented out by the initial "#":
#echo "l2gen,ifile='$dir/$name',geofile='$dir/$name2',ofile1='$dir/$name3',slat=34,elat=36,slon=138,elon=140,ctl_pt_incr=1,gas_opt=15,resolution=250,aer_opt=-9,aer_wave_short=1240,aer_wave_long=2130,wait=wait" >> $dir/load_seadas

Did you try the script I posted in the other thread?

Here's a version that should make both products:

cd /home/sarat/sarat/
rm -f load_seadas; touch load_seadas
set options1 = 'slat=32.25,elat=34.25,slon=130.75,elon=132.75,ctl_pt_incr=1,resolution=-1'
set options2 = 'slat=34,elat=36,slon=138,elon=140,ctl_pt_incr=1,gas_opt=15,resolution=250,aer_opt=-9,aer_wave_short=1240,aer_wave_long=2130'

foreach l1bfile (*.L1B_LAC)
  echo $l1bfile
  set gran = `echo $l1bfile | cut -d\. -f1`
  set geofile = $gran.GEO

  set l2file  = $gran.stdproc.L2
  set cmd = "l2gen,ifile='$l1bfile',geofile='$geofile',ofile1='$l2file',$options1,/wait"
  echo $cmd >> load_seadas

  set l2file  = $gran.L2_LAC
  set cmd = "l2gen,ifile='$l1bfile',geofile='$geofile',ofile1='$l2file',$options2,/wait"
  echo $cmd >> load_seadas

end
echo exit >> load_seadas
cat load_seadas
seadas -em -b load_seadas >&! seadas.log


HOWEVER,  according to the documentation your aer options aren't quite right:
   aer_opt (int) (default=99) = aerosol mode option
       -9: Multi-scattering with 2-band model selection using Wang & Shi
           turbidity index (1.30) to switch between SWIR and NIR. (MODIS only,
           requires aer_swir_short, aer_swir_long, aer_wave_short, aer_wave_long)

Did you mean to use aer_opt=-1?  Or to specify aer_swir_short and aer_swir_long?
Here's an example posted earlier.
By Sarat76 Date 2011-08-09 00:00
Thank you so much Gwyn.
Sorry to get back to you late.
Previous Next Up Topic SeaDAS / SeaDAS 6.x - General Questions / MODIS l1b to l2 processing problem (locked) (1760 hits)



Responsible NASA Official: Gene C. Feldman
Curator: OceanColor Webmaster
Authorized by: Gene C. Feldman
Updated: 27 November 2007
Privacy Policy and Important Notices NASA logo