Missing temperature data

From: Yann Bertrand <yann74_at_nyahnyahspammersnyahnyah>
Date: Tue Jun 19 2012 - 05:47:23 MDT

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