爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 1440|回复: 8

WRFDA+3dvar 同化 AHI data

[复制链接]

新浪微博达人勋

发表于 2023-11-7 10:40:01 | 显示全部楼层 |阅读模式

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

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

x
大神们可以看一下这是什么问题吗?我用3dvar同化AHI数据,rsl.out.0000出现下面问题
taskid: 0 hostname: changmh
Namelist logging not found in namelist.input. Using registry defaults for variables in logging.


*** VARIATIONAL ANALYSIS ***
    WRFDA V4.5.1                                                                    


Ntasks in X            1 , ntasks in Y            2
*************************************
Parent domain
ids,ide,jds,jde            1         150           1         130
ims,ime,jms,jme           -4         155          -4          72
ips,ipe,jps,jpe            1         150           1          65
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             533435800  bytes allocated
   Input data is acceptable to use: fg
hybrid_opt =  1
  
use_theta_m =  1
  
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE    149 JS      1 JE     65
WRF NUMBER OF TILES =   1
Domain mapping info:
map_proj =     1
cen_lat   =    0.330000E+02
cen_lon   =   -0.790000E+02
truelat1  =    0.300000E+02
truelat2  =    0.600000E+02
start_lat =    0.237806E+02
start_lon =   -0.895973E+02
pole_lat  =    0.900000E+02
dsm       =    0.150000E+02
  
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 1 1 1 0 0
  
Using PREPBUFR format observation input
  
PREPBUFR ob with quality marker <= 3 will be retained.
  
BUFR file date is: 2022010118
  
da_read_obs_bufr: num_report, num_outside_all, num_outside_time, num_thinned:   726057  675178       0       0
  
Using NCEP BUFR radiance 1b input
  
        3000
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
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/ahi_himawari08.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION:  8.02
   N_CHANNELS=10
