Appendix B
MM5 Evaluation using METSTAT and
MM5 Sample Files Appendix B
MM5 Evaluation using METSTAT
Figure B-1: The MM5 36-12-4 km horizontal grid structure
Figure B-2: The MM5 vertical structure and interpolation to CAM-X/CMAQ photochemical grid
Figure B-3 through B-8: METSTAT results for 12 km/4 km grids for June16-24, 2002
Figure B-9 through B-14: METSTAT results for 12 km/4 km grids for July 2-16, 2002
Figure B-15 through B-20: METSTAT results for 12 km/4 km grids for July29-Aug 4, 2002
Figure B-21 through B-23: METSTAT results for 12 km/4 km grids for July 2-16, 2002 (Base 1)
Sample configure.user file for August 1 – August 5, 2002
Sample mm5.deck file for August 1 – August 5, 2002
Changes to mm5.deck file for July 2002 Base 2 Runs to Correct Bias Figure B-1
MM5 Grid DimensionsFigure B-2
MM5 Vertical Structure and Conversion to CAM-X/CMAQ
*Collapsing of MM5 vertical layers to Photochemical Model layers. MM5 Vertical Layers = 34 MM5 Horizontal Structure: (Nx, Ny) 36 km = 165 x 129 12 km = 265 x 241 04 km = 271 x 235 Total Grid Cells = 5.0 mil CAMx/CMAQ Vertical Layers = 16 Horizontal Structure: (Nx,Ny) 36 km = 148 x 112 12 km = 203 x 200 04 km = 254 x 218 Total Grid Cells = 1.58 mil Bottom ofHeight(m)k(MM5)Depth(m)CAM-XDepth(m) Layer (m)1466234184116453514661128223314661135632122810127311062153611101269066309398127298437284287676517277041425336515581226652516025074553245693984235361315223982344822506 294221480246220367126342460209519266 1828182591142818261569171691400161661032913981235151631071141609318106991113158753127881557516751177 598107771535965219764458766151445369775294674514829422057414643743714610933733710973236236723613613636--Surface----Surface----Surface----Surface----Surface----Surface--Grid Conversion: MM5toCAMx METSTAT results for 12 km Grid for July 16-24, 2002
Figure B-3
June 16-24, 2002
12 km Grid Figure B-4
June 16-24, 2002
12 km Grid Figure B-5
Observed/Predicted Wind Direction
0
60
120
180
240
300
360
6/16 6/17 6/18 6/19 6/20 6/21 6/22 6/23 6/24
deg
ObsWndDir PrdWndDir
Bias Wind Direction
-90
-60
-30
0
30
60
90
6/16 6/17 6/18 6/19 6/20 6/21 6/22 6/23 6/24
deg
BiasWndDir
METSTAT results for 4 km
Grid for June 16-24, 2002
04 km Grid Figure B-86/20METSTAT results for 12 km Grid for July 2-16, 2002
Figure B-9
July 2-16, 2002
12 km Grid
Observed/Predicted Temperature
280
285
290
295
300
305
7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17
K
9Figure B-10
ly 2-16, 2002
y
7Figure B-11Observed/Predicted Wind Dir
0
60
120
180
240
300
360
7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17
deg
ection
9METSTAT results for 4 km Grid for July 2-16, 2002
Figure B-12
July 2-16, 2002
4 km Grid
Observed/Predicted Temperature
280
285
290
295
300
305
310
7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17
K
ObsTemp PrdTemp
Bias Temperature
-2
-1
0
1
2
3
4
7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17
K
BiasTemp
RMSE Temperature
0
0.5
1
1.5
2
2.5
3
3.5
4
7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17
K
RMSETemp RMSESTemp RMSEUTemp
IOA Temperature
0
0.2
0.4
0.6
0.8
1
7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17
IOATemp
Figure B-13
July 2-16, 2002
4 km Grid 4 km Grid DirectionMETSTAT results for 12 km Grid for July 29 – Aug 5, 2002
Figure B-15
12 km Grid
yFigure B-17
July 29- Aug 5, 2002
12 km Grid METSTAT results for 4 km Grid for July 29 – Aug 5, 2002
Figure B-18
Figure B-19Figure B-20Error! d DirectionMETSTAT results for 4 km Grid for July 2-16, 2002 (Base 1)
Figure B-21
STL ALL (Jul 2-Jul 16) 04km : Base1Observed/Predicted Temperature280285290295300305310 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17 ObsTemp PrdTemp Bias Temperature-4-3-2-10123456 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17 BiasTemp RMSE Temperature0123456 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17 RMSETemp RMSESTemp RMSEUTemp IOA Temperature00.10.20.30.40.50.60.70.80.91 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17 IOATemp Figure B-22
STL ALL (Jul 2-Jul 16) 04km: Base1Predicted/Observed Humidity02468101214161820 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17g/kgObsHum PrdHum Bias Humidity-1.5-1-0.500.511.52 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17g/kgBiasHum RMSE Humidity00.511.522.53 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17g/kgRMSEHum RMSESHum RMSEUHum IOA Humidity00.10.20.30.40.50.60.70.80.91 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/1
7Figure B-23
16) Observed/Predicted Windspeed0246 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17m/sObsWndSpd PrdWndS Bias Windspeed-1-0.500.511.5 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17m/sBiasWndS RMSE Windspeed0123 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17m/sRMSEWndSpdRMSESWndSpRMSEUW IOA Windspeed00.20.40.60.81 7/ 2 7/ 3 7/ 4 7/ 5 7/ 6 7/ 7 7/ 8 7/ 9 7/10 7/11 7/12 7/13 7/14 7/15 7/16 7/17 IOAWndS Observed/Predicted Wind Direction360 ObsWndDir PrdW# This configure.user.linux file is used to compile on PC running linux only. # For options to compile on Unix systems, please use configure.user # 1. System Variables # 2. User Variables # 3. Fortran options # 3i. PC_PGF77 (Linux) # 4. General commands # 5. Options for making "./include/parame.incl" # 6. Physics Options (memory related) #----------------------------------------------------------------------------# 1. System Variables #----------------------------------------------------------------------------SHELL = /bin/sh .SUFFIXES: .F .i .o .f #----------------------------------------------------------------------------# 2. User Variables #----------------------------------------------------------------------------# RUNTIME_SYSTEM - Currently supported systems. # PC_PGF77 RUNTIME_SYSTEM = "PC_PGF77" #----------------------------------------------------------------------------# 3. Fortran options #----------------------------------------------------------------------------LIBINCLUDE = $(DEVTOP)/include #----------------------------------------------------------------------------# 3i. PC running Linux and using pgf77 compiler # You may also need to unlimit stacksize by typing # limit stacksize unlimited # or # setenv MPSTKZ 8M # (or bigger.) # If your compiler supports both SGI and OpenMP parallel directives, # you need -Mnosgimp in your FCFLAGS. Otherwise remove it. # Add -Kieee option if you compile Gayno-Seaman PBL scheme (IBLTYP=6). # MM5 will fail without this option until the compiler bug is fixed. #----------------------------------------------------------------------------# 3i1. PC_PGF77 (LINUX/Portland Group Inc.) # pgf77 version 1.6 and above # May use pgf90 if the version is 3.1-4 #----------------------------------------------------------------------------FCFLAGS = -I$(LIBINCLUDE) -O2 -Mcray=pointer -tp p6 -p
#FCFLAGS = -
byteswapio -mp -Mnosgimp
CPP = /lib/cpp
CFLAGS = -O NCL
CPPFLAGS = -I$(LIBI
LDOPTIONS = -O2 -Mcra
Wl,-Bstatic y=poi
#LDOPTIONS = -O2 -Mcra
LOCAL_LIBRARIES =
MAKE = make -i -r -------------------------------------
#-
- # 4. General commands
#-----------------------
- AR = ar ru
RM = rm -f
RM_CMD = $(RM) *.CKP *.l
.emacs_* tags TAGS make.log MakeOut *.f ! GREP = grep -s
CC = cc
#-----------------------------------------------------------------------------
# 5. Options for making ./include/parame.incl ----------------
#---------------------------------------
-
#
# FDDAGD (integer) - "1" -> FDDA gridded run FDDAGD = 1
#
# FDDAOBS (integer) - "1" -> FDDA obs run FDDAOBS = 0
#
# MAXNES (integer) - Max Number of Domains in simulation MAXNES = 3
#
# MIX,MJX (integer) - Maximum Dimensions of
MIX = 275
MJX = 275
# MKX (integer)
MKX = 34 ---------------------
#-
-
# 6. Physics Options The first MAXNES values in the list will be used for the correspo
#
# model nests; the rest in the list can be used to compile other
options. # The exception is FRAD, of which only the first value is used in
model, # (i.e., only one radiation option is used for all nests). The rest
allow
# other options to be compiled.
# Compilation of Arakawa-Schubert cumulus sche
#-----------------------------------------------------------------------------
# IMPHYS - foIMPHYS = "5,5# - Dry,stable,warm rain,simple ice,miphase,
# - 1 ,2 ,3 ,4 ,5 #
#
MPHYSTBL = 0
#
# physics #
- 1=use look-up tables for moist phy# (currently only simple ice and mix
phase
#
are available) # # I
#
1,2,3,4,5,6,7,8 ICU
#
# IBLTYP - f
# - 0=no PBL fluxes,1=bulk,2=Blackak-Thompson,4=Eta M-Y,5=MR
# 3=Bur
#
IBLTYP =
# # F
# - Radiation
# 0=none,1=simple,2=cloud,3=ccm2,rrtm=4 FRA
#
# IPOLAR - (integer) for polar model used only if ISOIL=1
#
#
IPOLAR = 0
#
# I
# - 0=no,1=yes (only works with IBLTYP=2,
#
onl
# 3=Pleim-Xiu LSM (IBLTYP=7 only)
ISOIL = 3
#
# ISHALLO (array,integer) - Shallow Convection Option
#
1=shallow convection,0=No shallow convection ISH
#---------------------
-------------------------------------------------------
# Don't touch anything below this line #---------
-------------------------------------------------------------------- .c.o:
$(RM) $@ && \ $(
.F.o:
$(RM) $@
$(FC) -c $(FCFLAGS) $*.F .f.
$(RM) $@
$(FC) -c $(FCFLAGS) $*.f Sample “mm5.stl.base1.deck” (provided by Ken Anderson, Ameren)
#!/bin/sh
date > ./Run/.timing
START_YR="2002"
START_MO="08"
START_DY="01" START_HR="12"
END_YR="2002"
E
ND_MO="08"
END_DY="06" END_HR="12"
# Version 3 of mm5 job deck #
# The mm5 executable (mm5.exe) expects to find the following files
# in the Run/ directory: # MMINPUT_DOMAIN1 -|
#
BDYOUT_DOMAIN1 | --> o
# LOWBDY_DOMAIN1 -| # TERRAIN_DOMAIN[2,3..] if running nests
--> output from Terrain
# rt run:
# If it is a resta# RESTART_DOMAI
# SAVE_DOMAIN1[,2,3...]
# # If it is gridded FDDA run with surface analysis nudging:
# SFCF
DDA_DOMAIN1[2,3,...]
# #
If it is observational nudging run: # MM5OBS_DOMAIN1[,2,3..] --> user-created observation
files
# # Output from a MM5 run:
# If IFTAPE = 1 UT_DOMAIN1[,2,3...] --> one output for each domain
# MMO If IFS
# #
SAVE_DOMAIN1[,2,3...]
# # # temp files
should be accessible
umask 022 # # S
elect appropriate FDDAsw if doing gr
# FDDAsw=yes # gridded FDDA input switch #FD
DAs
#
# Sections # ptions for namelist ("mml
1. O
# 2. Running...
##---------------------------
----
- 1. O
#ptions for namelist ("mmli#--------------------------------------------------------
-
#
# The f
# identifier. Col 1 = Domain #1, Col 2 = Dom #2, etc.
#
# ara
cat > ./Run/op
&OPARAM
; *******
;
; .,
TIMAX = 7200
TISTEP = 108.
; ******
;
;
IFREST = .FALSE., ; whet IXTIMR = 0, ; res
IFSAVE = .TRUE., ; save data for restart nly save the last file for restart
SVLAST = .TRUE., ; T: o
multiple files
SAVFRQ = 720., ; how frequently to save data (in mdel output: 0,1
IFTAPE = 1, ; mo
TAPFRQ = 60., ; how frequently to output model results (in nutes)
mi
BUFFRQ = 1440., ;
minutes), ignored if < TAPFRQ for outputti
INCTAP = 1,1,1,1,1,1,1,1,1,1, ; multipliers of TAPFRQFSKIP = .FALSE., ; whether to skip input files - DO
I
restart RT_MO}-${START_DY}_${START_HR
CDATEST = '${START_YR}-${STATE for the starting file
DA
IFPRT = 0, ; sample print o
PRTFRQ = 720., ; Print frequency for sample output (in minuASCHK = 99999, ; mass conservation check (KTAU or no. of ti
M
IFTSOUT = .FALSE., ; o0.0,0.0,
TSLAT = 0.0,0.0,
negative) 0.0,0.0, ; longitudes of time series poi
TSLON = 0.0,0.0,0.0,
negative) END
&
EOF
cat > ./Run/lparam << EOF
&LPARAM
;
; 1. user-chosen options I
;
RADFRQ = 15., ;atmospheric radiation calculation f= 1, ;moist vertical diffusion in clouds
IMVDIF ly)
on
IVQADV =
linear - 1 ;vertical temperat
IVTADV = 1, linear - 1
0,
ITHADV = 1, ;advection of temperature uses potential temperature - 1, standard - 0
ITPDIF = 1, ;diffusion using p ICOR3D = 1, ;3D Coriolis forceIEXSI = 0, ;initial sea-ice - 0, 1(base on SST), 2(read in) (ISOIL=1 onl
IF
= 1, ;upper radiative boundary condition - 0, 1 ;
; 2. do not chan
;
IB
; (fixed, time-depe
;
;
; unless doing sensitivity runs
;
IFDRY = 0, ;fake-dry run (
0, 1 ;for IMPHYS = 2,3,4,5,6,7 (requires ICUPA = 1) TVAR= 1, ;varying SST in time -
ISS
0
IMO
used, ;1 - used w/o extra input, 2 - user w/ soil
IZ0TOPT= 0, ;thermal ro
2,5. ;0 - old option, 1 - Garratt, 2 - Zili
ISFMTHD= 1, ;method for calculation of 2m/10m diagnostics ;0
IFSNOW = 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, ;SNOW COVER EFFECTS - 0, 1, odel
no effect, 1 - with effect, 2 - simple snow m
ISFFLX = 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ;surface fluxes - 0, 1 ITGFLG =
1:yes, 3:no ISFPAR = 1, 1, 1, 1, 1, 1, 1, 1,
ICLOUD = 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ;cloud effects on radiation - 0, 1
; currently for IFRAD = 1,2 IEVAP =
;(currently for IMPHYS=3,4,5 only) ;MTJ The follwoing Three options pertain to PXLSM ONLY! NUDGE AND IFGROW
added by MTJ
ISMRD = 2, ;soil moisture initialization by P
LSM:
;=0, use moisture avail fr
LANDUSE.TBL ;=1, initialize soil moisture an
temperature from MMINPUT through INTERPX
NNRP
NUDG
E = 1,
requ
ires surface FDDA grid analys ROW = 2,
IFG
fraction between maximum and minimum values for each grid cell from monthly
VEGFRC data (from Terrain)
grid cell and use growth curves for crops. State of natural vegetation is estimated ac
requires an additional ASCII input file for planting dates (user provided).
;=2 same as IFGROW=1 except no planting dates are input. Emergence dates are estimated from VEGFRC. ; Next two sw
RDMAXALB=.FALSE. ;use climo maximum snow albe
landuse table) RDBRDALB=.FALSE. ;use monthly climo background albedo (not la
;
EOF
cat > ./Run/nparam << EOF
&NP
ARAM ;
; ************** NEST AND MOVING NEST OPTIONS **
;
LEVIDN = 0,1,2,1,1,1,1,1,1,1, ; level o
domain
NUM
NC = 1,1,2,1,1,1,1,1,1,1, ; ID of mother domain for eacnest NESTIX = 129, 241, 235, 1, 1, 1, 1, 1, 1,
i NE
j NESTI = 1, 20, 86, 1, 1, 1, 1, 1, 1, 1,
; start location i NESTJ = 1, 50, 102, 1, 1, 1, 1, 1, 1, 1
location i XSTNES = 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., ; dom
initiation XENNES = 7200., 7200., 7200., 0., 0., 0., 0., 0., 0., 0., ; domain
termination
IOVERW = 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, ; overwrite nest input
; 0=interpolate from coarse mesh (for nest domains) (mtj: must
1 for PX) ; 1=read in domain initial conditions;
; 2=read in nest terrain file IACTIV = 1, 1, 1, 0, 0,
; ; ******
; IMOVE = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; move domain 0,
IMOVCO = 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ; 1st move ; i
IMOVEI = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; I move #1 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, ; I move #2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; I move #3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; I move #4
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; I move #5
IMOVEJ = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; J move #
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; J move #2 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, ; J move #3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; J move #4
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; J move #5 IMOVET = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; time of move #1 0, 0, 0, 0
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; time of move 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; time of move #4 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; time of move #5 IFEED = 0, ; no feedback; 9-pt weighted average; 1-pt feedback w/o
smoothing / ; light smoothing / heavy smoothing - 0,1,2,3, and 4
&END
EOF cat > ./Run/ppar&PPARAM ;
;
;
; The values for the follo
; (i.
;
ZZLND = 0.1, ; roughness length over land in meters
ZZW
TR = 0.0001, ; roughness length over water in meters ALBLND = 0.15, ; albedo THINLD
XMAVA = 0.3, ; moisture availability over land as a decimal fracti
; CONF = 1.0, ; non-convective precipitation saturation threshold
(=1: 100%) &E
EOF cat > ./Run/
&FPARAM ;
; ************* 4DDA OPTIONS ********************** ;
; THE FIRST DIMENSION (COLUMN) IS THE DOMAIN IDENTIFIER: ; CO
; ; START TIME FOR FDDA (ANALYSIS OR OBS) FOR EACH DOMAIN
; (IN MINUTES RELATIVE TO MODEL INITIAL TIME) FDASTA=0.,
; ENDING TIME FOR FDDA (ANALYSIS OR OBS) FO
; (IN MINUTES RELATIVE TO MODEL IN
FDAEND=7200.,7200.,0.,0.,0.,0.,0.,0.,0.,
;
; **************** ANALYSIS NUDGING ******************
;
; THE FIRST DIMENSION (COLUMN) OF THE ARRAYS DENOTES THE
; DOMAIN IDENTIFIER:
; COLUMN 1 = DOMAIN #1, COLUMN 2 = DOMAIN #2, ETC.
; THE SECOND DIMENSION (ROW OR LINE) EITHER REFERS TO THE 3D VS
; SFC ANALYSIS OR WHICH VARIABLE IS ACCESSED:
; LINE 1 = 3D, LINE 2 = SFC OR
; LINE 1 = U, LINE 2 = V, LINE 3 = T, LINE 4 = Q
;
; IS THIS A GRID 4DDA RUN? 0 = NO; 1 = YES
I4D= 1,1,0,0,0,0,0,0,0,0,
1,1,0,0,0,0,0,0,0,0,
;
; SPECIFY THE TIME IN MINUTES BETWEEN THE INPUT (USUALLY
; FROM INTERP) USED FOR GRID FDDA
DIFTIM=180.,180.,180.,0.,0.,0.,0.,0.,0.,0., ; 3D ANALYSIS NUDGING
180.,180.,180.,0.,0.,0.,0.,0.,0.,0., ; SFC ANALYSIS NUDGING
;
; GRID NUDGE THE WIND FIELD? 0 = NO; 1 = YES
IWIND=1,1,0,0,0,0,0,0,0,0, ; 3D ANALYSIS NUDGING 1,1,0,0
,0,0,0,0,0,0, ; SFC ANALYSIS NUDGING ;
;
GV=2
.5E-4,1.0E-4,1.0E-4,0 2.5E-4,1.0E-4,1.0E-4,0.;
;
ITEMP=1,1,0,0,0,0,0,0,0,0, ; 3D ANALYSIS NUDGING
;
; NUDGING COEFFICIENT FOR TEMPERATURE ANALY
GT=
2.5E-4,1.0E-4,1.0E-4,0.,0.,0.,0.,0.,0.,0., ; SFC ANALY
;
IMOIS=1,1,0,0,0,0,0,0,0,0, ; 3
0,0,0,0,0,0,0,0,0,0, ; SFC ANALYSIS NUDGI
; ; NUDGING
GQ=
1.E-5,1.E-5,1.E-5,0.,0.,0.,0.,0.,0.,0., ; SFC ANALYSIS NUDGING ;
;
IROT
=0,0,0,0,0,0,0,0,0,0,;
; N5.E6,
GR=
;
;
; BOUNDARY LAYER FROM FDDA OF COARSE GRID THREE DIMENSI
; DATA (USUALLY FROM INTERP),
;
; 1 = YES, EXCLUDE BOUNDARY LAYER NUDGING
INONBL =0,0,0,0,0,0,0,0,0,0, ; U WIND
0,0,0,0,0,0,0,0,0,0, ;
1,1,1,1,1,1,1,1,1,1, ; TEMPERATURE
1,1,1,1,1,1,1,1,1,1, ; MIXING RATIO
; RADIUS OF INFLUENCE FOR SURFACE
;
; IF I4D(2,1)=1 OR I4D(2,2)=1, ETC, DEFINE RINBLW (KM)
;
; OF THE SURFACE-ANALYSIS NUDGING AS A FUNCTION OF SURF
; DATA DENSITY. OV
; ANALYSIS NUDGING IS LINEARLY DECREASED BY 80 PERCEN
; THOSE GRID POINTS GREATER THAN RINBLW FROM AN OBSERVATION
; TO ACCOUNT FOR DECREASED CONFIDENCE IN THE
; IN REGIONS NOT NEAR ANY OBSERVATI
RINBLW=250.,
;
; SET THE NUDGING PRINT FREQUENCY FOR SELEC
; PRINTS IN THE GRID
; TIMESTEPS)
NPF
;
; **************** OBSERVATION NUDGI
NG *************** ;
;
;
INDIVIDUAL OBSERVATION NUDGING. VARIABLES TH ; USE THE FIRST DIMENSION (COLUMN) AS THE DOMAIN
; COLUMN 1 = DOMAIN #1, COLUMN 2 = DOMAIN #2,
;
;
I4DI =0,0,0,0,0,0,0,0,0,0,
;
; OBS NUDGE THE WIND FIELD FROM STATION DATA? 0 = NO; 1 = YES ISW
;
; NUDGING COEFFICIENT FOR WINDS FROM STATION DA
GIV =4.E-4,4.E-4,0.,0.,0.,0.,0.,0.,0.,0.,
;
; OBS NUDGE THE TEMPERATURE FIELD FROM STATION D
ATA? 0 = NO; 1 = YES ISTEMP=1,0,0,0,0,0,0,0,0,0,
;
;
NUDGING COEFFICIENT FOR TEMPERATURES FROM ST GIT =4.E-4,4.E-4,0.,0.,0.,0.,0.,0.,0.,0.,
;
;
OBS NUDGE THE MIXING RATIO FIELD FROM STATION DATA ISMOIS=1,0,0,0,0,0,0,0,0,0,
;
; NUDGING COEFFICIENT FOR THE MIXING RATIO FROM STATION DATA
GIQ
;
; THE OBS NUDGING RADIUS OF INFLUENCE IN THE
;
; FUNCTIONS WHICH SPREAD THE OBS-NUDGING CORRECTION
; IN THE HORIZONTAL.
RIN
;
; THE OBS NUDGING RADIUS OF INFLUENCE IN THE
; VERTICAL IN SIGMA UNITS FOR CR
; WEIGHTED FUNCTIONS WHICH SPREAD THE OBS-NUDGING
; CORRECTION IN THE VERTICAL.
RINSIG=0.001,
;
; THE HALF-PERIOD OF THE TIME WINDOW, IN MINUT
; WHICH AN OBSERVATION WILL AFFECT THE FORECAST
;
; FROM TIMEOBS-TWINDO TO TIMEOBS+TWINDO. THE TEMPO
; WEIGHTING FUNCTION IS DEFINED SUCH THAT THE OBSERVATION
; IS APPLIED WITH FULL STRENGTH WITHIN TWINDO/2. MINUTES
; BEFORE OR AFTER THE OBSERVATION TIME, AND THEN LINEARLY
; DECREASES TO ZERO TWINDO MINUTES BEFORE OR AFTER THE
; OBSERVATION TIME.
TWINDO=40.0,
;
; THE NUDGING PRINT FREQUENCY FOR SELECTED OBS NUDGING CODE (IN CGM TIMESTE
; IN THE
NPF
I=20, ;
; FREQUENCY (IN CGM TIMESTEPS) TO COMPUTE OBS NUDGING W
IONF=2,
IDYNIN=0, ually
grad
; turn off the FDDA before the pure forecast, set i
n=0]
DTR
; nudging (obs nudging and analysis nudging) is ramped d
; from one to zero. Set dtramp negative if FDDA is to b
; down BEFORE the end-of-data time (DATEND), and p
;
FDDA ramp-down period extends beyond the end-of-data &END
EOF
#
#----------------------------------------------------------------- #
#
#
make mmlif
cd .
sed -f ../Util/no_comment.sed mmlif | grep [A-Z,a-z] > mmlif.tmp
mv mmlif.tmp mmlif
rm f
#
# copy gridded FDDA files
#
#if [ $FDDAsw = yes ]; then
# echo "Link grid fdda file"
#
# do
# Num=`echo $i | grep [1-9]$ | sed 's/.*\
#
# echo "ln -sf ../../INTERPF/MMINPUT_DOMAIN$Num M
MINPUT_DO# done
#fi
#---
---
log="log.mm5.out" ;test -f $log && rm $log
#---RUN MPP MM5---
#time mpirun -n
#----RUN 1CPU MM5----
/usr/bin/time -o .timecmd ./mm5.exe > $log 2>&1
#DATE >>.timing
#./logfmt --rm mmlif $log .timecmd .timing $log
#----------------------------------------------------------------
---
exit MM5Deck File Base2 Changes to 4DDA options
Changes
BOLD
************* 4DDA OPTIONS ************
;
; THE FIRST DIMENSION (COLUM
;
COLUMN 1 = DOMAIN ;
; START TIME FOR FDDA (AN
A ; (IN MINUTES RELATIVE TO M
FDASTA=
0.,0.,0.,0.,0.,0.,0.,0.,0.,0. ; ENDING TIME FOR FDDA (ANALYSIS OR OBS) FOR EACH
DOMAIN ; (IN MINUTES RELATIVE TO MODEL INITIAL TIME)
FDAEND=7200.,7200.,7200.,0.,0.,0.,0.,0.,0.,
;
;
;
;
THE FIRST DIMENSION (COLUMN) OF THE ARRAYS DENOTES THE ; DOMAIN IDENTIFIER: ;
COLUMN 1 = DOMAIN #1, COLUMN 2 = DOM ; THE SECOND DIMENSION (ROW OR LINE) EITHER
;
SFC ANALYSIS ; LINE 1 =
3D, LINE 2 = SFC OR ; LINE 1 = U, LINE 2 = V, LINE 3 = T, LINE 4 = Q
;
; IS THIS A GRID
4DDA RUN? 0 = NO; 1 = YE I4D= 1,1,1,0,0,0,0,0,0,0,
1,1,,0,0,0,0,0,0,0, ;
; SPECIFY THE TIME IN MINUTES BETWEEN THE IN
PUT (USUALLY ; FROM INTERP) USED FOR GRID FDDA DIFT
;
;
