Hi All

I was wondering if someone could please help me out with the steps involved to process MODIS L0 data to 250m res TSM data using a regional algorithm in place of tsm_clark. I have processed L0 to L1A then L1B and L2, and can display the 1km res L2 but not when I process LB to 250m res. L2 data; When attempting to display the L2 file I get the error msg- Unable to allocate memory: to array. Please consult supplier of application.

Also I am also not sure how to apply a regional TSM algorithm, or one more applicable to shallow, turbid coastal waters.

Any help would be much appreciated as clearly I am just a beginner. I am using Seadas v5.4.

Cheers
Sarah

> When attempting to display the L2 file I get the error msg- Unable to allocate memory: to array.

It sounds like you may not have enough physical RAM on your computer to display such large images. Before Loading the image you can set "Pixel Sample Rate" and "Line Sample Rate" to higher numbers to decimate the image so it's small enough to display.

Also, try extracting the MODIS L1A data file to create a smaller region before processing to L1B. You can do this using the SeaDAS GUI (Process->MODIS->l1aextract_modis) or the Unix command modis_L1A_extract.csh. If you can display smaller 250m images then it's extremely likely the problem is lack of RAM, and you might want to install more.

Thanks Mike. I tried the extraction and got 'error parsing input parameters' at the l2gen stage...

seadas@seadasva:~/mydata$ modis_L1A_extract.csh A2009217034500.L1A_LAC A2009217034500.GEO 150.9 -24.1 151.5 23.3 A2009217034500_sub.L1A_LAC A2009217034500_sub.GEO
lonlat2pixline A2009217034500.GEO 150.9 -24.1 151.5 23.3

lonlat2pixline 100% done. Reading scan: 1500 out of 1959
start/end line : 809/1960
start/end pixel: 451/754
l1aextract_modis A2009217034500.L1A_LAC 451 754 809 1960 A2009217034500_sub.L1A_LAC

***************************************************************************
MODIS L1A extraction complete. Generating new GEO file..
***************************************************************************

modis_L1A_to_GEO.csh A2009217034500_sub.L1A_LAC -o A2009217034500_sub.GEO
Level 1A input filename: A2009217034500_sub.L1A_LAC
Output geolocation file: A2009217034500_sub.GEO
Satellite: aqua
Year: 2009  Day: 217  Hour: 03  Minute: 45

*Terrain Correction Enabled*

Determining required definitive attitude and ephemeris files...
/home/seadas/seadas5.4/scripts/modis_definitive_atteph.csh aqua 2009217034500 2009217034949

Checking local hard disk for required definitive attitude and ephemeris files...
Present on hard disk: PM1ATTNR.P2009217.0200.003
Present on hard disk: PM1EPHND.P2009216.1200.003

All required definitive attitude and ephemeris files found on your hard disk.

att_file1: /home/seadas/seadas5.4/var/modis/atteph/2009/217/PM1ATTNR.P2009217.0200.003
att_file2: NULL
eph_file1: /home/seadas/seadas5.4/var/modis/atteph/2009/216/PM1EPHND.P2009216.1200.003
eph_file2: NULL

Thu Aug 27 20:06:09 EDT 2009
Creating MODIS geolocation file...
/home/seadas/seadas5.4/bin/geogen_modis
GEO version: 5.0.6 built on May 29 2009 (11:20:15)
scan: 0 out of 116 Thu Aug 27 20:06:14 2009
scan: 10 out of 116 Thu Aug 27 20:06:16 2009
scan: 20 out of 116 Thu Aug 27 20:06:16 2009
scan: 30 out of 116 Thu Aug 27 20:06:17 2009
scan: 40 out of 116 Thu Aug 27 20:06:18 2009
scan: 50 out of 116 Thu Aug 27 20:06:19 2009
scan: 60 out of 116 Thu Aug 27 20:06:20 2009
scan: 70 out of 116 Thu Aug 27 20:06:21 2009
scan: 80 out of 116 Thu Aug 27 20:06:22 2009
scan: 90 out of 116 Thu Aug 27 20:06:23 2009
scan: 100 out of 116 Thu Aug 27 20:06:23 2009
scan: 110 out of 116 Thu Aug 27 20:06:24 2009
geogen_modis exit status: 0
Thu Aug 27 20:06:25 EDT 2009

Running validation test on geolocation file...
/home/seadas/seadas5.4/scripts/modis_geocheck.csh ./A2009217034500_sub.GEO 95
Number of pixels with missing geolocation: 0
Percent of geolocated pixels (rounded): 100
Validation test passed for geolocation file A2009217034500_sub.GEO.

seadas@seadasva:~/mydata$ modis_L1A_to_L1B.csh A2009217034500_sub.L1A_LAC A2009217034500_sub.GEO

L1A filename: A2009217034500_sub.L1A_LAC
GEO filename: A2009217034500_sub.GEO
L1B 1KM filename: A2009217034500.L1B_LAC
L1B HKM filename: A2009217034500.L1B_HKM
L1B QKM filename: A2009217034500.L1B_QKM

