请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 266|回复: 1

WRF3DVAR+AHI同化遇到问题:同化成功,但是同化不进去数据,

[复制链接]

新浪微博达人勋

发表于 2024-2-23 10:58:00 | 显示全部楼层 |阅读模式
10金钱
本帖最后由 whuwjr 于 2024-2-26 12:34 编辑

请各位朋友帮忙指点指点,万分感谢
taskid: 0 hostname: s1804.para.bscc
Namelist logging not found in namelist.input. Using registry defaults for varia
bles in logging.

*** VARIATIONAL ANALYSIS ***
    WRFDA V4.5.2                                                                    

Ntasks in X            8 , ntasks in Y            8
*************************************
Parent domain
ids,ide,jds,jde            1         196           1         184
ims,ime,jms,jme           -4          32          -4          30
ips,ipe,jps,jpe            1          25           1          23
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,               48814276  bytes allocated
   Input data is acceptable to use: fg
hybrid_opt =  2

use_theta_m =  1

Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE     25 JS      1 JE     23
WRF NUMBER OF TILES =   1
Domain mapping info:
map_proj =     3
cen_lat   =    0.275000E+02
cen_lon   =    0.115500E+03
truelat1  =    0.300000E+02
truelat2  =    0.600000E+02
start_lat =    0.197139E+02
start_lon =    0.106433E+03
pole_lat  =    0.900000E+02
dsm       =    0.900000E+01

Set up observations (ob)

thin_conv_opt =  0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

Using ASCII format observation input

--------------------------- WARNING ---------------------------
WARNING FROM FILE:  da_scan_obs_ascii.inc  LINE:      69
Error   29 opening gts obs file ob.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE:  da_read_obs_ascii.inc  LINE:     107
Error   29 opening gts obs file ob.ascii
---------------------------------------------------------------
Using NCEP BUFR radiance 1b input

        3000
Reading extra data for AHI
Channel       BTLim      ca1     ca2     clearSkyBias
      1          0.000    0.000    0.000    0.000
      2          0.000    0.000    0.000    0.000
      3          0.000    0.000    0.000    0.000
      4          0.000    0.000    0.000    0.000
      5          0.000    0.000    0.000    0.000
      6          0.000    0.000    0.000    0.000
      7          0.000    0.000    0.000    0.000
      8          0.000    0.000    0.000    0.000
      9          0.000    0.000    0.000    0.000
     10          0.000    0.000    0.000    0.000
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: /public4/home/sc56185/assimilation/TEST4/crtm_coeffs/ahi_himawari08.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION:  8.01
   N_CHANNELS=10
Read_ODPS_Binary(INFORMATION) : FILE: /public4/home/sc56185/assimilation/TEST4/crtm_coeffs/ahi_himawari08.TauCoeff.bin;
   ODPS RELEASE.VERSION:  2.03  N_LAYERS=100  N_COMPONENTS=5  N_ABSORBERS=3  N_CHANNELS=10  N_COEFFS=42700
CloudCoeff_ReadFile(Binary)(INFORMATION) : FILE: /public4/home/sc56185/assimilation/TEST4/crtm_coeffs/CloudCoeff.bin;
CloudCoeff RELEASE.VERSION:  3.04  N_FREQUENCIES(MW)=  31  N_FREQUENCIES(IR)=  61  N_RADII(MW)=10  N_RADII(IR)=10  N_TEMPERATURES= 5  N_DENSITIES= 3  N_LEGENDRE_TERMS=38  N_PHASE_ELEMENTS= 1
AerosolCoeff_ReadFile(Binary)(INFORMATION) : FILE: /public4/home/sc56185/assimilation/TEST4/crtm_coeffs/AerosolCoeff.bin;
AerosolCoeff RELEASE.VERSION:  3.04  N_WAVELENGTHS=  61  N_RADII= 36  N_TYPES= 8  N_RH= 36  N_LEGENDRE_TERMS=37  N_PHASE_ELEMENTS= 1
SEcategory_ReadFile(INFORMATION) : FILE: /public4/home/sc56185/assimilation/TEST4/crtm_coeffs/IGBP.IRland.EmisCoeff.bin;
SEcategory RELEASE.VERSION:  3.01
   CLASSIFICATION: IGBP,  N_FREQUENCIES=20  N_SURFACE_TYPES=20
