Processing of MODIS L1B to L2 using near real time meterological data

Use this Forum to find information on, or ask a question about, NASA Earth Science data.
Post Reply
ivgbhavani
Posts: 16
Joined: Mon May 21, 2018 11:23 am America/New_York
Answers: 0

Processing of MODIS L1B to L2 using near real time meterological data

by ivgbhavani » Thu May 24, 2018 3:52 am America/New_York

Dear friends,
I could process MODIS L1B data if I use  $l2gen ifile=A2017358085500.L1B_LAC geofile=A2017358085500.GEO ofile=A2017358085500_L2.nc oformat=netCDF4
but when I am trying to process MODIS L1B data by giving meteorological file at met1 option in l2gen,I am getting below error. could anyone help in this regard?
Thanks in advance.

$ l2gen ifile=A2017358085500.L1B_LAC geofile=A2017358085500.GEO ofile=A2017358085500_L2.nc oformat=netCDF4  met1=N201735800_MET_NCEPR2_6h.hdf

Input file A2017358085500.L1B_LAC is MODIS Aqua Level-1B HDF-EOS product.
Loading default parameters from /home/bhavani/SEADAS74/ocssw/share/common/msl12_defaults.par
Loading default parameters for MODISA from /home/bhavani/SEADAS74/ocssw/share/modis/msl12_defaults.par
Loading default sub-sensor parameters for MODISA from /home/bhavani/SEADAS74/ocssw/share/modis/aqua/msl12_defaults.par
Loading parameters for suite OC from /home/bhavani/SEADAS74/ocssw/share/modis/msl12_defaults_OC.par
Loading command line parameters

Loading user parameters for MODISA

Loading characteristics for MODISA
Opening sensor information file /home/bhavani/SEADAS74/ocssw/share/modis/aqua/msl12_sensor_info.dat
  Bnd   Lam       Fo    Tau_r     k_oz    k_no2    t_co2   awhite       aw      bbw
    0  412.000  172.912 3.099e-01 1.987e-03 5.814e-19 1.000e+00 1.000e+00 2.030e-02 2.884e-03
    1  443.000  187.622 2.367e-01 3.189e-03 4.985e-19 1.000e+00 1.000e+00 1.110e-02 2.181e-03
    2  469.000  205.878 1.914e-01 8.745e-03 3.938e-19 1.000e+00 1.000e+00 1.034e-02 1.749e-03
    3  488.000  194.933 1.592e-01 2.032e-02 2.878e-19 1.000e+00 1.000e+00 1.676e-02 1.447e-03
    4  531.000  185.747 1.126e-01 6.838e-02 1.525e-19 1.000e+00 1.000e+00 4.494e-02 1.014e-03
    5  547.000  186.539 9.906e-02 8.622e-02 1.194e-19 1.000e+00 1.000e+00 5.491e-02 8.886e-04
    6  555.000  183.869 9.432e-02 9.553e-02 9.445e-20 1.000e+00 1.000e+00 5.839e-02 8.443e-04
    7  645.000  157.811 5.082e-02 7.382e-02 1.382e-20 1.000e+00 9.133e-01 3.458e-01 4.469e-04
    8  667.000  152.255 4.443e-02 4.890e-02 7.065e-21 1.000e+00 8.921e-01 4.923e-01 3.900e-04
    9  678.000  148.052 4.146e-02 3.787e-02 8.304e-21 1.000e+00 8.783e-01 5.116e-01 3.634e-04
   10  748.000  128.065 2.849e-02 1.235e-02 2.157e-21 1.000e+00 7.832e-01 2.800e+00 2.470e-04
   11  859.000   97.174 1.613e-02 2.347e-03 6.212e-23 1.000e+00 6.519e-01 4.396e+00 1.384e-04
   12  869.000   95.824 1.540e-02 1.936e-03 7.872e-23 1.000e+00 6.403e-01 4.727e+00 1.320e-04
   13 1240.000   45.467 3.617e-03 0.000e+00 0.000e+00 9.994e-01 0.000e+00 1.147e+02 3.021e-05
   14 1640.000   23.977 1.219e-03 0.000e+00 0.000e+00 9.890e-01 0.000e+00 6.689e+02 1.000e-05
   15 2130.000    9.885 4.286e-04 0.000e+00 0.000e+00 9.696e-01 0.000e+00 2.810e+03 3.465e-06

  Bnd   Lam
   16 3750.000
   17 3959.000
   18 4050.000
   19 6715.000
   20 7325.000
   21 8550.000
   22 11000.000
   23 12000.000

Internal data compression requested at compression level: 4
Opening filter file /home/bhavani/SEADAS74/ocssw/share/modis/msl12_filter.dat
Setting 7 x 5 straylight filter on HILT mask

Filter Kernel
1 1 1 1 1 1 1
1 1 1 1 1 1 1
1 1 1 1 1 1 1
1 1 1 1 1 1 1
1 1 1 1 1 1 1

Minimum fill set to 1 pixels

Setting 7 x 5 straylight filter on CLDICE mask

Filter Kernel
1 1 1 1 1 1 1
1 1 1 1 1 1 1
1 1 1 1 1 1 1
1 1 1 1 1 1 1
1 1 1 1 1 1 1

Minimum fill set to 1 pixels

