VIIRS SST Processing

[alakes]

Hi,
    I am trying to generate VIIRS SST product for Indian Ocean region.
used the following command to generate L2 files.
/bin/bash -c "l2gen ifile=$L1AFile geofile=$GEOFile suite=SST ofile=$L2File7 l2prod='sst,bias_sst,flags_sst,qual_sst,sstref,stdv_sst,l2_flags'" >$LogFile4

The Problem is, in L2 file qual_sst having value 3 in ocean and 4 in land for all over the Indian Ocean region.
Have a look on one of the LAC.
qual_sst
SST

As the qual_sst is always >=3 here so I am getting a l2bin error, total_filled_bins: 0  (as no bin is valid for the following command)

/bin/bash -c "l2bin  INFILE=SST_FileList.lst OFILE=Test_SST.nc RESOLVE=1 PRODTYPE=regional NOEXT=1 QUAL_PROD=qual_sst  NIGHT=0 FLAGUSE='LAND,HISOLZEN'  QUAL_MAX=2  L3BPROD='sst;-2.00000'"

Any suggestion why qual_sst is behaving like that. Or am I making any mistake?

Thanks in advance..

[OB.DAACx - SeanBailey]

Since you did not provide a reference SST (sstfile=), the code used an sst climatology and the algorithm explicitily categorizes data that use climatology for the reference as quality 3.

As you're not asking for anything that is not provided in the standard product, you could simply grab the one we produce... but if you want to process locally, call getanc.py first and provide the process with the NRT ancillary data (which will include the correct sstfile)

Sean

[alakes]

I am using the following code but same result :confused::confused:

$SeaDASPath/run/scripts/getanc.py $L1AFile --no2 --verbose

/bin/bash -c "l2gen ifile=$L1AFile geofile=$GEOFile suite=OC maskglint=on \
ofile=$L2File1 l2prod='chlor_a,Kd_490,pic,poc,aot_862,angstrom,par,ipar,l2_flags'"

/bin/bash -c "l2gen ifile=$L1AFile geofile=$GEOFile suite=SST \
ofile=$L2FileSST l2prod='sst,bias_sst,flags_sst,qual_sst,sstref,stdv_sst,l2_flags'"

as qual_sst is not working I am using flags_sst as follows:
l2bin  infile=$inFile.txt ofile=$l2binFile.$fileExtension sday=$Yyear$Yday eday=$Yyear$Yday resolve=$binRes prodtype='regional' noext=1 \
qual_prod=flags_sst  qual_max=0  night=0 flaguse='LAND,HISOLZEN'  l3bprod=${l3bProduct[$j]} oformat=$GACopFormat verbose=1 > $l2binFile.log

There is some problem in the overlapping region, I mean It can be seen clearly the overlap region of 2 LAC files referring the green region Northern Arabian Sea in the attached figure.
V-Mar2017-d31-4KM-IndiaSriLanka-SST.png

Is it due to the the flaguse?  But the binfile generated by OBPG using the same flag flaguse='LAND,HISOLZEN'..

[OB.DAACx - SeanBailey]

Running getanc.py will both print to the terminal and generate a <filename>.anc file with the ancillary parameters necessary for processing.
You need to either add the parameters to your own par file (or command line) or pass par=<filename>.anc (in your case $L1AFile.anc).

[liamgumley]

Hi,

I really like the SST image you generated named 'V-Mar2017-d31-4KM-IndiaSriLanka-SST.png'.
Would you mind sharing the recipe you used to create it?

Thanks,
Liam.