IRwaterCoeff_ReadFile(INFORMATION) : FILE: /public4/home/sc56185/assimilation/TEST4/crtm_coeffs/Nalli.IRwater.EmisCoeff.bin;
IRwaterCoeff RELEASE.VERSION:  3.02  N_ANGLES= 76  N_FREQUENCIES= 2223  N_WIND_SPEEDS= 11
SEcategory_ReadFile(INFORMATION) : FILE: /public4/home/sc56185/assimilation/TEST4/crtm_coeffs/NPOESS.IRsnow.EmisCoeff.bin;
SEcategory RELEASE.VERSION:  3.01
   CLASSIFICATION: NPOESS,  N_FREQUENCIES=20  N_SURFACE_TYPES=2
SEcategory_ReadFile(INFORMATION) : FILE: /public4/home/sc56185/assimilation/TEST4/crtm_coeffs/NPOESS.IRice.EmisCoeff.bin;
SEcategory RELEASE.VERSION:  3.01
   CLASSIFICATION: NPOESS,  N_FREQUENCIES=20  N_SURFACE_TYPES=1
Reading AHI data from JAXA NETCDF4
AHI nchan:           10
nlongitude,nlatitude:         1500        1000
lonstart,latstart:            1           1
observation date:         2019           6           5          18           0
           0
latitude,longitude(pixel=1,scan=1):    60.00000       80.00000   
tbb(pixel=1,scan=1,chan=           1 ):    268.8100   
tbb(pixel=1,scan=1,chan=           2 ):    223.9700   
tbb(pixel=1,scan=1,chan=           3 ):    229.5700   
tbb(pixel=1,scan=1,chan=           4 ):    237.7100   
tbb(pixel=1,scan=1,chan=           5 ):    259.8200   
tbb(pixel=1,scan=1,chan=           6 ):    237.8600   
tbb(pixel=1,scan=1,chan=           7 ):    266.3000   
tbb(pixel=1,scan=1,chan=           8 ):    266.2400   
tbb(pixel=1,scan=1,chan=           9 ):    262.8800   
tbb(pixel=1,scan=1,chan=          10 ):    244.9200   
satellite zenith angle(pixel=1,scan=1):    84.45000   
satellite azimuth angle(pixel=1,scan=1):    115.9500   
solar zenith angle(pixel=1,scan=1):    97.07000   
solar azimuth angle(pixel=1,scan=1):   -8.780000   
In file: L1AHITBR got num_ahi_file    :  222841
In file: L1AHITBR got num_ahi_global  :       0
In file: L1AHITBR got num_ahi_local   :       0
AHI data counts:
In file:  222841
Global :       0
Local  :       0
Used   :       0
Thinned:       0
Observation summary
   ob time  1 num_ob_indexes= 29

----------------------------------------------------------
[3.0] Set up background errors (be) for cv_option:  3
----------------------------------------------------------

    DA_Setup_Background_Errors: 3DVAR dry control variables are:
    psi, chi_u, t_u and psfc
    DA_Setup_Background_Errors: 3DVAR humidity control variable is q/qsg

Number of vertical level for stats =         42
Number of latitude           nlath =         96

Number of vertical level for WRFVar=    34

