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_seadasYou 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_seadasDid 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.logHOWEVER, 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.