I am batch processing some MODIS Aqua data (L1A to L3M, Rrs with the basic atmospheric correction) and I would like to include the szen and l2_flags.
I have 2 issues :
1) When I add szen, nothing happens and an error is returned saying szen does not exist. "-E- /home/seadas/ocssw/src/l2gen/l2_generic.c line 458: product szen not found." I can see different szen related to each waveband with the l2gen products window gui, but is there any general szen I can access to?
2) When I add l2_flags (note I cannot find it from the l2gen products window using the gui), I can see it and play with it using the L2 file created, but I cannot see it inside the L3M file at the end, only the Rrs and coordinates are created. See below an example of ncdump -h randomDay.L3m_DAY_OC with one Rrs processed:
I also did not ask for palette/eightbitcolor, but it is still returned, however I have no idea what it is.variables:
short Rrs_859(lat, lon) ;
Rrs_859:long_name = "Remote sensing reflectance at 859 nm" ;
Rrs_859:scale_factor = 2.e-06f ;
Rrs_859:add_offset = 0.05f ;
Rrs_859:units = "sr^-1" ;
Rrs_859:standard_name = "surface_ratio_of_upwelling_radiance_emerging_from_sea_water_to_downwelling_radiative_flux_in_air" ;
Rrs_859:_FillValue = -32767s ;
Rrs_859:valid_min = -30000s ;
Rrs_859:valid_max = 25000s ;
Rrs_859:display_scale = "linear" ;
Rrs_859:display_min = 0.f ;
Rrs_859:display_max = 0.015f ;
float lat(lat) ;
lat:long_name = "Latitude" ;
lat:units = "degrees_north" ;
lat:standard_name = "latitude" ;
lat:_FillValue = -999.f ;
lat:valid_min = -90.f ;
lat:valid_max = 90.f ;
float lon(lon) ;
lon:long_name = "Longitude" ;
lon:units = "degrees_east" ;
lon:standard_name = "longitude" ;
lon:_FillValue = -999.f ;
lon:valid_min = -180.f ;
lon:valid_max = 180.f ;
ubyte palette(rgb, eightbitcolor) ;
Here are the lines I ran, it worked when I processed Rrs/Rhos/Rhot in the past but it does not work for szen/l2_flags:
Code: Select all
##########################################################################
################ START THE PROCESSING FROM L1A to L3m ####################
##########################################################################
echo "Start processing from L1A to L3m.."
cd $fullpath # Go to the folder where you need to process
################ START THE PROCESSING FROM L1A to GEO #####################
# Process the L1A subscene files to GEO
echo "modis_GEO : Processing L1A to GEO.."
modis_GEO.py $FILE -o $GEOFILE
echo "modis_GEO Done!"
################ START THE PROCESSING FROM L1A to L1B #####################
echo "modis_L1B : Processing Level 1A to Level 1B.."
# Process the L1A/GEO subscene files to L1B
modis_L1B.py $FILE $GEOFILE -o $L1BFILE --del-hkm --del-qkm
# Determine ancillary data
# Create a file in l2gen's par file format called L1BFILE.anc
getanc.py $L1BFILE
echo "modis_L1B Done!"
################ START PROCESSING FROM L1B to L2A #####################
echo "l2gen : Processing $L1BFILE to Level 2.."
l2gen ifile=$L1BFILE geofile=$GEOFILE ofile=$L2FILE \ l2prod1='Rrs_859,l2_flags' \
echo "l2gen Done!"
################ START PROCESSING FROM L2A to L3BIN #####################
echo "l2bin : Processing $L3BINFILE to Level 3 binned.."
l2bin ifile=$L2FILE ofile=$L3BINFILE \
resolve=1 \
flaguse=ATMFAIL \
echo "l2bin Done!"
################ START PROCESSING FROM L3BIN to L3MAP #####################
echo "l3mapgen : Processing $L3BINFILE to Level 3 mapped.."
l3mapgen ifile=$L3BINFILE ofile=$L3MAPFILE \
resolution=1km \
interp=area \
north=65 \
south=48 \
east=13 \
west=-25 \
fudge=3 \
Am I missing a detail ?
Thanks