【求助】WRF 模拟城市扩张，结果和预期相反

wz123

之前做城市化方案，将土地利用类型分为33种，城市类型为31,32,33.模拟城市扩张后温度状况，结果出来后发现：1城市区域的温度低于周边郊区的温度，2城市扩张后温度反而降低了！不知道是什么地方设置出了问题，请各位多指教！
case1是模拟城市现状的：
case3是城市扩张后的：
下面是我的namelist.input

1. 
   
2. <P> &time_control
   
3. run_days                            = 3,
   
4. run_hours                           = 0,
   
5. run_minutes                         = 0,
   
6. run_seconds                         = 0,
   
7. start_year                          = 2012, 2012, 2012,
   
8. start_month                         = 01,   01,   01,
   
9. start_day                           = 07,   07,   07,
   
10. start_hour                          = 18,   18,   18,
    
11. start_minute                        = 00,   00,   00,
    
12. start_second                        = 00,   00,   00,
    
13. end_year                            = 2012, 2012, 2012,
    
14. end_month                           = 01,   01,   01,
    
15. end_day                             = 10,   10,   10,
    
16. end_hour                            = 18,   18,   18,
    
17. end_minute                          = 00,   00,   00,
    
18. end_second                          = 00,   00,   00,
    
19. interval_seconds                    = 21600
    
20. input_from_file                     = .true.,.true.,.true.,
    
21. history_interval                    = 60,  60,   60,
    
22. frames_per_outfile                  = 1000, 1000, 1000,
    
23. restart                             = .false.,
    
24. restart_interval                    = 1440,
    
25. write_hist_at_0h_rst                = .true.
    
26. io_form_history                     = 2
    
27. io_form_restart                     = 2
    
28. io_form_input                       = 2
    
29. io_form_boundary                    = 2
    
30. debug_level                         = 0
    
31. /</P>
    
32. <P> &domains
    
33. time_step                           = 27,
    
34. time_step_fract_num                 = 0,
    
35. time_step_fract_den                 = 1,
    
36. max_dom                             = 3,
    
37. e_we                                = 103,    103,   103,
    
38. e_sn                                = 121,    121,    121,
    
39. s_vert                              = 1,     1,     1,
    
40. e_vert                              = 35,    35,    35,
    
41. ztop                                = 12000, 12000,
    
42. p_top_requested                     = 5000,
    
43. num_metgrid_levels                  = 27,
    
44. num_metgrid_soil_levels             = 4,
    
45. dx                                  = 4500, 1500,  500,
    
46. dy                                  = 4500, 1500,  500,
    
47. grid_id                             = 1,     2,     3,
    
48. parent_id                           = 1,     1,     2,
    
49. i_parent_start                      = 1,     35,    35,
    
50. j_parent_start                      = 1,     41,    41,
    
51. parent_grid_ratio                   = 1,     3,     3,
    
52. parent_time_step_ratio              = 1,     3,     3,
    
53. feedback                            = 1,
    
54. smooth_option                       = 0
    
55. /</P>
    
56. <P> &physics
    
57. mp_physics                          = 8,     8,     8,
    
58. ra_lw_physics                       = 1,     1,     1,
    
59. ra_sw_physics                       = 2,     2,     2,
    
60. radt                                = 4.5,   4.5,   4.5,
    
61. sf_sfclay_physics                   = 2,     2,     2,
    
62. sf_surface_physics                  = 2,     2,     2,
    
63. bl_pbl_physics                      = 2,     2,     2,
    
64. bldt                                = 0,     0,     0,
    
65. cu_physics                          = 0,     0,     0,
    
66. cudt                                = 5,     5,     5,
    
67. isfflx                              = 1,
    
68. ifsnow                              = 0,
    
69. icloud                              = 1,
    
70. surface_input_source                = 3,
    
71. num_land_cat                         = 33,
    
72. num_soil_layers                     = 4,
    
73. sf_urban_physics                    = 1,     1,     1,
    
74. maxiens                             = 1,
    
