Reprojection of multiple of SST4 images

[arnabp893]

I want to reproject multiple SST4 L2 image to projected SST4 L2 data with WGS84 mercator projection.

I have this .par file given below,
[main]
ifile=/Volumes/My_Passport/ifile_text/new2023L1Asst.txt
odir=/Volumes/My_Passport
use_ancillary=1
overwrite=1

resolution=1km
projection=Mercator
product=SST4



# ANCILLARY INPUTS Default = climatology (select 'Get Ancillary' to download ancillary files)

The error is
(base) Euricos-iMac:~ gulfcarbon$ /Users/gulfcarbon/SeaDAS/ocssw/bin/mapgen par=/Users/gulfcarbon/Downloads/MODIS_L1A/mlp_aqua_l2gen_rep.par
usage: mapgen [-h] [--parfile PARFILE] [--ifile IFILE] [--geofile GEOFILE]
[--ofile OFILE] [--logfile LOGFILE] [--use_rgb]
[--product PRODUCT | --product_rgb PRODUCT_RGB]
[--resolution RESOLUTION] [--oformat {netcdf4,png,ppm,tiff}]
[--use_transparency] [--north NORTH] [--south SOUTH]
[--east EAST] [--west WEST] [--projection PROJECTION]
[--central_meridian CENTRAL_MERIDIAN] [--palfile PALFILE]
[--fudge FUDGE] [--datamin DATAMIN] [--datamax DATAMAX]
[--scale_type {linear,log,arctan}] [--threshold THRESHOLD]
[--trimNSEW] [--write_projtext] [--keep-intermediates]
[--verbose]
mapgen: error: unrecognized arguments: par=/Users/gulfcarbon/Downloads/MODIS_L1A/mlp_aqua_l2gen_rep.par



the new2023L1Asst.txt consists file names with path.

Any help will be appreciated. Thank you.

[OB ODPS - towens]

From the usage message, mapgen is expecting --parfile PARFILE or -p PARFILE. It is not handling par=PARFILE.

Tommy

[arnabp893]

I did and the result is

(base) Euricos-iMac:ocssw gulfcarbon$ /Users/gulfcarbon/SeaDAS/ocssw/bin/mapgen -p /Users/gulfcarbon/Downloads/MODIS_L1A/mlp_aqua_l2gen_rep.par
--Error-- : Product not found = sst, L2 file = /Volumes/Extreme_SSD/SST_MODIS_4u/2023/AQUA_MODIS.20230102T080001.L2.SST4.nc
Exiting. See /Users/dshea/ocssw/src/libl2/expand3D.cpp:327

Exiting: /Users/gulfcarbon/SeaDAS/ocssw/bin/l2bin suite=SST flaguse=LAND,BOWTIEDEL ifile=/Volumes/My_Passport/ifile_text/new2023L1Asst.txt ofile=mapgen_tmp.l2bin.nc area_weighting=1 prodtype=regional resolution=H returned a status 1

I changed the par file variables from sst4 to sst

but problem remained

[OB SeaDAS - xuanyang02]

I used multilevel_processor to bin two sst4 files together and then projected to a nc file.

Code: Select all

multilevel_processor par=/Users/*/Scenes/MODIS_AQUA/forum/multilevel_processor-ParFile.par

Here is the content of the multilevel_processor-ParFile.par file -
[main]
ifile=/Users/*/Scenes/MODIS_AQUA/forum/l2_sst_files.txt
odir=/Users/*/Scenes/MODIS_AQUA/forum
overwrite=1

[l2bin]
l3bprod=sst4

[l3mapgen]
product=sst4
projection=mercator
resolution=1km

l2_sst_files.txt is a text file with the content -
/Users/*/Scenes/MODIS_AQUA/forum/AQUA_MODIS.20101010T181508.L2.SST4.nc
/Users/*/Scenes/MODIS_AQUA/forum/AQUA_MODIS.20101010T182008.L2.SST4.nc

[arnabp893]

Thank you
It solved the problem. Also I had to add flaguse=NONE under [l2bin] otherwise it was showing error as,

(base) Euricos-iMac:~ gulfcarbon$ /Users/gulfcarbon/SeaDAS/ocssw/bin/multilevel_processor par=myfolder/mlp_aqua_l2gen_rep.par
Running l2bin on file /Volumes/Extreme_SSD/SST_MODIS_4u/2023/AQUA_MODIS.20230123T081000.L2.SST4.nc.
-I- The l2bin program produced a bin file with no data. No further processing will be done.

[OB SeaDAS - xuanyang02]

Thanks for the tip on flaguse= NONE for l2bin

[arnabp893]

Although I am a little confused.

In the par file .txt I took one file as sample and it worked and I got the file name as filename.L3m.DAY.1km.nc.
but when I started with all the images of a single year (>100) it gave me filename(range of time).L3m.CU.1km.nc
Is it making annual average climatological composite of the available SST4 images?
can I create monthly mean sst4 images if I only put images of one month?


Thank you

[OB SeaDAS - xuanyang02]

When you do l2bin, the output file names will be DAY, 8D(means there is a 7 day difference between the start day and stop day), CU(custom day range), such as L3b.DAY or L3b.CU, or whatever name you specify.

When you do l3mapgen on the L3b file, the file name will be L3m.DAY, L3m.8D or L3m.CU accordingly

multilevel_processor just does l2bin and l3mapgen according to your par file. You can certainly place a month's file names in the txt file, and multilevel_processor can do l2bin and l3mapgen for you

[arnabp893]

Thank you for the clarification.