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

whuwjr · 发表于 2024-2-23 10:58:00

本帖最后由 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 ***

任性小王子 · 发表于 2024-4-12 19:36:08

本帖最后由任性小王子于 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吗？

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