教你转Noah-MP陆面模式，获取陆面多种变量

novelbean · 发表于 2019-11-3 19:39:22

学习学习，好久没有整陆面了

zxc521122 · 发表于 2019-11-5 15:57:05

您好，我按照步驟進行梨縣的Noah-MP模擬
使用的是GLDAS的資料當初始場條件，
將前置處理都處理好後跑出LDASIN files在LDASIN的資料夾裡
然後運行./hrldas.exe
跑出了此行error message 如下
NetCDF: Attribute not found

***** READ_HRLDAS_HDRINFO: Problems finding global attribute 'GRID_ID' or 'grid_id'

請問出了什麼問題啊.....

andrewsoong · 发表于 2019-11-5 21:24:02

zxc521122 发表于 2019-11-5 15:57
您好，我按照步驟進行梨縣的Noah-MP模擬
使用的是GLDAS的資料當初始場條件，
將前置處理都處理好後跑出LD ...

你的文件输入量没有包含属性

andrewsoong · 发表于 2019-11-5 21:24:10

zxc521122 发表于 2019-11-5 15:57
您好，我按照步驟進行梨縣的Noah-MP模擬
使用的是GLDAS的資料當初始場條件，
將前置處理都處理好後跑出LD ...

你的文件输入量没有包含属性

zxc521122 · 发表于 2019-11-7 17:06:08

很不好意思的我又來發問了
想請問FULL_IC_FRQ    : set to 0 if you want initial conditions in the init file, -1 otherwise
當我設0的時候會出現下面的error
[raymon34@air-cluster HRLDAS_forcing]$ ./create_forcing.exe namelist.input.GLDAS
  11 | 1  | 0  |    | T2D    | K       | Temperature    at 2 m             |  0  |  0  |  0  | 103 |
  51 | 1  | 0  |    | Q2D    | kg/kg    | Specific Humidity at 2 m             |  0  |  1  |  0  | 103 |
  33 | 1  | 0  |    | U2D    | m/s    | U                at 10 m             |  0  |  2  |  2  | 103 |
  34 | 1  | 0  |    | V2D    | m/s    | V                at 10 m             |  0  |  2  |  3  | 103 |
1 | 1  | 0  |    | PSFC    | Pa       | Surface Pressure                      |  0  |  3  |  0  | 1 |
  59 | 1  | 0  |    | RAINRATE | kg/m^2/s  | Precipitation Rate                   |  0  |  1  |  8  | 1 |
204 | 1  | 0  |    | SWDOWN | W/m^2    | Downward short-wave radiation flux    |  0  |  4  | 192 | 1 |
205 | 1  | 0  |    | LWDOWN | W/m^2    | Downward long-wave radiation flux    |  0  |  5  | 192 | 1 |
7 | 1  | 0  |    | TERRAIN  | m       | Terrain field of source analysis       |  2  |  0  |  7  | 1 |
138 | 1  | 0  |    | TSK    | K       | Skin temperature                      |  0  |  0  |  0  | 1 |
  65 | 1  | 0  |    | SNOW    | kg/m^2 | Water equivalent snow depth          |  0  |  1  | 13  | 1 |
  71 | 1  | 0  |    | CANWAT | kg/m^2 | Plant Canopy Surface Water             |  2  |  0  | 196 | 1 |
  86 | 112  | 0  | 4  | SMOIS_1  | gldas    | Soil Moist 0-10 cm below grn layer (Up) |  2  |  0  | 192 | 106 |
  86 | 112  | 0  | 3  | SMOIS_2  | gldas    | Soil Moist 10-40 cm below grn layer    |  2  |  0  | 192 | 106 |
  86 | 112  | 0  | 2  | SMOIS_3  | gldas    | Soil Moist 40-100 cm below grn layer |  2  |  0  | 192 | 106 |
  86 | 112  | 0  | 1  | SMOIS_4  | gldas    | Soil Moist 100-200 cm below gr layer |  2  |  0  | 192 | 106 |
  85 | 112  | 0  | 4  | STEMP_1  | K       | T 0-10 cm below ground layer (Upper) |  2  |  0  |  2  | 106 |
  85 | 112  | 0  | 3  | STEMP_2  | K       | T 10-40 cm below ground layer (Upper) |  2  |  0  |  2  | 106 |
  85 | 112  | 0  | 2  | STEMP_3  | K       | T 40-100 cm below ground layer (Upper)  |  2  |  0  |  2  | 106 |
  85 | 112  | 0  | 1  | STEMP_4  | K       | T 100-200 cm below ground layer (Bottom)|  2  |  0  |  2  | 106 |
