WRF and rcm2rgrid

Hi,

I am having problems when trying to regrid my WRF output to a gaussian
grid (to compare to a reanalysis). The regridded data is mostly, but
not all, missing values. The WRF domain is pan-arctic. I have a
feeling it is due to the ordering of the 2D lat/lon arrays. I have
enclosed both my script and output.

Thanks
Dave

  Copyright (C) 1995-2009 - All Rights Reserved
  University Corporation for Atmospheric Research
  NCAR Command Language Version 5.1.1
  The use of this software is governed by a License Agreement.
  See http://www.ncl.ucar.edu/ for more details.

Variable: fi
Type: float
Total Size: 1804000 bytes
             451000 values
Number of Dimensions: 3
Dimensions and sizes: [Time | 8] x [south_north | 205] x [west_east
| 275]
Coordinates:
Number Of Attributes: 6
   FieldType : 104
   MemoryOrder : XY
   description : TOA OUTGOING LONG WAVE
   units : W m-2
   stagger :
   coordinates : XLONG XLAT
(0)
(0) TOA OUTGOING LONG WAVE: min=95.4118 max=304.182

Variable: lat2d
Type: float
Total Size: 225500 bytes
             56375 values
Number of Dimensions: 2
Dimensions and sizes: [south_north | 205] x [west_east | 275]
Coordinates:
Number Of Attributes: 5
   FieldType : 104
   MemoryOrder : XY
   description : LATITUDE, SOUTH IS NEGATIVE
   units : degrees_north
   stagger :

Variable: lon2d
Type: float
Total Size: 225500 bytes
             56375 values
Number of Dimensions: 2
Dimensions and sizes: [south_north | 205] x [west_east | 275]
Coordinates:
Number Of Attributes: 5
   FieldType : 104
   MemoryOrder : XY
   description : LONGITUDE, WEST IS NEGATIVE
   units : degrees_east
   stagger :

Variable: fir
Type: float
Total Size: 1638400 bytes
             409600 values
Number of Dimensions: 3
Dimensions and sizes: [Time | 8] x [lat | 160] x [lon | 320]
Coordinates:
             lat: [-89.14152..89.14152]
             lon: [ 0..358.875]
Number Of Attributes: 8
   _FillValue : -999
   FieldType : 104
   MemoryOrder : XY
   description : TOA OUTGOING LONG WAVE
   units : W m-2
   stagger :
   coordinates : XLONG XLAT
   ncl : rcm2rgrid used for interpolation

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
;load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "/home/porterdf/src/NCL/WRFUserARW.ncl"

begin
   f = addfile("/data4/porterdf/wrfout20001231.nc","r")
   fi = f->OLR
   lat2d = f->XLAT(0,:,:)
   lon2d = f->XLONG(0,:,:)
   lat2d_at_units = "degrees_north"
   lon2d_at_units = "degrees_east"
           printVarSummary(fi)
           printMinMax(fi, True )
           printVarSummary(lat2d)
           printVarSummary(lon2d)

   jlat = latGau(160, "lat", "latitude", "degrees_north")
   jlon = lonGlobeF(320, "lon", "longitude", "degrees_east")

   fir = rcm2rgrid_Wrap(lat2d,lon2d,fi,jlat,jlon,1)
     printVarSummary(fir)

;************************************************
; create plot
;************************************************
   plDir = "./" ; output dir for plots
   plFile = "test"
   type = "x11" ; create variable
to hold info
   type_at_wkColorModel = "cmyk" ; make output cmyk
   wks = gsn_open_wks(type,plDir+plFile) ; open a ps
file
   gsn_define_colormap(wks,"testcmap") ; choose colormap

   res = True ; plot mods desired
   res_at_gsnPolar = "NH" ; select the N. Hem.
   res_at_mpMinLatF = 60 ; minimum lat to plot
   res_at_mpFillOn = False
   res_at_cnFillOn = True ; color fill
   res_at_cnLinesOn = False ; no contour lines

   res_at_gsnSpreadColors = True ; use full colormap
   res_at_gsnSpreadColorStart = 2
   res_at_gsnSpreadColorEnd = -3

     plot = gsn_csm_contour_map_polar(wks,fir(0,:,:), res)

end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Oct 14 2009 - 17:37:57 MDT