Dear all,

I have some problems when using seadas batch mode to l2gen the L1B datasets  to L2, I do step by step to L2 from L1A, want to produce all the nlw and rrs products.
1. i tested in Linux fedora10 using seadas 6.2up 2b, using the script attach below. I use command line: chmod a+x l1b_2_l2.csh
                                                                                                                   again command line: l1b_2_l2.csh
all the process work, producing all the products in the one L2 files(all the nlw and rrs and other parameters i want included), switching from SWIR to NIR according to wang and shi turbidity inidex 1.3.

then i copy same L1B datasets to Mac 10.5 seadas 6.2up 2b, I  only change the cut -c 1-15(in linux system in the command script) to cut -c 1-14([/b]in Mac system in the command script, if i do not change, the script cannot work) ,then i ran the command in the bash mode:
                                                                                                                           chmod a+x l1b_2_l2.csh
                                                                                                                           ./l1b_2_l2.csh
however, there is just Rrs and chl_a products , I was very curious almost the script no difference, would produce such different products? what's meaning of cut -c 1-4 and cut -c 1-15? i find here is no one couting difference. how to get all the products as same as in linux?
Do i miss something?

Thanks in advance!"
joe

Your script does not define l2prod for the l2gen program, therefore, l2gen will output the sensor default products
which do NOT include nLw.

Sean

thanks very much, sean!
i was curious that why Linux can produce all the nlw and rrs products, even there is default?
mac just output default rrs using same script.

joe

dear sean,
can you do me some help to check the script from L1A to L2?
previously i can do it from one level to another level, but here i find another script in the forum, but there some error i cannot solved.
I have the L1A.x.hdf file in /home/sarat/sarat/ECS/L1A/
I search some you mentioned that the seadas no longer input L1B_QKM file,  here, i just want to creat the L1B and L2, why  appear message below?

here is the message:
**************************************************************************
Processing                                                                       
Begin processing. Filebasename:                                                  
A2002187041500                                                                   
*** /home/sarat/seadas62/run/bin3/ncdump: ncopen failed on /home/sarat/sarat/ECS/L1A/A2002187041500.L1A_LAC.x.hd����i                                                                  
modis_GEO.csh: ERROR: Unable to determine start time from L1A file.                        

modis_L1A_extract.csh: ERROR: File /home/sarat/sarat/ECS/L1A/A2002187041500.GEO does not exist. Exiting.                                                                               

modis_L1B.csh: ERROR: File /home/sarat/sarat/ECS/L1A/A2002187041500_sub.L1A does not exist.

Usage: modis_L1B.csh MODIS_L1A_file MODIS_GEO_file [OPTIONS]

Options:

  -b                        Use base of input filename for output filenames
                            (without -b, output filenames are automatically determined)
  -o L1B_1KM_filename       Output MODIS L1B 1KM HDF filename                         
  -h L1B_HKM_filename       Output MODIS L1B HKM HDF filename                         
  -q L1B_QKM_filename       Output MODIS L1B QKM HDF filename                         
  -rlut REFL_LUT_filepath   Full path and filename of Reflective LUT file             
  -elut EMIS_LUT_filepath   Full path and filename of Emissive LUT file               
  -qlut QA_LUT_filepath     Full path and filename of QA LUT file                     
  -delete-1km               Delete 1km resolution L1B file                            
  -delete-hkm               Delete 500m resolution L1B file                           
  -delete-qkm               Delete 250m resolution L1B file                           
  -save-log                 Save processing log file(s)                               

Getting ancillary data
Determining pass start and end times..

Input file : /home/sarat/sarat/ECS/L1B/A2002187041500.L1B_QKM
Start time : 2002187041500                                  
End time   : 2002187041959                                  

Created 'A2002187041500.L1B_QKM.anc' l2gen parameter text file:

met1=/home/sarat/seadas62/run/var/anc/2002/187/N200218700_MET_NCEPN_6h.hdf
met2=/home/sarat/seadas62/run/var/anc/2002/187/N200218706_MET_NCEPN_6h.hdf
met3=/home/sarat/seadas62/run/var/anc/2002/187/N200218706_MET_NCEPN_6h.hdf
ozone1=/home/sarat/seadas62/run/var/anc/2002/186/N200218600_O3_TOMSOMI_24h.hdf
ozone2=/home/sarat/seadas62/run/var/anc/2002/187/N200218700_O3_TOMSOMI_24h.hdf
ozone3=/home/sarat/seadas62/run/var/anc/2002/187/N200218700_O3_TOMSOMI_24h.hdf
sstfile=/home/sarat/seadas62/run/var/anc/2002/187/N2002187_SST_OIV2AVAM_24h.nc
icefile=/home/sarat/seadas62/run/var/anc/2002/187/N200218700_SEAICE_NSIDC_24h.hdf

- All optimal ancillary data files were determined and found on local disk. -

Processing  to Level 2..
./l1a_to_l1b_2_l2.sh: line 59: l2gen,ifile=/home/sarat/sarat/ECS/L1B/A2002187041500.L1B_QKM,geofile=/home/sarat/sarat/ECS/L1A/A2002187041500_sub.GEO,ofile1=/home/sarat/sarat/ECS/L2/A2002187041500.L2_LAC: No such file or directory

remove the commas on the l2gen line

Sarat -

We've seen this kind of geolocation failure before.  One workaround is to run from the same directory as the geolocation file - or you could upgrade to SeaDAS version 6.3.  The new python scripts haven't had that problem.

Gwyn

Dear sean, gwyn,
Thanks for your previous answer!
here i came across another problem,

I use the same script from l1A to L2 that attached, i just changed the l2gen as you said and switch NIR and SWIR, all the process is good, but the atmospheric correction is
still using NIR, didnot switch, i was so curious about this, because i defin aer_pot=-9, why still select the NIR? can you help to solve this? looking forward!
**************************************************************************************
l2gen ifile=$L1B1KMFILE geofile=$GEOFILE_sub ofile1=$L2FILE \
#l2prod1='Rrs_412,Rrs_443,Rrs_469,Rrs_488,Rrs_531,Rrs_547,Rrs_555,Rrs_645,Rrs_667,Rrs_678,l2_flags' \
aer_opt=-9 \
ctl_pt_incr=1 \
gas_opt=15 \
aer_wave_short=748 \
aer_wave_long=869 \
aer_wave_short=1240 \
aer_wave_long=2130 \

here is the message:
*******************************************************************
Using  869.0 nm channel for cloud flagging over water.
Using  412.0 nm channel for cloud flagging over land.

Processing scan #     0 (1 of 2030) after      2 seconds

Loading SSES table from /Users/ishizaka_lab/seadas6.2/run/data/modisa/cal/sst_sses_modisa_v5.hdf

Loading SSES table from /Users/ishizaka_lab/seadas6.2/run/data/modisa/cal/sst4_sses_modisa_v5.hdf
Loading SST4 coefficients from /Users/ishizaka_lab/seadas6.2/run/data/modisa/cal/sst4_coef_modisa_v5.dat:
2002001 2002212  0.474  1.030  0.497  1.489  0.000
Loading SST coefficients from /Users/ishizaka_lab/seadas6.2/run/data/modisa/cal/sst_coef_modisa_v5.dat:
2002001 2002212  1.223  0.958  0.140  2.211
2002001 2002212  2.052  0.926  0.126  1.285

Aerosol selection bands 748 and 869
NIR correction enabled.

Hi, even i set define aer_opt=-1 to fix the SWIR, the script still choose NIR only
I was confused by this, looking forwar to hear your advice!

**********************************************
Limiting aerosol models based on RH.

Using Gordon & Wang aerosol model selection
  and NIR correction with up to 10 iterations
Using bands at 748.0 and 869.0 nm for model selection
Extrapolating from 869.0 nm
80 aerosol models: 8 humidities x 10 size fractions
model 0, rh=30.000000, sd=16, alpha=2.157241, name=r30f95v01
model 1, rh=30.000000, sd=17, alpha=2.046257, name=r30f80v01
model 2, rh=30.000000, sd=18, alpha=1.716705, name=r30f50v01
model 3, rh=30.000000, sd=19, alpha=1.337765, name=r30f30v01

It's my mistake , now fix it!