ACCoeff_ReadFile(Binary)(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/amsua_n15.SpcCoeff.bin;
ACCoeff RELEASE.VERSION:  1.04
   N_FOVS=30  N_CHANNELS=15
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/amsua_n15.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION:  8.01
   N_CHANNELS=15
Read_ODPS_Binary(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/ahi_himawari08.TauCoeff.bin;
   ODPS RELEASE.VERSION:  2.01  N_LAYERS=100  N_COMPONENTS=5  N_ABSORBERS=3  N_CHANNELS=10  N_COEFFS=42700
Read_ODPS_Binary(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/amsua_n15.TauCoeff.bin;
   ODPS RELEASE.VERSION:  2.01  N_LAYERS=100  N_COMPONENTS=2  N_ABSORBERS=1  N_CHANNELS=15  N_COEFFS=21600
CloudCoeff_ReadFile(Binary)(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/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: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/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: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/USGS.IRland.EmisCoeff.bin;
SEcategory RELEASE.VERSION:  3.01
   CLASSIFICATION: USGS,  N_FREQUENCIES=20  N_SURFACE_TYPES=27
IRwaterCoeff_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/Nalli.IRwater.EmisCoeff.bin;
IRwaterCoeff RELEASE.VERSION:  3.02  N_ANGLES= 76  N_FREQUENCIES= 2223  N_WIND_SPEEDS= 11
SEcategory_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/NPOESS.IRsnow.EmisCoeff.bin;
SEcategory RELEASE.VERSION:  3.01
   CLASSIFICATION: NPOESS,  N_FREQUENCIES=20  N_SURFACE_TYPES=2
SEcategory_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/NPOESS.IRice.EmisCoeff.bin;
SEcategory RELEASE.VERSION:  3.01
   CLASSIFICATION: NPOESS,  N_FREQUENCIES=20  N_SURFACE_TYPES=1
FitCoeff_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 2
FitCoeff_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 6
FitCoeff_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 3, 7, 2
FitCoeff_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 10, 4, 3
FitCoeff_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 8
FitCoeff_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 3, 6, 2
MWwaterLUT_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/FASTEM5.MWwater.EmisCoeff.bin;
MWwaterLUT RELEASE.VERSION:  1.01
   N_ANGLES=15  N_FREQUENCIES=18  N_TEMPERATURES=12  N_WIND_SPEEDS=16
MWwaterCoeff_ReadFile(INFORMATION) : FILE: /mnt/e/Models/Build_WRFDA/crtm_coeffs_2.3.0/FASTEM5.MWwater.EmisCoeff.bin;
MWwaterCoeff RELEASE.VERSION: 1.5
Reading AHI data from JAXA NETCDF4
AHI nchan:           10
nlongitude,nlatitude:         2401        2401
lonstart,latstart:            1           1
observation date:         2022           1           1           0           0           0
latitude,longitude(pixel=1,scan=1):    60.0000000       80.0000000   
tbb(pixel=1,scan=1,chan=           1 ):    236.839996   
tbb(pixel=1,scan=1,chan=           2 ):    212.940002   
tbb(pixel=1,scan=1,chan=           3 ):    221.839996   
tbb(pixel=1,scan=1,chan=           4 ):    224.970001   
tbb(pixel=1,scan=1,chan=           5 ):    232.850006   
tbb(pixel=1,scan=1,chan=           6 ):    216.369995   
tbb(pixel=1,scan=1,chan=           7 ):    236.779999   
tbb(pixel=1,scan=1,chan=           8 ):    234.979996   
tbb(pixel=1,scan=1,chan=           9 ):    235.139999   
tbb(pixel=1,scan=1,chan=          10 ):    224.350006   
satellite zenith angle(pixel=1,scan=1):    84.4599991   
satellite azimuth angle(pixel=1,scan=1):    115.959999   
solar zenith angle(pixel=1,scan=1):    115.059998   
solar azimuth angle(pixel=1,scan=1):    93.0599976   
In file: L1AHITBR got num_ahi_file    : 4315723
In file: L1AHITBR got num_ahi_global  :       0
In file: L1AHITBR got num_ahi_local   :       0
AHI data counts:
In file: 4315723
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=    44
  
Calculate innovation vector(iv)
  
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
  
convert T to THM when use_theta_m =  1
  
*** WRF-Var completed successfully ***



我的namelist文件添加到附件了:
希望大神指点一下。

namelist.input

2.24 KB, 下载次数: 12, 下载积分: 金钱 -5

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

新浪微博达人勋

发表于 2024-2-23 11:08:29 | 显示全部楼层
解决了一部分,已经能读到数据了,需要修改ahi_info文件。如果你直接下载的是葵花的nc文件,要把ahi_info中第二行的数字修改一下,这个数字表示你输入的数据格式。然后我觉得还应该根据同化的区域,修改第六行的area information,这个我还没太搞清楚,在尝试修改中。
1.png
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2024-2-21 19:27:48 | 显示全部楼层
嗨你好,我也遇到了这个问题,请问你现在解决了这个问题吗(求助
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-2-21 20:59:55 | 显示全部楼层
Anya 发表于 2024-2-21 19:27
嗨你好,我也遇到了这个问题,请问你现在解决了这个问题吗(求助

现在没解决,我只能再看源代码
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2024-2-22 09:30:52 | 显示全部楼层
任性小王子 发表于 2024-2-21 20:59
现在没解决,我只能再看源代码

好的好的,谢谢啦
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2024-2-22 11:57:44 | 显示全部楼层
Anya 发表于 2024-2-22 09:30
好的好的,谢谢啦

请问您解决了吗 我也遇到同样的问题了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2024-2-22 12:06:20 | 显示全部楼层
遇到同样的问题,帮顶
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2024-2-23 11:50:29 | 显示全部楼层
Anya 发表于 2024-2-23 11:08
解决了一部分,已经能读到数据了,需要修改ahi_info文件。如果你直接下载的是葵花的nc文件,要把ahi_info中 ...

我的ahi_info 第二行是3,是没错的,第六行你可以改为2401试试,我现在还是同样的问题
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-2-23 17:13:55 | 显示全部楼层
whuwjr 发表于 2024-2-23 11:50
我的ahi_info 第二行是3,是没错的,第六行你可以改为2401试试,我现在还是同样的问题

源代码上再看一下,v4.5.2 和v4.4版本同化模块有差异
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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