修改USGS静态地形数据方法探讨

yaohll

辛苦楼主了 正好用到唉 不过还是有点难看懂

yaohll

如何改地面类型呢 我现在要把某一地区的地面全部换成snow and ice   求助

wenzerg

可以直接用read_wrf_nc修改，wrf-manual里有使用说明，不过还是羡慕直接写程序的牛人！

天目神眉

wenzerg 发表于 2013-7-20 14:09

                               
登录/注册后可看大图

可以直接用read_wrf_nc修改，wrf-manual里有使用说明，不过还是羡慕直接写程序的牛人！

是的，那个程序确实很方面！！！

lgg00828

program geogrid_clc
implicit none

integer :: i,j
integer :: isigned, endian, wordsize
integer :: nx, ny, nz
real :: scalefactor
real*8 :: xllcorner, yllcorner, cellsize, missvalue
character :: head12
real, allocatable :: rarray(:,:), iarray(:,:)

isigned = 1
endian = 0
wordsize = 2
scalefactor = 1.0
nz = 1


! read in the ascii new landuse data
open (10, file = 'file.asc')

!read in the header
read(10,*) head12, nx
read(10,*) head12, ny
read(10,*) head12, xllcorner
read(10,*) head12, yllcorner
read(10,*) head12, cellsize
read(10,*) head12, missvalue

allocate(rarray(nx,ny))
allocate(iarray(nx,ny))

!read in the data
do j = 1,ny
read(10,*) iarray(:,j)
end do

! reverse the data so that it begins at the lower-left corner
do j = 1,ny
rarray(:,j) = iarray(:,ny-(j-1))
enddo

!set the missing values
do j = 1, ny
do i = 1, nx
if ( rarray(i,j) < 0 ) then
rarray(i,j) = -99 ! set negative terrain to be zero since those are near coastal or river
banks
end if
end do
end do


call write_geogrid(rarray, nx, ny, nz, isigned, endian, scalefactor, wordsize)

end program
到这一步是怎么运行呢，把write_geogrid.c这个程序接在end program后面就好了吗，还是怎么弄得，程序write_geogrid.c里看不到输出到哪个目录，请指导，最好给我发一个脚本就好了！

清晨的雨

非常感谢

sywyx

随缘 发表于 2012-4-9 10:54

                               
登录/注册后可看大图

改地形有实际意义吗？好多人在质疑改地形好假，改变地形的那段高度的气压怎么补偿？

最近模拟一次暴雨过程,用了多个微物理方案,单参数和双参数的,报的过程降水结果,个人认为全部脱靶,
问题在哪里?
直觉是WRF模式对地形考虑的不到位.

适当的修改地形高度,是不是就不那么假了呢?也就不用考虑要素补偿问题了呢?

sywyx

直接修改USGS静态地形数据，怕是要有后遗症。

damien1987

看到楼主此文真的如获至宝，但是看完又没看懂，请教楼主一下：
其实我是学遥感的，不是做数值预报，但是要为别人更新一个区域的静态数据也就是landuse，目前我得到两个文件一个是土地利用现状还有dem文件，调用什么模块可以把这二者写成可输入的二进制文件？write_geogrid.c对无编程基础的人来说可以胜任吗？
ps：我还在尝试安装WRF的阶段

andrewsoong

sywyx 发表于 2014-3-23 19:14
最近模拟一次暴雨过程,用了多个微物理方案,单参数和双参数的,报的过程降水结果,个人认为全部脱靶,
问题 ...

你是觉得原有的地形高度太低了？应该增大地形高度？这样结果如何？