geo_em_flnm = "/work/GROUP/raymon/WRF3.7.1/WPS3.9.1/geo_em.d02.nc"
flnm = /work/GROUP/raymon/WRF3.7.1/WPS3.9.1/geo_em.d02.nc
Done with subroutine read_geo_em_file
Date = 2014-10-01_00  ihour =          0
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/Tair/GLDAS_Tair.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/Tair/GLDAS_Tair.2014100100.grb
T2D       GRIB Editon 1
         :  Checking for file '/work/GROUP/raymon/hrldas-release-3.8/HRLDAS/HRLDAS_forcing/run/examples/GLDAS/GLDAS_ELEVATION.grb'
         :  Found file /work/GROUP/raymon/hrldas-release-3.8/HRLDAS/HRLDAS_forcing/run/examples/GLDAS/GLDAS_ELEVATION.grb
TERRAIN GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/Qair/GLDAS_Qair.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/Qair/GLDAS_Qair.2014100100.grb
Q2D       GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/Psurf/GLDAS_Psurf.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/Psurf/GLDAS_Psurf.2014100100.grb
PSFC    GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/U/GLDAS_U.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/U/GLDAS_U.2014100100.grb
U2D       GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/V/GLDAS_V.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/V/GLDAS_V.2014100100.grb
V2D       GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/LWdown/GLDAS_LWdown.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/LWdown/GLDAS_LWdown.2014100100.grb
LWDOWN    GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/SWdown24/GLDAS_SWdown.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/SWdown24/GLDAS_SWdown.2014100100.grb
SWDOWN    GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/SWdown24/GLDAS_SWdown.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/SWdown24/GLDAS_SWdown.2014100100.grb
SWDOWN    GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/Precip/GLDAS_Precip.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/Precip/GLDAS_Precip.2014100100.grb
RAINRATE GRIB Editon 1
         :  Checking for file '/backup2/raymon/offline_mp/2014_extracted/Precip/GLDAS_Precip.2014100100.grb'
         :  Found file /backup2/raymon/offline_mp/2014_extracted/Precip/GLDAS_Precip.2014100100.grb
RAINRATE GRIB Editon 1
Warning: ieee_inexact is signaling
String is all blanks.

好像並沒有完全地做完初始場，導致我後面的步驟沒辦法繼續下去，請問該如何解決
那如果設成-1又會有什麼影響呢?

andrewsoong · 发表于 2019-11-7 17:24:13

zxc521122 发表于 2019-11-7 17:06
很不好意思的我又來發問了
想請問FULL_IC_FRQ : set to 0 if you want initial conditions in the i ...

有没有更详细的报错啊？你模式是不是用的PGI编译器？

zxc521122 · 发表于 2019-11-7 23:25:46

andrewsoong 发表于 2019-11-7 17:24
有没有更详细的报错啊？你模式是不是用的PGI编译器？

我是用PGI編譯器沒錯
更詳細的error資料就只有如下

geo_em_flnm = "/raid/raymon/raw_GLDASv1/2014_extracted/geo_em.d02.nc"
flnm = /raid/raymon/raw_GLDASv1/2014_extracted/geo_em.d02.nc
Done with subroutine read_geo_em_file
Date = 2014-10-01_00  ihour =          0
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/Tair/GLDAS_Tair.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/Tair/GLDAS_Tair.2014100100.grb
T2D       GRIB Editon 1
         :  Checking for file '/raid/raymon/hrldas-release-3.8/HRLDAS/HRLDAS_forcing/run/examples/GLDAS/GLDAS_ELEVATION.grb'
         :  Found file /raid/raymon/hrldas-release-3.8/HRLDAS/HRLDAS_forcing/run/examples/GLDAS/GLDAS_ELEVATION.grb
TERRAIN GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/Qair/GLDAS_Qair.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/Qair/GLDAS_Qair.2014100100.grb
Q2D       GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/Psurf/GLDAS_Psurf.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/Psurf/GLDAS_Psurf.2014100100.grb
PSFC    GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/U/GLDAS_U.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/U/GLDAS_U.2014100100.grb
U2D       GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/V/GLDAS_V.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/V/GLDAS_V.2014100100.grb
V2D       GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/LWdown/GLDAS_LWdown.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/LWdown/GLDAS_LWdown.2014100100.grb
LWDOWN    GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/SWdown24/GLDAS_SWdown.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/SWdown24/GLDAS_SWdown.2014100100.grb
SWDOWN    GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/SWdown24/GLDAS_SWdown.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/SWdown24/GLDAS_SWdown.2014100100.grb
SWDOWN    GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/Precip/GLDAS_Precip.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/Precip/GLDAS_Precip.2014100100.grb
RAINRATE GRIB Editon 1
         :  Checking for file '/raid/raymon/raw_GLDASv1/2014_extracted/Precip/GLDAS_Precip.2014100100.grb'
         :  Found file /raid/raymon/raw_GLDASv1/2014_extracted/Precip/GLDAS_Precip.2014100100.grb
