WRF Utilities

Document technical information

Format pdf
Size 4.5 MB
First found Nov 13, 2015

Document content analisys

Category Also themed
Language
English
Type
not defined
Concepts
no text concepts found

Organizations

Places

Transcript

WRF Utilities
Cindy Bruyère
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
1
Utilities
 
Graphics
 
Designing a model domain
 
Data
UG: Chapter 3 & 10
Input
Intermediate
Output
grib 1&2
intermediate format
netcdf
netcdf
 
netCDF tools
 
MET
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
2
Graphics : ImageMagick
 
Convert graphical files from one format to another
o 
Many options available (rotate frames, trim white space, etc.)
o 
Can be used for files with single or multiple frames
o 
Cannot deal with .ncgm files
o 
http://www.imagemagick.org convert
convert
convert
convert
file.pdf
file.png
file.pdf
file.ras
file.png
file.bmp file.gif
file.png
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
3
Making Movies
•  Run graphical package
•  Create individual frames for each
image
–  Either directly from graphical
package;
–  Or with a tool like
ImageMagick
•  Use a movie making tool to
create movie
–  {GIF Movie Gear ; Windows ;
commercial software }
•  convert –delay 20 *png movie.gif
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
4
WRF Model Domain Design
[email protected] =
(/"background","DeepSkyBlue",
"ForestGreen","DeepSkyBlue",
"transparent"/)"
"
[email protected]
= 45"
"
[email protected]
= (/ "white",
"Red" , "Red" , "Blue" /)"
"
[email protected]
“NoBoundaries” ; “Geophysical”
“National”
; “USStates”
“GeophysicalAndUSStates”
“AllBoundaries”
•  NCL
–  WPS/util/plotgrids.ncl
–  plotgrids_new.ncl (NCL 6.1.0+)
–  reads namelist information to
generate plot
–  create plot x11/png/pdf
pares = True
[email protected] = "White”
[email protected] = 16
[email protected] = 0.01
gsn_polymarker(wks,mp,-77.26,38.56,
pmres)
"
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
5
WRF Model Domain Design
DOMS = 1
DX = 36.
MAP = "mercator"
LAT1 = (/ -35.0, -45., -27. /)
LAT2 = (/
0., -20., -23. /)
LON1 = (/ 131., 121., 125./)
LON2 = (/ 171., 159., 131./)
parent_id = (/ 0, 1, 2 /)
parent_grid_ratio = (/ 1, 3, 3 /)
ncl design_grids.ncl
Suggested namelist options
parent_id = 0,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 123,
e_sn = 107,
dx = 36000,
dy = 36000,
map_proj = 'mercator',
ref_lat = -17.50,
ref_lon = 151.00,
truelat1 = -17.00,
truelat2 =
0.00,
stand_lon = 151.00,
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
6
obsgrid.exe
geogrid.exe
ungrib.exe
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
7
Fields in geo_em files
 
All invariant (static) data
o 
XLAT & XLONG
o 
MAPFAC
o 
LANDMASK
o 
HGT
o 
GREENFRAC
o 
LANDUSEF & LU_INDEX
o 
ALBEDO12M & SNOALB
o 
SOILTEMP; SOILCTOP & SOILCBOT
Data availability: http://www.mmm.ucar.edu/wrf/users/
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
8
Observational Data
 
http://rda.ucar.edu/datasets/ds353.4
 
http://rda.ucar.edu/datasets/ds464.0  
http://rda.ucar.edu/datasets/ds351.0 BUFRdecode_ADPuprair_littlr.tar is available at: http://rda.ucar.edu/datasets/ds351.0/software
/glade/p/rda/data/ds351.0/little_r
 
ADP or ON29
NCEP dump-bufr
http://rda.ucar.edu/datasets/ds461.0 BUFRdecode_ADPsfc_littlr.tar is available at:
http://rda.ucar.edu/datasets/ds461.0/software
/glade/p/rda/data/ds461.0/little_r
OBSGRID/util/get_rda_data.f
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
9
Observational Data
 
