Ocean Color Forum - Not logged in
Forum Ocean Color Home Help Search Login
Previous Next Up Topic SeaDAS / SeaDAS 6.x - General Questions / Missing algorithm coefficients (locked) (5606 hits)
By rogers Date 2006-05-18 09:22
After updating to SeaDas 4.9.4 I have problem processing MODIS L2 data.
The msl1,2 gives this error message, even when I uncheck the CHL_OC2 and CHL_OC3 products from the lists:
Loading SST coefficients from /seadas4.9/data/modisa/cal/sst_modisa.dat:
2006091 2006120  1.192  0.945  0.171  1.794
2006091 2006120  2.078  0.925  0.125  1.239

chl_oc2: algorithm coefficients not provided for this sensor.


Roger
By mark Date 2006-05-18 19:37
Roger,

I'm not sure how that could be happening if you de-selected the chl_oc2 product.  My suspicion is that it might somehow still be trying to create a chl_oc2 product.  I would like to see the msl12.par file that gets generated when you do this processing to see what products it is trying to create.  Another possibility is that you might not have installed all of the updates in order.

Mark
By rogers Date 2006-07-19 12:50
Mark,
I have tried to re-run all updated again in successiv order, just to make sure.....
This is the output from the msl12.par file:
ifile=/data/import/PO_TERRA1_L1B_20060718124308997/MOD021KM.01.20060718095006000.000.hdf
geofile=/data/import/PO_TERRA1_L1B_20060718124308997/MOD03.01.20060718095006000.000.hdf
ofile1=level2
l2prod1=chl_oc2,chl_oc3,l2_flags,K_490,sst,Rrs_412,Rrs_443,Rrs_488,Rrs_531,Rrs_551,Rrs_667,Rrs_678,Rrs_748,Rrs_869
polfile=$MSL12_DATA/modist/cal/polcor_modist
spixl=       1
epixl=      -1
dpixl=       1
sline=       1
eline=      -1
dline=       1
ctl_pt_incr=       8
proc_ocean=       1
proc_land=       0
atmocor=       1
proc_sst=       1
pol_opt=       3
aer_opt=      -3
aer_iter_max=      10
brdf_opt=       7
iop_opt=       0
qaa_opt=       0
glint_opt=       1
outband_opt=       2
filter_opt=       1
filter_file=$MSL12_DATA/modisa/msl12_filter.dat
met1=/seadas4.9/data/common/metoz/2006/199/S200619906_NCEP.MET
met2=/seadas4.9/data/common/metoz/2006/199/S200619912_NCEP.MET
met3=/seadas4.9/data/common/metoz/2006/199/S200619912_NCEP.MET
ozone1=$MSL12_DATA/common/S19891991_TOMS.OZONE
sstfile=$MSL12_DATA/common/sst_climatology.hdf
land=$MSL12_DATA/common/landmask.dat
water=$MSL12_DATA/common/watermask.dat
icefile=$MSL12_DATA/common/ice_mask.hdf
alphafile=$MSL12_DATA/common/alpha510_climatology.hdf
tauafile=$MSL12_DATA/common/taua865_climatology.hdf
gain=[0.983,0.9805,0.9796,0.9789,0.9866,0.9863,0.9863,0.992,1.0]
offset=[0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0]
albedo=    0.0270000
rhoamin=  1.00000e-04
qaa_s=    0.0150000
chloc2_wave=[]
chloc2_coef=[]
chloc3_wave=[443,489,550]
chloc3_coef=[0.283,-2.753, 1.457, 0.659,-1.403]
chloc4_wave=[]
chloc4_coef=[]
chlclark_wave=[443,488,551]
chlclark_coef=[0.789273,-3.925523,11.637764,-27.157997,27.936958,-10.398587]
nlwmin=     0.150000
wsmax=      8.00000
tauamax=     0.300000
epsmin=     0.850000
epsmax=      1.35000
glint=   0.00500000
windspeed=   -1000
windangle=   -1000
pressure=   -1000
ozone=   -1000
watervapor=   -1000
relhumid=   -1000
sunzen=      70.0000
satzen=      60.0000
maskland=       1
maskcloud=       1
maskglint=       0
maskbath=       0
masksunzen=       0
masksatzen=       0
maskhilt=       1
maskstlight=       0

Thank you,
Roger
By mark Date 2006-07-19 14:35
Roger,

Well, this is confusing isn't it?  Not only is it specifying the chl_oc2 and chl_oc3 products in the parameter file, but you are also processing Terra data and it's trying to load in Aqua coefficients.  Are you doing all of this through the gui or are you using a command script?  In any case, I need to see the entire terminal output when you try to process this, from beginning to end.  And I need to know if you are doing this by loading in a parameter file.

Mark
By rogers Date 2006-07-20 07:17
Mark,
At this poing I am only using the GUI.
Deselcting chl_oc2 and chl_oc3 and selecting only chlor_a I am able to process the data.
However, selecting all of the products results in this terminal output:

Created MSL12 parameter file: msl12.par
/bin/bash -c "msl12 par=msl12.par"
Input file /data/import/PO_TERRA1_L1B_20060718124308997/MOD021KM.01.20060718095006000.000.hdf is MODIS Terra Level-1B HDF-EOS product.
Loading default parameters for MODIST from /seadas4.9/data/modist/msl12_defaults.par

Loading user parameters for MODIST

Opening filter file /seadas4.9/data/modisa/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

