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 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