-
Notifications
You must be signed in to change notification settings - Fork 2
/
Registry.EM_COMMON
3389 lines (3106 loc) · 391 KB
/
Registry.EM_COMMON
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
# At the present time this file is managed manually and edited by hand.
#
################################################################################
# Dimension specifications
#
# This section of the Registry file is used to specify the dimensions
# that will be used to define arrays. Dim is the one-letter name of the
# dimension. How defined can either be "standard_domain", which means
# that the dimension (1) is one of the three spatial dimensions and (2)
# it will be set using the standard namelist mechanism and domain data
# structure dimension fields (e.g. sd31,ed31,sd32...).
#
# Order refers to which of the three sets of just-mentioned internal
# dimension variables the dimension is referred to by in the driver.
# That is, is it the first, second, or third dimension. The registry
# infers the mapping of its internal dimensions according to the
# combination of Order and Coord-axis that are specified in this table.
# Note that it is all right to more than one dimension name for, say, the
# x dimension. However, the Order and Coord-axis relationship must be
# consistent throughout.
#
# Note: these entries do not enforce storage order on a particular field.
# That is determined by the dimension strings for each field. But it does
# relate the dimspec to the internal data structures that the driver uses
# to maintain the three physical domain dimensions.
#
# "How defined" can also specify the name of a namelist variable from which
# the definition for the dimension will come; this is specified as
# "namelist=<variable name>". The namelist variable must have been
# defined as an integer and with only one entry in the rconfig table. Or
# a constant can be specified. The coordinate axis for the dimension is
# either X, Y, Z, or C (for "not a spatial dimension"). The Dimname is
# the descriptive name of the dimension that will be included in the
# metadata in data sets. Note that the b, f, and t modifiers that appear
# as the last characters of dimension strings used # in state and # i1
# registry definitions are not dimensions and do not need to be declared
# here.
#
################################################################################
################################################################################
################################################################################
# Lines that start with the word 'state' form a table that is
# used by the script use_registry to generate module_state_descript.F
# and other files. Also see documentation in use_registry.
#
# table entries are of the form
#<Table> <Type> <Sym> <Dims> <Use> <NumTLev> <Stagger> <IO> <DNAME> <DESCRIP> <UNITS>
#
state real XLAT ij misc 1 - i0123rh0156{22}{23}du=(copy_fcnm) "XLAT" "LATITUDE, SOUTH IS NEGATIVE" "degree_north"
state real XLONG ij misc 1 - i0123rh0156{22}{23}d=(interp_fcn_blint_ll:xlat,input_from_file)u=(copy_fcnm) "XLONG" "LONGITUDE, WEST IS NEGATIVE" "degree_east"
# It is required that LU_INDEX appears before any variable that is
# interpolated with a mask, as lu_index supplies that mask.
# this next 1 is for the HFSoLE/PET demo; writing these to auxhist1 output over MCEL for coupling
# with wave model, only if compiled with -DMCELIO, JM 2003/05/29
state real LU_INDEX ij misc 1 - i012rh01d=(interp_fcnm_lu:xlat,xlong,dx,grid_id)u=(copy_fcnm) "LU_INDEX" "LAND USE CATEGORY" ""
state real LU_MASK ij misc 1 - i3h1 "LU_MASK" "0 land 1 water" ""
# znw, znu, dzs, and zs must be listed before any 3-d fields
# in order for the grib output module to work correctly. The grib output
# module retrieves the vertical levels from these parameters. If znw, znu
# dzs, and zs are not listed first, vertical level will not be encoded at
# time 0.
state real znu k dyn_em 1 - irh "znu" "eta values on half (mass) levels" ""
state real znw k dyn_em 1 Z i0rh "znw" "eta values on full (w) levels" ""
state real ZS l misc - Z irh "ZS" "DEPTHS OF CENTERS OF SOIL LAYERS" "m"
state real DZS l misc - Z irh "DZS" "THICKNESSES OF SOIL LAYERS" "m"
# cyl : variables for trajectory calculation
state real traj_i {ntraj} dyn_em 1 Z irh "traj_i" "grid number of trajectory " ""
state real traj_j {ntraj} dyn_em 1 Z irh "traj_j" "grid number of trajectory " ""
state real traj_k {ntraj} dyn_em 1 Z irh "traj_k" "vertical eta level of trajectory " ""
state real traj_long {ntraj} dyn_em 1 Z irh "traj_long" "longitude of trajectory " ""
state real traj_lat {ntraj} dyn_em 1 Z irh "traj_lat" "latitude of trajectory " ""
#
# Variables from WPS
#
state real u_gc igj dyn_em 1 XZ i1 "UU" "x-wind component" "m s-1"
state real v_gc igj dyn_em 1 YZ i1 "VV" "y-wind component" "m s-1"
state real t_gc igj dyn_em 1 Z i1 "TT" "temperature" "K"
state real rh_gc igj dyn_em 1 Z i1 "RH" "relative humidity" "%"
state real ght_gc igj dyn_em 1 Z i1 "GHT" "geopotential height" "m"
state real p_gc igj dyn_em 1 Z i1 "PRES" "pressure" "Pa"
state real prho_gc igj dyn_em 1 Z i1 "PRHO" "for UM data, pressure of U and V" "Pa"
state real xlat_gc ij dyn_em 1 - i1 "XLAT_M" "latitude, positive north" "degrees"
state real xlong_gc ij dyn_em 1 - i1 "XLONG_M" "longitude, positive east" "degrees"
state real ht_gc ij dyn_em 1 - i1 "HGT_M" "topography elevation" "m"
state real var_sso ij dyn_em 1 - i012hrd "var_sso" "variance of subgrid-scale orography" "m2"
state real lap_hgt ij dyn_em 1 - r "lap_hgt" "Laplacian of orography" "m"
state real tsk_gc ij dyn_em 1 - i1 "SKINTEMP" "skin temperature" "K"
state real tavgsfc ij dyn_em 1 - i1 "TAVGSFC" "daily mean of surface air temperature" "K"
state real tmn_gc ij dyn_em 1 - i1 "SOILTEMP" "annual mean deep soil temperature" "K"
state real pslv_gc ij dyn_em 1 - i1 "PMSL" "sea level pressure" "Pa"
state real sct_dom_gc ij dyn_em 1 - i1 "SCT_DOM" "Dominant soil (top) category from GEOGRID" "cat"
state real scb_dom_gc ij dyn_em 1 - i1 "SCB_DOM" "Dominant soil (bottom) category from GEOGRID" "cat"
state real greenfrac imj dyn_em 1 Z i1 "GREENFRAC" "monthly greenness fraction" "0 - 1 fraction"
state real albedo12m imj dyn_em 1 Z i1 "ALBEDO12M" "background albedo" "0 - 1 fraction"
state real lai12m imj dyn_em 1 Z i1 "LAI12M" "monthly LAI" "m2/m2"
state real pd_gc igj dyn_em 1 Z - "PD" "dry pressure" "Pa"
state real pdrho_gc igj dyn_em 1 Z - "PDRHO" "dry pressure for UM data for the variables U and V" "Pa"
state real psfc_gc ij dyn_em 1 - - "PSFC_GC" "surface pressure" "Pa"
state real intq_gc ij dyn_em 1 - - "INTQ" "integrated mixing ratio" "Pa"
state real pdhs ij dyn_em 1 - - "PDHS" "hydrostatic dry surface pressure" "Pa"
state real qv_gc igj dyn_em 1 Z i1 "QV" "mixing ratio" "kg kg-1"
state real sh_gc igj dyn_em 1 Z i1 "SPECHUMD" "Specific humidity" "kg kg-1"
state real cl_gc igj dyn_em 1 Z i1 "SPECCLDL" "Cloud water content, liquid" "kg kg-1"
state real cf_gc igj dyn_em 1 Z i1 "SPECCLDF" "Cloud water content, frozen" "kg kg-1"
state real icefrac_gc ij dyn_em 1 - i1 "ICEFRAC" "Sea ice fraction" "0 - 1 fraction"
state real icepct ij dyn_em 1 - i1 "ICEPCT" "Sea ice percent" "%"
state real qr_gc igj dyn_em 1 Z i1 "QR" "rain water mixing ratio" "kg kg-1"
state real qc_gc igj dyn_em 1 Z i1 "QC" "cloud water mixing ratio" "kg kg-1"
state real qs_gc igj dyn_em 1 Z i1 "QS" "snow mixing ratio" "kg kg-1"
state real qi_gc igj dyn_em 1 Z i1 "QI" "cloud ice mixing ratio" "kg kg-1"
state real qg_gc igj dyn_em 1 Z i1 "QG" "graupel mixing ratio" "kg kg-1"
state real qh_gc igj dyn_em 1 Z i1 "QH" "hail mixing ratio" "kg kg-1"
state real qni_gc igj dyn_em 1 Z i1 "QNI" "ice num concentration" "# kg-1"
state real qnc_gc igj dyn_em 1 Z i1 "QNC" "cloud water num concentration" "# kg-1"
state real qnr_gc igj dyn_em 1 Z i1 "QNR" "rain num concentration" "# kg-1"
state real qns_gc igj dyn_em 1 Z i1 "QNS" "snow num concentration" "# kg-1"
state real qng_gc igj dyn_em 1 Z i1 "QNG" "graupel num concentration" "# kg-1"
state real qnh_gc igj dyn_em 1 Z i1 "QNH" "hail num concentration" "# kg-1"
state real qnwfa_gc igj dyn_em 1 Z i1 "QNWFA" "water-friendly aerosol num concentration" "# kg-1"
state real qnifa_gc igj dyn_em 1 Z i1 "QNIFA" "water-friendly aerosol num concentration" "# kg-1"
state real qnwfa_now igj dyn_em 1 Z - "QNWFA_NOW" "num water-friendly aerosol Now" "# kg-1"
state real qnwfa_jan igj dyn_em 1 Z i1 "QNWFA_JAN" "num water-friendly aerosol Jan" "# kg-1"
state real qnwfa_feb igj dyn_em 1 Z i1 "QNWFA_FEB" "num water-friendly aerosol Feb" "# kg-1"
state real qnwfa_mar igj dyn_em 1 Z i1 "QNWFA_MAR" "num water-friendly aerosol Mar" "# kg-1"
state real qnwfa_apr igj dyn_em 1 Z i1 "QNWFA_APR" "num water-friendly aerosol Apr" "# kg-1"
state real qnwfa_may igj dyn_em 1 Z i1 "QNWFA_MAY" "num water-friendly aerosol May" "# kg-1"
state real qnwfa_jun igj dyn_em 1 Z i1 "QNWFA_JUN" "num water-friendly aerosol Jun" "# kg-1"
state real qnwfa_jul igj dyn_em 1 Z i1 "QNWFA_JUL" "num water-friendly aerosol Jul" "# kg-1"
state real qnwfa_aug igj dyn_em 1 Z i1 "QNWFA_AUG" "num water-friendly aerosol Aug" "# kg-1"
state real qnwfa_sep igj dyn_em 1 Z i1 "QNWFA_SEP" "num water-friendly aerosol Sep" "# kg-1"
state real qnwfa_oct igj dyn_em 1 Z i1 "QNWFA_OCT" "num water-friendly aerosol Oct" "# kg-1"
state real qnwfa_nov igj dyn_em 1 Z i1 "QNWFA_NOV" "num water-friendly aerosol Nov" "# kg-1"
state real qnwfa_dec igj dyn_em 1 Z i1 "QNWFA_DEC" "num water-friendly aerosol Dec" "# kg-1"
state real qnifa_now igj dyn_em 1 Z - "QNIFA_NOW" "num ice-friendly aerosol Now" "# kg-1"
state real qnifa_jan igj dyn_em 1 Z i1 "QNIFA_JAN" "num ice-friendly aerosol Jan" "# kg-1"
state real qnifa_feb igj dyn_em 1 Z i1 "QNIFA_FEB" "num ice-friendly aerosol Feb" "# kg-1"
state real qnifa_mar igj dyn_em 1 Z i1 "QNIFA_MAR" "num ice-friendly aerosol Mar" "# kg-1"
state real qnifa_apr igj dyn_em 1 Z i1 "QNIFA_APR" "num ice-friendly aerosol Apr" "# kg-1"
state real qnifa_may igj dyn_em 1 Z i1 "QNIFA_MAY" "num ice-friendly aerosol May" "# kg-1"
state real qnifa_jun igj dyn_em 1 Z i1 "QNIFA_JUN" "num ice-friendly aerosol Jun" "# kg-1"
state real qnifa_jul igj dyn_em 1 Z i1 "QNIFA_JUL" "num ice-friendly aerosol Jul" "# kg-1"
state real qnifa_aug igj dyn_em 1 Z i1 "QNIFA_AUG" "num ice-friendly aerosol Aug" "# kg-1"
state real qnifa_sep igj dyn_em 1 Z i1 "QNIFA_SEP" "num ice-friendly aerosol Sep" "# kg-1"
state real qnifa_oct igj dyn_em 1 Z i1 "QNIFA_OCT" "num ice-friendly aerosol Oct" "# kg-1"
state real qnifa_nov igj dyn_em 1 Z i1 "QNIFA_NOV" "num ice-friendly aerosol Nov" "# kg-1"
state real qnifa_dec igj dyn_em 1 Z i1 "QNIFA_DEC" "num ice-friendly aerosol Dec" "# kg-1"
state real qntemp imj dyn_em 1 Z - "QNTEMP" "temporary var for time interp" ""
state real qntemp2 ij dyn_em 1 - - "QNTEMP2" "temporary var2D for time interp" ""
state real t_max_p ij dyn_em 1 - i0d "T_MAX_P" "temperature at max pressure" "K"
state real ght_max_p ij dyn_em 1 - i0d "GHT_MAX_P" "geopotential height at max pressure" "m"
state real max_p ij dyn_em 1 - i0d "MAX_P" "max pressure " "Pa"
state real t_min_p ij dyn_em 1 - i0d "T_MIN_P" "temperature at min pressure" "K"
state real ght_min_p ij dyn_em 1 - i0d "GHT_MIN_P" "geopotential height at min pressure" "m"
state real min_p ij dyn_em 1 - i0d "MIN_P" "min pressure " "Pa"
state real hgtmaxw ij dyn_em 1 - i1 "HGTMAXW" "Height of the max wind speed" "m"
state real hgttrop ij dyn_em 1 - i1 "HGTTROP" "Height of the tropopause" "m"
state real pmaxw ij dyn_em 1 - i1 "PMAXW" "Pressure of the max wind speed" "Pa"
state real pmaxwnn ij dyn_em 1 - i1 "PMAXWNN" "PMAXW, nearest neighbor interp" "Pa"
state real ptrop ij dyn_em 1 - i1 "PTROP" "Pressure of the tropopause" "Pa"
state real ptropnn ij dyn_em 1 - i1 "PTROPNN" "PTROP, nearest neighbor interp" "Pa"
state real tmaxw ij dyn_em 1 - i1 "TMAXW" "Temperature of the max wind speed" "K"
state real ttrop ij dyn_em 1 - i1 "TTROP" "Temperature of the tropopause" "K"
state real umaxw ij dyn_em 1 X i1 "UMAXW" "U-component of the max wind speed" "m s-1"
state real utrop ij dyn_em 1 X i1 "UTROP" "U-component of the tropopause wind" "m s-1"
state real vmaxw ij dyn_em 1 Y i1 "VMAXW" "V-component of the max wind speed" "m s-1"
state real vtrop ij dyn_em 1 Y i1 "VTROP" "V-component of the tropopause wind" "m s-1"
state real erod ij. misc 1 - i012rd "EROD" "fraction of erodible surface in each grid cell (0-1)" "none"
#-----------------------------------------------------------------------------------------------------------------------------------------------------------------
#
# Variables for Eulerian mass coordinate dynamics
#
# Velocities
#
# U Vel
state real u ikjb dyn_em 2 X \
i0rhusdf=(bdy_interp:dt) "U" "x-wind component" "m s-1"
state real ru ikj dyn_em 1 X - "MU_U" "mu-coupled u" "Pa m s-1"
state real ru_m ikj dyn_em 1 X - "ru_m" "" ""
state real ru_tend ikj dyn_em 1 X - "ru_tend" "" ""
i1 real ru_tendf ikj dyn_em 1 X
state real u_save ikj dyn_em 1 X - "u_save"
state real z_force | dyn_em 1 - i3rh "Z_FORCE" "height of forcing input" "m"
state real z_force_tend | dyn_em 1 - i3rh "Z_FORCE_TEND" "tendency height of forcing input" "m"
state real u_g | dyn_em 1 - i3rh "U_G" "x-direction geostrophic wind" "m s-1"
state real u_g_tend | dyn_em 1 - i3rh "U_G_TEND" "tendency x-direction geostrophic wind" "m s-1"
#
# V Vel
state real v ikjb dyn_em 2 Y \
i0rhusdf=(bdy_interp:dt) "V" "y-wind component" "m s-1"
state real rv ikj dyn_em 1 Y - "MU_V" "mu-coupled v" "Pa m s-1"
state real rv_m ikj dyn_em 1 Y - "rv_m"
state real rv_tend ikj dyn_em 1 Y - "rv_tend"
i1 real rv_tendf ikj dyn_em 1 Y
state real v_save ikj dyn_em 1 Y - "v_save"
state real v_g | dyn_em 1 - i3rh "V_G" "y-direction geostrophic wind" "m s-1"
state real v_g_tend | dyn_em 1 - i3rh "V_G_TEND" "tendency y-direction geostrophic wind" "m s-1"
#
# Vertical Vel
state real w ikjb dyn_em 2 Z \
irhusdf=(bdy_interp:dt) "w" "z-wind component" "m s-1"
state real ww ikj dyn_em 1 Z r "ww" "mu-coupled eta-dot" "Pa s-1"
state real rw ikj dyn_em 1 Z - "rw" "mu-coupled w" "Pa m s-1"
i1 real ww1 ikj dyn_em 1 Z
state real ww_m ikj dyn_em 1 Z r "ww_m" "time-avg mu-coupled eta-dot" "Pa s-1"
i1 real wwp ikj dyn_em 1 Z
i1 real rw_tend ikj dyn_em 1 Z
i1 real rw_tendf ikj dyn_em 1 Z
i1 real w_save ikj dyn_em 1 Z
state real w_subs | dyn_em 1 - i3rh "W_SUBS" "large-scale vertical velocity" "m s-1"
state real w_subs_tend | dyn_em 1 - i3rh "W_SUBS_TEND" "tendency large-scale vertical velocity" "m s-1"
# Geopotential
state real ph ikjb dyn_em 2 Z \
irhusdf=(bdy_interp:dt) "ph" "perturbation geopotential" "m2 s-2"
state real phb ikj dyn_em 1 Z irhdus "phb" "base-state geopotential" "m2 s-2"
state real phb_fine ikj dyn_em 1 Z - "phb_fine" "for nesting, temp holding interpolated coarse grid phb" "m2 s-2"
state real ph0 ikj dyn_em 1 Z r "ph0" "initial geopotential" "m2 s-2"
state real php ikj dyn_em 1 - r "php" "geopotential" "m2 s-2"
i1 real ph_tend ikj dyn_em 1 Z
i1 real ph_tendf ikj dyn_em 1 Z
i1 real ph_save ikj dyn_em 1 Z
# Potential Temperature
state real th_phy_m_t0 ikj dyn_em 1 - irhd "t" "perturbation potential temperature theta-t0" "K"
state real t ikjb dyn_em 2 - \
i0rhusdf=(bdy_interp:dt) "thm" "either 1) pert moist pot temp=(1+Rv/Rd Qv)*(theta)-T0, or 2) pert dry pot temp=theta-T0; based on use_theta_m setting" "K"
state real t_init ikj dyn_em 1 - ird "t_init" "initial potential temperature" "K"
i1 real t_tend ikj dyn_em 1 -
i1 real t_tendf ikj dyn_em 1 -
i1 real t_2save ikj dyn_em 1 -
state real t_save ikj dyn_em 1 - - "t_save"
state real th_upstream_x | dyn_em 1 - i3rh "TH_UPSTREAM_X" "upstream theta x-advection" "K s-1"
state real th_upstream_x_tend | dyn_em 1 - i3rh "TH_UPSTREAM_X_TEND" "tendency upstream theta x-advection" "K s-2"
state real th_upstream_y | dyn_em 1 - i3rh "TH_UPSTREAM_Y" "upstream theta y-advection" "K s-1"
state real th_upstream_y_tend | dyn_em 1 - i3rh "TH_UPSTREAM_Y_TEND" "tendency upstream theta y-advection" "K s-2"
state real qv_upstream_x | dyn_em 1 - i3rh "QV_UPSTREAM_X" "upstream qv x-advection" "kg kg-1 s-1"
state real qv_upstream_x_tend | dyn_em 1 - i3rh "QV_UPSTREAM_X_TEND" "tendency upstream qv x-advection" "kg kg-1 s-2"
state real qv_upstream_y | dyn_em 1 - i3rh "QV_UPSTREAM_Y" "upstream qv y-advection" "kg kg-1 s-1"
state real qv_upstream_y_tend | dyn_em 1 - i3rh "QV_UPSTREAM_Y_TEND" "tendency upstream qv y-advection" "kg kg-1 s-2"
state real ql_upstream_x | dyn_em 1 - i3rh "QL_UPSTREAM_X" "upstream ql x-advection" "kg kg-1 s-1"
state real ql_upstream_x_tend | dyn_em 1 - i3rh "QL_UPSTREAM_X_TEND" "tendency upstream ql x-advection" "kg kg-1 s-2"
state real ql_upstream_y | dyn_em 1 - i3rh "QL_UPSTREAM_Y" "upstream ql y-advection" "kg kg-1 s-1"
state real ql_upstream_y_tend | dyn_em 1 - i3rh "QL_UPSTREAM_Y_TEND" "tendency upstream ql y-advection" "kg kg-1 s-2"
state real u_upstream_x | dyn_em 1 - i3rh "U_UPSTREAM_X" "upstream u x-advection" "m s-2"
state real u_upstream_x_tend | dyn_em 1 - i3rh "U_UPSTREAM_X_TEND" "tendency upstream u x-advection" "m s-3"
state real u_upstream_y | dyn_em 1 - i3rh "U_UPSTREAM_Y" "upstream u y-advection" "m s-2"
state real u_upstream_y_tend | dyn_em 1 - i3rh "U_UPSTREAM_Y_TEND" "tendency upstream u y-advection" "m s-3"
state real v_upstream_x | dyn_em 1 - i3rh "V_UPSTREAM_X" "upstream v x-advection" "m s-2"
state real v_upstream_x_tend | dyn_em 1 - i3rh "V_UPSTREAM_X_TEND" "tendency upstream v x-advection" "m s-3"
state real v_upstream_y | dyn_em 1 - i3rh "V_UPSTREAM_Y" "upstream v y-advection" "m s-2"
state real v_upstream_y_tend | dyn_em 1 - i3rh "V_UPSTREAM_Y_TEND" "tendency upstream v y-advection" "m s-3"
state real th_t_tend | dyn_em 1 - i3rh "TH_T_TEND" "tendency theta time" "K s-2"
state real qv_t_tend | dyn_em 1 - i3rh "QV_T_TEND" "tendency qv time" "kg kg-1 s-1"
state real th_largescale | dyn_em 1 - i3rh "TH_LARGESCALE" "SCM largescale theta" "K"
state real th_largescale_tend | dyn_em 1 - i3rh "TH_LARGESCALE_TEND" "SCM tendency largescale theta" "K s-1"
state real qv_largescale | dyn_em 1 - i3rh "QV_LARGESCALE" "SCM largescale qv" "kg kg-1 s-1"
state real qv_largescale_tend | dyn_em 1 - i3rh "QV_LARGESCALE_TEND" "SCM tendency largescale qv" "kg kg-1 s-2"
state real ql_largescale | dyn_em 1 - i3rh "QL_LARGESCALE" "SCM largescale ql" "kg kg-1 s-1"
state real ql_largescale_tend | dyn_em 1 - i3rh "QL_LARGESCALE_TEND" "SCM tendency largescale ql" "kg kg-1 s-2"
state real u_largescale | dyn_em 1 - i3rh "U_LARGESCALE" "SCM largescale u" "m s-2"
state real u_largescale_tend | dyn_em 1 - i3rh "U_LARGESCALE_TEND" "SCM tendency largescale u" "m s-3"
state real v_largescale | dyn_em 1 - i3rh "V_LARGESCALE" "SCM largescale v" "m s-2"
state real v_largescale_tend | dyn_em 1 - i3rh "V_LARGESCALE_TEND" "SCM tendency largescale v" "m s-3"
state real tau_largescale | dyn_em 1 - i3rh "TAU_LARGESCALE" "SCM largescale timescale" "s"
state real tau_largescale_tend | dyn_em 1 - i3rh "TAU_LARGESCALE_TEND" "SCM tendency largescale timescale" ""
state real tau_x | dyn_em 1 - i3rh "TAU_X" "X-direction advective timescale" "s"
state real tau_x_tend | dyn_em 1 - i3rh "TAU_X_TEND" "tendency X-direction advective timescale" ""
state real tau_y | dyn_em 1 - i3rh "TAU_Y" "Y-direction advective timescale" "s"
state real tau_y_tend | dyn_em 1 - i3rh "TAU_Y_TEND" "tendency Y-direction advective timescale" ""
dimspec fslay 2 namelist=num_force_soil_layers z force_soil_layers
state real t_soil_forcing_val {fslay} dyn_em 1 - i3rh "T_SOIL_FORCING_VAL" "Soil temp value for SCM forcing" "K"
state real t_soil_forcing_tend {fslay} dyn_em 1 - i3rh "T_SOIL_FORCING_TEND" "tendency soil temp for SCM forcing" "K s-1"
state real q_soil_forcing_val {fslay} dyn_em 1 - i3rh "Q_SOIL_FORCING_VAL" "Soil moisture value for SCM forcing" "1"
state real q_soil_forcing_tend {fslay} dyn_em 1 - i3rh "Q_SOIL_FORCING_TEND" "tendency soil moisture for SCM forcing" "s-1"
state real tau_soil {fslay} dyn_em 1 - i3rh "TAU_SOIL" "SCM soil forcing timescale" "s"
state real soil_depth_force {fslay} dyn_em 1 - i3rh "SOIL_DEPTH_FORCE" "SCM depth at center of soil layers in forcing file" "1"
state real hfx_force - dyn_em 1 - i3rh "HFX_FORCE" "SCM ideal surface sensible heat flux" "W m-2"
state real lh_force - dyn_em 1 - i3rh "LH_FORCE" "SCM ideal surface latent heat flux" "W m-2"
state real tsk_force - dyn_em 1 - i3rh "TSK_FORCE" "SCM ideal surface skin temperature" "W m-2"
state real hfx_force_tend - dyn_em 1 - i3rh "HFX_FORCE_TEND" "SCM ideal surface sensible heat flux tendency" "W m-2 s-1"
state real lh_force_tend - dyn_em 1 - i3rh "LH_FORCE_TEND" "SCM ideal surface latent heat flux tendency" "W m-2 s-1"
state real tsk_force_tend - dyn_em 1 - i3rh "TSK_FORCE_TEND" "SCM ideal surface skin temperature tendency" "W m-2 s-1"
# Mass
state real mu ijb dyn_em 2 - \
irh01usdf=(bdy_interp:dt) "mu" "perturbation dry air mass in column" "Pa"
state real mub ij dyn_em 1 - irhdus "mub" "base state dry air mass in column" "Pa"
state real mub_fine ij dyn_em 1 - - "mub_fine" "nest temp, holds interpolated coarse grid mub" "Pa"
state real mub_save ij dyn_em 1 - - "mub_save" "nest temp, holds orig fine grid mub" "Pa"
state real mu0 ij dyn_em 1 - i1 "mu0" "initial dry mass in column" "Pa"
state real mudf ij dyn_em 1 - - "mudf" "" ""
state real muu ij dyn_em 1 X - "muu"
state real muus ij dyn_em 1 X - "muus"
state real muv ij dyn_em 1 Y - "muv"
state real muvs ij dyn_em 1 Y - "muvs"
state real mut ij dyn_em 1 - - "mut"
state real muts ij dyn_em 1 - - "muts"
i1 real muave ij dyn_em 1 -
i1 real mu_save ij dyn_em 1 -
i1 real mu_tend ij dyn_em 1 -
i1 real mu_tendf ij dyn_em 1 -
#diagnostic for looking at nest position in output. A mungy version of terrain height.
state real nest_pos ij misc 1 - rhu=(mark_domain) "NEST_POS"
state real nest_mask ij misc 1 - ru=(mark_domain) "NEST_MASK" "LOCATION OF NEST IF ANY"
state real ht_coarse ij misc 1 - r - "STORAGE FOR LOW-RES TERRAIN"
# TKE
state real tke ikj dyn_em 2 - r "tke" "TURBULENCE KINETIC ENERGY" "m2 s-2"
i1 real tke_tend ikj dyn_em 1 -
# Pressure and Density
state real p ikj dyn_em 1 - irh "p" "perturbation pressure" "Pa"
state real al ikj dyn_em 1 - ir "al" "inverse perturbation density" "m3 kg-1"
state real alt ikj dyn_em 1 - r "alt" "inverse density" "m3 kg-1"
state real alb ikj dyn_em 1 - irdus "alb" "inverse base density" "m3 kg-1"
state real zx ikj dyn_em 1 XZ - " " " " " "
state real zy ikj dyn_em 1 YZ - " " " " " "
state real rdz ikj dyn_em 1 Z - " " " " " "
state real rdzw ikj dyn_em 1 Z - " " " " " "
state real pb ikj dyn_em 1 - irhdus "pb" "BASE STATE PRESSURE " "Pa"
#
# Other dyn
#
i1 real advect_tend ikj dyn_em 1 -
i1 real alpha ikj dyn_em 1 -
i1 real a ikj dyn_em 1 -
i1 real gamma ikj dyn_em 1 -
i1 real c2a ikj dyn_em 1 - -
state real rho ikj misc 1 - r "RHO" "DENSITY" "Kg m-3"
i1 real phm ikj dyn_em 1 - -
i1 real cqu ikj dyn_em 1 - -
i1 real cqv ikj dyn_em 1 - -
i1 real cqw ikj dyn_em 1 - -
i1 real pm1 ikj dyn_em 1 - -
state real fnm k dyn_em 1 - irh "fnm" "upper weight for vertical stretching" ""
state real fnp k dyn_em 1 - irh "fnp" "lower weight for vertical stretching" ""
state real rdnw k dyn_em 1 - irh "rdnw" "inverse d(eta) values between full (w) levels" ""
state real rdn k dyn_em 1 - irh "rdn" "inverse d(eta) values between half (mass) levels" ""
state real dnw k dyn_em 1 - irh "dnw" "d(eta) values between full (w) levels" ""
state real dn k dyn_em 1 - irh "dn " "d(eta) values between half (mass) levels" ""
state real t_base k dyn_em 1 - ir "t_base" "BASE STATE T IN IDEALIZED CASES" "K"
state real z ikj dyn_em 1 - - " " " " " "
state real z_at_w ikj dyn_em 1 Z
state real cfn - misc - - irh "cfn" "extrapolation constant" ""
state real cfn1 - misc - - irh "cfn1" "extrapolation constant" ""
state integer step_number - misc - - ir "step_number" ""
# Idealized run
state logical this_is_an_ideal_run - misc - - irh "this_is_an_ideal_run" "T/F flag: this is an ARW ideal simulation"
# For the adaptive timestep restart
state logical stepping_to_time - misc - - r "stepping_to_time" ""
state integer last_step_updated - misc - - r "last_step_updated" ""
state logical adapt_step_using_child - misc - - r "adapt_step_using_child" ""
state integer last_dt_sec - misc - - r "last_dt_sec" "Whole seconds for last timestep" "sec"
state integer last_dt_sec_num - misc - - r "last_dt_sec_num" "Fractional secs, numerator" "sec"
state integer last_dt_sec_den - misc - - r "last_dt_sec_den" "Fractional secs, denominator" "sec"
state integer last_dt_yr - misc - - r "last_dt_yr" "Relative year" "years"
state integer last_dt_mm - misc - - r "last_dt_mm" "Relative month" "months"
# hydrostatic pressure vars
state real p_hyd ikj dyn_em 1 - irh "p_hyd" "hydrostatic pressure" "Pa"
state real p_hyd_w ikj dyn_em 1 Z r "p_hyd_w" "hydrostatic pressure at full levels" "Pa"
# 2m and 10m output diagnostics
state real Q2 ij misc 1 - irh01{22}{23}du "Q2" "QV at 2 M" "kg kg-1"
state real T2 ij misc 1 - i01rh01{22}{23}du "T2" "TEMP at 2 M" "K"
state real TH2 ij misc 1 - irhdu "TH2" "POT TEMP at 2 M" "K"
state real PSFC ij misc 1 - i01rh01du "PSFC" "SFC PRESSURE" "Pa"
# these next 2 are for the HFSoLE/PET demo; writing these to auxhist1 output over MCEL for coupling
# with wave model, only if compiled with -DMCELIO, JM 2003/05/29
state real U10 ij misc 1 - irh01{22}{23}du "U10" "U at 10 M" "m s-1"
state real V10 ij misc 1 - irh01{22}{23}du "V10" "V at 10 M" "m s-1"
# LPI
state real LPI ij misc 1 - rhdu "LPI" "Lightning Potential Index" "m^2 s-2"
# these next 4 are for observational nudging
state real uratx ij misc 1 - r "URATX" "Ratio of U over U10 on mass points " "dimensionless"
state real vratx ij misc 1 - r "VRATX" "Ratio of V over V10 on mass points " "dimensionless"
state real tratx ij misc 1 - r "TRATX" "Ratio of T over TH2 on mass points " "dimensionless"
state real obs_savwt hikj dyn_em 1 X - "OBS_SAVWT" "Internal space holding weights of each ob, for each i,j,k" "dimensionless"
# Other
state real rdx - misc - - irh "rdx" "INVERSE X GRID LENGTH" ""
state real rdy - misc - - irh "rdy" "INVERSE Y GRID LENGTH" ""
state real dts - misc - - ir "dts" "SMALL TIMESTEP" ""
state real dtseps - misc - - ir "dtseps" "TIME WEIGHT CONSTANT FOR SMALL STEPS" ""
state real resm - misc - - irh "resm" "TIME WEIGHT CONSTANT FOR SMALL STEPS" ""
state real zetatop - misc - - irh "zetatop" "ZETA AT MODEL TOP" ""
state real cf1 - misc - - irh "cf1" "2nd order extrapolation constant" ""
state real cf2 - misc - - irh "cf2" "2nd order extrapolation constant" ""
state real cf3 - misc - - irh "cf3" "2nd order extrapolation constant" ""
state integer number_at_same_level - - - - - "number_at_same_level" "" ""
state real radtacttime - - - - r "radtacttime" "RADTACTTIME" "LW SW ACTIVATION TIME in s"
state real bldtacttime - - - - r "bldtacttime" "BLDTACTTIME" "PBL ACTIVATION TIME in s"
state real cudtacttime - - - - r "cudtacttime" "CUDTACTTIME" "CPS ACTIVATION TIME in s"
state real ltngacttime - - - - r "ltngacttime" "LTNGACTTIME" "LTNG ACTIVATION TIME in s"
state real power ij misc 1 - irh "Power" "Power production" "W"
# State for derived time quantities.
state integer itimestep - - - - rh "itimestep" "" ""
state real xtime - - - - rh "xtime" "minutes since YYYY-MM-DD hh:mm:ss" "minutes since simulation start"
state real julian - - - - - "julian" "day of year, 0.0 at 0Z on 1 Jan." "days"
# input file descriptor for lbcs on parent domain
state integer lbc_fid - - - - - "lbc_fid" "" ""
# indicates if tiling has been computed
state logical tiled - - - - - "tiled" "" ""
# indicates if patches have been computed
state logical patched - - - - - "patched" "" ""
# indicates whether to read input from file or generate
#state logical input_from_file - - - - - "input_from_file" "" ""
# indicates whether to recompute mu
state logical press_adj - - - - - "press_adj" "T/F flag adjust mu" ""
# Mask for moving nest interpolations
state integer imask_nostag ij misc 1 -
state integer imask_xstag ij misc 1 X
state integer imask_ystag ij misc 1 Y
state integer imask_xystag ij misc 1 XY
# vortex center indices; need for restarts of moving nests
state real xi - misc - - r
state real xj - misc - - r
state real vc_i - misc - - r
state real vc_j - misc - - r
#-----------------------------------------------------------------------------------------------------------------------------------------------------------------
# Scalar (4D) arrays
# Moist Scalars
#
# The first line ensures that there will be identifiers named moist and
# moist_tend even if there are not any moist scalars (so the essentially
# dry code will will still link properly)
#
state real - ikjftb moist 1 - - -
state real qv ikjftb moist 1 - \
i0rh01usdf=(bdy_interp:dt) "QVAPOR" "Water vapor mixing ratio" "kg kg-1"
state real qc ikjftb moist 1 - \
i0rhusdf=(bdy_interp:dt) "QCLOUD" "Cloud water mixing ratio" "kg kg-1"
state real qr ikjftb moist 1 - \
i0rhusdf=(bdy_interp:dt) "QRAIN" "Rain water mixing ratio" "kg kg-1"
state real qi ikjftb moist 1 - \
i0rhusdf=(bdy_interp:dt) "QICE" "Ice mixing ratio" "kg kg-1"
state real qi2 ikjftb moist 1 - \
i0rhusdf=(bdy_interp:dt) "QICE2" "Ice mixing ratio cat 2" "kg kg-1"
state real qi3 ikjftb moist 1 - \
i0rhusdf=(bdy_interp:dt) "QICE3" "Ice mixing ratio cat 3" "kg kg-1"
state real qs ikjftb moist 1 - \
i0rhusdf=(bdy_interp:dt) "QSNOW" "Snow mixing ratio" "kg kg-1"
state real qg ikjftb moist 1 - \
i0rhusdf=(bdy_interp:dt) "QGRAUP" "Graupel mixing ratio" "kg kg-1"
state real qh ikjftb moist 1 - \
i0rhusdf=(bdy_interp:dt) "QHAIL" "Hail mixing ratio" "kg kg-1"
state real - ikjftb dfi_moist 1 - - -
state real dfi_qv ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QVAPOR" "Water vapor mixing ratio" "kg kg-1"
state real dfi_qc ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QCLOUD" "Cloud water mixing ratio" "kg kg-1"
state real dfi_qr ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QRAIN" "Rain water mixing ratio" "kg kg-1"
state real dfi_qi ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QICE" "Ice mixing ratio" "kg kg-1"
state real dfi_qi2 ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QICE2" "Ice mixing ratio cat 2" "kg kg-1"
state real dfi_qi3 ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QICE3" "Ice mixing ratio cat 3" "kg kg-1"
state real dfi_qs ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QSNOW" "Snow mixing ratio" "kg kg-1"
state real dfi_qg ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QGRAUP" "Graupel mixing ratio" "kg kg-1"
state real dfi_qh ikjftb dfi_moist 1 - \
rusdf=(bdy_interp:dt) "DFI_QHAIL" "Hail mixing ratio" "kg kg-1"
state real qvold ikj misc 1 - rdu "QVOLD" "Water vapor mixing ratio, old time step" "kg kg-1"
state real rimi ikj misc 1 - irh "RIMI" "riming intensity" "fraction"
state real qnwfa2d ij misc 1 - i014rhdu "QNWFA2D" "Surface aerosol number conc emission" "kg-1 s-1"
state real qnifa2d ij misc 1 - i014rhdu "QNIFA2D" "Surface dust number conc emission" "kg-1 s-1"
state real re_cloud ikj misc 1 - r "RE_CLOUD" "Effective radius cloud water" "m"
state real re_ice ikj misc 1 - r "RE_ICE" "Effective radius cloud ice" "m"
state real re_snow ikj misc 1 - r "RE_SNOW" "Effective radius snow" "m"
state real re_cloud_gsfc ikj misc 1 - rh "RE_CLOUD_GSFC" "Cloud Water effective radius" "micron"
state real re_rain_gsfc ikj misc 1 - rh "RE_RAIN_GSFC" "Rain Water effective radius" "micron"
state real re_ice_gsfc ikj misc 1 - rh "RE_ICE_GSFC" "Cloud Ice effective radius" "micron"
state real re_snow_gsfc ikj misc 1 - rh "RE_SNOW_GSFC" "Snow effective radius" "micron"
state real re_graupel_gsfc ikj misc 1 - rh "RE_GRAUPEL_GSFC" "Graupel Water effective radius" "micron"
state real re_hail_gsfc ikj misc 1 - rh "RE_HAIL_GSFC" "Hail Water effective radius" "micron"
state real dfi_re_cloud ikj misc 1 - - "DFI_RE_CLOUD" "DFI Effective radius cloud water" "m"
state real dfi_re_ice ikj misc 1 - - "DFI_RE_ICE" "DFI Effective radius cloud ice" "m"
state real dfi_re_snow ikj misc 1 - - "DFI_RE_SNOW" "DFI Effective radius snow" "m"
state real dfi_re_cloud_gsfc ikj misc 1 - - "DFI_RE_CLOUD_GSFC" "DFI Cloud Water effective radius" "micron"
state real dfi_re_rain_gsfc ikj misc 1 - - "DFI_RE_RAIN_GSFC" "DFI Rain Water effective radius" "micron"
state real dfi_re_ice_gsfc ikj misc 1 - - "DFI_RE_ICE_GSFC" "DFI Cloud Ice effective radius" "micron"
state real dfi_re_snow_gsfc ikj misc 1 - - "DFI_RE_SNOW_GSFC" "DFI Snow effective radius" "micron"
state real dfi_re_graupel_gsfc ikj misc 1 - - "DFI_RE_GRAUPEL_GSFC" "DFI Graupel Water effective radius" "micron"
state real dfi_re_hail_gsfc ikj misc 1 - - "DFI_RE_HAIL_GSFC" "DFI Hail Water effective radius" "micron"
state integer has_reqc - misc 1 - r "has_reqc" "Flag for having effective radius cloud water" ""
state integer has_reqi - misc 1 - r "has_reqi" "Flag for having effective radius cloud ice" ""
state integer has_reqs - misc 1 - r "has_reqs" "Flag for having effective radius snow" ""
# Other Scalars
state real - ikjftb scalar 1 - - -
state real qndrop ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNDROP" "Droplet number mixing ratio" "# kg-1"
state real qni ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNICE" "Ice Number concentration" "# kg-1"
state real qni2 ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNICE2" "Ice Number concentration cat 2" "# kg-1"
state real qni3 ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNICE3" "Ice Number concentration cat 3" "# kg-1"
state real qt ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "CWM" "Total condensate mixing ratio" "kg kg-1"
state real qns ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNSNOW" "Snow Number concentration" "# kg(-1)"
state real qnr ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNRAIN" "Rain Number concentration" "# kg(-1)"
state real qng ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNGRAUPEL" "Graupel Number concentration" "# kg(-1)"
state real qnh ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNHAIL" "Hail Number concentration" "# kg(-1)"
state real qnn ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNCCN" "CCN Number concentration" "# kg(-1)"
state real qnc ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QNCLOUD" "cloud water Number concentration" "# kg(-1)"
state real qvolg ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QVGRAUPEL" "Graupel Particle Volume" "m(3) kg(-1)"
state real qvolh ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QVHAIL" "Hail Particle Volume" "m(3) kg(-1)"
state real qrimef ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QRIMEF" "rime factor * qi" "kg kg-1"
state real qir ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QIR" "Rime ice mass-1 mixing ratio" "kg kg(-1)"
state real qib ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QIB" "Rime ice volume-1 mixing ratio" "m(3) kg(-1)"
state real qir2 ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QIR2" "Rime ice mass-2 mixing ratio" "kg kg(-1)"
state real qib2 ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QIB2" "Rime ice volume-2 mixing ratio" "m(3) kg(-1)"
state real qvoli ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QVOLI" "Ice volume-1 mixing ratio" "m(3) kg(-1)"
state real qaoli ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QAOLI" "Ice volume-1 times aspect ratio mixing ratio" "m(3) kg(-1)"
state real qvoli2 ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QVOLI2" "Ice volume-2 mixing ratio" "m(3) kg(-1)"
state real qaoli2 ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QAOLI2" "Ice volume-2 times aspect ratio mixing ratio" "m(3) kg(-1)"
state real qvoli3 ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QVOLI3" "Ice volume-3 mixing ratio" "m(3) kg(-1)"
state real qaoli3 ikjftb scalar 1 - \
i0rhusdf=(bdy_interp:dt) "QAOLI3" "Ice volume-3 times aspect ratio mixing ratio" "m(3) kg(-1)"
state real - ikjftb dfi_scalar 1 - - -
state real dfi_qndrop ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNDROP" "DFI Droplet number mixing ratio" "# kg-1"
state real dfi_qni ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNICE" "DFI Ice Number concentration" "# kg-1"
state real dfi_qni2 ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNICE2" "DFI Ice Number concentration cat 2" "# kg-1"
state real dfi_qni3 ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNICE3" "DFI Ice Number concentration cat 3" "# kg-1"
state real dfi_qt ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_CWM" "DFI Total condensate mixing ratio" "kg kg-1"
state real dfi_qns ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNSNOW" "DFI Snow Number concentration" "# kg(-1)"
state real dfi_qnr ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNRAIN" "DFI Rain Number concentration" "# kg(-1)"
state real dfi_qng ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNGRAUPEL" "DFI Graupel Number concentration" "# kg(-1)"
state real dfi_qnh ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNHAIL" "DFI Hail Number concentration" "# kg(-1)"
state real dfi_qnn ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNCC" "DFI CNN Number concentration" "# kg(-1)"
state real dfi_qnc ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "DFI_QNCLOUD" "DFI Cloud Number concentration" "# kg(-1)"
state real dfi_qvolg ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QVGRAUPEL" "DFI Graupel Particle Volume" "m(3) kg(-1)"
state real dfi_qvolh ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QVHAIL" "DFI Hail Particle Volume" "m(3) kg(-1)"
state real dfi_qir ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QIR" "DFI Rime ice mass-1 mixing ratio" "kg kg(-1)"
state real dfi_qib ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QIB" "DFI Rime ice volume-1 mixing ratio" "m(3) kg(-1)"
state real dfi_qir2 ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QIR2" "DFI Rime ice mass-2 mixing ratio" "kg kg(-1)"
state real dfi_qib2 ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QIB2" "DFI Rime ice volume-2 mixing ratio" "m(3) kg(-1)"
state real dfi_qvoli ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QVOLI" "DFI Ice volume-1 mixing ratio" "kg kg(-1)"
state real dfi_qaoli ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QAOLI" "DFI Ice volume-1 times aspect ratio mixing ratio" "kg kg(-1)"
state real dfi_qvoli2 ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QVOLI2" "DFI Ice volume-2 mixing ratio" "kg kg(-1)"
state real dfi_qaoli2 ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QAOLI2" "DFI Ice volume-2 times aspect ratio mixing ratio" "kg kg(-1)"
state real dfi_qvoli3 ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QVOLI3" "DFI Ice volume-3 mixing ratio" "kg kg(-1)"
state real dfi_qaoli3 ikjftb dfi_scalar 1 - \
rhusdf=(bdy_interp:dt) "DFI_QAOLI3" "DFI Ice volume-3 times aspect ratio mixing ratio" "kg kg(-1)"
state real dfi_qke_adv ikjftb dfi_scalar 1 - \
rusdf=(bdy_interp:dt) "dfi_qke_adv" "DFI twice TKE from MYNN" "m2 s-2"
#-----------------------------------------------------------------------------------------------------------------------------------------------------------------
# Arrays for Specified LBCs (lbc arrays REMOVED; Boundary arrays are now specified with the state array; see above, 20050413 JM )
state real fcx w misc - - ir "fcx" "RELAXATION TERM FOR BOUNDARY ZONE" ""
state real gcx w misc - - ir "gcx" "2ND RELAXATION TERM FOR BOUNDARY ZONE" ""
state real dtbc - misc - - r "dtbc" "TIME SINCE BOUNDARY READ" ""
#-------------------------------------------------------------------------------------------------------------------------------------------
#-------------------------------------------------------------------------------------------------------------------------------------------
# Physics Related State Varibles
#-------------------------------------------------------------------------------------------------------------------------------------------
# SI - start variables from netCDF format from Standard Initialization, most eventually for use in LSM schemes
#-------------------------------------------------------------------------------------------------------------------------------------------
state real soil_layers i{lin}j misc 1 Z i1 "SOIL_LAYERS" "SOIL LAYERS" "cm"
state real soil_levels i{lin}j misc 1 Z i1 "SOIL_LEVELS" "SOIL LEVELS" "cm"
state real st i{lin}j misc 1 Z i1 "ST" "SOIL TEMPERATURES" "K"
state real sm i{lin}j misc 1 Z i1 "SM" "SOIL MOISTURES" "m3 m-3"
state real sw i{lin}j misc 1 Z i1 "SW" "SOIL LIQUIDS" "m3 m-3"
state real soilt i{lin}j misc 1 Z i1 "SOILT" "RUC SOIL TEMPERATURES" "K"
state real soilm i{lin}j misc 1 Z i1 "SOILM" "RUC SOIL MOISTURES" "m3 m-3"
state real sm000007 ij misc 1 - i1 "SM000007" "LAYER SOIL MOISTURE" "m3 m-3"
state real sm007028 ij misc 1 - i1 "SM007028" "LAYER SOIL MOISTURE" "m3 m-3"
state real sm028100 ij misc 1 - i1 "SM028100" "LAYER SOIL MOISTURE" "m3 m-3"
state real sm100255 ij misc 1 - i1 "SM100255" "LAYER SOIL MOISTURE" "m3 m-3"
state real st000007 ij misc 1 - i1 "ST000007" "LAYER SOIL TEMPERATURE" "K"
state real st007028 ij misc 1 - i1 "ST007028" "LAYER SOIL TEMPERATURE" "K"
state real st028100 ij misc 1 - i1 "ST028100" "LAYER SOIL TEMPERATURE" "K"
state real st100255 ij misc 1 - i1 "ST100255" "LAYER SOIL TEMPERATURE" "K"
state real sm000010 ij misc 1 - i1 "SM000010" "LAYER SOIL MOISTURE" "m3 m-3"
state real sm010040 ij misc 1 - i1 "SM010040 " "LAYER SOIL MOISTURE" "m3 m-3"
state real sm040100 ij misc 1 - i1 "SM040100 " "LAYER SOIL MOISTURE" "m3 m-3"
state real sm100200 ij misc 1 - i1 "SM100200 " "LAYER SOIL MOISTURE" "m3 m-3"
state real sm010200 ij misc 1 - i1 "SM010200" "LAYER SOIL MOISTURE" "m3 m-3"
state real soilm000 ij misc 1 - i1 "SOILM000" "LAYER SOIL MOISTURE" "m3 m-3"
state real soilm005 ij misc 1 - i1 "SOILM005" "LAYER SOIL MOISTURE" "m3 m-3"
state real soilm020 ij misc 1 - i1 "SOILM020" "LAYER SOIL MOISTURE" "m3 m-3"
state real soilm040 ij misc 1 - i1 "SOILM040" "LAYER SOIL MOISTURE" "m3 m-3"
state real soilm160 ij misc 1 - i1 "SOILM160" "LAYER SOIL MOISTURE" "m3 m-3"
state real soilm300 ij misc 1 - i1 "SOILM300" "LAYER SOIL MOISTURE" "m3 m-3"
state real sw000010 ij misc 1 - i1 "SW000010" "LAYER SOIL LIQUID" "m3 m-3"
state real sw010040 ij misc 1 - i1 "SW010040" "LAYER SOIL LIQUID" "m3 m-3"
state real sw040100 ij misc 1 - i1 "SW040100" "LAYER SOIL LIQUID" "m3 m-3"
state real sw100200 ij misc 1 - i1 "SW100200" "LAYER SOIL LIQUID" "m3 m-3"
state real sw010200 ij misc 1 - i1 "SW010200" "LAYER SOIL LIQUID" "m3 m-3"
state real soilw000 ij misc 1 - i1 "SOILW000" "LAYER SOIL LIQUID" "m3 m-3"
state real soilw005 ij misc 1 - i1 "SOILW005" "LAYER SOIL LIQUID" "m3 m-3"
state real soilw020 ij misc 1 - i1 "SOILW020" "LAYER SOIL LIQUID" "m3 m-3"
state real soilw040 ij misc 1 - i1 "SOILW040" "LAYER SOIL LIQUID" "m3 m-3"
state real soilw160 ij misc 1 - i1 "SOILW160" "LAYER SOIL LIQUID" "m3 m-3"
state real soilw300 ij misc 1 - i1 "SOILW300" "LAYER SOIL LIQUID" "m3 m-3"
state real st000010 ij misc 1 - i1 "ST000010" "LAYER SOIL TEMPERATURE" "K"
state real st010040 ij misc 1 - i1 "ST010040" "LAYER SOIL TEMPERATURE" "K"
state real st040100 ij misc 1 - i1 "ST040100" "LAYER SOIL TEMPERATURE" "K"
state real st100200 ij misc 1 - i1 "ST100200" "LAYER SOIL TEMPERATURE" "K"
state real st010200 ij misc 1 - i1 "ST010200" "LAYER SOIL TEMPERATURE" "K"
state real soilt000 ij misc 1 - i1 "SOILT000" "LAYER SOIL TEMPERATURE" "K"
state real soilt005 ij misc 1 - i1 "SOILT005" "LAYER SOIL TEMPERATURE" "K"
state real soilt020 ij misc 1 - i1 "SOILT020" "LAYER SOIL TEMPERATURE" "K"
state real soilt040 ij misc 1 - i1 "SOILT040" "LAYER SOIL TEMPERATURE" "K"
state real soilt160 ij misc 1 - i1 "SOILT160" "LAYER SOIL TEMPERATURE" "K"
state real soilt300 ij misc 1 - i1 "SOILT300" "LAYER SOIL TEMPERATURE" "K"
state real topostdv ij misc 1 - i12 "TOPOSTDV" "ELEVATION STD DEV" "m"
state real toposlpx ij misc 1 - i012rdu "TOPOSLPX" "ELEVATION X SLOPE" ""
state real toposlpy ij misc 1 - i012rdu "TOPOSLPY" "ELEVATION Y SLOPE" ""
state real slope ij misc 1 - rdu "SLOPE" "ELEVATION SLOPE" ""
state real slp_azi ij misc 1 - rdu "SLP_AZI" "ELEVATION SLOPE AZIMUTH" "rad"
state real shdmax ij misc 1 - i012rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "SHDMAX" "ANNUAL MAX VEG FRACTION" ""
state real shdmin ij misc 1 - i012rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "SHDMIN" "ANNUAL MIN VEG FRACTION" ""
state real snoalb ij misc 1 - i012rhd "SNOALB" "ANNUAL MAX SNOW ALBEDO IN FRACTION" ""
state real toposoil ij misc 1 - i12 "SOILHGT" "ELEVATION OF LSM DATA" "m"
state real landusef iuj misc 1 Z i012rdu "LANDUSEF" "LANDUSE FRACTION BY CATEGORY" ""
state real soilctop isj misc 1 Z i012rdu "SOILCTOP" "SOIL CAT FRACTION (TOP)" ""
state real soilcbot isj misc 1 Z i012rdu "SOILCBOT" "SOIL CAT FRACTION (BOTTOM)" ""
state real soilcat ij misc 1 - i12 "SOILCAT" "SOIL CAT DOMINANT TYPE" ""
state real vegcat ij misc 1 - i12 "VEGCAT" "VEGETATION CAT DOMINANT TYPE" ""
#---------------------------------------------------------------------------------------------------------------------------------------
# SI - end variables from netCDF format from Standard Initialization
#---------------------------------------------------------------------------------------------------------------------------------------
# soil model variables (Note that they are marked as staggered in the vertical dimension
# because they are "fully dimensioned" -- they use every element in that dim
state real TSLB ilj misc 1 Z i02rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TSLB" "SOIL TEMPERATURE" "K"
# Time series variables
state real ts_hour ?! misc - - - "TS_HOUR" "Model integration time, hours"
state real ts_u ?! misc - - - "TS_U" "Surface wind U-component, earth-relative"
state real ts_v ?! misc - - - "TS_V" "Surface wind V-component, earth-relative"
state real ts_q ?! misc - - - "TS_Q" "Surface mixing ratio"
state real ts_t ?! misc - - - "TS_T" "Surface temperature"
state real ts_psfc ?! misc - - - "TS_PSFC" "Surface pressure"
state real ts_glw ?! misc - - - "TS_GLW" "Downward long wave flux at surface"
state real ts_gsw ?! misc - - - "TS_GSW" "Net short wave flux at surface"
state real ts_hfx ?! misc - - - "TS_HFX" "Upward heat flux at surface"
state real ts_lh ?! misc - - - "TS_LH" "Upward moisture flux at surface"
state real ts_tsk ?! misc - - - "TS_TSK" "Skin temperature"
state real ts_tslb ?! misc - - - "TS_TSLB" "Soil temperature"
state real ts_clw ?! misc - - - "TS_CLW" "Column integrated cloud water"
state real ts_rainc ?! misc - - - "TS_RAINC" "Cumulus precip"
state real ts_rainnc ?! misc - - - "TS_RAINNC" "Grid-scale precip"
# Time series of vertical profile of U, V, GHT, THETA, QVAPOR, PRESSURE
state real ts_u_profile ?!k misc - - - "TS_U_PROFILE" "Wind u-Component, earth relative, vertical profile" "m/s"
state real ts_v_profile ?!k misc - - - "TS_V_PROFILE" "Wind v-Component, earth relative, vertical profile" "m/s"
state real ts_w_profile ?!k misc - - - "TS_W_PROFILE" "Wind w-Component, vertical profile" "m/s"
state real ts_gph_profile ?!k misc - - - "TS_GPH_PROFILE" "Total geopotential height, vertical profile" "m"
state real ts_th_profile ?!k misc - - - "TS_TH_PROFILE" "Total potential temperature, vertical profile" "K"
state real ts_qv_profile ?!k misc - - - "TS_QV_PROFILE" "Water vapor mixing ratio, vertical profile" "kg/kg"
state real ts_p_profile ?!k misc - - - "TS_P_PROFILE" "Pressure, vertical profile" "Pa"
# urban model variables
state real DZR l em - Z r "DZR" "THICKNESSES OF ROOF LAYERS" "m"
state real DZB l em - Z r "DZB" "THICKNESSES OF WALL LAYERS" "m"
state real DZG l em - Z r "DZG" "THICKNESSES OF ROAD LAYERS" "m"
state real URB_PARAM i{urb}j misc 1 - i1 "URB_PARAM" "NUDAPT_NBSD Urban Parameters" "parameter"
state real LP_URB2D ij misc 1 - i0 "BUILD_AREA_FRACTION" "BUILDING PLAN AREA DENSITY" "dimensionless"
state real HI_URB2D i{uhi}j misc 1 Z i0 "HEIGHT_HISTOGRAMS" "DISTRIBUTION OF BUILDING HEIGHTS" "dimensionless"
state real LB_URB2D ij misc 1 - i0 "BUILD_SURF_RATIO" "BUILDING SURFACE AREA TO PLAN AREA RATIO" "dimensionless"
state real HGT_URB2D ij misc 1 - i0 "BUILD_HEIGHT" "AVERAGE BUILDING HEIGHT WEIGHTED BY BUILDING PLAN AREA" "m"
state real MH_URB2D ij misc 1 - i0 "MH_URB2D" "Mean Building Height" "m"
state real STDH_URB2D ij misc 1 - i0 "STDH_URB2D" "Standard Deviation of Building Height" "m2"
state real LF_URB2D i{udr}j misc 1 Z i0 "LF_URB2D" "Frontal Area Index" "dimensionless"
# lsm State Variables
state real SMOIS ilj - 1 Z i02rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "SMOIS" "SOIL MOISTURE" "m3 m-3"
state real SH2O ilj - 1 Z i02rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "SH2O" "SOIL LIQUID WATER" "m3 m-3"
state real SMCREL ilj - 1 Z i02rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "SMCREL" "RELATIVE SOIL MOISTURE" ""
state real XICE ij misc 1 - i0124rhd=(interp_mask_field:lu_index,isice)u=(copy_fcnm) "SEAICE" "SEA ICE FLAG" ""
state real ICEDEPTH ij misc 1 - i0124rhd=(interp_mask_field:lu_index,isice)u=(copy_fcnm) "ICEDEPTH" "SEA ICE THICKNESS" "m"
state real XICEM ij misc 1 - rhd=(interp_mask_field:lu_index,isice)u=(copy_fcnm) "XICEM" "SEA ICE FLAG (PREVIOUS STEP)" ""
state real ALBSI ij misc 1 - i0124rhd=(interp_mask_field:lu_index,isice)u=(copy_fcnm) "ALBSI" "SEA ICE ALBEDO" ""
state real SNOWSI ij misc 1 - i0124rhd=(interp_mask_field:lu_index,isice)u=(copy_fcnm) "SNOWSI" "SNOW DEPTH ON SEA ICE" "m"
state real SMSTAV ij misc 1 - rd=(interp_mask_field:lu_index,iswater) "SMSTAV" "MOISTURE AVAILABILITY" ""
state real SMSTOT ij misc 1 - r "SMSTOT" "TOTAL SOIL MOISTURE" "m3 m-3"
state real SOLDRAIN ij misc 1 - r "SOLDRAIN" "soil column drainage" "mm"
state real SFCHEADRT ij misc 1 - r "SFCHEADRT" "surface water depth" "mm"
state real INFXSRT ij misc 1 - r "INFXSRT" "time step infiltration excess" "mm"
state real SFCRUNOFF ij misc 1 - rhd=(interp_mask_field:lu_index,iswater) "SFROFF" "SURFACE RUNOFF" "mm"
state real UDRUNOFF ij misc 1 - rhd=(interp_mask_field:lu_index,iswater) "UDROFF" "UNDERGROUND RUNOFF" "mm"
state integer IVGTYP ij misc 1 - i02rhd=(interp_fcni)u=(copy_fcni) "IVGTYP" "DOMINANT VEGETATION CATEGORY" ""
state integer ISLTYP ij misc 1 - i02rhd=(interp_mask_soil:lu_index)u=(copy_fcni) "ISLTYP" "DOMINANT SOIL CATEGORY" ""
state real VEGFRA ij misc 1 - i024rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "VEGFRA" "VEGETATION FRACTION" ""
state real SFCEVP ij misc 1 - r "SFCEVP" "ACCUMULATED SURFACE EVAPORATION" "kg m-2"
state real GRDFLX ij misc 1 - rh "GRDFLX" "GROUND HEAT FLUX" "W m-2"
state real ACGRDFLX ij misc 1 - rhdu "ACGRDFLX" "ACCUMULATED GROUND HEAT FLUX" "J m-2"
state real SFCEXC ij misc 1 - r "SFCEXC " "SURFACE EXCHANGE COEFFICIENT" "m s-1"
state real ACSNOW ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "ACSNOW" "ACCUMULATED SNOW" "kg m-2"
state real ACRUNOFF ij misc 1 - irhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "ACRUNOFF" "ACCUMULATED RUNOFF" "kg m-2"
state real ACSNOM ij misc 1 - rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "ACSNOM" "ACCUMULATED MELTED SNOW" "kg m-2"
state real SNOW ij misc 1 - i012rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "SNOW" "SNOW WATER EQUIVALENT" "kg m-2"
state real SNOWH ij misc 1 - i012rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "SNOWH" "PHYSICAL SNOW DEPTH" "m"
#state real RHOSN ij misc 1 - i02rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "RHOSN" " SNOW DENSITY" "kg m-3"
state real CANWAT ij misc 1 - i012rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "CANWAT" "CANOPY WATER" "kg m-2"
state integer IFNDSNOWH - misc 1 - i "FNDSNOWH" "SNOWH_LOGICAL"
state integer IFNDSOILW - misc 1 - i "FNDSOILW" "SOILW_LOGICAL"
state integer IFNDALBSI - misc 1 - ir "FNDALBSI" "ALBSI_LOGICAL"
state integer IFNDSNOWSI - misc 1 - ir "FNDSNOWSI" "SNOWSI_LOGICAL"
state integer IFNDICEDEPTH - misc 1 - ir "FNDICEDEPTH" "ICEDEPTH_LOGICAL"
state real XLAIDYN ij misc 1 - - "XLAIDYN" "Noah Dynamic LEAF AREA INDEX" "m-2/m-2"
# SKIN SST
state real SSTSK ij misc 1 - rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "SSTSK" "SKIN SEA SURFACE TEMPERATURE" "K"
state real lake_depth ij misc 1 - i012rd=(interp_mask_water_field:lu_index,islake) "lake_depth" "lake depth" "m"
state real DTW ij misc 1 - r "DTW" "WARM LAYER TEMP DIFF" "C"
# Ocean surface currents
state real UOCE ij misc 1 - i0124rd=(interp_mask_water_field:lu_index,iswater) "UOCE" "SEA SURFACE ZONAL CURRENTS" "m s-1"
state real VOCE ij misc 1 - i0124rd=(interp_mask_water_field:lu_index,iswater) "VOCE" "SEA SURFACE MERIDIONAL CURRENTS" "m s-1"
# DFI variables
state real hcoeff {ndfi} misc 1 - - "HCOEFF" "initialization weights"
state real hcoeff_tot - misc 1 - - "HCOEFF_TOT" "initialization weights"
state real dfi_p ikj misc 1 - r "P_DFI" "perturbation pressure" "Pa"
state real dfi_al ikj misc 1 - r "AL_DFI" "inverse perturbation density" "m3 kg-1"
state real dfi_mu ij misc 1 - r "MU_DFI" "perturbation dry air mass in column" "Pa"
state real dfi_phb ikj misc 1 Z r "PHB_DFI" "base-state geopotential" "m2 s-2"
state real dfi_ph0 ikj misc 1 Z r "PH0_DFI" "initial geopotential" "m2 s-2"
state real dfi_php ikj misc 1 Z r "PHP_DFI" "geopotential" "m2 s-2"
state real dfi_u ikj misc 1 - r "U_DFI" "u accumulation array" " "
state real dfi_v ikj misc 1 - r "V_DFI" "v accumulation array" " "
state real dfi_w ikj misc 1 - r "W_DFI" "w accumulation array" " "
state real dfi_ww ikj misc 1 Z r "WW_DFI" "mu-coupled eta-dot" "Pa s-1"
state real dfi_t ikj misc 1 - r "TT_DFI" "t accumulation array" " "
state real dfi_rh ikj misc 1 - r "RH_DFI" "initial relative humidity" " "
state real dfi_ph ikj misc 1 - r "PH_DFI" "p accumulation array" " "
state real dfi_pb ikj misc 1 - r "PB_DFI" "pb accumulation array" " "
state real dfi_alt ikj misc 1 - r "ALT_DFI" "1/rho accumulation array" " "
state real dfi_tke ikj misc 1 - r "TKE_DFI" "TURBULENCE KINETIC ENERGY" "m2 s-2"
state real dfi_tten_rad ikj misc 1 - ir "RAD_TTEN_DFI" "RADAR POT. TEMP. TENDENCY" "K s-1"
state real dfi_TSLB ilj misc 1 Z r "TSLB_dfi" "SOIL TEMPERATURE" "K"
state real dfi_SMOIS ilj - 1 Z r "SMOIS_dfi" "SOIL MOISTURE" "m3 m-3"
state real dfi_SNOW ij misc 1 - r "SNOW_dfi" "SNOW WATER EQUIVALENT" "kg m-2"
state real dfi_SNOWH ij misc 1 - r "SNOWH_dfi" "PHYSICAL SNOW DEPTH" "m"
state real dfi_CANWAT ij misc 1 - r "CANWAT_dfi" "CANOPY WATER" "kg m-2"
state real dfi_SMFR3D ilj misc 1 Z r "SMFR3D_dfi" "SOIL ICE" ""
state real dfi_KEEPFR3DFLAG ilj misc 1 Z r "KEEPFR3DFLAG_dfi" "FLAG - 1. FROZEN SOIL YES, 0 - NO" ""
# urban state variables
state real TSK_RURAL ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TSK_RURAL" "TSK for rural fraction" "K"
state real TR_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TR_URB" "URBAN ROOF SKIN TEMPERATURE" "K"
state real TGR_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TGR_URB" "URBAN GREEN ROOF SKIN TEMPERATURE" "K"
state real TB_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TB_URB" "URBAN WALL SKIN TEMPERATURE" "K"
state real TG_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TG_URB" "URBAN ROAD SKIN TEMPERATURE" "K"
state real TC_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TC_URB" "URBAN CANOPY TEMPERATURE" "K"
state real QC_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "QC_URB" "URBAN CANOPY HUMIDITY" "kg kg{-1}"
state real UC_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "UC_URB" "URBAN CANOPY WIND" "m s{-1}"
state real XXXR_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "XXXR_URB" "M-O LENGTH ABOVE URBAN ROOF" "dimensionless"
state real XXXB_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "XXXB_URB" "M-O LENGTH ABOVE URBAN WALL" "dimensionless"
state real XXXG_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "XXXG_URB" "M-O LENGTH ABOVE URBAN ROAD" "dimensionless"
state real XXXC_URB2D ij misc 1 - rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "XXXC_URB" "M-O LENGTH ABOVE URBAN CANOPY" "dimensionless"
state real CMCR_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "CMCR_URB" "GREEN ROOF CANOPY INTERCAPTED WATER" "m"
state real DRELR_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "DRELR_URB" "WATER HOLDING DEPTH ON ROOF IMPERVIOUS SURFACE" "m"
state real DRELB_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "DRELB_URB" "WATER HOLDING DEPTH ON WALL IMPERVIOUS SURFACE" "m"
state real DRELG_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "DRELG_URB" "WATER HOLDING DEPTH ON ROAD IMPERVIOUS SURFACE" "m"
state real FLXHUMR_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "FLXHUMR_URB" "WATER FLUX ON ROOF IMPERVIOUS SURFACE" "m/s"
state real FLXHUMB_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "FLXHUMB_URB" "WATER FLUX ON WALL IMPERVIOUS SURFACE" "m/s"
state real FLXHUMG_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "FLXHUMG_URB" "WATER FLUX ON ROAD IMPERVIOUS SURFACE" "m/s"
state real TGRL_URB3D ilj misc 1 Z rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TGRL_URB" "GREEN ROOF LAYER TEMPERATURE"
state real SMR_URB3D ilj misc 1 Z rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "SMR_URB" "GREEN ROOF LAYER SOIL MOISTURE"
state real TRL_URB3D ilj misc 1 Z rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TRL_URB" "ROOF LAYER TEMPERATURE" "K"
state real TBL_URB3D ilj misc 1 Z rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TBL_URB" "WALL LAYER TEMPERATURE" "K"
state real TGL_URB3D ilj misc 1 Z rd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "TGL_URB" "ROAD LAYER TEMPERATURE" "K"
state real SH_URB2D ij misc 1 - rd=(interp_fcnm)u=(copy_fcnm) "SH_URB" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real LH_URB2D ij misc 1 - rd=(interp_fcnm)u=(copy_fcnm) "LH_URB" "LATENT HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real G_URB2D ij misc 1 - rd=(interp_fcnm)u=(copy_fcnm) "G_URB" "GROUND HEAT FLUX INTO URBAN" "W m{-2}"
state real RN_URB2D ij misc 1 - rd=(interp_fcnm)u=(copy_fcnm) "RN_URB" "NET RADIATION ON URBAN SFC" "W m{-2}"
state real TS_URB2D ij misc 1 - rd=(interp_fcnm)u=(copy_fcnm) "TS_URB" "SKIN TEMPERATURE" "K"
state real FRC_URB2D ij misc 1 - i012rd=(interp_fcnm)u=(copy_fcnm) "FRC_URB2D" "URBAN FRACTION" "dimensionless"
state integer UTYPE_URB2D ij misc 1 - rd=(interp_fcnm)u=(copy_fcnm) "UTYPE_URB" "URBAN TYPE" "dimensionless"
state real TRB_URB4D i{umap1}j misc 1 Z r "TRB_URB4D" "ROOF LAYER TEMPERATURE" "K"
state real TW1_URB4D i{umap2}j misc 1 Z r "TW1_URB4D" "WALL LAYER TEMPERATURE" "K"
state real TW2_URB4D i{umap2}j misc 1 Z r "TW2_URB4D" "WALL LAYER TEMPERATURE" "K"
state real TGB_URB4D i{umap3}j misc 1 Z r "TGB_URB4D" "ROAD LAYER TEMPERATURE" "K"
state real TLEV_URB3D i{umap6}j misc 1 Z r "TLEV_URB3D" "INDOOR TEMPERATURE" "K"
state real QLEV_URB3D i{umap6}j misc 1 Z r "QLEV_URB3D" "SPECIFIC HUMIDITY" "dimensionless"
state real TW1LEV_URB3D i{umap7}j misc 1 Z r "TW1LEV_URB3D" "WINDOW TEMPERATURE" "K"
state real TW2LEV_URB3D i{umap7}j misc 1 Z r "TW2LEV_URB3D" "WINDOW TEMPERATURE" "K"
state real TGLEV_URB3D i{umap8}j misc 1 Z r "TGLEV_URB3D" "GROUND TEMPERATURE BELOW A BUILDING" "K"
state real TFLEV_URB3D i{umap9}j misc 1 Z r "TFLEV_URB3D" "FLOOR TEMPERATURE" "K"
state real SF_AC_URB3D ij misc 1 - r "SF_AC_URB3D" "SENSIBLE HEAT FLUX FROM THE AIR COND." "W m{-2}"
state real LF_AC_URB3D ij misc 1 - r "LF_AC_URB3D" "LATENT HEAT FLUX FROM THE AIR COND." "W m{-2}"
state real CM_AC_URB3D ij misc 1 - r "CM_AC_URB3D" "CONSUMPTION OF THE AIR COND." "W m{-2}"
state real SFVENT_URB3D ij misc 1 - r "SFVENT_URB3D" "SENSIBLE HEAT FLUX FROM URBAN VENTILATION" "W m{-2}"
state real LFVENT_URB3D ij misc 1 - r "LFVENT_URB3D" "LATENT HEAT FLUX FROM URBAN VENTILATION" "W m{-2}"
state real SFWIN1_URB3D i{umap7}j misc 1 Z r "SFWIN1_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC WINDOW" "W m{-2}"
state real SFWIN2_URB3D i{umap7}j misc 1 Z r "SFWIN2_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC WINDOW" "W m{-2}"
state real SFW1_URB3D i{umap4}j misc 1 Z r "SFW1_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real SFW2_URB3D i{umap4}j misc 1 Z r "SFW2_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real SFR_URB3D i{umap5}j misc 1 Z r "SFR_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real SFG_URB3D i{umap0}j misc 1 Z r "SFG_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real CMR_SFCDIF ij misc 1 - r "CMR_SFCDIF" "" ""
state real CHR_SFCDIF ij misc 1 - r "CHR_SFCDIF" "" ""
state real CMC_SFCDIF ij misc 1 - r "CMC_SFCDIF" "" ""
state real CHC_SFCDIF ij misc 1 - r "CHC_SFCDIF" "" ""
state real CMGR_SFCDIF ij misc 1 - r "CMGR_SFCDIF" "" ""
state real CHGR_SFCDIF ij misc 1 - r "CHGR_SFCDIF" "" ""
# solar location variables from radiation driver
state real COSZEN ij misc 1 - rh "COSZEN" "COS of SOLAR ZENITH ANGLE" "dimensionless"
state real HRANG ij misc 1 - r "HRANG" "SOLAR HOUR ANGLE" "radians"
state real DECLIN - misc 1 - r "DECLIN" "SOLAR DECLINATION" "radians"
state real SOLCON - misc 1 - r "SOLCON" "SOLAR CONSTANT" "W m-2"
# RUC LSM
state real RHOSNF ij misc 1 - irh "RHOSNF" "DENSITY OF FROZEN PRECIP" "kg/m^3"
state real SNOWFALLAC ij misc 1 - irh "SNOWFALLAC" "RUN-TOTAL ACCUMULATED SNOWFALL [mm]" ""
state real PRECIPFR ij misc 1 - - "PRECIPFR" "TIME-STEP FROZEN PRECIP [mm]" ""
state real SMFR3D ilj misc 1 Z r "SMFR3D" "SOIL ICE" ""
state real KEEPFR3DFLAG ilj misc 1 Z r "KEEPFR3DFLAG" "FLAG - 1. FROZEN SOIL YES, 0 - NO" ""
state real SWVISDIR ij misc 1 Z r "SWVISDIR" "SWR VIS DIR component" ""
state real SWVISDIF ij misc 1 Z r "SWVISDIF" "SWR VIS DIF component" ""
state real SWNIRDIR ij misc 1 Z r "SWNIRDIR" "SWR NIR DIR component" ""
state real SWNIRDIF ij misc 1 Z r "SWNIRDIF" "SWR NIR DIF component" ""
state real ALSWVISDIR ij misc 1 Z r "ALSWVISDIR" "ALB VIS DIR component" ""
state real ALSWVISDIF ij misc 1 Z r "ALSWVISDIF" "ALB VIS DIF component" ""
state real ALSWNIRDIR ij misc 1 Z r "ALSWNIRDIR" "ALB NIR DIR component" ""
state real ALSWNIRDIF ij misc 1 Z r "ALSWNIRDIF" "ALB NIR DIF component" ""
# From P-X PBL and LSM; RS (RC) and RA also from Noah and NoahMP
state real RA ij misc 1 - r "RA" "AERODYNAMIC RESISTANCE" "s m-1"
state real RS ij misc 1 - r "RS" "SURFACE RESISTANCE" "s m-1"
state real LAI ij misc 1 - i0124rh "LAI" "LEAF AREA INDEX" "m-2/m-2"
state real VEGF_PX ij misc 1 - rh "VEGF_PX" "Vegetation Fraction for PX LSM" "area/area"
state real T2OBS ij misc 1 - r "T2OBS" "2-m temperature from analysis " "K"
state real Q2OBS ij misc 1 - r "Q2OBS" "2-m mixing ratio from analysis " "kg/kg"
state real IMPERV ij misc 1 - i012rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "IMPERV" "Impervious surface fraction NLCD" "percent"
state real CANFRA ij misc 1 - i012rhd=(interp_mask_field:lu_index,iswater)u=(copy_fcnm) "CANFRA" "Satellite canopy fraction" "percent"
state real LAI_PX ij misc 1 - rh "LAI_PX" "LAI used for PX LSM" "m2/m2"
state real WWLT_PX ij misc 1 - rh "WWLT_PX" "Soil wilting point for PX LSM" "m3/m3"
state real WFC_PX ij misc 1 - rh "WFC_PX" "Soil field capacity for PX LSM" "m3/m3"
state real WSAT_PX ij misc 1 - rh "WSAT_PX" "Soil saturation for PX LSM" "m3/m3"
state real CLAY_PX ij misc 1 - rh "CLAY_PX" "Clay perent for PX LSM" "fraction"
state real CSAND_PX ij misc 1 - rh "CSAND_PX" "Coarse sand perent for PX LSM" "fraction"
state real FMSAND_PX ij misc 1 - rh "FMSAND_PX" "Fine-medium sand perent for PX LSM" "fraction"
# sfclay PBL variables
i1 real PSIM ij misc 1 - - "PSIM" "SIMILARITY FUNCTION FOR MOMENTUM" ""
i1 real PSIH ij misc 1 - - "PSIH" "SIMILARITY FUNCTION FOR HEAT" ""
state real FM ij misc 1 - - "FM" "INTEGRATED FUNCTION FOR MOMENTUM" ""
state real FH ij misc 1 - - "FH" "INTEGRATED FUNCTION FOR HEAT" ""
state real WSPD ij misc 1 - - "WSPD" "Wind Speed At Lowest Model Level (may contain vconv)" "m s-1"
i1 real GZ1OZ0 ij misc 1 - - "GZ1OZ0" "LOG OF Z1 over Z0" ""
state real BR ij misc 1 - - "BR" "Bulk Richardson" ""
state real ZOL ij misc 1 - - "ZOL" "z/L" ""
# ysupbl variables for grims shallow convection
state real WSTAR_YSU ij misc 1 - - "WSTAR_YSU" "mixed-layer velocity scale from ysupbl" "m/s"
state real DELTA_YSU ij misc 1 - - "DELTA_YSU" "entrainment layer depth from ysupbl" "m"
# MYJ PBL variables; GBM PBL: EXCH_H, EXCH_M
state real EXCH_H ikj misc 1 Z r "EXCH_H" "SCALAR EXCHANGE COEFFICIENTS " "m2 s-1"
state real EXCH_M ikj misc 1 Z r "EXCH_M" "EXCHANGE COEFFICIENTS " "m2 s-1"
state real CT ij misc 1 - r "CT" "COUNTERGRADIENT TERM" "K"
state real THZ0 ij misc 1 - r "THZ0" "POTENTIAL TEMPERATURE AT ZNT" "K"
state real Z0 ij misc 1 - r "Z0" "Background ROUGHNESS LENGTH" "m"
state real QZ0 ij misc 1 - r "QZ0" "SPECIFIC HUMIDITY AT ZNT" "kg kg-1"
state real UZ0 ij misc 1 - r "UZ0" "U WIND COMPONENT AT ZNT" "m s-1"
state real VZ0 ij misc 1 - r "VZ0" "V WIND COMPONENT AT ZNT" "m s-1"
state real QSFC ij misc 1 - r "QSFC" "SPECIFIC HUMIDITY AT LOWER BOUNDARY" "kg kg-1"
state real AKHS ij misc 1 - r "AKHS" "SFC EXCH COEFF FOR HEAT" "m s-1"
state real AKMS ij misc 1 - r "AKMS" "SFC EXCH COEFF FOR MOMENTUM" "m s-1"
state integer KPBL ij misc 1 - r "KPBL" "LEVEL OF PBL TOP" ""
#BSINGH - For CuP
#~wig: added real pbl level index for output tests
state real AKPBL ij misc 1 - r "AKPBL" "LEVEL OF PBL TOP" ""
#BSINGH - ENDS
state real TSHLTR ij misc 1 - r "TSHLTR" "SHELTER THETA FROM MYJ" "K"
state real QSHLTR ij misc 1 - r "QSHLTR" "SHELTER SPECIFIC HUMIDITY FROM MYJ" "kg kg-1"
state real PSHLTR ij misc 1 - r "PSHLTR" "SHELTER PRESSURE FROM MYJ" "Pa"
state real TH10 ij misc 1 - r "TH10" "10-M THETA FROM MYJ" "K"
state real Q10 ij misc 1 - r "Q10" "10-M SPECIFIC HUMIDITY FROM MYJ" "kg kg-1"
i1 real CHKLOWQ ij misc 1 - - "CHKLOWQ" "SURFACE SATURATION FLAG" ""
# MF_SHCONV variables
state real massflux_EDKF ikj misc 1 - - "MASS_FLUX" "MASS FLUX FROM EDKF" "Kg m s-1"
state real entr_EDKF ikj misc 1 - - "ENTR" "ENTRAINMENT FROM EDKF" "m-1"
state real detr_EDKF ikj misc 1 - - "DETR" "ENTRAINMENT FROM EDKF" "m-1"
state real thl_up ikj misc 1 - - "THL_UP" "THL OF UPDRAFT FROM EDKF" "K"
state real thv_up ikj misc 1 - - "THV_UP" "THL OF UPDRAFT FROM EDKF" "K"
state real rv_up ikj misc 1 - - "RV_UP" "RV OF UPDRAFT FROM EDKF" "kq/kg"
state real rt_up ikj misc 1 - - "RT_UP" "RT OF UPDRAFT FROM EDKF" "kq/kg"
state real rc_up ikj misc 1 - - "RC_UP" "RC OF UPDRAFT FROM EDKF" "kq/kg"
state real u_up ikj misc 1 - - "U_UP" "U OF UPDRAFT FROM EDKF" "m/s"
state real v_up ikj misc 1 - - "V_UP" "U OF UPDRAFT FROM EDKF" "m/s"
state real frac_up ikj misc 1 - - "FRAC_UP" "FRACTION OF UPDRAFT FROM EDKF" ""
state real rc_mf ikj misc 1 - r "RC_MF" "RC IN THE GRID COMPUTED BY EDKF" "kg/kg"
# TEMF PBL variables
state real te_temf ikj misc 1 - rh "te_temf" "Total energy from TEMF PBL scheme" "m2 s-2"
state real kh_temf ikj misc 1 - rh "kh_temf" "Diffusion coefficient for heat from TEMF PBL"
state real km_temf ikj misc 1 - rh "km_temf" "Diffusion coefficient for momentum from TEMF PBL"
state real shf_temf ikj misc 1 - rh "shf_temf" "Sensible heat flux from TEMF PBL" "K m s-1"
state real qf_temf ikj misc 1 - rh "qf_temf" "Sensible heat flux from TEMF PBL" "kg/kg m s-1"
state real uw_temf ikj misc 1 - rh "uw_temf" "U momentum flux from TEMF PBL" "m2 s-2"
state real vw_temf ikj misc 1 - rh "vw_temf" "V momentum flux from TEMF PBL" "m2 s-2"
state real wupd_temf ikj misc 1 - rh "wupd_temf" "Updraft velocity from TEMF PBL" "m s-1"
state real mf_temf ikj misc 1 - rh "mf_temf" "Mass flux from TEMF PBL" "m s-1"
state real thup_temf ikj misc 1 - rh "thup_temf" "Updraft thetal from TEMF PBL" "K"
state real qtup_temf ikj misc 1 - rh "qtup_temf" "Updraft qt from TEMF PBL" "1"
state real qlup_temf ikj misc 1 - rh "qlup_temf" "Updraft ql (liquid water) from TEMF PBL" "1"