WRFDA4.0 Bias Correction

prj

Radiance Data资料同化，使用的版本是WRFDA4.0，偏差订正没有效果，完全没有起作用，3.8版本也使用过，依然是不起作用
OMB with BC 和OMB no BC完全一样
然后迭代开始的时候，又 gradient=NaN，然后报错

1. &wrfvar1
   
2. var4d=true,
   
3. var4d_lbc=false,
   
4. var4d_bin=3600,
   
5. var4d_bin_rain=21600,
   
6. print_detail_obs=.true.,
   
7. print_detail_grad=.true.,
   
8. /
   
9. &wrfvar2
   
10. /
    
11. &wrfvar3
    
12. num_fgat_time=7,
    
13. /
    
14. &wrfvar4
    
15. use_amsuaobs=.true.
    
16. /
    
17. &wrfvar5
    
18. check_max_iv=true,
    
19. /
    
20. &wrfvar6
    
21. max_ext_its=1,
    
22. ntmax=1,
    
23. orthonorm_gradient=true,
    
24. /
    
25. &wrfvar7
    
26. cv_options=3
    
27. as1(3) = 0.1,0.4,0.6
    
28. as2(3) = 0.1,0.4,0.6
    
29. as3(3) = 0.1,0.4,0.6
    
30. as4(3) = 0.1,0.4,0.6
    
31. as5(3) = 0.1,0.4,0.6
    
32. /
    
33. &wrfvar8
    
34. /
    
35. &wrfvar9
    
36. /
    
37. &wrfvar10
    
38. test_transforms=false,
    
39. test_gradient=false,
    
40. /
    
41. &wrfvar11
    
42. /
    
43. &wrfvar12
    
44. /
    
45. &wrfvar13
    
46. /
    
47. &wrfvar14
    
48. rtminit_nsensor=1,
    
49. rtminit_platform=1
    
50. rtminit_satid=19
    
51. rtminit_sensor=3
    
52. thinning=true,
    
53. thinning_mesh=60.0,
    
54. qc_rad=true,
    
55. write_iv_rad_ascii=true,
    
56. write_oa_rad_ascii=true,
    
57. rtm_option=2,
    
58. only_sea_rad=false,
    
59. crtm_cloud=true,
    
60. use_varbc=true,
    
61. use_crtm_kmatrix =true,
    
62. freeze_varbc=false,
    
63. varbc_factor=1,
    
64. /
    
65. &wrfvar15
    
66. /
    
67. &wrfvar16
    
68. /
    
69. &wrfvar17
    
70. /
    
71. &wrfvar18
    
72. analysis_date=2015-08-03_06:00:00.0000,
    
73. /
    
74. &wrfvar19
    
75. /
    
76. &wrfvar20
    
77. /
    
78. &wrfvar21
    
79. time_window_min=2015-08-03_06:00:00.0000,
    
80. /
    
81. &wrfvar22
    
82. time_window_max=2015-08-03_12:00:00.0000,
    
83. /
    
84. &wrfvar23
    
85. /
    
86. &time_control
    
87. run_hours=6,
    
88. start_year=2015,
    
89. start_month=08,
    
90. start_day=03,
    
91. start_hour=06,
    
92. end_year=2015,
    
93. end_month=08,
    
94. end_day=03,
    
95. end_hour=12,
    
96. interval_seconds=21600,
    
97. /
    
98. &fdda
    
99. /
    
100. &domains
     
101. time_step=180,
     
102. e_we=121,
     
103. e_sn=121,
     
104. e_vert=42,
     
105. num_metgrid_levels=32,
     
106. p_top_requested=5000,
     
107. interp_type=2,
     
108. t_extrap_type=2,
     
109. dx=30000,
     
110. dy=30000,
     
111. /
     
112. &dfi_control
     
113. /
     
114. &tc
     
115. /
     
116. &physics
     
117. mp_physics=2,
     
118. ra_lw_physics=1,
     
119. ra_sw_physics=1,
     
120. radt=30,
     
121. sf_sfclay_physics=1,
     
122. sf_surface_physics=2,
     
123. bl_pbl_physics=1,
     
124. cu_physics=1,
     
125. cudt=0,
     
126. num_soil_layers=4,
     
127. mp_zero_out=2,
     
128. co2tf=0,
     
129. /
     
130. &scm
     
131. /
     
132. &dynamics
     
133. hybrid_opt=0,
     
134. etac=0.0,
     
135. use_theta_m=0,
     
136. w_damping=1,
     
137. diff_opt=1,
     
138. km_opt=4,
     
139. dampcoef=0.2,
     
140. time_step_sound=4,
     
141. base_temp=290.0,
     
142. use_baseparam_fr_nml=true,
     
143. /
     
144. &bdy_control
     
145. spec_bdy_width                      = 5,
     
146. spec_zone                           = 1,
     
147. relax_zone                          = 4,
     
148. specified                           = .true.,
     
149. nested                              = .false.,
     
150. /
     
151. &grib2
     
152. /
     
153. &fire
     
154. /
     
155. &namelist_quilt
     
156. /
     
157. &perturbation
     
158. trajectory_io=true,
     
159. enable_identity=false,
     
160. jcdfi_use=false,
     
161. jcdfi_diag=1,
     
162. jcdfi_penalty=1000.0,
     
163. /

复制代码

ljchen1989

同问，我也偶尔遇到gradient=NaN
然后就投降了

prj

ljchen1989 发表于 2020-7-22 16:11
同问，我也偶尔遇到gradient=NaN
然后就投降了

你做的也是四维变分的辐射资料同化吗？偏差订正有效果吗

ljchen1989

prj 发表于 2020-7-24 08:38
你做的也是四维变分的辐射资料同化吗？偏差订正有效果吗

没有nan的都有效果
nan的跑不通

prj

ljchen1989 发表于 2020-7-24 10:36
没有nan的都有效果
nan的跑不通

你好，我用3.8.1版本的就算跑通也没有效果，namelist的设置和4.0的几乎一样，可以的话能看看你的namelist设置吗

任性小王子

prj 发表于 2020-7-24 15:03
你好，我用3.8.1版本的就算跑通也没有效果，namelist的设置和4.0的几乎一样，可以的话能看看你的namelist ...

在上面已经贴出来了，楼主找到原因了吗？可以分享一下经验吗？