城市冠层模型UCM与城市土地利用的修改

Kecheng

pigzero527 发表于 2025-7-27 22:34
LANDUSE.TBL里要进行相应的更改，比如找到modis的类型，然后给对应的31、32、33分配上low,middle,high
31, ...

您好我还有个问题想请教一下，num_land_cat是否需要修改呢？依然还是使用21吗？

pigzero527

Kecheng 发表于 2025-7-27 22:51
您好我还有个问题想请教一下，num_land_cat是否需要修改呢？依然还是使用21吗？

应该是不用修改，你试试看吧

Kecheng

                               
登录/注册后可看大图

Kecheng

pigzero527 发表于 2025-7-27 23:06
应该是不用修改，你试试看吧

您好，因为回复里好像贴不了图片，所以我把修改之后的LANDUSE.TBL和VEGPARM.TBL的部分截图贴在第13层楼了。想请问一下我这样修改对吗？还是修改完LANDUSE.TBL和VEGPARM.TBL需要把wrf重新编译呀？我目前这样修改完之后，然后用python修改了geo_em的类别，met_em依然会改成31、32、33，但是运行完real，wrfinput又变成13了

Kecheng

pigzero527 发表于 2025-7-27 23:06
应该是不用修改，你试试看吧

您好，想再请问一下：
我把surface_input_source改成3之后，wrfinput目前已经改成了31、32、33
不过我看到运行wrf的输出里有

USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY

和

USING DEFAULT URBAN MORPHOLOGY
WARNING, FRC_URB2D = 0 BUT IVGTYP IS URBAN
WARNING, THE URBAN FRACTION WILL BE READ FROM URBPARM.TBL
USING DEFAULT URBAN MORPHOLOGY
WARNING, FRC_URB2D = 0 BUT IVGTYP IS URBAN
WARNING, THE URBAN FRACTION WILL BE READ FROM URBPARM.TBL
USING DEFAULT URBAN MORPHOLOGY
WARNING, FRC_URB2D = 0 BUT IVGTYP IS URBAN
WARNING, THE URBAN FRACTION WILL BE READ FROM URBPARM.TBL

这样两段文字
想请问一下，为什么它会读取默认的形态参数？不应该都读取URBPARM.TBL吗？

pigzero527

我是直接改的wrfinput里的LU_INDEX和LANDUSEF，再运行wrf，wrfout里的lu_index就保留了31，32，33，改geo_em里的我没试过
“USING DEFAULT URBAN MORPHOLOGY
WARNING, FRC_URB2D = 0 BUT IVGTYP IS URBAN
WARNING, THE URBAN FRACTION WILL BE READ FROM URBPARM.TBL”
这里面warning的意思是你的地表文件里没有urban fraction，但是IVGTYP显示的是城市，那就读取URBPARM.TBL里设定的数值，31、32、33分别对应0.5，0.9，0.95

Kecheng

pigzero527 发表于 2025-7-29 21:12
我是直接改的wrfinput里的LU_INDEX和LANDUSEF，再运行wrf，wrfout里的lu_index就保留了31，32，33，改geo_e ...

好的，谢谢您的帮助，我现在已经成功修改并运行啦。

不过根据我今天下午在NCAR论坛上的交流，貌似虽然rsl.error会输出WARNING, THE URBAN FRACTION WILL BE READ FROM URBPARM.TBL，但实际上会直接设定为0.9而非读取URAPARM.TBL里设定的值

在WRF/dyn_em/ module_initialize_real.F中有：
IF ( grid%FRC_URB2D(i,j) == 0. ) THEN
                 IF ( (MMINLU == 'NLCD40' .OR. MMINLU == 'MODIFIED_IGBP_MODIS_NOAH') .AND. &
                   (grid%ivgtyp(i,j)==24 .OR. grid%ivgtyp(i,j)==25 .OR. grid%ivgtyp(i,j)==26 .OR. grid%ivgtyp(i,j)==13) ) grid%FRC_URB2D(i,j) = 0.9
                 IF ( MMINLU == 'USGS' .AND. grid%ivgtyp(i,j)==1 ) grid%FRC_URB2D(i,j) = 0.9
              ENDIF
貌似需要注释掉这几行才会从Table里读取

tbag

Kecheng 发表于 2025-7-30 00:06
好的，谢谢您的帮助，我现在已经成功修改并运行啦。

不过根据我今天下午在NCAR论坛上的交流，貌似虽然 ...

兄弟，能分享一下如何将urban built-up按照人口分布分解成low/high/commercial 的方法或者论文吗，多谢。