Re: Missing temperature data

Yann,

In the future, it would help if you could provide a sample image that shows the problem.

One guess is that you may need to set:

   res@gsnAddCyclic = True

This is necessary if you have global data *and* it doesn't quite wrap all the way around in the globe in the longitude direction.
This resource will cause an additional longitude point to be added at the end of your array, which simply a duplicate of the
value at the first longitude in your array.

--Mary

On Jun 19, 2012, at 5:47 AM, Yann Bertrand wrote:

> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 19 08:17:19 2012