VARBC: Reading VARBC.in file
VARBC: no matching for platform/satid/sensor   1  15   3
VARBC: no matching for platform/satid/sensor   1  16   3
VARBC: no matching for platform/satid/sensor   1  18   3
VARBC: no matching for platform/satid/sensor   1  19   3
VARBC: no matching for platform/satid/sensor  23   1  40
VARBC: no matching for platform/satid/sensor  23   2  40
VARBC: no matching for platform/satid/sensor  10   1   3
VARBC: no matching for platform/satid/sensor  10   2   3
VARBC: no matching for platform/satid/sensor   1  15   4
VARBC: no matching for platform/satid/sensor   1  16   4
VARBC: no matching for platform/satid/sensor   1  17   4
VARBC: no matching for platform/satid/sensor   1  18  15
VARBC: no matching for platform/satid/sensor   1  19  15
VARBC: no matching for platform/satid/sensor  23   1  41
VARBC: no matching for platform/satid/sensor  23   2  41
VARBC: no matching for platform/satid/sensor  10   1  15
VARBC: no matching for platform/satid/sensor  10   2  15
VARBC: no matching for platform/satid/sensor  17   0  19
VARBC: no matching for platform/satid/sensor   2  16  10
VARBC: no matching for platform/satid/sensor   9   2   3
VARBC: no matching for platform/satid/sensor   9   2  11
VARBC: no matching for platform/satid/sensor  19   0   0
VARBC: no matching for platform/satid/sensor   1  11   0
VARBC: no matching for platform/satid/sensor   1  12   0
VARBC: no matching for platform/satid/sensor   1  14   0
VARBC: no matching for platform/satid/sensor  19   0   1
VARBC: no matching for platform/satid/sensor   1  11   1
VARBC: no matching for platform/satid/sensor   1  12   1
VARBC: no matching for platform/satid/sensor   1  14   1
VARBC: no matching for platform/satid/sensor  10   1  16
VARBC: no matching for platform/satid/sensor  10   2  16
VARBC: no matching for platform/satid/sensor  12   2  21
VARBC: no matching for platform/satid/sensor  12   3  21
VARBC: no matching for platform/satid/sensor  29   1  63
Calculate innovation vector(iv)

VARBC: Applying bias correction for himawari-8-ahi
VARBC: Estimate Hessian for preconditioning
Minimize cost function using CG method


Diagnostics
   Final cost function J       =                           0.00

   Total number of obs.        =        0
   Final value of J            =                        0.00000
   Final value of Jo           =                        0.00000
   Final value of Jd           =         0.00000
   Final value of Jb           =         0.00000
   Final value of Jc           =         0.00000
   Final value of Je           =         0.00000
   Final value of Jp           =         0.00000
   Final value of Jl           =         0.00000
   Final value of Jm           =         0.00000
   Final value of Jt           =         0.00000
   Jb factor used              =         1.00000
   Je factor used              =         1.00000
   VarBC factor used           =         1.00000

   Total number of radiances    =        0
   Cost function for radiances  =         0.00000

Writing radiance OMA ascii file

VARBC: Updating bias parameters
VARBC: Writing information in VARBC.out file
convert T to THM when use_theta_m =  1

*** WRF-Var completed successfully ***


密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 6 天前 | 显示全部楼层
本帖最后由 任性小王子 于 2024-4-12 19:37 编辑

请问楼主解决了吗?但是有个问题Reading extra data for AHI
Channel       BTLim      ca1     ca2     clearSkyBias
      1          0.000    0.000    0.000    0.000
      2          0.000    0.000    0.000    0.000
      3          0.000    0.000    0.000    0.000
      4          0.000    0.000    0.000    0.000
      5          0.000    0.000    0.000    0.000
      6          0.000    0.000    0.000    0.000
      7          0.000    0.000    0.000    0.000
      8          0.000    0.000    0.000    0.000
      9          0.000    0.000    0.000    0.000
     10          0.000    0.000    0.000    0.000
你的这些数据怎么都是0呐。但是我运行的结果 Reading extra data for AHI
Channel       BTLim      ca1     ca2     clearSkyBias
      7        273.300    0.500   14.500   -1.023
      8        238.000    1.000   25.000    0.000
      9        230.350    0.500   14.500   -1.221
     10        220.720    0.500   14.500   -1.621
     11        220.720    0.500   14.500   -1.621
     12        220.720    0.500   14.500   -1.621
     13        220.720    0.500   14.500   -1.621
     14        220.720    0.500   14.500   -1.621
     15        220.720    0.500   14.500   -1.621
     16        220.720    0.500   14.500   -1.621

这些都是有值的呐?是namelist设置的问题吗?可以贴一下你的namelist吗?
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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