Loading characteristics for MODIST
Opening sensor information file /seadas4.9/data/modist/msl12_sensor_info.dat
  Bnd Lam       Fo    Tau_r     k_oz   awhite       aw      bbw
    0 412  170.367 3.167e-01 1.480e-03 1.000e+00 1.804e-02 3.298e-03
    1 443  186.503 2.377e-01 3.780e-03 1.000e+00 1.709e-02 2.454e-03
    2 488  191.820 1.610e-01 2.220e-02 1.000e+00 2.004e-02 1.634e-03
    3 531  188.567 1.135e-01 6.570e-02 1.000e+00 4.965e-02 1.134e-03
    4 551  187.159 9.990e-02 8.320e-02 1.000e+00 5.720e-02 9.922e-04
    5 667  154.151 4.460e-02 4.870e-02 8.892e-01 5.649e-01 4.239e-04
    6 678  149.877 4.170e-02 3.990e-02 8.892e-01 6.184e-01 3.947e-04
    7 748  128.075 2.860e-02 1.200e-02 7.600e-01 2.813e+00 2.664e-04
    8 869   97.297 1.560e-02 3.750e-03 6.449e-01 4.705e+00 1.408e-04

Calibration to be applied:
Calfile:
Vicarious Gain[1]=0.983000, Offset[1]=0.000000
Vicarious Gain[2]=0.980500, Offset[2]=0.000000
Vicarious Gain[3]=0.979600, Offset[3]=0.000000
Vicarious Gain[4]=0.978900, Offset[4]=0.000000
Vicarious Gain[5]=0.986600, Offset[5]=0.000000
Vicarious Gain[6]=0.986300, Offset[6]=0.000000
Vicarious Gain[7]=0.986300, Offset[7]=0.000000
Vicarious Gain[8]=0.992000, Offset[8]=0.000000
Vicarious Gain[9]=1.000000, Offset[9]=0.000000

Input file /data/import/PO_TERRA1_L1B_20060718124308997/MOD021KM.01.20060718095006000.000.hdf is MODIS Terra Level-1B HDF-EOS product.

Opening: 200607181243_level2

Begin MSl12 Version 5.4.2 Processing
Sensor is MODIST
Sensor ID is 5
Sensor has 9 visible 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 1710 by 1
Ocean processing enabled
Land processing disabled
Atmospheric correction enabled
Using multiple scattering aerosol model with 7/8 algorithm
  and NIR iteration with up to 10 iterations

Begin MSl12 processing at 2006201092921000

Radiance to brightness temperature loaded from /seadas4.9/data/modist/cal/bt_modist.hdf
Loading land mask file from /seadas4.9/data/common/landmask.dat
Loading bathymetry mask file from /seadas4.9/data/common/watermask.dat
Loading ice mask file from /seadas4.9/data/common/ice_mask.hdf
Loading climatology file /seadas4.9/data/common/sst_climatology.hdf

Opening meteorological files.
  met1   = /seadas4.9/data/common/metoz/2006/199/S200619906_NCEP.MET
  met2   = /seadas4.9/data/common/metoz/2006/199/S200619912_NCEP.MET
  met3   = /seadas4.9/data/common/metoz/2006/199/S200619912_NCEP.MET
  ozone1 = /seadas4.9/data/common/metoz/2006/198/S20061980019823_TOAST.OZONE
  ozone2 = /seadas4.9/data/common/metoz/2006/199/S20061990019923_TOAST.OZONE
  ozone3 = /seadas4.9/data/common/metoz/2006/199/S20061990019923_TOAST.OZONE

Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_412_iqu.hdf
Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_443_iqu.hdf
Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_488_iqu.hdf
Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_531_iqu.hdf
Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_551_iqu.hdf
Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_667_iqu.hdf
Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_678_iqu.hdf
Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_748_iqu.hdf
Opening Rayleigh file located at
/seadas4.9/data/modist/rayleigh/rayleigh_modist_869_iqu.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_412.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_443.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_488.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_531.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_551.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_667.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_678.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_748.hdf
Loading polarization file /seadas4.9/data/modist/cal/polcor_modist_869.hdf
Processing scan #     0 (1 of 1710) after      2 seconds

Loading SSES table for sst from /seadas4.9/data/modist/cal/sst_sses_modist.hdf

Loading SSES table for sst4 from /seadas4.9/data/modist/cal/sst4_sses_modist.hdfLoading SST4 coefficients from /seadas4.9/data/modist/cal/sst4_modist.dat:
2006182 2006212 -0.124  1.036  0.683  1.100

Loading SST coefficients from /seadas4.9/data/modist/cal/sst_modist.dat:
2006182 2006212  1.128  0.964  0.151  1.771
2006182 2006212  2.054  0.925  0.128  1.327

Loading mirror-side SST coefficients from /seadas4.9/data/modist/cal/sst_mside_modist.dat.
Using mirror-side SST correction of 0.120000 from 2006335.
chl_oc2: algorithm coefficients not provided for this sensor.

Thank you,
Roger
By @sean Date 2006-07-20 12:15
Roger,

OC2 is not a valid algorithm for MODIS, the coefficients were derived for SeaWiFS
bands.  If you really want it to work for MODIS, you'll have to copy the chloc2_coef 
from the $SEADAS/data/seawifs/msl12_defaults.par and add chloc2_wave=[488,551]. 
I don't recommend doing this, but unless you remove chl_oc2 from your parameter file
for MODIS processing, this is the only way to make it work.

Sean
Previous Next Up Topic SeaDAS / SeaDAS 6.x - General Questions / Missing algorithm coefficients (locked) (5606 hits)



Responsible NASA Official: Gene C. Feldman
Curator: OceanColor Webmaster
Authorized by: Gene C. Feldman
Updated: 27 November 2007
Privacy Policy and Important Notices NASA logo