-
Notifications
You must be signed in to change notification settings - Fork 663
/
README.namelist
1936 lines (1781 loc) · 159 KB
/
README.namelist
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
Description of namelist variables
---------------------------------
Note: variables followed by (max_dom) indicate that this variable needs to
be defined for the nests when max_dom > 1.
&time_control
run_days = 0, ! run time in days
run_hours = 0, ! run time in hours
Note: if it is more than 1 day, one may use both run_days and run_hours
or just run_hours. e.g. if the total run length is 36 hrs, you may
set run_days = 1, and run_hours = 12, or run_days = 0, and run_hours = 36.
run_minutes = 0, ! run time in minutes
run_seconds = 0, ! run time in seconds
start_year (max_dom) = 2001, ! four digit year of starting time
start_month (max_dom) = 06, ! two digit month of starting time
start_day (max_dom) = 11, ! two digit day of starting time
start_hour (max_dom) = 12, ! two digit hour of starting time
start_minute (max_dom) = 00, ! two digit minute of starting time
start_second (max_dom) = 00, ! two digit second of starting time
Note: the start time is used to name the first wrfout file.
It also controls the start time for nest domains, and the time to restart
end_year (max_dom) = 2001, ! four digit year of ending time
end_month (max_dom) = 06, ! two digit month of ending time
end_day (max_dom) = 12, ! two digit day of ending time
end_hour (max_dom) = 12, ! two digit hour of ending time
end_minute (max_dom) = 00, ! two digit minute of ending time
end_second (max_dom) = 00, ! two digit second of ending time
It also controls when the nest domain integrations end
All start and end times are used by real.exe.
Note that one may use either run_days/run_hours etc. or
end_year/month/day/hour etc. to control the length of
model integration. But run_days/run_hours
takes precedence over the end times.
Program real.exe uses start and end times only.
interval_seconds = 10800, ! time interval between incoming real data, which will be the interval
between the lateral boundary condition file (in seconds)
input_from_file (max_dom) = T, ! whether nested run will have input files for domains other than 1
fine_input_stream (max_dom) = 0, ! field selection from nest input for its initialization
0: all fields are used! 2: only static and time-varying, masked land
surface fields are used. This requires the use of
io_form_auxinput2
history_interval (max_dom) = 60, ! history output file interval in minutes
frames_per_outfile (max_dom) = 1, ! number of output times per history output file,
used to split output into multiple files
into smaller pieces
restart = F, ! whether this run is a restart run
cycling = F, ! whether this run is a cycling run, if so, initializes look-up table for Thompson schemes only
restart_interval = 1440, ! restart output file interval in minutes
reset_simulation_start = F, ! whether to overwrite simulation_start_date with forecast start time
io_form_history = 2, ! 2 = netCDF
io_form_restart = 2, ! 2 = netCDF
io_form_input = 2, ! 2 = netCDF
io_form_boundary = 2, ! netCDF format
= 4, ! PHD5 format
= 5, ! GRIB1 format
= 10, ! GRIB2 format
= 11, ! pnetCDF format
ncd_nofill = .true., ! only a single write, not the write/read/write sequence
debug_level = 0, ! 50,100,200,300 values give increasing prints
diag_print = 0, ! print out time series of model diagnostics
0 = no print
1 = domain averaged 3-hourly hydrostatic surface pressure tendency
(Dpsfc/Dt), and dry-hydrostatic column pressure tendency (Dmu/Dt)
will appear in stdout file
2 = in addition to those above, domain averaged rainfall,
surface evaporation, and sensible and latent heat fluxes will be output
all_ic_times = .false., ! whether to write out wrfinput for all processing times
adjust_output_times = .false., ! adjust output times to the nearest hour
override_restart_timers = .false., ! whether to change the alarms from what is previously set
write_hist_at_0h_rst = .false., ! whether to output history file at the start of restart run
write_restart_at_0h = .false., ! whether to output restart file at the start of restart run
output_ready_flag = .true., ! asks the model to write-out an empty file with the name
'wrfoutReady_d<domain>_<date>.
Useful in production runs so that post-processing code can check on the
completeness of this file
force_use_old_data = .false., ! if set to .true., allow WRF Version 3 input data
=.false., stop when WRF model detects Version 3 input data
To choose between SI and WPS input to real for EM core:
auxinput1_inname = "met_em.d<domain>.<date>" ! Input to real from WPS (default since 3.0)
= "wrf_real_input_em.d<domain>.<date>" ! Input to real from SI
Other output options: Note all auxhist[1-24], auxinput[2-24] interval variables are domain dependent
auxhist9_outname = "rainfall" ! file name for extra output! if not specified,
auxhist9_d<domain>_<date> will be used
also note that to write variables in output other
than the history file requires Registry.EM file change
auxhist9_interval (max_dom) = 10, ! interval in minutes
io_form_auxhist9 = 2, ! output in netCDF
frames_per_auxhist9 = 1000, ! number of output times in this file
For SST updating (used only with sst_update=1):
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval (max_dom) = 360 ! minutes generally matches time given by interval_seconds
io_form_auxinput4 = 2 ! IO format
nwp_diagnostics = 0 ! set to = 1 to add 7 history-interval max diagnostic fields
For additional regional climate surface fields
output_diagnostics = 0 ! set to = 1 to add 36 surface diagnostic arrays (max/min/mean/std)
auxhist3_outname = 'wrfxtrm_d<domain>_<date>' ! file name for added diagnostics
io_form_auxhist3 = 2 ! netcdf
auxhist3_interval (max_dom) = 1440 ! minutes between outputs (1440 gives daily max/min)
frames_per_auxhist3 = 1 ! output times per file
Note: do restart only at multiple of auxhist3_intervals
For observation nudging:
auxinput11_interval = 10 ! interval in minutes for observation data. It should be
set as or more frequently as obs_ionf (with unit of
coarse domain time step).
auxinput11_end_h = 6 ! end of observation time in hours.
Options for run-time IO:
iofields_filename (max_dom) = "my_iofields_list.txt",
(example: +:h:21:rainc, rainnc, rthcuten)
ignore_iofields_warning = .true., ! what to do when encountering an error in the user-specified files
.false., ! abort when encountering an error in iofields_filename file
nocolons = .false., ! whether to write file date with colons
.true., e.g. 2020-05-20_12_00_00
Additional settings when running WRFVAR:
write_input = t, ! write input-formatted data as output
inputout_interval (max_dom) = 180, ! interval in minutes when writing input-formatted data
input_outname = 'wrfinput_d<domain>_<date>' ! you may change the output file name
inputout_begin_y (max_dom) = 0
inputout_begin_m = 0
inputout_begin_d (max_dom) = 0
inputout_begin_h (max_dom) = 3
inputout_begin_m (max_dom) = 0
inputout_begin_s (max_dom) = 0
inputout_end_y (max_dom) = 0
inputout_end_m = 0
inputout_end_d (max_dom) = 0
inputout_end_h (max_dom) = 12
inputout_end_m (max_dom) = 0
inputout_end_s (max_dom) = 0 ! the above shows that the input-formatted data are output
starting from hour 3 to hour 12 in 180 min interval.
For automatic moving nests: requires special input data, and environment variable TERRAIN_AND_LANDUSE set at compile time
(This option will overwrite input_from_file for nest domains)
input_from_hires (max_dom) = .true.,
rsmas_data_path = "path-to-terrain-and-landuse-dataset"
&domains
time_step = 60, ! time step for integration in integer seconds
recommend 6*dx (in km) for typical real-data cases
time_step_fract_num = 0, ! numerator for fractional time step
time_step_fract_den = 1, ! denominator for fractional time step
Example, if you want to use 60.3 sec as your time step,
set time_step = 60, time_step_fract_num = 3, and
time_step_fract_den = 10
time_step_dfi = 60, ! time step for DFI, may be different from regular time_step
reasonable_time_step_ratio = 6., ! Any d01, real-data case with a time step ratio larger than this is stopped. Except for specific circumstances (e.g., using IEVA), this value should be no larger than 6 (default).
max_dom = 1, ! number of domains - set it to > 1 if it is a nested run
s_we (max_dom) = 1, ! start index in x (west-east) direction (leave as is)
e_we (max_dom) = 91, ! end index in x (west-east) direction (staggered dimension)
s_sn (max_dom) = 1, ! start index in y (south-north) direction (leave as is)
e_sn (max_dom) = 82, ! end index in y (south-north) direction (staggered dimension)
s_vert (max_dom) = 1, ! start index in z (vertical) direction (leave as is)
e_vert (max_dom) = 30, ! end index in z (vertical) direction (staggered dimension)
Note: this refers to full levels including surface and top
vertical dimensions need to be the same for all nests
Note: most variables are unstaggered (= staggered dim - 1)
dx (max_dom) = 10000, ! grid length in x direction; unit in meters.
dy (max_dom) = 10000, ! grid length in y direction; unit in meters.
ztop (max_dom) = 19000. ! used in mass model for idealized cases
grid_id (max_dom) = 1, ! domain identifier
parent_id (max_dom) = 0, ! id of the parent domain
i_parent_start (max_dom) = 0, ! starting LLC I-indices from the parent domain
j_parent_start (max_dom) = 0, ! starting LLC J-indices from the parent domain
parent_grid_ratio (max_dom) = 1, ! parent-to-nest domain grid size ratio: for real-data cases
the ratio has to be odd! for idealized cases,
the ratio can be even if feedback is set to 0.
parent_time_step_ratio (max_dom) = 1, ! parent-to-nest time step ratio! it can be different
from the parent_grid_ratio
feedback = 1, ! feedback from nest to its parent domain; 0 = no feedback
smooth_option = 2 ! smoothing option for parent domain, used only with feedback
option on. 0: no smoothing; 1: 1-2-1 smoothing; 2: smoothing-desmoothing (default)
Namelist variables specifically for the WPS input for real:
num_metgrid_soil_levels = 4 ! number of vertical soil levels or layers input
from WPS metgrid program
num_metgrid_levels = 27 ! number of vertical levels of 3d meteorological fields coming
from WPS metgrid program
interp_type = 2 ! vertical interpolation
1 = linear in pressure
2 = linear in log(pressure)
extrap_type = 2 ! vertical extrapolation of non-temperature fields
1 = extrapolate using the two lowest levels
2 = use lowest level as constant below ground
t_extrap_type = 2 ! vertical extrapolation for potential temperature
1 = isothermal
2 = -6.5 K/km lapse rate for temperature
3 = constant theta
use_levels_below_ground = .true. ! in vertical interpolation, use levels below input surface level
T = use input isobaric levels below input surface
F = extrapolate when WRF location is below input surface value
use_surface = .true. ! use the input surface level data in the vertical interp and extrap
T = use the input surface data
F = do not use the input surface data
lagrange_order = 2 ! vertical interpolation order
1 = linear
2 = quadratic
9 = cubic spline
zap_close_levels = 500 ! ignore isobaric level above surface if delta p (Pa) < zap_close_levels
lowest_lev_from_sfc = .false. ! place the surface value into the lowest eta location
! T = use surface value as lowest eta (u,v,t,q)
! F = use traditional interpolation
force_sfc_in_vinterp = 1 ! use the surface level as the lower boundary when interpolating
through this many eta levels
0 = perform traditional trapping interpolation
n = first n eta levels directly use surface level
maxw_horiz_pres_diff = 5000 ! Pressure threshold (Pa). For using the level of max winds, when the
pressure difference between neighboring values exceeds this maximum,
the variable is NOT inserted into the column for vertical interpolation.
trop_horiz_pres_diff = 5000 ! Pressure threshold (Pa). For using the tropopause level, when the
pressure difference between neighboring values exceeds this maximum,
the variable is NOT inserted into the column for vertical interpolation.
maxw_above_this_level = 30000 ! Minimum height (actually it is pressure in Pa) to allow using the
level of max wind information in real. With a value of 300 hPa, then
a max wind value at 500 hPa will be ignored.
use_maxw_level 0=do not use max wind speed level in vertical interpolation inside
of the real program, 1 = use level
use_trop_level as above, with tropopause level data
sfcp_to_sfcp = .false. ! Optional method to compute model's surface pressure when incoming
data only has surface pressure and terrain, but not SLP
smooth_cg_topo = .false. ! Smooth the outer rows and columns of domain 1's topography w.r.t.
the input data
use_tavg_for_tsk = .false. ! whether to use diurnally averaged surface temp as skin temp. The
diurnally averaged surface temp can be computed using WPS utility
avg_tsfc.exe. May use this option when SKINTEMP is not present.
aggregate_lu = .false. ! whether to aggregate the grass, shrubs, trees in dominant landuse;
default is false.
rh2qv_wrt_liquid = .true., ! whether to compute RH with respect to water (true) or ice (false)
rh2qv_method = 1, ! which method to use to computer mixing ratio from RH:
default is option 1, the old MM5 method; option 2 uses a WMO
recommended method (WMO-No. 49, corrigendum, August 2000) -
use_sh_qv = .false., ! whether to use specific humidity or mixing ratio data from input
recommended if input data has high vertical resolution
interp_theta = .false. ! If set to .false., it will vertically interpolate temperature
instead of potential temperature, which may reduce bias when
compared with input data
hypsometric_opt = 2, ! = 1: default method
= 2: it uses an alternative way (less biased
when compared against input data) to compute height in program
real and pressure in model.
wif_input_opt = 0 ! = 1: option to process the Water Ice Friendly Aerosol input from metgrid for use with mp_physics=28
= 2: since V4.4, option to use black carbon aerosol category with mp_physics=28, as well as its radiative effect. Must include
file QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat during WPS
num_wif_levels = 30 ! number of levels in the Thompson Water Ice Friendly aerosols (mp_physic=28)
p_top_requested = 5000 ! p_top (Pa) to use in the model
vert_refine_fact = 1 ! vertical refinement factor for ndown, not used for concurrent vertical grid refinement
vert_refine_method (max_dom) = 0 ! vertical refinement method
0: no vertical refinement
1: integer vertical refinement (must set rebalance=1)
2: use specified or computed eta levels for vertical refinement (must set rebalance=1)
rebalance = 0 ! must be set to =1 if vertical nesting is used.
ts_buf_size = 200 ! size of time series buffer
max_ts_locs = 5, ! maximum number of time series locations
max_ts_level = 15, ! highest model level for time series output
tslist_unstagger_winds = .false., ! whether to interpolate 3D u and v to cell centers
Users may explicitly define full eta levels. Given are two distributions for 28 and 35 levels. The number
of levels must agree with the number of eta surfaces allocated (e_vert). Users may alternatively request
only the number of levels (with e_vert), and the real program will compute values. There are
two methods selected with auto_levels_opt = 1 (old) or 2 (new). The old computation assumes
a known first several layers, then generates equi-height spaced levels up to the top of the model.
The new method has surface and upper stretching factors (dz_stretch_s and dz_stretch_u) to stretch levels according to
log p up to where it reaches the maximum thickness (max_dz) and starting from thickness dzbot.
The stretching transitions from dzstretch_s to dzstretch_u by the time the thickness reaches max_dz/2.
max_dz = 1000. ! maximum level thickness allowed (m)
auto_levels_opt = 1 ! old
= 2 ! new default (also set dzstretch_s, dzstretch_u, dzbot, max_dz)
dzbot = 50. ! thickness of lowest layer (m) for auto_levels_opt=2
dzstretch_s = 1.3 ! surface stretch factor for auto_levels_opt=2
dzstretch_u = 1.1 ! upper stretch factor for auto_levels_opt=2
eta_levels = 1.000, 0.990, 0.978, 0.964, 0.946,
0.922, 0.894, 0.860, 0.817, 0.766,
0.707, 0.644, 0.576, 0.507, 0.444,
0.380, 0.324, 0.273, 0.228, 0.188,
0.152, 0.121, 0.093, 0.069, 0.048,
0.029, 0.014, 0.000,
eta_levels = 1.000, 0.993, 0.983, 0.970, 0.954,
0.934, 0.909, 0.880, 0.845, 0.807,
0.765, 0.719, 0.672, 0.622, 0.571,
0.520, 0.468, 0.420, 0.376, 0.335,
0.298, 0.263, 0.231, 0.202, 0.175,
0.150, 0.127, 0.106, 0.088, 0.070,
0.055, 0.040, 0.026, 0.013, 0.000
= 0,2, ! this allows vertical nesting in the nest domain
Note that with vertical nesting one can only use RRTM and RRTMG radiation physics
An example to define vertical nested levels (in program real):
e_vert = 35, 45,
eta_levels(1:35) = 1., 0.993, 0.983, 0.97, 0.954, 0.934, 0.909, 0.88, 0.8406663, 0.8013327,
0.761999, 0.7226653, 0.6525755, 0.5877361, 0.5278192, 0.472514,
0.4215262, 0.3745775, 0.3314044, 0.2917579, 0.2554026, 0.2221162,
0.1916888, 0.1639222, 0.1386297, 0.1156351, 0.09525016, 0.07733481,
0.06158983, 0.04775231, 0.03559115, 0.02490328, 0.0155102, 0.007255059, 0.
eta_levels(36:81) = 1.0000, 0.9946, 0.9875, 0.9789, 0.9685, 0.9562, 0.9413, 0.9238, 0.9037, 0.8813, 0.8514,
0.8210, 0.7906, 0.7602, 0.7298, 0.6812, 0.6290, 0.5796, 0.5333, 0.4901, 0.4493, 0.4109,
0.3746, 0.3412, 0.3098, 0.2802, 0.2524, 0.2267, 0.2028, 0.1803, 0.1593, 0.1398, 0.1219,
0.1054, 0.0904, 0.0766, 0.0645, 0.0534, 0.0433, 0.0341, 0.0259, 0.0185, 0.0118, 0.0056, 0.
ideal_init_method method to compute alb in idealized cases in start_em
= 1, alb from phb (default); = 2, alb from t_init
Horizontal interpolation options, coarse grid to fine grid. The default is to use
the Smolarkiewicz "SINT" method. However, this is known to break with the
implementation inside of WRF for large refinement ratios (such as 15:1). For those
extreme (and quite rare occurrences), other schemes are available. For options
1, 3, 4, and 12, the FG lateral boundaries use the same horizontal scheme for the
lateral BC computations.
interp_method_type = 1 ! bi-linear interpolation
= 2 ! SINT, (default)
= 3 ! nearest neighbor - only to be used for
testing purposes
= 4 ! overlapping quadratic
=12 ! again for testing, uses SINT horizontal
interpolation, and same scheme for
computation of FG lateral boundaries
Variables specifically for the 3d ocean initialization with a single profile. Set
the ocean physics option to #2. Specify a number of levels. For each of those levels,
provide a depth (m) below the surface. At each depth provide a temperature (K) and
a salinity (ppt). The default is not to use the 3d ocean. Even when the 3d ocean is
activated, the user must specify a reasonable ocean. Currently, this is the only way
available to run the 3d ocean option.
&physics
sf_ocean_physics = activate ocean model (0=no, 1=1d mixed layer; 2=3D PWP, no bathymetry)
&domains
ocean_levels = 30,
ocean_z ; vertical profile of layer depths for ocean (in meters), e.g.:
= 5., 15., 25., 35., 45., 55.,
65., 75., 85., 95., 105., 115.,
125., 135., 145., 155., 165., 175.,
185., 195., 210., 230., 250., 270.,
290., 310., 330., 350., 370., 390.
ocean_t ; vertical profile of ocean temps, e.g.:
= 302.3493, 302.3493, 302.3493, 302.1055, 301.9763, 301.6818,
301.2220, 300.7531, 300.1200, 299.4778, 298.7443, 297.9194,
297.0883, 296.1443, 295.1941, 294.1979, 293.1558, 292.1136,
291.0714, 290.0293, 288.7377, 287.1967, 285.6557, 284.8503,
284.0450, 283.4316, 283.0102, 282.5888, 282.1674, 281.7461
ocean_s ; vertical profile of salinity, e.g.:
= 34.0127, 34.0127, 34.0127, 34.3217, 34.2624, 34.2632,
34.3240, 34.3824, 34.3980, 34.4113, 34.4220, 34.4303,
34.6173, 34.6409, 34.6535, 34.6550, 34.6565, 34.6527,
34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247,
34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247
Namelist variables for controlling the specified moving nest:
Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
to the ARCHFLAGS. The maximum number of moves, max_moves, is set to 50
but can be modified in source code file frame/module_driver_constants.F.
num_moves = 0 ! total number of moves
move_id(max_moves) = 2,2,2,2, ! a list of nest domain id's, one per move
move_interval(max_moves) = 60,120,150,180, ! time in minutes since the start of this domain
move_cd_x(max_moves) = 1,1,0,-1, ! the number of parent domain grid cells to move in i direction
move_cd_y(max_moves) = 1,0,-1,1, ! the number of parent domain grid cells to move in j direction
positive is to move in increasing i and j direction, and
negative is to move in decreasing i and j direction.
0 means no move. The limitation now is to move only 1 grid cell
at each move.
Namelist variables for controlling the automatic moving nest:
Note that this moving nest option needs to be activated at the compile time by adding -DMOVE_NESTS
and -DVORTEX_CENTER to the ARCHFLAGS. This option uses a mid-level vortex following algorithm to
determine the nest move. This option is experimental.
vortex_interval(max_dom) = 15 ! how often the new vortex position is computed
max_vortex_speed(max_dom) = 40 ! used to compute the search radius for the new vortex position
corral_dist(max_dom) = 8 ! how many coarse grid cells the moving nest is allowed to get
near the mother domain boundary
track_level = 50000 ! pressure value in Pa where the vortex is tracked
time_to_move(max_dom) = 0. ! time (in minutes) to start the moving nests
tile_sz_x = 0, ! number of points in tile x direction
tile_sz_y = 0, ! number of points in tile y direction
can be determined automatically
numtiles = 1, ! number of tiles per patch (alternative to above two items)
nproc_x = -1, ! number of processors in x for decomposition
nproc_y = -1, ! number of processors in y for decomposition
-1: code will do automatic decomposition
>1: for both: will be used for decomposition
Namelist variables for controlling the adaptive time step option:
use_adaptive_time_step = .false. ! T/F use adaptive time stepping.
step_to_output_time = .true. ! if adaptive time stepping, T/F modify the
time steps so that the exact history time is reached
target_cfl(max_dom) = 1.2,1.2 ! vertical and horizontal CFL <= to this value implies
no reason to reduce the time step, and to increase it
target_hcfl(max_dom) = .84,.84 ! horizontal CFL <= to this value implies
max_step_increase_pct(max_dom) = 5,51 ! percentage of previous time step to increase, if the
max(vert cfl, horiz cfl) <= target_cfl, then the time
will increase by max_step_increase_pct. Use something
large for nests (51% suggested)
starting_time_step(max_dom) = -1,-1 ! flag = -1 implies use 4*dx (defined in start_em),
starting_time_step = 100 means the starting time step
for the coarse grid is 100 s
max_time_step(max_dom) = -1,-1 ! flag = -1 implies max time step is 8*dx,
max_time_step = 100 means that the time step will not
exceed 100 s
min_time_step(max_dom) = -1,-1 ! flag = -1 implies max time step is 3*dx,
min_time_step = 100 means that the time step will not
be less than 100 s
adaptation_domain = 1 ! default, all fine grid domains adaptive dt driven by coarse-grid
2 = Fine grid domain #2 determines the fundamental adaptive dt.
&dfi_control
dfi_opt = 0 ! which DFI option to use (3 is recommended)
0 = no digital filter initialization
1 = digital filter launch (DFL)
2 = diabatic DFI (DDFI)
3 = twice DFI (TDFI)
dfi_nfilter = 7 ! digital filter type to use (7 is recommended)
0 = uniform
1 = Lanczos
2 = Hamming
3 = Blackman
4 = Kaiser
5 = Potter
6 = Dolph window
7 = Dolph
8 = recursive high-order
dfi_write_filtered_input = .true. ! whether to write wrfinput file with filtered
model state before beginning forecast
dfi_write_dfi_history = .false. ! whether to write wrfout files during filtering integration
dfi_cutoff_seconds = 3600 ! cutoff period, in seconds, for the filter
dfi_time_dim = 1000 ! maximum number of time steps for filtering period
this value can be larger than necessary
dfi_bckstop_year = 2004 ! four-digit year of stop time for backward DFI integration
dfi_bckstop_month = 03 ! two-digit month of stop time for backward DFI integration
dfi_bckstop_day = 14 ! two-digit day of stop time for backward DFI integration
dfi_bckstop_hour = 12 ! two-digit hour of stop time for backward DFI integration
dfi_bckstop_minute = 00 ! two-digit minute of stop time for backward DFI integration
dfi_bckstop_second = 00 ! two-digit second of stop time for backward DFI integration
dfi_fwdstop_year = 2004 ! four-digit year of stop time for forward DFI integration
dfi_fwdstop_month = 03 ! two-digit month of stop time for forward DFI integration
dfi_fwdstop_day = 13 ! two-digit month of stop time for forward DFI integration
dfi_fwdstop_hour = 12 ! two-digit month of stop time for forward DFI integration
dfi_fwdstop_minute = 00 ! two-digit month of stop time for forward DFI integration
dfi_fwdstop_second = 00 ! two-digit month of stop time for forward DFI integration
dfi_savehydmeteors = 0 ! option for radar DA: 0: set hydrometeors to 0 before DFI and
let them spin up in DFI; 1: keep them unchanged in DFI
&physics
Note: even the physics options can be different in different nest domains,
caution must be used as what options are sensible to use
chem_opt (max_dom) = 0, ! chemistry option - use WRF-Chem
mp_physics (max_dom) microphysics option
= 0, no microphysics
= 1, Kessler scheme
= 2, Lin et al. scheme
= 3, WSM 3-class simple ice scheme
= 4, WSM 5-class scheme
= 5, Ferrier (new Eta) microphysics, operational High-Resolution Window version
= 6, WSM 6-class graupel scheme
= 7, Goddard 4-ice scheme
= 8, Thompson scheme
= 9, Milbrandt-Yau 2-moment scheme
= 10, Morrison (2 moments)
= 11, CAM 5.1 microphysics
= 13, SBU_YLIN scheme
= 14, WDM 5-class scheme
= 16, WDM 6-class scheme
= 18, NSSL 2-moment 4-ice scheme with predicted (unactivated) CCN (or activated CCN)
to change global CCN value, use
nssl_cccn = 0.7e9 ; CCN (#/m^3 at sea level pressure) for NSSL scheme (18) or nssl_ccn_on=1
Also sets ccn_conc for mp_physics=18
For NSSL 1-moment schemes, intercept and particle densities can be set for snow,
graupel, hail, and rain. For the 1- and 2-moment schemes, the shape parameters
for graupel and hail can be set.
PLEASE SEE doc/README.NSSLmp for options affecting the NSSL scheme
= 17, 19, 21, 22: Legacy NSSL-MP options: see README.NSSLmp for equivalent settings with 18
= 24, WSM 7-class scheme (separate hail and graupel categories)
= 26, WDM 7-class scheme (separate hail and graupel categories)
= 28, aerosol-aware Thompson scheme with water- and ice-friendly aerosol climatology
This option has two climatological aerosol input options:
use_aero_icbc = .F. : use constant values
use_aero_icbc = .T. : use climatological aerosol input from WPS
use_rap_aero_icbc = .false. ! Set to .true. to ingest real-time data containing aerosols (new in 4.4)
qna_update = 0 ! set to 1 to update time-varying sfc aerosol emission from climatology or
real-time data
with mp_physics = 28. Use with input file ‘wrfqnainp_d0*’
(must set auxinput17_interval and io_form_auxinput17; new in 4.4)
wif_fire_emit = .false. ! set to .true. to include biomass burning organic and black carbon
aerosols with mp_physics = 28 (new in 4.4)
wif_fire_inj = 1 ! (default) vertically distribute biomass burning emissions
in mp_physics = 28 (new in 4.4)
= 30, HUJI (Hebrew University of Jerusalem, Israel) spectral bin microphysics,
fast version
= 32, HUJI spectral bin microphysics, full version
= 38, Thompson scheme with 2-moment graupel/hail (new in 4.5)
= 40, Morrison (2 moments) with consideration of CESM-NCSU RCP4.5 climatological
aerosol
= 50, P3 1-ice category, 1-moment cloud water
= 51, P3 1-ice category plus double-moment cloud water
= 52, P3 2-ice categories plus double-moment cloud water
= 53, P3 1-ice category, 3-moment ice, plus double-moment cloud water
= 55, Jensen-ISHMAEL (Ice-Spheroids Habit Model with Aspect-ratio Evolution) scheme
predicting qc, qr, and three ice species with four predictive equations for each ice
= 56, NTU multi-moment scheme ! for ntu3m
ccnty = 1 ! for aerosol background type, marine, for ntu3m
2 ! continental clean (default), for ntu3m
3 ! continental average, for ntu3m
4 ! continental urban, for ntu3m
= 95, Ferrier (old Eta) microphysics
= 96, Madwrf
= 97, Goddard GCE scheme (also uses gsfcgce_hail, gsfcgce_2ice)
For non-zero mp_physics options, to keep Qv .GE. 0, and to set the other moisture
fields .LT. a critical value to zero
mp_zero_out = 0, ! no action taken, no adjustment to any moist field
= 1, ! except for Qv, all other moist arrays are set to zero
if they fall below a critical value ('moist' array only)
= 2, ! Qv is .GE. 0, all other moist arrays are set to zero
if they fall below a critical value ('moist' array only)
mp_zero_out_thresh = 1.e-8 ! critical value for moist array threshold, below which
moist arrays (except for Qv) are set to zero (kg/kg)
mp_zero_out_all = 0, ! if =1 and mp_zero_out>0, then reproduce old behavior and
apply threshold to scalar, chem, and tracer arrays
gsfcgce_hail = 0 ! for running gsfcgce microphysics with graupel
= 1 ! for running gsfcgce microphysics with hail
default value = 0
gsfcgce_2ice = 0 ! for running with snow, ice and graupel/hail
= 1 ! for running with only ice and snow
= 2 ! for running with only ice and graupel
(only used in very extreme situation)
default value = 0
gsfcgce_hail is ignored if gsfcgce_2ice is set to 1 or 2.
hail_opt = 0 ! hail switch for WSM6 and WDM6 : 0 - off, 1 - on
morr_rimed_ice = 1 ! hail switch for Morrison schemes (mp_physics=10 and 40: 0 - off,
1 - on (new in 4.0)
clean_atm_diag = 0 ! if set to =1, turns on clean sky diagnostics (for chem); default is 0=off
acc_phy_tend = 0 ! set to =1 to output 16 accumulated physics tendencies for potential temp,
water vaopr mixing ratio, and U/V wind components; default is 0=off (new in 4.4)
progn (max_dom) = 0 ! switch to use mix-activate scheme (Only for Morrison, WDM6, WDM5,
and NSSL_2MOM)
ccn_conc = 1.E8 ! CCN concentration, used by WDM schemes (set automatically for NSSL_2MOM using nssl_cccn)
no_mp_heating = 0 ! normal
= 1 ! turn off latent heating from a microphysics scheme
use_mp_re = 1 ! whether to use effective radii computed in mp schemes in RRTMG
0: do not use; 1: use effective radii
(The mp schemes that compute effective radii are 3,4,6,7,8,10,14,16,18,24,26,28,50-53,55)
force_read_thompson = .false. ! whether to read tables for mp_physics = 8,28
write_thompson_tables = .true. ! whether to read or compute tables for mp_phyiscs = 8,28
write_thompson_mp38table = .false. ! whether to read table (qr_acr_qg_mp38V1.dat) for mp_physics = 38
ra_lw_physics (max_dom) longwave radiation option
= 0, no longwave radiation
= 1, rrtm scheme
(Default values for GHG in V4.2: co2vmr=(280. + 90.*exp(0.02*(yr-2000)))*1.e-6
n2ovmr=319.e-9, ch4vmr=1774.e-9
Values used in previous versions: co2vmr=330.e-6, n2ovmr=0., ch4vmr=0.)
= 3, cam scheme
also requires levsiz, paerlev, cam_abs_dim1/2 (see below)
= 4, rrtmg scheme
(Default values for GHG in V4.2: co2vmr=(280. + 90.*exp(0.02*(yr-2000)))*1.e-6
n2ovmr=319.e-9, ch4vmr=1774.e-9,
cfc11=0.251e-9, cfc12=0.538e-9,
= 14, rrtmg-k scheme from KIAPS
= 24, fast rrtmg scheme for GPU and MIC (since 3.7)
(Default values for GHG in V4.2: co2vmr=(280. + 90.*exp(0.02*(yr-2000)))*1.e-6
n2ovmr=319.e-9, ch4vmr=1774.e-9,
cfc11=0.251e-9, cfc12=0.538e-9
= 5, Goddard longwave scheme
= 7, FLG (UCLA) scheme
= 31, Earth Held-Suarez forcing
= 99, GFDL (Eta) longwave (semi-supported)
also must use co2tf = 1 for ARW
ra_sw_physics (max_dom) shortwave radiation option
= 0, no shortwave radiation
= 1, Dudhia scheme
= 2, Goddard short wave
= 3, cam scheme
also must set levsiz, paerlev, cam_abs_dim1/2 (see below)
= 4, rrtmg scheme
= 14, rrtmg-k scheme from KIAPS
= 24, fast rrtmg scheme for GPU and MIC (since 3.7)
= 5, Goddard shortwave scheme
= 7, FLG (UCLA) scheme
= 99, GFDL (Eta) longwave (semi-supported)
also must use co2tf = 1 for ARW
radt (max_dom) = 30, ! minutes between radiation physics calls
recommend 1 min per km of dx (e.g. 10 for 10 km);
use the same value for all nests.
cldovrlp = 2, ! cloud overlapping option for RRTMG only. 1=random, 2=maximum-random (default),
3=maximum, 4=exponential, 5=exponential-random
idcor = 0, ! decorrelation length flag for cldovrlp=4 or 5
0 = constant decorrelation length, 2500 m;
1 = latitude-varying decorrelation length
ra_sw_eclipse = 0, ! eclipse effect on shortwave radiation. 0: off, 1: on.
Applies to RRTMG, Goddard, old Goddard and Dudhia schemes
ghg_input = 1, ! Option to read CAMtr_volume_mixing_ratio files of green house gas values,
as of v4.4, the default is SSP 2 with RCP 4.5 => SSP245
Used for CAM LW and SW, RRTM, RRTMG LW and SW, RRTMG_fast LW and SW
0 = do not read in the annual data
1 = read in time dependent data for CO2, N2O, CH4, CFC11, CFC12
co2tf CO2 transmission function flag only for GFDL radiation
= 0, ! read CO2 function data from pre-generated file
= 1, ! generate CO2 functions internally in the forecast
ra_call_offset radiation call offset
= 0 ! (=0, no offset; =-1, old offset)
swint_opt Interpolation of short-wave radiation based on the updated solar zenith angle
between SW call
= 0, ! no interpolation
= 1, ! use interpolation
= 2, ! Calls the Fast All-Sky Radiation Model for Solar applications (FARMS)
! every model time step.
couple_farms = .false. ! True) uses FARMS SW radiation to drive the LSM, False) Uses the SW radiation from rad_sw_physics
cam_abs_freq_s = 21600 ! default CAM clearsky longwave absorption calculation frequency
(recommended minimum value to speed scheme up)
levsiz = 59 ! for CAM radiation input ozone levels, set automatically
paerlev = 29 ! for CAM radiation input aerosol levels, set automatically
cam_abs_dim1 = 4 ! for CAM absorption save array, set automatically
cam_abs_dim2 = value of e_vert for CAM 2nd absorption save array, set automatically
o3input = ozone input option for radiation (currently rrtmg only)
= 0, ! using profile inside the code
= 2, ! using CAM ozone data (ozone.formatted)
aer_opt = aerosol input option for radiation (currently rrtmg only)
= 0, ! none
= 1, ! using Tegen (1997) data,
= 2, ! using J. A. Ruiz-Arias method (see other aer_* options)
= 3, ! using G. Thompson's water/ice friendly climatological aerosol
alevsiz = 12 ! for Tegen aerosol input levels, set automatically
no_src_types = 6 ! for Tegen aerosols: organic and black carbon, sea salt, sulfalte, dust,
and stratospheric aerosol (volcanic ashes - currently 0), set automatically
The following aerosol options allow RRTMG and new Goddard radiation schemes to see it, but the aerosols are
constant during the model integration.
aer_aod550_opt (max_dom) = [1,2]
1 = input constant value for AOD at 550 nm from namelist.
In this case, the value is read from aer_aod550_val;
2 = input value from auxiliary input 15. It is a time-varying 2D grid in netcdf
wrf-compatible format. The default is aer_aod550_opt=1 and aer_aod550_val=0.12
aer_aod550_val (max_dom) = 0.12
aer_angexp_opt (max_dom) = [1,2,3] :
1 = input constant value for Angstrom exponent from namelist. In this case,
the value is read from aer_angexp_val;
2 = input value from auxiliary input 15, as in aer_aod550_opt;
3 = Angstrom exponent value estimated from the aerosol type defined in aer_type, and modulated
with the RH in WRF. Default operation is aer_angexp_opt = 1, and aer_angexp_val=1.3.
aer_angexp_val (max_dom) = 1.3
aer_ssa_opt (max_dom) = [1,2,3] similar to aer_angexp_opt.
aer_ssa_val (max_dom) = 0.85
aer_asy_opt (max_dom) = [1,2,3] similar to aer_angexp_opt.
aer_asy_val (max_dom) = 0.9
aer_type (max_dom) = [1,2,3] : 1 for rural, 2 is urban and 3 is maritime.
sf_sfclay_physics (max_dom) surface-layer option (old bl_sfclay_physics option)
= 0, no surface-layer
= 1, Revised MM5 Monin-Obukhov scheme (Jimenez, renamed in v3.6)
= 2, Monin-Obukhov (Janjic) scheme
= 3, NCEP Global Forecast System scheme
= 4, QNSE surface layer
= 5, MYNN surface layer
= 7, Pleim-Xiu surface layer
= 10, TEMF surface layer
= 91, Old MM5 scheme (previously option 1)
sf_surface_physics (max_dom) land-surface option (old bl_surface_physics option)
= 0, no surface temp prediction
= 1, thermal diffusion scheme
= 2, Unified Noah land-surface model
= 3, RUC land-surface model
= 4, Noah-MP land-surface model (see additional &noah_mp namelist)
= 5, Community Land Model version 4 (CLM4), adapted from CAM
= 6, Community Terrestrial Systems Model (CTSM). See doc/README.CTSM for info on running this option.
= 7, Pleim-Xiu LSM
= 8, Simplified Simple Biosphere Model (SSiB)
- can be used with Dudhia/RRTM, CAM or RRTMG radiation options
sf_urban_physics(max_dom) = 0, ! activate urban canopy model (in Noah and Noah-MP LSMs only)
= 0: no
= 1: Single-layer, UCM
= 2: Multi-layer, Building Environment Parameterization (BEP) scheme
(works only with YSU, MYJ and BouLac PBL)
= 3: Multi-layer, Building Environment Model (BEM) scheme
(works only with YSU, MYJ and BouLac PBL)
num_urban_ndm = 1! (= 2 if BEP or BEM active) maximum number of street dimensions (ndm in BEP or BEM header)
num_urban_ng = 1! (= 10 if BEP or BEM active) number of grid levels in the ground (ng_u in BEP or BEM header)
num_urban_nwr = 1! (= 10 if BEP or BEM active) number of grid levels in the walls or roof (nwr_u in BEP or BEM header)
num_urban_nz = 1! (= 18 if BEP or BEM active) maximum number of vertical levels in the urban grid (nz_um in BEP or BEM header)
num_urban_ngb = 1! (= 10 if BEM active) number of grid levels in the ground below building (ngb_u in BEM header)
num_urban_nf = 1! (= 10 if BEM active) number of grid levels in the floors (nf_u in BEM header)
num_urban_nbui = 1! (= 15 if BEM active) maximum number of types of buildings in an urban class (nbui_max in BEM header)
sf_surf_irr_scheme (max_dom) = 0: no irrigation scheme
= 1: channel scheme (in Noah only)
= 2: drip scheme
= 3: sprinkler scheme
irr_daily_amount (max_dom) = daily equivalent irrigation water amount in mm (float, e.g. 5.7)
irr_start_hour (max_dom) = UTC start hour for irrigation (int, e.g. 5)
irr_num_hours (max_dom) = number of consecutive hours for irrigation (int, e.g. 3)
irr_start_julianday (max_dom) = Julian start day of irrigation (int, e.g. 135)
irr_end_julianday (max_dom) = Julian end day of irrigation (int, e.g. 255)
irr_ph (max_dom) = 0: all the static IRRIGATION field is activated
= 1: pseudo random activation field: ensures repeatability across compilers
= 2: random activation field: results might depend on the compiler fortran RANDOM_SEED
irr_freq (max_dom) = irrigation frequency in days (int, e.g. 3)
bl_pbl_physics (max_dom) boundary-layer option
= 0, no boundary-layer
= 1, YSU scheme
= 2, Mellor-Yamada-Janjic TKE scheme
= 3, Hybrid EDMF GFS scheme
= 4, Eddy-diffusivity Mass Flux, Quasi-Normal Scale Elimination PBL
= 5, MYNN 2.5 level TKE scheme, works with
sf_sfclay_physics=1 or 2 as well as 5
(option 6 removed in 4.5, replaced by bl_mynn_closure options)
bl_mynn_closure = 2.5: level 2.5
2.6: level 2.6
3.0: level 3.0
= 7, ACM2 (Pleim) PBL
= 8, Bougeault and Lacarrere (BouLac) PBL
= 9, UW boundary layer scheme from CAM5 (CESM 1_0_1)
= 10, TEMF (Total Energy Mass Flux) scheme
sf_sfclay_physics=10
= 11, Shin-Hong 'scale-aware' PBL scheme
= 12, Grenier-Bretherton-McCaa scheme
= 16, TKE+TKE dissipation rate (epsilon) scheme
works with surface layer options 1,91,2,5.
= 17, TKE+TKE dissipation rate+TPE (temparature variance) scheme
= 99, MRF scheme
bldt (max_dom) = 0, ! minutes between boundary-layer physics calls
grav_settling (max_dom) gravitational settling of fog/cloud droplets (Now works for any PBL scheme)
= 0, No settling of cloud droplets
= 1, Settling from Dyunkerke 1991 (in atmos and at surface)
= 2, Fogdes (vegetation & wind speed dependent; Katata et al. 2008) at surface
and Dyunkerke in the atmos.
ysu_topdown_pblmix = 0, ! whether to turn on top-down, radiation-driven mixing (1=yes)
mfshconv (max_dom) = 1, ! whether to turn on new day-time EDMF QNSE (0=no)
topo_wind (max_dom) = 0, ! turn off,
= 1, turn on topographic surface wind correction from Jimenez
(YSU PBL only, and require extra input from geogrid)
= 2, turn on topographic surface wind correction from Mass (YSU PBL only)
tke_budget (max_dom) = 0, ! default off; = 1 adds MYNN tke budget terms to output (new in 4.5)
(replacing bl_mynn_tkebudget in prior versions)
bl_mynn_tkeadvect (max_dom) = .false., ! default off; = .true. do MYNN tke advection
icloud_bl option to couple the subgrid-scale clouds from the PBL scheme (MYNN only)
to radiation schemes
0: no coupling; 1: activate coupling to radiation (default)
bl_mynn_cloudmix (max_dom) = 0 ! default off; =1 activates mixing of qc and qi in MYNN
0: no mixing of qc & qi; 1: mixing activated (default).
Note qnc and qni are mixed when scalar_pblmix =1.
bl_mynn_mixlength option to change mixing length formulation in MYNN
0: original as in Nakanishi and Niino 2009,
1: RAP/HRRR (including BouLac in free atmosphere),
2: experimental (default; includes cloud-specific mixing length and a scale-aware mixing
length, following Ito et al. 2015, BLM). Option 2 has been well tested with
the edmf options.
bl_mynn_cloudpdf option to switch to different cloud PDFs to represent subgrid clouds
0: original (Sommeria and Deardorf 1977);
1: Kuwano et al 2010, similar to option 0, but uses resolved scale gradients
as opposed to higher order moments ;
2: from Chaboureau and Bechtold (2002, JAS, with mods, default)
bl_mynn_edmf (max_dom) = 1 ! activate mass-flux option in MYNN, 0: mass-flux option off
Related (hidden) options:
bl_mynn_edmf_mom (max_dom) = 1 ! - activates momentum transport in MYNN mass-flux scheme
(assuming bl_mynn_edmf > 0); 0=off (default)
0 no momentum transport; 1: momentum transport activated (default)
bl_mynn_edmf_tke (max_dom) = 0, ! default; = 1 activates TKE transport in MYNN mass-flux scheme
(assuming bl_mynn_edmf > 0)
0 no TKE transport (default);1: activate TKE transport
bl_mynn_output (max_dom) = 0, ! do not output extra arrays
1 allocate and output extra 3D arrays from MYNN PBL
bl_mynn_mixscalars (max_dom) = 0 ! off, 1: activate mixing of scalars
bl_mynn_mixqt (max_dom) = 0 ! mixing moisture species separately, 1: mixing total water
scalar_pblmix (max_dom) = 1 ! mix scalar fields consistent with PBL option (exch_h)
tracer_pblmix (max_dom) = 1 ! mix tracer fields consistent with PBL option (exch_h)
shinhong_tke_diag (max_dom) = 0 ! diagnostic TKE and mixing length from Shin-Hong PBL
opt_thcnd option to treat thermal conductivity in Noah LSM
= 1, original (default)
= 2, McCumber and Pielke for silt loam and sandy loam
sf_surface_mosaic option to mosaic landuse categories for Noah LSM
= 0 ! default; use dominant category only
= 1 ! use mosaic landuse categories
mosaic_cat = 3 ! number of mosaic landuse categories in a grid cell
mosaic_lu = 0, ! default ; set to = 1 to use mosaic landuse categories in RUC
mosaic_soil = 0, ! default ; set to = 1 to use mosaic soil categories in RUC
flag_sm_adj = 0, ! default ; set to = 1 to adjust soil moisture for RUC in real.exe
(e.g. in case RUC is initialized from GFS LSM data)
cu_physics (max_dom) cumulus option
= 0, no cumulus
= 1, Kain-Fritsch (new Eta) scheme
= 2, Betts-Miller-Janjic scheme
= 3, Grell-Freitas ensemble scheme
= 4, Scale-aware GFS Simplified Arakawa-Schubert (SAS) scheme
= 5, Grell 3D ensemble scheme
= 6, Modifed Tiedtke scheme
= 7, Zhang-McFarlane scheme from CAM5 (CESM 1_0_1)
= 10, Modified Kain-Fritsch scheme with trigger function based on PDFs
= 11, Multi-scale Kain-Fritsch scheme
= 14, KIM Simplified Arakawa-Schubert scheme (KSAS) across gray-zone resolutions
= 16, A newer Tiedtke scheme
= 94, 2015 GFS Simplified Arakawa-Schubert scheme (from HWRF, use with caution)
= 95, Previous GFS Simplified Arakawa-Schubert scheme (from HWRF, use with caution)
= 96, Previous NEW GFS simplified Arakawa-Schubert scheme from YSU
= 93, Grell-Devenyi ensemble scheme
= 99, previous Kain-Fritsch scheme
shcu_physics (max_dom) independent shallow cumulus option (not tied to deep convection)
= 0, no independent shallow cumulus
= 1, Grell 3D ensemble scheme (use with cu_physics=93 or 5) (PLACEHOLDER: SWITCH NOT YET IMPLEMENTED--use ishallow)
= 2, Park and Bretherton shallow cumulus from CAM5 (CESM 1_0_1)
= 3, GRIMS shallow cumulus from YSU group
= 4, NSAS shallow cululus scheme. Must be used with KIAPS SAS cumulus scheme (cu_physics = 14)
= 5, Deng cumulus scheme (including both shallow and deep convections) from PSU and WRF-Solar.
However the deep part isn't very active. Not recommended to use alone for deep convection
case. Could work well for grid sizes 3-9 km.
This scheme only works with MYJ PBL scheme, and should not be combined with icloud=3
This scheme can also work with MYNN PBL scheme, but one should turn off EDMF (bl_mynn_edmf=0)
ishallow = 0, ! = 1 turns on shallow convection, used with Grell 3D ensemble schemes (cu_physics = 3 or 5)
clos_choice = 0, ! closure choice (place holder only)
cu_diag (max_dom) = 0, ! additional t-averaged stuff for cu physics (cu_phys = 3, 5, 10, and 93 only)
kf_edrates (max_dom) = 0, ! Add entrainment/detrainment rates and convective timescale output variables for KF-based
cumulus schemes (cu_phys = 1, 11 and 99 only)
= 0, ! no output; = 1, additional output
convtrans_avglen_m = 30, ! averaging time for variables used by convective transport (call cu_phys options) and radiation routines (only cu_phys=3,5 and 93) (minutes)
cu_rad_feedback (max_dom) = .false. ! sub-grid cloud effect to the optical depth in radiation
currently it works only for GF, G3, GD and KF scheme
One also needs to set cu_diag = 1 for GF, G3, KF-CuP, and GD schemes
bmj_rad_feedback (max_dom) = .false. ! BMJ convective precipitation-derived sub-grid cloud effect to radiation
cudt (max_dom) = 0, ! minutes between cumulus physics calls
kfeta_trigger KF trigger option (cu_physics=1 only):
= 1, default option
= 2, moisture-advection based trigger (Ma and Tan [2009])
= 3, RH-dependent additional perturbation to option 1 (JMA)
cugd_avedx ; number of grid boxes over which subsidence is spread.
= 1, default, for large grid distances
= 3, for small grid distances (DX < 5 km)
nsas_dx_factor = 0, ! default option
= 1, NSAS grid-distance dependent option
For KF-CuP scheme: recommended to use with cu_rad_feedback
shallowcu_forced_ra(max_dom) radiative impact of shallow Cu by a prescribed maximum cloud fraction
= .false., option off, default
= .true., radiative impact of shallow cu with a cloud fraction value of 0.36
numbins(max_dom) number of perturbations for potential temperature and mixing ratio in the CuP PDF,
should be an odd number (21 is a recommended value)
thBinSize(max_dom) bin size of potential temperature perturbation increment (0.01 K)
rBinSize(max_dom) bin size of mixing ratio perturbation increment (1.0e-4 kg/kg)
minDeepFreq(max_dom) minimum frequency required before deep convection is allowed (0.333)
minShallowFreq(max_dom) minimum frequency required before shallow convection is allowed (1.0e-2)
shcu_aerosols_opt(max_dom) whether aerosols in shcu: 0=none, 2=prognostic (run with WRF-Chem),
For MSKF scheme:
aercu_opt option to control aerosol interaction in MSKF and Morrison microphysics (mp_physics=40)
= 0: no aerosol interaction
= 1: aerosol interaction with MSKF only
= 2: aerosol interaction with both MSKF and Morrison
! aercu_opt = 1 and 2 both require the "CESM_RCP4.5_Aerosol_Data.dat" static runtime data
file. This can be downloaded from:
http:https://www2.mmm.ucar.edu/wrf/src/wrf_files/CESM_RCP4.5_Aerosol_Data.tar.gz
Once unpacked link/copy either of the two files to CESM_RCP4.5_Aerosol_Data.dat in
your run directory prior to running wrf.exe
aercu_fct factor to multiply with aerosol amount
= 1. ! default value
no_src_types_cu = 1 ! number of aerosol species in global aerosol data: 10 for CESM input,
set automatically
alevsiz_cu = 1 ! number of levels in global aerosol data: 30 for CESM input,
set automatically
isfflx = 1, ! heat and moisture fluxes from the surface
(only works for sf_sfclay_physics = 1,5,7,11)
1 = with fluxes from the surface
0 = no flux from the surface
with bl_pbl_physics=0 this uses tke_drag_coefficient
and tke_heat_flux in vertical diffusion
2 = use drag from sf_sfclay_physics and heat flux from
tke_heat_flux with bl_pbl_physics=0
ideal_xland = 1, ! sets XLAND (1=land,2=water) for ideal cases with no input land-use
run-time switch for wrf.exe physics_init (default 1 as before)
ifsnow = 1, ! snow-cover effects
(only works for sf_surface_physics = 1)
1 = with snow-cover effect
0 = without snow-cover effect
icloud = 1, ! cloud effect to the optical depth in radiation
(only works for ra_sw_physics = 1,4 and ra_lw_physics = 1,4)
Since 3.6, this also controls the cloud fraction options
1 = with cloud effect, and use cloud fraction option 1
(Xu-Randall method)
0 = without cloud effect
2 = with cloud effect, and use cloud fraction option 2 (0/1 based
on threshold
3 = with cloud effect, and use cloud fraction option 3, based on
Sundqvist et al. (1989) (since 3.7)
insert_init_cloud = .false., ! Default
= .true., ! For use with cold start (zero cloud fields), this option will
use the icloud=3 scheme to produce initial cloud water and cloud
ice fields as well as ensure initial water vapor field matches
saturation, which retains clouds beyond the first few timesteps.
swrad_scat = 1. ! scattering tuning parameter (default 1. is 1.e-5 m2/kg)
(works for ra_sw_physics = 1 option only)
surface_input_source = 3, ! where landuse and soil category data come from:
1 = WPS/geogrid but with dominant categories recomputed
2 = GRIB data from another model (only possible
(VEGCAT/SOILCAT are in met_em files from WPS)
3 = use dominant land and soil categories from WPS/geogrid (default since 3.8)
num_soil_layers = 5, ! number of soil layers in land surface model
= 5: thermal diffusion scheme
= 4: Noah landsurface model
= 6 or 9: RUC landsurface model
= 10: CLM4 landsurface model
= 2: Pleim-Xu landsurface model
= 3: SSiB landsurface model
num_land_cat = 21, ! number of land categories in input data.
24 - for USGS (default); 20 for MODIS
28 - for USGS if including lake category
21 - for MODIS if including lake category (default since 3.8)
40 - for NCLD
use_wudapt_lcz = 0, ! Option to use WUDAPT LCZ urban landuse categories, along with standard
urban classes (31-33).
= 0: use the traditional 33 landuse classes (31-33 for urban)
= 1: use the WUDAPT LCZ urban landuse categories
* Note: If the number of urban category in the input files is
inconsistent with the namelist option, error messages will occur.
The method to create the LCZ data is described here: http:https://www.wudapt.org/
slucm_distributed_drag = .false. ! option to use spatially varying 2-D urban Zero-plane Displacement, Roughness length for momentum, Frontal area index
! currently does not work with LCZ, only works with single-layer urban physics (urban_physics=1)
! need additional aforementioned 3 input variables in wrfinput file
distributed_ahe_opt = 0, ! option to handle anthropogenic surface heat flux (need additional input in wrfinput file)
= 0: no anthropogenic surface heat flux from input data
= 1: add to first level temperature tendency
= 2: add to surface sensible heat flux
num_soil_cat = 16, ! number of soil categories in input data
pxlsm_smois_init(max_dom) = 1 ! PXLSM Soil moisture initialization option
0 - From analysis, 1 - From moisture availability
or SLMO in LANDUSE.TBL
pxlsm_modis_veg = 1, ! PX LSM LAI and VEGFRA
0 - Old PX method that uses PX landuse look-up table
1 - Use VEGFRA and LAI in wrflowinp_d0* file
Note: Values used are called VEGF_PX and LAI_PX in output.
maxiens = 1, ! Grell-Devenyi only
maxens = 3, ! G-D only
maxens2 = 3, ! G-D only
maxens3 = 16 ! G-D only
ensdim = 144 ! G-D only
These are recommended numbers. If you would like to use
any other number, consult the code, know what you are doing.
seaice_threshold = 100. ! tsk < seaice_threshold, if water point and 5-layer slab
scheme, set to land point and permanent ice; if water point
and Noah scheme, set to land point, permanent ice, set temps
from 2 m to surface, and set smois and sh2o. The default value has changed
from 271 to 100 K in v3.5.1 to avoid mixed-up use with fractional seaice input
Used by land model option 1,2,3,4 and 8
sst_update = 0 ! time-varying sea-surface temp (0=no, 1=yes). If selected real
puts SST, XICE, ALBEDO and VEGFRA in wrflowinp_d01 file, and wrf updates
these from it at same interval as boundary file. Also requires
namelists in &time_control: auxinput4_interval, auxinput4_end_h,
auxinput4_inname = "wrflowinp_d<domain>", and io_form_auxinput4
usemonalb = .false. ! use monthly albedo map instead of table value
(recommended for sst_update=1)
rdmaxalb = .true. ! use snow albedo from geogrid; false means using values from table
rdlai2d = .false. ! use LAI from input; false means using values from table
if sst_update=1, LAI will also be in wrflowinp file
dust_emis = 0 ! Enable (0=no, 1=yes) surface dust emission scheme to enter mp_physics=28 QNIFA (ice-friendly aerosol variable)
erosion_dim = 3 ! In conjunction with dust_emis=1, this value can only be set equal to 3 (erodibility information)
bucket_mm = -1. ! bucket reset value for water accumulations (value in mm, -1.=inactive)
bucket_J = -1. ! bucket reset value for energy accumulations (value in J, -1.=inactive)
tmn_update = 0 ! update deep soil temperature (1, yes; 0, no)