wrf-hydro运行问题

小白猪

请教一下各位大佬，运行wrf-hydro模型时，设置的起始时间是2011年1月1日00时，但是运行的时候，读取forcing数据时显示读取的是1月1日03时的数据，请问这是什么原因呢？

read forcing data at 2011-01-01_00:00:00./FORCING/2011010103.LDASIN_DOMAIN1

小白猪

在运行的时候，还是读取下一个时次的FORCING文件，日志中显示报错信息如下
emitted longwave <0; skin T may be wrong due to inconsistent
input of SHDFAC with LAI
         245         203 SHDFAC=  0.0000000E+00 VAI=  0.0000000E+00 TV=
   266.0056     TG=  -10.11275
LWDN=  -9999.000     FIRA=   9699.031     SNOWH=  0.0000000E+00
FATAL ERROR: STOP in Noah-MP

请教各位大佬这是什么原因造成的呢

1300214974

你好，不知道问题是否已经解决：
“emitted longwave <0; skin T may be wrong due to inconsistent
input of SHDFAC with LAI
         245         203 SHDFAC=  0.0000000E+00 VAI=  0.0000000E+00 TV=
   266.0056     TG=  -10.11275
LWDN=  -9999.000     FIRA=   9699.031     SNOWH=  0.0000000E+00
FATAL ERROR: STOP in Noah-MP”

这个报错应该是由于驱动数据在陆地表面为缺测导致的，只需要对驱动数据当中，landmask==1（陆地）的部分对缺测值进行填充即可。附可能用到的缺测值填充函数：

def replace_missing_values(array, missing_value):
    array = np.nan_to_num(array, nan=missing_value)
    rows, cols = array.shape
    for i in range(rows):
        for j in range(cols):
            if array[i, j] == missing_value:
                # Find valid neighbors
                neighbors = []
                if i > 0 and array[i-1, j] != missing_value:
                    neighbors.append(array[i-1, j])
                if i < rows-1 and array[i+1, j] != missing_value:
                    neighbors.append(array[i+1, j])
                if j > 0 and array[i, j-1] != missing_value:
                    neighbors.append(array[i, j-1])
                if j < cols-1 and array[i, j+1] != missing_value:
                    neighbors.append(array[i, j+1])
               
                # Replace missing value with a random valid neighbor
                if neighbors:
                    array[i, j] = np.random.choice(neighbors)
                else:
                    # If no valid neighbors found, replace with a random value
                    # Adjust the range as needed
                    array[i, j] = np.random.randint(1, 10)
   
    return array

qwelzx

小白猪 发表于 2023-10-7 22:07
在运行的时候，还是读取下一个时次的FORCING文件，日志中显示报错信息如下
emitted longwave

请问大佬这个问题你怎么解决的呢？我也遇到相同问题了
emitted longwave <0; skin T may be wrong due to inconsistent
input of SHDFAC with LAI
          61           1 SHDFAC=  0.895097554     VAI=   4.58524609     TV=   296.341431     TG=   329.706451   
LWDN=   367.595581     FIRA=  -417.277618     SNOWH=   0.00000000   
FATAL ERROR: STOP in Noah-MP