75. maxens                              = 3,
    
76. maxens2                             = 3,
    
77. maxens3                             = 16,
    
78. ensdim                              = 144,
    
79. /</P>
    
80. <P> &fdda
    
81. /</P>
    
82. <P> &dynamics
    
83. w_damping                           = 0,
    
84. diff_opt                            = 1,      1,      1,
    
85. km_opt                              = 4,      4,      4,
    
86. diff_6th_opt                        = 0,      0,      0,
    
87. diff_6th_factor                     = 0.12,   0.12,   0.12,
    
88. base_temp                           = 290.
    
89. damp_opt                            = 0,
    
90. zdamp                               = 5000.,  5000.,  5000.,
    
91. dampcoef                            = 0.2,    0.2,    0.2
    
92. khdif                               = 0,      0,      0,
    
93. kvdif                               = 0,      0,      0,
    
94. non_hydrostatic                     = .true., .true., .true.,
    
95. moist_adv_opt                       = 1,      1,      1,     
    
96. scalar_adv_opt                      = 1,      1,      1,     
    
97. /</P>
    
98. <P> &bdy_control
    
99. spec_bdy_width                      = 5,
    
100. spec_zone                           = 1,
     
101. relax_zone                          = 4,
     
102. specified                           = .true., .false.,.false.,
     
103. nested                              = .false., .true., .true.,
     
104. /</P>
     
105. <P> &grib2
     
106. /</P>
     
107. <P> &namelist_quilt
     
108. nio_tasks_per_group = 0,
     
109. nio_groups = 1,
     
110. / </P>

复制代码

刘振

这只是一个时次的温度是么，其实我想问热岛的话，是早上和晚上都热。城市夸张确认只是修改landuse么，是不是还有land fraction变量啥的，修改土地利用主要是通过albedo和 emissivity吧，你可以查看这些相关的变量是否有变化，仅仅是我瞎扯，你参考下

andrewsoong

不应该是城市热到么效应么？

xyan88

请问run_days为什么设置为3，不应该是 end - start 么  我看你多这个是从1月份到8月份  

xyan88

请问run_days为什么设置为3，不应该是 end - start 么  我看你的这个是从1月份到8月份  

wz123

xyan88 发表于 2016-1-6 19:28
请问run_days为什么设置为3，不应该是 end - start 么  我看你的这个是从1月份到8月份

不好意思，笔误，应该是1月7日至1月10日共3天

wz123

刘振 发表于 2016-1-6 16:49
这只是一个时次的温度是么，其实我想问热岛的话，是早上和晚上都热。城市夸张确认只是修改landuse么，是不 ...

谢谢您的回复，
1其实我是做了三天的，只贴出了一个时次的结果，每个时次都是这种情况，城市温度低于周围郊区温度，不知为何。
2城市扩张我修改了geo_em文件的LU_INDEX，LANDUSEF,LANDMASK,其他城市化参数在URBPARM.TBL中修改，不知道程序有没有问题？
3您提到了反照率辐射率参数是应该在哪儿修改呢？

刘振

wz123 发表于 2016-1-7 21:46
谢谢您的回复，
1其实我是做了三天的，只贴出了一个时次的结果，每个时次都是这种情况，城市温度低于周 ...

额，先不管你的城市扩张实验，你的case1是什么都没改的情况？如果什么都没改，而且你的模拟没错且物理机制正确的话，那你得城市扩张实验就不矛盾了。
说到你的原始模拟，我知道的是，地表温度主要是白天受反照率（albedo）的影响吸收短波辐射影响温度，晚上是受发射率emissivity放出长波辐射，温度降低。你先看白天的吧，albedo应该是有输出的。而且你每步都要检查你程序是否出错的话，看下wrfinput的landuse，land_index等变量是不是你预期的结果。至于反照率的大小应该是根据landuse这些确定的吧，你得具体查下。
热岛效应是白天晚上都存在的？根据你读的文章？

鹰隼

感谢楼主分享

xinyao

建议统一图例再看