http://rda.ucar.edu/datasets/ds353.4
 
http://rda.ucar.edu/datasets/ds464.0  
http://rda.ucar.edu/datasets/ds351.0 BUFRdecode_ADPuprair_littlr.tar is available at: http://rda.ucar.edu/datasets/ds351.0/software
/glade/p/rda/data/ds351.0/little_r
 
ADP or ON29
NCEP dump-bufr
http://rda.ucar.edu/datasets/ds461.0 BUFRdecode_ADPsfc_littlr.tar is available at:
http://rda.ucar.edu/datasets/ds461.0/software
/glade/p/rda/data/ds461.0/little_r
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
10
Fields needed for IM files – 6hourly
3D Data (e.g. data on pressure levels)
Temperature ; U and V components of Wind
Geopotential Height; Relative Humidity
2D Data
Surface Pressure; Mean Sea Level Pressure
Skin Temperature
2 meter Temperature and Relative Humidity
10 meter U and V components of wind
Recommended
LANDSEA
Soil data (temperature and moisture) & soil height
SST (climate runs – requires)
Water equivalent snow depth & SEAICE
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
11
ungrib : External Data Sources
 
 
 
 
 
NCEP/NCAR Global Reanalysis (NNRP / R1) (~ 2.5o Global; Jan 1948 to present)
http://rda.ucar.edu/datasets/ds090.0 NCEP / DOE Reanalysis II (~ 2.5o Global; Jan 1979 to present)
http://rda.ucar.edu/datasets/ds091.0 GFS 0.5o Global data (Dec 2002 to present)
http://rda.ucar.edu/datasets/ds335.0 ERA Interim Data (~ 0.7o Global; Jan 1979 to present)
http://rda.ucar.edu/datasets/ds627.0 Final Analysis (~ 1.0o Global; FNL; Aug 1999 to present)
http://rda.ucar.edu/datasets/ds083.2
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
12
ungrib : External Data Sources
 
 
 
 
NCEP Climate Forecast System Reanalysis (CFSR)
(~38km, Global; Jan 1979 to Dec 2010)
http://rda.ucar.edu/datasets/ds093.0
NCEP Climate Forecast System Reanalysis (CFSR2)
(0.2o, Global; Jan 2011 to present)
http://rda.ucar.edu/datasets/ds094.0
SST data
http://polar.ncep.noaa.gov/sst (Weekly to Daily; 1o to 1/12o)
http://nomad3.ncep.noaa.gov/pub/sst (Monthly to Weekly)
http://rda.ucar.edu/datasets/277.0 (1894-present; Weekly)
NCAR CESM – CMIP5 data (netCDF)
http://rda.ucar.edu/datasets/ds316.0/
 
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
13
ungrib : External Data Sources
 
 
 
 
GFS Real-time 1o Global
http://www.emc.ncep.noaa.gov
NAM Real-time 32/12km North America
http://www.emc.ncep.noaa.gov GCIP NCEP Eta model output 40km North America
http://rda.ucar.edu/datasets/ds609.2 NCEP North American Regional Reanalysis (NARR) 32km North America
http://rda.ucar.edu/datasets/ds608.2
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
14
http://nomads.ncdc.noaa.gov
•  NAM
•  GFS
•  RUC
•  CFS
•  NARR
•  R1/R2
•  SST
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
15
GRIB
 
Documents and decoders:
o 
 
g1print.exe & g2print.exe
o 
 
Show data available in GRIB1 & GRIB2 files
Available from util/ directory in WPS
grib2ctl.pl
o 
 
