Hi,
i just begin to use NCL (6.1.0 beta),
i try to draw a temperature map about sea level,
my problem is that on zones where altitude is quite high (like mountains) nothing appears
on my map (white gaps).
How can i display this data?
Thanks.
NCL SCRIPT:
load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "region"; define MNLAT, MNLON, MXLAT MXLON.
;***************************************
begin
;;;;
;;; read data from grib file
;;;;
outfile="zzz"
grib_file="gfs.grib2" ; real name = gfs.t06z.pgrb2bf09
fin = addfile(grib_file,"r")
;;; get variables: temperature, longitude, latitude
t2m = fin->TMP_P0_L102_GLL0
bis = t2m(0,:,:)
bis = bis-273.15
printVarSummary(t2m)
wks = gsn_open_wks("ps",outfile) ; open wk station
gsn_define_colormap(wks,"BkBlAqGrYeOrReViWh200")
;;; Font Style
res1 = True ; plots modification on
res1@txFont = "Helvetica"
res1@txFontQuality = "High"
;;; Title
res1@tiMainString = "yann" ; add title
res1@gsnLeftString = "Temp 2M"
res1@gsnRightString = ""
;;; Map
res1@mpDataBaseVersion = "Ncarg4_1" ; choose more recent database
res1@mpDataSetName = "Earth..4" ; high resolution
res1@mpOutlineBoundarySets = "National" ; National borders
res1@mpNationalLineColor = "Black" ; National borders color
res1@mpFillOn = False
;;; T2m
res1@cnFillOn = True ; turn on color fill
res1@cnLinesOn = False ; turn of contour lines
;;; color fill tresholds and colors
res1@cnLevelSelectionMode = "ExplicitLevels" ; set explicit contour levels
res1@cnLevels = (/ -26.,-24.,-22.,-20.,-18.,-16.,-14.,-12.,-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.,12.,14.,16.,18.,20.,22.,24.,26.,28.,30.,32.,34.,36.,38.,40./) ; set levels
res1@cnFillColors = (/ 10,12,14,15,17,19,22,25,31,35,39,43,48,53,61,64,67,71,89,95,99,104,108,111,115,123,130,135,143,148,150,155,159,163,165/) ; set the colors to be used
res1@mpProjection = "LambertConformal"
res1@mpMinLonF = MNLON
res1@mpMaxLonF = MXLON
res1@mpMinLatF = MNLAT
res1@mpMaxLatF = MXLAT
res1@mpLimitMode="LatLon"
plot = gsn_csm_contour_map(wks,bis,res1) ; create plot
Variable: t2m
Type: float
Total Size: 5198400 bytes
1299600 values
Number of Dimensions: 3
Dimensions and sizes: [lv_AMSL1 | 5] x [lat_0 | 361] x [lon_0 | 720]
Coordinates:
lv_AMSL1: [305..4572]
lat_0: [90..-90]
lon_0: [ 0..359.5]
Number Of Attributes: 12
center : US National Weather Service - NCEP (WMC)
production_status : Operational products
long_name : Temperature
units : K
_FillValue : 1e+20
grid_type : Latitude/longitude
parameter_discipline_and_category : Meteorological products, Temperature
parameter_template_discipline_category_number : ( 0, 0, 0, 0 )
level_type : Specific altitude above mean sea level (m)
forecast_time : 9
forecast_time_units : hours
initial_time : 06/19/2012 (06:00)
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 19 05:47:32 2012
This archive was generated by hypermail 2.1.8 : Mon Jun 25 2012 - 09:57:23 MDT