MODIS data plot: vertical cross section

Hi NCL help,

I am trying to plot vertical cross section of MODIS data (MOD07_L2) for
temparature and humidity for a particular lattude.

I am not able to fix lalitude and vary longitude & vertical levels. Scrtipts
is below:

=================================
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;*********************************
begin
;*********************************
; Read in four MODIS HDF files
;*********************************
  f = addfile("MOD07_L2.A2010277.0535.005.2010277185110.hdf","r")

;*************************************************************
; Read in cloud temperature
;*************************************************************
  wv1s = f->Retrieved_Temperature_Profile(0,:,:)

; Apply scale and offset and convert to double
  wv1 = wv1s@scale_factor*1.d * (wv1s - wv1s@add_offset)

  time = doubletoint(f->Scan_Start_Time(0,0))
gdate = ut_inv_calendar(1993,1,1,0,0,time(0),"years since 1993-01-01",0)
  date2d = ut_calendar(gdate,0)
  date = date2d(0,:)
  lev = (/300,400,500,600,700,800,900,1000/)

  lat2d1 = f->Latitude
  lon2d1 = f->Longitude

;*************************************************************
; Start the graphics
;*************************************************************
  wks = gsn_open_wks("ps", "hdf4sds")
  gsn_define_colormap(wks,"rainbow+gray") ; choose colormap

;*************************************************************
; Set some resources that will apply to the base
; contour/map plot that we are going to use to
; overlay the other contour plots on.
;*************************************************************
  res = True

  res@gsnMaximize = True ; maximize pot in frame
  res@gsnFrame = False ; don't advance frame
  res@gsnDraw = False ; don't draw plot

  res@cnFillOn = True ; color Fill
  res@cnFillMode = "CellFill" ; Raster Mode
  res@cnRasterSmoothingOn = True
  res@cnRasterMinCellSizeF = 0.0005
  res@cnLinesOn = False ; Turn off contour lines
  res@cnLineLabelsOn = False ; Turn off contour lines
  res@cnMaxLevelCount = 100

  res@gsnSpreadColors = True ; use full colormap
res@gsnSpreadColorStart = 5
  res@gsnSpreadColorEnd = -2

  res@trGridType = "TriangularMesh"

  res@gsnAddCyclic = False ; Data is not cyclic

;******************************************************************
; Make a copy of the resources at this point, because these are
; the resources we want to apply to the rest of the contour plots
; we're going to create later.
;******************************************************************
  res1 = res

;
;******************************************************************
; Set the rest of the resources that we only want to apply to
; the base map/contour plot.
;
;******************************************************************
  res@lbTitleString = wv1s@long_name
  res@lbTitleFontHeightF = 0.01 ; Make font smaller
  res@lbLabelAutoStride = True ; Nice stride for labels
  res@lbLabelFontHeightF = 0.012 ; Make labels bigger
  res@lbBoxLinesOn = False

  res@mpDataBaseVersion = "MediumRes" ; Higher res coastline
  res@mpProjection = "Orthographic"

; Set limits of map, based on the min/max of all four datasets.
  res@mpLimitMode = "LatLon"
  res@mpMinLatF = min((/lat2d1/))
  res@mpMaxLatF = max((/lat2d1/))
  res@mpMinLonF = min((/lon2d1/))
  res@mpMaxLonF = max((/lon2d1/))
res@mpGridAndLimbOn = True
  res@mpCenterLonF = 80 ; change map center
  res@mpCenterLatF = 20

  res@pmTickMarkDisplayMode = "Always"

  res@tiMainString = "Four MODIS swaths : " + date(0) + "/"
+ \
                                     date(1) + "/" + date(2) + " " + \
                                     date(3) + ":" + date(4) + " - 18:45"

  res@sfXArray = lon2d1
  res@sfYArray = lat2d1

 res = True
  res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
  res@cnMinLevelValF = 200
  res@cnMaxLevelValF = 300
  res@cnLevelSpacingF = 5
  plot = gsn_csm_contour_map(wks,wv1, res)
draw(plot)
  frame(wks)
end
=================================

Kindly help me to sort the problems.

Thanking you,

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Dec 29 05:25:14 2010