爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 1731|回复: 2

【求助】用ncl处理grib1数据来求水汽通量,变量高度不一致该怎么处理?

[复制链接]
发表于 2021-9-7 17:26:56 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
$ ncl_filedump -v /home/lenovo/fnl_20000101_00_00.grib1
ncl_filedump: warning: no variable(s) specified.
Copyright (C) 1995-2017 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.4.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
Variable: f
Type: file
filename:       fnl_20000101_00_00
path:   /home/lenovo/fnl_20000101_00_00.grib1
   file global attributes:
   dimensions:
      lat_3 = 181
      lon_3 = 360
      lv_ISBL2 = 26
      lv_GPML3 = 3
      lv_DBLY4 = 2
      lv_ISBL5 = 2
      lv_ISBL6 = 21
      lv_SIGY7 = 4
      lv_ISBL8 = 6
   variables:
      float PRES_3_SFC ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Pressure
         units :        Pa
         _FillValue :   1e+20
         level_indicator :      1
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     1
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float PRES_3_MWSL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Pressure
         units :        Pa
         _FillValue :   1e+20
         level_indicator :      6
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     1
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float PRES_3_TRO ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Pressure
         units :        Pa
         _FillValue :   1e+20
         level_indicator :      7
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     1
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float PRES_3_CCBL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Pressure
         units :        Pa
         _FillValue :   1e+20
         level_indicator :      242
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     1
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float PRES_3_CCTL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Pressure
         units :        Pa
         _FillValue :   1e+20
         level_indicator :      243
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     1
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float PRMSL_3_MSL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Pressure reduced to MSL
         units :        Pa
         _FillValue :   1e+20
         level_indicator :      102
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     2
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float HGT_3_SFC ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Geopotential height
         units :        gpm
         _FillValue :   1e+20
         level_indicator :      1
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     7
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float HGT_3_MWSL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Geopotential height
         units :        gpm
         _FillValue :   1e+20
         level_indicator :      6
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     7
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float HGT_3_TRO ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Geopotential height
         units :        gpm
         _FillValue :   1e+20
         level_indicator :      7
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     7
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float HGT_3_ISBL ( lv_ISBL2, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Geopotential height
         units :        gpm
         _FillValue :   1e+20
         level_indicator :      100
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     7
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TOZONE_3_EATM ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Total ozone
         units :        Dobson
         _FillValue :   1e+20
         level_indicator :      200
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     10
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        0
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_SFC ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      1
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_MWSL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      6
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_TRO ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      7
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_ISBL ( lv_ISBL2, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      100
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_GPML ( lv_GPML3, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      103
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_HTGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      105
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        2
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_SIGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      107
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        9950
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_DBLY ( lv_DBLY4, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      112
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float TMP_3_SPDY ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      116
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     11
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        ( 30, 0 )
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float POT_3_SIGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Potential temperature
         units :        K
         _FillValue :   1e+20
         level_indicator :      107
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     13
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        9950
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float GPA_3_ISBL ( lv_ISBL5, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Geopotential height anomaly
         units :        gpm
         _FillValue :   1e+20
         level_indicator :      100
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     27
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float U_GRD_3_MWSL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    u-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      6
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     33
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float U_GRD_3_TRO ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    u-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      7
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     33
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float U_GRD_3_ISBL ( lv_ISBL2, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    u-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      100
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     33
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float U_GRD_3_GPML ( lv_GPML3, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    u-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      103
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     33
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float U_GRD_3_HTGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    u-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      105
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     33
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        10
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float U_GRD_3_SIGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    u-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      107
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     33
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        9950
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float U_GRD_3_SPDY ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    u-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      116
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     33
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        ( 30, 0 )
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_GRD_3_MWSL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    v-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      6
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     34
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_GRD_3_TRO ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    v-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      7
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     34
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_GRD_3_ISBL ( lv_ISBL2, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    v-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      100
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     34
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_GRD_3_GPML ( lv_GPML3, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    v-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      103
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     34
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_GRD_3_HTGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    v-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      105
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     34
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        10
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_GRD_3_SIGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    v-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      107
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     34
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        9950
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_GRD_3_SPDY ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    v-component of wind
         units :        m/s
         _FillValue :   1e+20
         level_indicator :      116
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     34
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        ( 30, 0 )
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_VEL_3_ISBL ( lv_ISBL6, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Pressure vertical velocity
         units :        Pa/s
         _FillValue :   1e+20
         level_indicator :      100
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     39
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float V_VEL_3_SIGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Pressure vertical velocity
         units :        Pa/s
         _FillValue :   1e+20
         level_indicator :      107
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     39
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        9950
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float ABS_V_3_ISBL ( lv_ISBL2, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Absolute vorticity
         units :        /s
         _FillValue :   1e+20
         level_indicator :      100
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     41
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float SPF_H_3_HTGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Specific humidity
         units :        kg/kg
         _FillValue :   1e+20
         level_indicator :      105
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     51
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        2
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float SPF_H_3_SPDY ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Specific humidity
         units :        kg/kg
         _FillValue :   1e+20
         level_indicator :      116
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     51
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        ( 30, 0 )
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float R_H_3_ISBL ( lv_ISBL6, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Relative humidity
         units :        %
         _FillValue :   1e+20
         level_indicator :      100
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     52
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float R_H_3_HTGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Relative humidity
         units :        %
         _FillValue :   1e+20
         level_indicator :      105
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     52
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        2
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float R_H_3_SIGL ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Relative humidity
         units :        %
         _FillValue :   1e+20
         level_indicator :      107
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     52
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        9950
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float R_H_3_SIGY ( lv_SIGY7, lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Relative humidity
         units :        %
         _FillValue :   1e+20
         level_indicator :      108
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     52
         model :        Analysis from GDAS (Global Data Assimilation System)
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float R_H_3_SPDY ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Relative humidity
         units :        %
         _FillValue :   1e+20
         level_indicator :      116
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     52
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        ( 30, 0 )
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00)
      float R_H_3_EATM ( lat_3, lon_3 )
         center :       US National Weather Service - NCEP (WMC)
         long_name :    Relative humidity
         units :        %
         _FillValue :   1e+20
         level_indicator :      200
         grid_number :  3
         parameter_table_version :      2
         parameter_number :     52
         model :        Analysis from GDAS (Global Data Assimilation System)
         level :        0
         forecast_time :        0
         forecast_time_units :  hours
         initial_time : 01/01/2000 (00:00) 搜狗截图20210907170921.jpg

搜狗截图20210907172717.jpg
密码修改失败请联系微信:mofangbao
发表于 2021-9-7 17:53:32 | 显示全部楼层
和U和V提取为比湿相同层
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-9-7 17:57:14 | 显示全部楼层
werewolf 发表于 2021-9-7 17:53
和U和V提取为比湿相同层

老师我提取了某个层的u/v/rh,但是不知道为什么总是出问题,难道求平均不能像这样弄吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表