GRIB1 and GRIB2
wgrib; wgrib2; unpackgrib2.c; grib2to1.c
http://rda.ucar.edu/#!GRIB
Create .ctl and .idx files, so one can plot GRIB files with GrADS
http://www.cpc.ncep.noaa.gov/products/wesley/grib2ctl.html
ncl_convert2nc
o 
http://www.ncl.ucar.edu/Document/Tools/ncl_convert2nc.shtml
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
16
WPS Intermediate Files
 
Output format of ungrid
 
WPS util/ directory
 
o 
plotfmt.ncl (graphical interface to view intermediate file)
o 
rd_intermediate
Create your own intermediate files
o 
example if you have input data in netCDF format
o 
http://www.mmm.ucar.edu/wrf/OnLineTutorial/WPS/IM_files.htm
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
17
Utility: rd_intermediate
The rd_intermediate lists information about the fields found in an
intermediate-format file
================================================
FIELD = TT
UNITS = K DESCRIPTION = Temperature
DATE = 2000-01-24_12:00:00 FCST = 0.000000
SOURCE = unknown model from NCEP GRID 212
LEVEL = 200100.000000
I,J DIMS = 185, 129
IPROJ = 1
REF_X, REF_Y = 1.000000, 1.000000
REF_LAT, REF_LON = 12.190000, -133.459000
DX, DY = 40.635250, 40.635250
TRUELAT1 = 25.000002
DATA(1,1)=295.910950
================================================
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
18
Utility: plotfmt
The plotfmt program plots the fields in the ungrib intermediateformatted files
ncl plotfmt.ncl 'filename=”FNL:2007-09-15_00”'
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
19
Intermediate Files in netCDF Format
Utility: int2nc.exe
Plot: plotfmt_nc.ncl
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
20
netCDF
 
 
netCDF stands for network Common Data Form
netCDF is one of the current supported data formats chosen
for WRF I/O API
o 
o 
 
 
WRF I/O supports netCDF (not fully CF compliant - Climate and
Forecast Metadata Convention)/ binary/GRIB/HDF
Most support graphical packages currently only support netCDF
file format
http://www.unidata.ucar.edu (documentation)
http://www.unidata.ucar.edu/software/netcdf/docs/netcdf-f77.pdf
http://www.unidata.ucar.edu/software/netcdf/docs/netcdf-f90.pdf (writing Fortran programs to read/write netCDF files)
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
21
NCO tools
http://nco.sourceforge.net/
 
ncdiff
o 
 
Write specified variables / times to a new file
ncrcat -v RAINNC wrfout* -o RAINNC.nc
ncrcat -d Time,0,231 -v RAINNC wrfout* -o RAINNC.nc
ncra (nc average)
o 
 
-o output.nc
ncrcat (nc cat)
o 
 
Difference two file
ncdiff input1.nc input2.nc
Average variables and write to a new file
ncra -v OLR wrfout* -o OLR.nc
ncks (nc kitchen sink)
o 
Combination of NCO tools all in one (handy: one tool for multiple
operations)
Specifically handy to split files
ncks -d Time,1,1 wrfout -o wrfout1.nc
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
22
NCO tools
http://nco.sourceforge.net/
 
ncap2 Arithmetic Processor
 
ncatted ATTribute Editor
 
ncbo Binary Operator (includes ncadd, ncsubtract, ncmultiply, ncdivide)
 
ncea Ensemble Averager
 
ncecat Ensemble conCATenator
 
ncflint FiLe INTerpolator
 
ncks Kitchen Sink
 
ncpdq Permute Dimensions Quickly, Pack Data Quietly
 
ncra Record Averager
 
ncrcat Record conCATenator
 
ncrename RENAMEer
 
ncwa Weighted Averager
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
23
netCDF : Utilities
 