Reflective LUT file: /home/seadas/seadas5.4/data/modisa/cal/MYD02_Reflective_LUTs.V5.0.35.12_OC1.hdf
Emissive LUT file: /home/seadas/seadas5.4/data/modisa/cal/MYD02_Emissive_LUTs.V5.0.35.12_OC1.hdf
QA LUT file: /home/seadas/seadas5.4/data/modisa/cal/MYD02_QA_LUTs.V5.0.35.12_OC1.hdf

Satellite: aqua
Year: 2009  Day: 217  Hour: 03  Minute: 45

Thu Aug 27 20:14:29 EDT 2009
Processing MODIS L1A file to L1B...
/home/seadas/seadas5.4/bin/l1bgen_modisa
MYD_PR02 version 5.0.7 built on May 29 2009, at 11:19:50
scan: 0 out of 116 Thu Aug 27 20:14:34 2009
scan: 10 out of 116 Thu Aug 27 20:14:51 2009
scan: 20 out of 116 Thu Aug 27 20:15:00 2009
scan: 30 out of 116 Thu Aug 27 20:15:12 2009
scan: 40 out of 116 Thu Aug 27 20:15:23 2009
scan: 50 out of 116 Thu Aug 27 20:15:34 2009
scan: 60 out of 116 Thu Aug 27 20:15:47 2009
scan: 70 out of 116 Thu Aug 27 20:15:59 2009
scan: 80 out of 116 Thu Aug 27 20:16:11 2009
scan: 90 out of 116 Thu Aug 27 20:16:22 2009
scan: 100 out of 116 Thu Aug 27 20:16:33 2009
scan: 110 out of 116 Thu Aug 27 20:16:45 2009
l1bgen_modisa exit status: 0
Thu Aug 27 20:16:53 EDT 2009

MODIS L1B processing complete.
seadas@seadasva:~/mydata$ ms_met.csh A2009217034500.L1B_LAC
Determining pass start and end times..

Input file : A2009217034500.L1B_LAC
Start time : 2009217034510
End time   : 2009217034949

- None of the required MET files were detected in the LUT archive file -
- Checking local disk for files that are -/+ 6 hours of pass start/end time -

All optimal MET files found on hard disk:

MET1: S200921700_NCEP.MET (optimal)
MET2: S200921706_NCEP.MET (optimal)
MET3: S200921706_NCEP.MET (optimal)

seadas@seadasva:~/mydata$ ms_ozone.csh A2009217034500.L1B_LAC
Determining pass start and end times..

Input file : A2009217034500.L1B_LAC
Start time : 2009217034510
End time   : 2009217034949

Parsing OZONE LUT file /home/seadas/seadas5.4/var/metoz/OZONE_EPTOMS_LUT.txt
*No appropriate files found in LUT*
Parsing OZONE LUT file /home/seadas/seadas5.4/var/metoz/OZONE_TOAST_LUT.txt
*No appropriate files found in LUT*
Parsing OZONE LUT file /home/seadas/seadas5.4/var/metoz/OZONE_TOVS_LUT.txt
*No appropriate files found in LUT*

- None of the required ozone files were detected in the LUT archive files -
- Checking local disk for files that are -/+ 36 hours of pass start/end time -
- Querying ftp server for files that are -/+ 36 hours of pass start/end time -

*** No EPTOMS files found in remote directories:
*** ftp://oceans.gsfc.nasa.gov/METOZ/2009/215/
*** ftp://oceans.gsfc.nasa.gov/METOZ/2009/216/
*** ftp://oceans.gsfc.nasa.gov/METOZ/2009/217/
*** ftp://oceans.gsfc.nasa.gov/METOZ/2009/218/
*** ftp://oceans.gsfc.nasa.gov/METOZ/2009/219/
- No EPTOMS files found on hard disk. Now searching for TOAST -

All optimal OZONE (TOAST) files found on hard disk:

OZONE1: S20092160021623_TOAST.OZONE (optimal)
OZONE2: S20092170021723_TOAST.OZONE (optimal)
OZONE3: S20092170021723_TOAST.OZONE (optimal)

seadas@seadasva:~/mydata$ l2gen ifile=A2009217034500.L1B_LAC geofile=A2009217034500_sub.GEO ofile1=A2009217034500.L2 par=hires.par par=A2009217034500.L1B_LAC.met_list par=A2009217034500.L1B_LAC.ozone_list
Invalid argument "nLw_531,nLw_551,nLw_555,nLw_645,nLw_667,nLw_678,l2_flags "
-E- msl12_input.c Line 367: Error reading input parameters.
-E- l2gen: Error parsing input parameters.

Sorry, disregard that- just realised it was my mistake in the hires.par file.

Would you have any advice about how to apply a tsm algorithm?

Some previous discussions related to TSM, turbid water/shallow water
in this forum may be good references for you:

http://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?pid=8416;hl=tsm#pid8416

http://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?pid=7580;hl=tsm#pid7580

http://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?pid=5469;hl=tsm#pid5469

http://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?pid=1971;hl=turbit#pid1971

http://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?pid=3020;hl=shallow#pid3020

Long

If you wish to manipulate loaded band data (MODIS/SeaWiFS data)
with SeaDAS and wish to have a regional TSM algorithm approach,
you may try to study sample scripts in SeaDAS's User Defined Command
Operation function:

Utilities -> Data Manipulation-> User Defined Operations -> Example 3.

This example may provide a useful tool to calculate and visualize
new L2 product with your modified algorithm.

Long