dtrend_msg

Hi

I'm trying to plot the slopes of two 3d data sets (time,lat,lon). It is
27 years of seasonal temperature data. In one case reasonable results
appear in my plots for all seasons in another case three of four plots
look reasonable but the last is very different. I would like to know if
you see anything in the way I'm using dtrend_msg or how I act on the
results of using this function that would produced the inconsistent
output plots.

Attached is the case with the inconsistent set of plots.

thanks
Larry

-- 
Larry R. McDaniel mcdaniel_at_ucar.edu Phone: 303 497 8128
Institute for the Study of Society and Environment
NCAR P.O. Box 3000
Boulder, CO  80307

;*************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
 begin
;************************************************
; create pointer to file and read in temperature
;************************************************
;cdo seasavg /home/mcdaniel/home/lmcd/UDEL/for-larry/cru/from-seth/cru/sub/temp_CRU_all.nc ofile.nc
inf = addfile("ofile.nc","r")
i = 0
do while(i.le.3)
  n1 = i * 27
  n2 = n1 + 26
  dataA = inf->temp(n1:n2,:,:)
  LAT = inf->lat(:)
  LON = inf->lon(:)
  dataA!0 = "time"
  dataA!1 = "lat"
  dataA!2 = "lon"
  dataA&lat_at_units="degrees_north"
  dataA&lon_at_units="degrees_east"
  x_reorder = dataA(lat|:,lon|:,time|:)
  dataADtrend = dtrend_msg(x_reorder&time,x_reorder(lat|:,lon|:,time|:),False,True)
  slope2D = new((/121,261/),double,-999.)
  slope2D!0 = "lat"
  slope2D!1 = "lon"
  LAT1D= inf->dataA!1
  LON1D= inf->dataA!2
  slope2D&lat = dataA&lat
  slope2D&lon = dataA&lon
  slope2D = onedtond(dataADtrend_at_slope, (/121,261/) )
  delete (dataADtrend_at_slope)
  delete (dataADtrend_at_y_intercept)
  res = True
  res_at_gsnAddCyclic = False
  res_at_mpProjection = "Mercator"
  res_at_mpLambertMeridianF = -97.5
  res_at_mpLambertParallel1F = 30.
  res_at_mpLambertParallel2F = 60.
  res_at_tfDoNDCOverlay = False ;; because data is not on native grid
  res_at_mpLimitMode = "Corners"
  res_at_mpLeftCornerLatF = LAT(120)
  res_at_mpLeftCornerLonF = LON(0)
  res_at_mpRightCornerLatF = LAT(0)
  res_at_mpRightCornerLonF = LON(260)
  res_at_mpFillOn = False
  res_at_gsnMaximize = True
  res_at_cnFillOn = True
  res_at_cnLinesOn = False
  res_at_lbLabelAutoStride = True
  res_at_lbBoxLinesOn = False
  res_at_cnFillOn = True ; color Fill
  res_at_cnFillMode = "RasterFill" ; Raster Mode
  res_at_cnLinesOn = False ; Turn off contour lines
  res_at_cnLevelSelectionMode = "ExplicitLevels"
  res_at_cnLevels = (/-1.-.0004,-.0003,-.00025,-.0002,-.00015,-.0001,0.,.0001,.0002,.0003,.0004,1./)

  if (i.eq.0) then
     wks = gsn_open_wks("ps", "Winter_cru-trends")
     res_at_tiMainString = "Winter CRU Trends"
  end if
  if (i.eq.1) then
     wks = gsn_open_wks("ps", "Spring_cru-trends")
     res_at_tiMainString = "Spring CRU Trends"
  end if
  if (i.eq.2) then
     wks = gsn_open_wks("ps", "Summer_cru-trends")
     res_at_tiMainString = "Summer CRU Trends"
  end if
  if (i.eq.3) then
     wks = gsn_open_wks("ps", "Fall_cru-trends")
     res_at_tiMainString = "Fall CRU Trends"
  end if
  gsn_define_colormap(wks,"BlueGreen14")
  plot = gsn_csm_contour_map(wks, slope2D, res)

  i=i+1
end do

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk