Fwd: mask on 1D (X Y Z) satellite data

From: Lorraine Manlay <lorraine.manlay_at_nyahnyahspammersnyahnyah>
Date: Sun, 12 Apr 2009 22:47:14 +0200

Hi ncl users,
I try to do some plot of temperature from satellite data.

The data are 1D : lon lat T (temperature).

I use sfXArray and sfYArray to generate something that I can plot after
with gsn_csm_contour_map.
The result is fine, you can see it in the attach file. And the script is
paste below.

My question is :
Is there a possibility to mask the zone were the contour is not "physical",
i.e where there is no point ... In fact I just want the zone that are
covered by the satellite and I want to mask the ones that are not. Is it
possible ? Is there a simple solution using mask ?

Thanks ncl users.

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/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

;;;; Load data



  print("Reading file: "+infile)

  T=a->$fldname$ ; load Temperature data


  indice_col= 0 ;



;;; Just 1 level : indice_level, make a mask






;;;; Plot

  wks = gsn_open_wks("eps","tempe_"+lev(indice_level)*0.01)

  res = True
  res_at_gsnFrame = False
  res_at_gsnMaximize = True

  res_at_gsnSpreadColors = True
  res_at_gsnSpreadColorStart = 2
  res_at_gsnSpreadColorEnd = 90

  res_at_mpMaxLatF = -40.
  res_at_gsnPolar = "SH"
  res_at_cnRasterModeOn = True

  res_at_trGridType = "TriangularMesh"

  res_at_cnLineLabelPlacementMode = "Constant"

  res_at_cnFillOn = True
  res_at_cnFillMode = "RasterFill"
  res_at_cnLinesOn = False
  res_at_cnLineLabelsOn = False

  res_at_lbLabelStride = 2

  res_at_sfXArray = lon_mask
  res_at_sfYArray = lat_mask

  res_at_tiMainString = "Temperature at
"+lev(indice_level)*0.01+" hPa"

  map = gsn_csm_contour_map(wks,T_mask,res)

; Draw markers on the plot in the lat/lon locations.

  mkres = True
  mkres_at_gsMarkerIndex = 17
  mkres_at_gsMarkerSizeF = 0.0005




ncl-talk mailing list
List instructions, subscriber options, unsubscribe:

Received on Sun Apr 12 2009 - 14:47:14 MDT

This archive was generated by hypermail 2.2.0 : Thu Apr 16 2009 - 10:05:49 MDT