drawing from binary(grads) format file

From: Ȳ <hyj1004_at_nyahnyahspammersnyahnyah>
Date: Mon, 4 Dec 2006 15:19:29 +0900

Hi ,
I'm NCL beginner.
I need help .

data pediod : 212
lat : 121 (20N-90N, 0.5 interval) grids
lon : 241 (190-310, 0.5 interval ) grids
Data value for North America except ocean is non-missing.

I drew contour plot for 1st time.
But the figure using NCL was different from using Grads.

I use the following NCL script.

I want to draw exactly.

Help me~
Thank you ~

***********************************************************************************
  latS = 20.
  latN = 80.
  lonL = 190.
  lonR = 310.

  file1="na_ano_cru_ts2_19502002_djfm.dat"
  SLP = fbindirread(file1,0,(/212,121,241/),"float")

  SLP!0 = "time"
  SLP!1 = "lat"
  SLP!2 = "lon"
  SLP@_FillValue = -999.999

  lat = fspan(20,80,121)
  lon = fspan(190,310,241)
  time = ispan(1950,2161,1)

  SLP&lat = lat
  SLP&lon = lon
  SLP&time = time

  res = True

  res_at_mpMinLatF = latS
  res_at_mpMaxLatF = latN
  res_at_mpMinLonF = lonL
  res_at_mpMaxLonF = lonR

  res_at_cnFillOn = True
  res_at_cnLineLabelsOn = False
  res_at_gsnSpreadColors = True

  wks = gsn_open_wks("X11","SLP")
  gsn_define_colormap(wks,"gui_default") ; choose colormap
  contour = gsn_csm_contour_map(wks,SLP(0,:,:),res)

***********************************************************************************

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Dec 03 2006 - 23:19:29 MST

This archive was generated by hypermail 2.2.0 : Tue Dec 05 2006 - 13:29:24 MST