Reading Thuillier_F0.dat
Processing at 1000 meter resolution.
    1000-meter file: A2017358085500.L1B_LAC

Allocated 2343838 bytes in L1 record.
Allocated 920720 bytes in L2 record.

Opening: A2017358085500_L2.nc

The following products will be included in A2017358085500_L2.nc.
0 aot_869
1 angstrom
2 Rrs_412
3 Rrs_443
4 Rrs_469
5 Rrs_488
6 Rrs_531
7 Rrs_547
8 Rrs_555
9 Rrs_645
10 Rrs_667
11 Rrs_678
12 chlor_a
13 chl_ocx
14 Kd_490
15 pic
16 poc
17 ipar
18 nflh
19 par
20 l2_flags

Begin l2gen Version 9.3.0--128-NOTFOUND Processing
Sensor is MODISA
Sensor ID is 7
Sensor has 16 reflective bands
Sensor has 8 emissive bands
Number of along-track detectors per band is 10
Number of input pixels per scan is 1354
Processing pixels 1 to 1354 by 1
Processing scans 1 to 2040 by 1
Ocean processing enabled
Land processing enabled
Atmospheric correction enabled

Begin MSl12 processing at 2018144131304000

Allocated 2343838 bytes in L1 record.
Allocated 2343838 bytes in L1 record.
Allocated 2343838 bytes in L1 record.
Allocated 2343838 bytes in L1 record.
Allocated 2343838 bytes in L1 record.
Loading radiance to brightness temperature from /home/bhavani/SEADAS74/ocssw/share/modis/aqua/cal/bt_modisa.hdf
Loading land mask file from /home/bhavani/SEADAS74/ocssw/share/common/landmask_GMT15ARC.nc
Loading bathymetry mask file from /home/bhavani/SEADAS74/ocssw/share/common/watermask.dat
Loading ice mask file from /home/bhavani/SEADAS74/ocssw/share/common/ice_climatology.hdf
Loaded monthly NSIDC ice climatology HDF file.
Loading elevation file from /home/bhavani/SEADAS74/ocssw/share/common/ETOPO1_ocssw.nc
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_412.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_443.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_469.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_488.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_531.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_547.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_555.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_645.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_667.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_678.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_748.hdf
Loading XCAL rvs and polarization sensitivities from /home/bhavani/SEADAS74/ocssw/var/modisa/xcal/OPER/xcal_modisa_axc_oc_v1.12d_869.hdf
Loading DEM info from /home/bhavani/SEADAS74/ocssw/share/common/ETOPO1_ocssw.nc
Loading climatology file /home/bhavani/SEADAS74/ocssw/share/common/sst_climatology.hdf
Loading SSS reference from Climatology file: /home/bhavani/SEADAS74/ocssw/share/common/sss_climatology_woa2009.hdf

Opening meteorological files.
  met1   = N201735800_MET_NCEPR2_6h.hdf
  met2   =
  met3   =
  ozone1 = /home/bhavani/SEADAS74/ocssw/share/common/ozone_climatology_v2014.hdf
  ozone2 =
  ozone3 =
  no2    = /home/bhavani/SEADAS74/ocssw/share/common/no2_climatology_v2013.hdf

-I- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/anc_acq.c 194: Ancillary file: N201735800_MET_NCEPR2_6h.hdf is not netcdf, thus a std anc file
-I- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/anc_acq.c 194: Ancillary file: /home/bhavani/SEADAS74/ocssw/share/common/ozone_climatology_v2014.hdf is not netcdf, thus a std anc file

****ERROR: get_ancillary: get_clim_data: Vfind failed to find vgroup December
-E- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/setanc.c 652: Error loading relative humidity ancillary data. /home/bhavani/SEADAS74/ocssw/share/common/anc_cor_file_28jan2014.nc

****ERROR: get_ancillary: get_clim_data: Vfind failed to find vgroup December
-E- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/setanc.c 667: Error loading Zonal wind speed ancillary data.

****ERROR: get_ancillary: get_clim_data: Vfind failed to find vgroup December
-E- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/setanc.c 678: Error loading Meridional wind speed ancillary data.

****ERROR: get_ancillary: get_clim_data: Vfind failed to find vgroup December
-E- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/setanc.c 723: Error loading surface pressure ancillary data.

****ERROR: get_ancillary: get_clim_data: Vfind failed to find vgroup December
-E- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/setanc.c 757: Error loading precipitable water ancillary data.

Opening ozone file /home/bhavani/SEADAS74/ocssw/share/common/ozone_climatology_v2014.hdf

Opening NO2 file /home/bhavani/SEADAS74/ocssw/share/common/no2_climatology_v2013.hdf

Opening NO2 frac file /home/bhavani/SEADAS74/ocssw/share/common/trop_f_no2_200m.hdf

-E- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/setanc.c 832: Error loading ancillary data.
-E- /home/bhavani/SEADAS74/ocssw/ocssw-src/src/l2gen/getl1rec.c 153: Error reading A2017358085500.L1B_LAC at scan 0.

Tags:

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

Processing of MODIS L1B to L2 using near real time meterological data

by OB.DAAC - SeanBailey » Fri May 25, 2018 7:18 pm America/New_York

The ancillary input processing is peculiar in that it needs more than one... unless that one is a climatology.  Use getanc.py to retrieve the required set of ancillary inputs.

Sean

Post Reply