HICO L-2 Processing

Use this Forum to find information on, or ask a question about, NASA Earth Science data.
avmehta
Posts: 184
Joined: Mon Feb 03, 2020 10:27 am America/New_York
Answers: 0
Been thanked: 1 time

HICO L-2 Processing

by avmehta » Wed Jul 08, 2020 6:10 pm America/New_York

Hello,
   I am using l2gen to get HICO L2 data.  My question is about 'Derived Geophysical Parameters'. Are all the options valid for HICO?
Also, I see multiple algorithm options for certain parameters -- e.g. Kd 490, sst, Zsd (Sechhi Depth).  Is there any default recommended?
I also  have the same question for OLCI, MSI, and OLI! So far I have calculated cholr_a for all of them but wondering if all the geophysical parameters are derived for these sensors.
Also,
When I try to run l2gen for HICO, I get the following message:
  " Reading geometry limiting angles of ATREM
     ATREM Angle_Limit File=/pathname/atren_angle)limit.nc
     ncopen: filename abscf_gas.nc: No such file or directory"

Thank you.
Amita

Tags:

OB.DAAC - SeanBailey
User Services
User Services
Posts: 1464
Joined: Wed Sep 18, 2019 6:15 pm America/New_York
Answers: 1
Been thanked: 5 times

HICO L-2 Processing

by OB.DAAC - SeanBailey » Thu Jul 09, 2020 10:29 am America/New_York

Amita,

Most algorithms can be computed with HICO...but not all.  SST for example.  That requires  thermal IR bands that HICO does not have.  If you ask l2gen to create a product that the sensor cannot support, it'll let you know :grin:

There are indeed multiple algorithms for a number of products.  There are a set of algorithms that we deem "default", and if you run l2gen without setting the products to generate, the defaults will be made.  These are listed in the $OCDATAROOT/<sensor>[/platform]/msl12_defaults.par file (or in the suite default file defined in that file, usually OC, so msl12_defaults_OC.par).  For HICO, they are: l2prod=aot_868 angstrom Rrs_vvv chlor_a chl_ocx Kd_490 pic poc

If it's not in that list, the choice is up to the user as to which to produce.

As for the ncopen: filename abscf_gas.nc: No such file or directory" issue.  That file is under $OCDATAROOT/hico/.  Make sure you've installed the HICO data tables.  You can do so via:
install_ocssw.py -b v7.5 --hico

However, it seems the code in current SeaDAS release expects the file to be in the working directory. 
We'll see about fixing that in the next release.

Until then, you can make it work by putting a symbolic link to the file, e.g.:

$ cd <directory where l2gen is to be run>
$ ln -s $OCDATAROOT/hico/abscf_gas.nc .


Sean

avmehta
Posts: 184
Joined: Mon Feb 03, 2020 10:27 am America/New_York
Answers: 0
Been thanked: 1 time

HICO L-2 Processing

by avmehta » Thu Jul 09, 2020 4:29 pm America/New_York

Thanks much Sean.
I am  running l2gen for HICO  on  SeaDAS GUI. I see the absc_gas.nc file on my Mac in /seadas-7.5.3/ocssw/share/hico directory.
I moved the file to /seadas-7.5.3/ocssw/share/common where atrem_angle_limit.nc is and just got read before I get the 'No such file or directory' message for the absc_gas.nc file. Still l2gen can't find it! Not sure where it is looking for the file.

Thanks again.
Amita

OB.DAAC - SeanBailey
User Services
User Services
Posts: 1464
Joined: Wed Sep 18, 2019 6:15 pm America/New_York
Answers: 1
Been thanked: 5 times

HICO L-2 Processing

by OB.DAAC - SeanBailey » Fri Jul 10, 2020 7:41 am America/New_York

In the working directory...the one from which you launch l2gen.  No need to move it, just use a symbolic link.

avmehta
Posts: 184
Joined: Mon Feb 03, 2020 10:27 am America/New_York
Answers: 0
Been thanked: 1 time

HICO L-2 Processing

by avmehta » Wed Sep 09, 2020 2:24 pm America/New_York

Hi Sean,
     Getting back to HICO processing on GUI again after weeks! I tried what you suggested -- I found l2g3n in /ocssw/bin/. In this directory I used
      ln -s $OCDATAROOT/hico/abscf_gas.nc .
      Then I opened the GUI and tried to run l2gen with a HICO file. But I still get the same error:
    
      ncopen: filename abscf_gas.nc: No such file or directory

  Not sure how to resolve this!
  Thank you.
  Amita

jomoga
Posts: 1
Joined: Fri Jan 07, 2005 4:14 pm America/New_York
Answers: 0

HICO L-2 Processing

by jomoga » Wed Sep 09, 2020 3:31 pm America/New_York

Amita,

I have looked at the code that opens this file and it is looking for it in the $OCDATAROOT/common directory.  Try moving the file there.

Joel

avmehta
Posts: 184
Joined: Mon Feb 03, 2020 10:27 am America/New_York
Answers: 0
Been thanked: 1 time

HICO L-2 Processing

by avmehta » Wed Sep 09, 2020 3:53 pm America/New_York

Hi Joel,
  That is the first place I had copied this file to (/seadas-7.5.3/ocssw/share/common)! Still I get the same error.
Thanks.
Amita

gnwiii
Posts: 713
Joined: Fri Jan 29, 2021 5:51 pm America/New_York
Answers: 2
Has thanked: 1 time

HICO L-2 Processing

by gnwiii » Wed Sep 09, 2020 7:44 pm America/New_York

The "working directory" is generally not the "bin" directory that contains l2gen.  When you run l2gen using the GUI, the environment is set using the $OCSSWROOT/scripts/ocssw_runner.sh.  You can edit this script to add a line such as echo "0 -- working directory is $PWD" > ~/ocssw_runner_pwd.text just before the final line ($*).

avmehta
Posts: 184
Joined: Mon Feb 03, 2020 10:27 am America/New_York
Answers: 0
Been thanked: 1 time

HICO L-2 Processing

by avmehta » Thu Sep 10, 2020 9:21 am America/New_York

Thanks  for the info George.
In ocssw/scripts I only see python scripts!  When I grep  ocssw_runner I get the following:

install_ocssw.py:        commandStr = os.path.join(installDir, 'scripts', 'ocssw_runner')
install_ocssw.py:        commandStr = os.path.join(installDir, 'run', 'scripts', 'ocssw_runner')
install_ocssw.py:    commandStr = os.path.join(scriptsDir, 'ocssw_runner')
install_ocssw.py:        commandStr = os.path.join(scriptsDir, 'ocssw_runner')

I do not know python!

Once I am able to run l2gen for HICO to get L1 to L2 on GUI on my Mac and check results, I will do bulk processing on my Linux server.

Also, I downloaded  HICO files (.bz2) using direct ftp and use .ISS file in the GUI. If I go through OC L1&L2 Data Browser I get .nc files which GUI does not recognize it seems!

One more question: is it still true that HICO bands >720 nm are not processed for Rrs? Can we get Rhot (and may be Rhos) for all bands?

Sorry for my ignorance/confusion and so many questions as I have not worked with HICO data before.

Thanks again.
Amita

OB.DAAC - SeanBailey
User Services
User Services
Posts: 1464
Joined: Wed Sep 18, 2019 6:15 pm America/New_York
Answers: 1
Been thanked: 5 times

HICO L-2 Processing

by OB.DAAC - SeanBailey » Fri Sep 11, 2020 9:38 am America/New_York

Amita,
The ocssw_runner script is part of SeaDAS and is a bash wrapper to ensure the environment is properly set up when the programs are called by the GUI.

The GUI does recognize the netCDF version of the HICO files we distribute.  By "direct ftp", I assume you mean our Direct Data Access page - (which BTW, is HTTPS, not FTP).  You should preferentially choose the .nc files, as these have an improved navigation over the older HDF5 files (the ISS.bz2).  A small fraction of the data were not able to be processed as netCDF, so we've left the older files accessible.

The l2gen code will only output Rrs in the 'visible' spectrum, which we define a max wavelength of  720nm (above that we consider it NIR...until it's SWIR... then TIR :grin: ).

In the version of l2gen currently distributed with SeaDAS (9.5.0-V2019.3), the abscf_gas.nc needs to be linked to the working directory - which is the directory you launch l2gen from, not the directory l2gen resides in.

Sean

Post Reply