Page 1 of 1
local_planet error???
Posted: Wed Apr 16, 2025 1:05 pm America/New_York
by 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.
Re: local_planet error???
Posted: Mon Apr 21, 2025 11:13 am America/New_York
by alaroy
I'm bumping this thread (I think it got lost over the holiday weekend)
Re: local_planet error???
Posted: Tue Apr 22, 2025 2:31 pm America/New_York
by 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
Re: local_planet error???
Posted: Fri Apr 25, 2025 2:33 pm America/New_York
by 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.
Re: local_planet error???
Posted: Fri Apr 25, 2025 2:53 pm America/New_York
by 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
Re: local_planet error???
Posted: Fri Apr 25, 2025 3:04 pm America/New_York
by alaroy
That worked great! I figured it was something obvious to someone who is more of an expert than me.
Thank you Sean!!!!!