local_planet error???

[alaroy]

Greetings Earth-data-lings!

We’re consistently getting an error when we try to process NPP L1A to L2 data using l2gen.
We’re trying to run a command of the following form:

Code: Select all

l2gen maskland=0 proc_cloud=0  maskcloud=on cloud_thresh=0.1  maskhilt=0 maskstlight=0 resolution=1000 gas_opt=8 l2prod=rhos_NNNN,rhos_MMMM,rhos_PPPP ifile=S20250414154200-L1A-NPP-OBPG.nc geofile=S20250414154200-GEO-NPP-OBPG.nc ofile=S20250414154200-L2-NPP-OBPG.nc.tmp

It always fails with a bunch of errors of the form:

Code: Select all

….
Allocated 4080040 bytes in L1 record.
Allocated 4080040 bytes in L1 record.
Allocated 4080040 bytes in L1 record.
VcstCmnGeo::NOVAS-C function local_planet returned error code.
VcstCmnGeo::Common Geolocation lunar vectorAtSat() failure at scan 0
VcstCmnGeo::NOVAS-C function local_planet returned error code.
VcstCmnGeo::Common Geolocation lunar vectorAtSat() failure at scan 1
VcstCmnGeo::NOVAS-C function local_planet returned error code.
VcstCmnGeo::Common Geolocation lunar vectorAtSat() failure at scan 2
VcstCmnGeo::NOVAS-C function local_planet returned error code.
VcstCmnGeo::Common Geolocation lunar vectorAtSat() failure at scan 3
…

In case it helps we have l2gen version 9.10.0-V2025.0 (Jan 16 2025 23:07:44)

Thank you in advance for any pointers you can give me on this one. Hopefully I don’t have to change my local planet in order for it to work !

-Andrew L.

[alaroy]

I'm bumping this thread (I think it got lost over the holiday weekend)

[OB.DAAC - SeanBailey]

Andrew,

There does seem to be an issue with l2gen calling the VIIRS calibration code. We have never done this in operations and in fact plan to depreceate the ability to process L1A files directly. In production we process the L1A to L1B using calibrate_viirs then use that L1B as input to l2gen. I confirmed that works as expected with V2025.1

Sean

[alaroy]

Thank you Sean!

I tried running the command

Code: Select all

calibrate_viirs ifile=S20250423143600-L1A-NPP-OBPG.nc l1bfile_img=S20250423143600-L1B-NPP-OBPG.nc

Then the usual

Code: Select all

l2gen maskland=0 proc_cloud=0  maskcloud=on cloud_thresh=0.1  maskhilt=0 maskstlight=0 resolution=1000 gas_opt=8 l2prod=rhos_MMMM,rhos_NNNN,rhos_PPP ifile=S20250423143600-L1B-NPP-OBPG.nc geofile=S20250423143600-GEO-NPP-OBPG.nc ofile=/data/NC/L2/FAI/2025-04-23/S20250423143600-FAI-L2-NPP-OBPG.nc.tmp

l2gen failed with the error:

Code: Select all

-E- /home/dshea/ocssw/oel_hdf4/libl1/l1_nc_generic_read.c line 142: nc_inq_ncid failed for file, S20250423143600-L1B-NPP-OBPG.nc  group, geophysical_data.

I'm pretty sure I'm missing something "obvious" here but its not so obvious to me right now. Any hints on what I'm doing wrong, please?

ETA: In case it makes a difference this was with l2gen Version 9.11.0-V2025.1

Thank you in advance,
-Andrew L.

[OB.DAAC - SeanBailey]

Andrew,

It's a wee little thing: l1bfile_mod instead of l1bfile_img
l2gen requires the moderate resolution bands (750m), and cannot run on the imagery resolution (375m). You can output both, if you like, just pass l2gen the MOD file.

Sean

[alaroy]

That worked great! I figured it was something obvious to someone who is more of an expert than me.

Thank you Sean!!!!!