ncdump
o 
reads a netCDF dataset and prints information from the dataset
o 
ncdump –h file
o 
print header (inc. list of variables in the file) ncdump –v VAR file
print data of the variable VAR
ncdump –v Times file WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
24
netCDF : ncdump -v Times
netcdf wrfout_d01_2000-01-24_12:00:00 {
dimensions:
Time = UNLIMITED ; // (3 currently)
DateStrLen = 19 ;
west_east = 73 ;
south_north = 60 ;
west_east_stag = 74 ;
bottom_top = 27 ;
south_north_stag = 61 ;
bottom_top_stag = 28 ;
variables:
char Times(Time, DateStrLen) ;
float LU_INDEX(Time, south_north, west_east) ;
LU_INDEX:FieldType = 104 ;
LU_INDEX:MemoryOrder = "XY " ;
LU_INDEX:description = "LAND USE CATEGORY" ;
LU_INDEX:units = "" ;
LU_INDEX:stagger = "" ;
………
………
global attributes:
………
………
:TITLE = " OUTPUT FROM WRF V3.4.1 MODEL";
:START_DATE = "2000-01-24_12:00:00" ;
:WEST-EAST_GRID_DIMENSION = 74 ;
:SOUTH-NORTH_GRID_DIMENSION = 61 ;
:BOTTOM-TOP_GRID_DIMENSION = 28 ;
:DX = 30000.f ;
:DY = 30000.f ;
data:
Times =
"2000-01-24_12:00:00",
"2000-01-24_18:00:00",
"2000-01-25_00:00:00"
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
25
Wrfout output fields (ncdump -h)
ALBBCK
CF3
DNW
F
GRDFLX
ISTEP
MAPFAC_M
MAPFAC_UY
MAX_MSTFY
OLR
PH
Q2
QVAPOR
RDNW
SEAICE
SNOPCX
SOILTB
T2
TSLB
V
XICEM
XLONG
ZNW
ALBEDO
CFN
DZS
FNM
HFX
IVGTYP
MAPFAC_MX
MAPFAC_V
MF_VX_INV
P_TOP
PHB
QCLOUD
RAINC
RDX
SFROFF
SNOW
SR
TH2
U
V10
XLAND
XLONG_U
ZS
CANWAT
CFN1
E
FNP
HGT
LANDMASK
MAPFAC_MY
MAPFAC_VX
MU
P
POTEVP
QFX
RAINCV
RDY
SH2O
SNOWC
SST
Times
U10
VEGFRA
XLAT
XLONG_V
Total Geopotential, staggered (PH+PHB)
Wind components, grid relative,
staggered (U & V)
10m wind components, grid relative,
mass points (U10 & V10)
CF1
COSALPHA
EDT_OUT
GLW
HGT_SHAD
LH
MAPFAC_U
MAPFAC_VY
MUB
PB
PRATEC
QNDROPSOURCE
RAINNC
RESM
SINALPHA
SNOWH
SWDOWN
TMN
UDROFF
W
XLAT_U
ZETATOP
CF2
DN
EMISS
GRAUPELNC
ISLTYP
LU_INDEX
MAPFAC_UX
MAX_MSTFX
NEST_POS
PBLH
PSFC
QRAIN
RDN
RHOSN
SMOIS
SNOWNC
T
TSK
UST
X
XLAT_V
ZNU
Total Pressure in Pa (P+PB)
Total Potential Temperature (T+300)
Surface temperature in K (T2)
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
26
ncview
http://meteora.ucsd.edu/~pierce/ncview_home_page.html
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
27
ncview
http://meteora.ucsd.edu/~pierce/ncview_home_page.html
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
28
ncview
http://meteora.ucsd.edu/~pierce/ncview_home_page.html
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
29
ncBrowse
http://www.epic.noaa.gov/java/ncBrowse/
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
30
Panoply
http://www.giss.nasa.gov/tools/panoply/
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
31
MET verification software
 
Model Evaluation Tools
 
All the basics (e.g. RMSE, bias, skill scores)
 
Plus
o 
advanced spatial methods (wavelets, objects)
o 
confidence intervals
 
Get it here: http://www.dtcenter.org/met/users/downloads/
 
Get help from [email protected] or the documentation
WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
32
×

Report this document