Skip to content. | Skip to navigation

Personal tools
You are here: GES DISC Home Services Data Tools simap example5.bat

example5.bat

;This is the batch file for simap.pro
;The sample is for Ocean L2 Group 1 gallery



.full_reset_session
.compile simap

sel={name:"",chans:STRARR(38),numch:0}
sel=REPLICATE(sel,20)


inpath="/ftp/data/modis/samples/oceans/L2/"
outpath="/var/scratch"



hdfname="MODOCL2A.A2000338.1810.002.2001033042045.hdf"



;====================FLAGS DEFINITIONS============================

intrpl=1 ;Set to 1 activate filling
;up of holes in the images. This may slow
;the process, but will clean the images.
;Set to 0 to deactivate.


extgeo=1 ;set to 1 to use EXTERNAL geolocations
;set to 0 to use INTERNAL geolocations

igeo=1 ;set to 1 to interpolate geolocations
;It must be 1 when using internal geolocations
;and L1B MODIS files

sdsprint=0 ;Flag to list SDS
;0=no list
;1=short list, SDS names only, no attributes
;2=long list, SDS names, attributes and their
;values


True_col=1 ;Flag to generate image after all procedures
;are completed.
;Set to 3 to save RGB image composite from the
;first 3 channels selected.
;Set to 1 to save grey/color scale separate
;images of all selected channels.

flcnt=1 ;set to 1 to fill continents
;set to 0 not to fill continents


automap=0 ;Set to 1 to automaticaly
;set map limits.
; Set to 0 to use the default
; map limits.


autointen=0 ;Set to 1 to automatically set image intensity.
;Set to 0 to use default image intensity
;ranges "intens"


zlog=1 ;Set to different than 0 to show image
;in log scale

gallery=1 ;Flag to generate galleries of image files.
;0=no gallery, will do global map
;1=gif or jpeg gallery of gray/color scale,
;of true color, RGB, image files. If
;True_col=1 then grey/color scale.
;If True_col=3 then RGB


debug=-1 ;debug mode: -1 will print to monitor
; 4 will print to log file
;============================================================




;============PARAMETERS DEFINITIONS===========================


;---------SUBSET-------
xstart="" ;leave as NULL string for full size
xend="" ;or give numerical value

ystart=""
yend=""


;---------Default map settings-----------------
latmin=21.
latmax=31.
lonmin=-118. ;westboundcoordinate
lonmax=-106. ;eastboundcoordinate

map_limit=[latmin,lonmin,latmax,lonmax]
center_lat=0;latmin+(latmax-latmin)/2.
center_lon=lonmin+ (lonmax-lonmin)/2.
center_lon=atan((sin(lonmin*!dtor)+sin(lonmax*!dtor)),$
(cos(lonmin*!dtor)+cos(lonmax*!dtor)))*!radeg


mappos=[.05,.08,.95,.98] ;map position and size,
; relative to the graphic window

proj=1 ;1 - cylindrical
;2 - stereo
;--------------------------------------------

;IMAGE INTENSITY
;default ranges for data planes:
intens=FLTARR(2,20)

slice=3 ;For 3D arrays, this is the dimension
;to fix. In IDL, dimensions are
;e.g. array[x,y,z]. Set slice=1 or 2 or 3 to fix
; x or y or z dimension, correspondingly.
;For most SDS's slice=3, i.e. they are stacked
;over the z-dimension. However, for example MODIS
;Atmosphere L2 Quality_Assurance needs slice=1 !!!
;I.e. it is stacked over the x-dimension


NXM=800 ;Request in pxl on the size of map
NYM=800 ;to be produced.


ctab=33 ;IDL color table

backgnd=70 ;background inside map
;and ouside the granule(s)

;this will be the size of the graphic window to
;hold the map:
NX=NXM/(mappos[2]-mappos[0])
NY=NYM/(mappos[3]-mappos[1])
;================================================================


sel[0:4].name=['chlor_fluor_ht','chlor_fluor_base','chlor_fluor_effic',$
'phycoeryth_conc','phycou_conc']
sel[0:4].chans=['1']

zlog=0
intens[*,0]=[0.,0.5] ;chlor_fluor_ht
intens[*,1]=[0.,5.] ;chlor_fluor_base
intens[*,2]=[0,0.4] ;chlor_fluor_effic
intens[*,3]=[0,1] ;phycoeryth_conc
intens[*,4]=[0,0.1] ;phycou_conc


simap,hdfname




sel[0:7].name=['CZCS_pigment','chlor_MODIS','pigment_c1_total',$
'susp_solids_conc','cocco_pigmnt_conc','calcite_conc','K_490','cocco_conc_detach']
sel[0:7].chans=['1']

zlog=1
intens[*,0]=[.1,10.] ;CZCS
intens[*,1]=[.1,10.] ;chlor_MODIS
intens[*,2]=[.1,10.] ;pigment_c1
intens[*,3]=[1.,1000.] ;susp_solids
intens[*,4]=[.1,10] ;cocco_pigmnt
intens[*,5]=[.1,10] ;calcite_conc
intens[*,6]=[0.01,10] ;K_490
intens[*,7]=[10.^9,10.^11];cocco_conc_detach


simap,hdfname

Document Actions
NASA Logo - nasa.gov
NASA Privacy Policy and Important Notices
Last updated: Sep 10, 2010 02:20 PM ET
Top