RAINRATE GRIB Editon 1
         2          2          2          2
         16          16          16          16
         14          14          14          14
Warning: ieee_inexact is signaling
String is all blanks.

另外附上我的namelist.input.GLDAS
&files
STARTDATE       = "2014-10-01_00"
ENDDATE          = "2014-12-31_23"
DataDir          = "/raid/raymon/raw_GLDASv1/2014_extracted"
OutputDir       = "/raid/raymon/raw_GLDASv1/2014_LDASIN/"
FULL_IC_FRQ       = 0
RAINFALL_INTERP = 0
RESCALE_SHORTWAVE  = .FALSE.
UPDATE_SNOW       = .FALSE.
FORCING_HEIGHT_2D  = .FALSE.
TRUNCATE_SW       = .FALSE.
EXPAND_LOOP       = 8
INIT_LAI          = .TRUE.
VARY_LAI          = .TRUE.
MASK_WATER       = .TRUE.

geo_em_flnm       = "/raid/raymon/raw_GLDASv1/2014_extracted/geo_em.d02.nc"

Zfile_template    = "/raid/raymon/hrldas-release-3.8/HRLDAS/HRLDAS_forcing/run/examples/GLDAS/GLDAS_ELEVATION.grb"

Tfile_template    = "<DataDir>/Tair/GLDAS_Tair.<date>.grb"
Ufile_template    = "<DataDir>/U/GLDAS_U.<date>.grb",
Vfile_template    = "<DataDir>/V/GLDAS_V.<date>.grb",
Pfile_template    = "<DataDir>/Psurf/GLDAS_Psurf.<date>.grb",
Qfile_template    = "<DataDir>/Qair/GLDAS_Qair.<date>.grb",
LWfile_template = "<DataDir>/LWdown/GLDAS_LWdown.<date>.grb",
SWfile_primary    = "<DataDir>/SWdown24/GLDAS_SWdown.<date>.grb",
SWfile_secondary = "<DataDir>/SWdown24/GLDAS_SWdown.<date>.grb",
PCPfile_primary = "<DataDir>/Precip/GLDAS_Precip.<date>.grb",
PCPfile_secondary  = "<DataDir>/Precip/GLDAS_Precip.<date>.grb",

CANWTfile_template = "<DataDir>/INIT/GLDAS_Canopint.<date>.grb",
SKINTfile_template = "<DataDir>/INIT/GLDAS_AvgSurfT.<date>.grb",

STfile_template = "<DataDir>/INIT/GLDAS_TSoil_000-010.<date>.grb",
                  "<DataDir>/INIT/GLDAS_TSoil_010-040.<date>.grb",
                  "<DataDir>/INIT/GLDAS_TSoil_040-100.<date>.grb",
                  "<DataDir>/INIT/GLDAS_TSoil_100-200.<date>.grb",

SMfile_template = "<DataDir>/INIT/GLDAS_SoilM_000-010.<date>.grb",
                  "<DataDir>/INIT/GLDAS_SoilM_010-040.<date>.grb",
                  "<DataDir>/INIT/GLDAS_SoilM_040-100.<date>.grb",
                  "<DataDir>/INIT/GLDAS_SoilM_100-200.<date>.grb",

andrewsoong · 发表于 2019-11-8 00:08:22

zxc521122 发表于 2019-11-7 23:25
我是用PGI編譯器沒錯
更詳細的error資料就只有如下

pgi编译器的确会出现这个问题

zxc521122 · 发表于 2019-11-8 09:33:11

本帖最后由 zxc521122 于 2019-11-10 18:21 编辑

andrewsoong 发表于 2019-11-8 00:08
pgi编译器的确会出现这个问题

哇............
那請問我該怎麼辦..........我的pgi編譯器版本為18.10
很抱歉我積分不足無法用私訊功能，或許可以考慮加個微信或信箱來討論相關問題可以嗎?

zxc521122 · 发表于 2019-11-10 18:19:15

請問用PGI編譯器無法順利做好初始場，該怎麼解決

		自动登录	找回密码
密码			立即注册