WRFDA 3dVAR循环模式模拟值出错

lirui12139

大家好，我在做WRFDA的cycle模式的时候，分别在0点、6点、12点和18点进行资料同化，转到6点先根据手册更新WRFOUT的下边界，然后做DA的同化，再更新测边界，再运行下一个6小时的数据；但是我现在的模拟结果有些问题，我感觉在我循环的时刻，模式重新启动了，而不是接着上一个小时继续转；例如下面是个某个站点的数据值，在0、6、12、18时刻都像是模式重启了；time,pm25(ug/m3),pm10(ug/m3),o3(ug/m3),no2(ug/m3),co(mg/m3),so2(ug/m3)
2024-04-22 00,0.20244463,0.20244463,64.30935,0.033279184,0.1,0.2857143
2024-04-22 01,6.6517353,8.997257,64.78694,19.273441,0.1760601,1.0644994
2024-04-22 02,8.402133,11.223334,73.285675,23.432602,0.198681,1.3076808
2024-04-22 03,7.945503,10.757943,84.93985,22.55281,0.20321059,1.471113
2024-04-22 04,7.1627665,9.568329,93.98237,17.53452,0.19261995,1.6030457
2024-04-22 05,6.4791307,8.586087,100.89378,13.56076,0.18418005,1.7519838
2024-04-22 06,0.19254985,0.19254985,64.30935,0.033279184,0.1,0.2857143
2024-04-22 07,2.9175656,4.042787,80.17402,8.309897,0.12740694,0.6784844
2024-04-22 08,4.810759,6.532378,84.790245,15.817197,0.14632753,0.9595211
2024-04-22 09,7.012762,9.677109,82.01494,28.049314,0.17808314,1.4064645
2024-04-22 10,20.627956,26.851313,73.77171,87.235825,0.335143,3.28203
2024-04-22 11,11.917521,18.29122,72.67229,91.260666,0.33285907,2.8657918
2024-04-22 12,0.19944471,0.19944471,64.30935,0.033279184,0.1,0.2857143
2024-04-22 13,10.289372,13.328677,54.274044,56.029385,0.2509546,1.7263021
2024-04-22 14,14.368815,18.041348,50.667812,68.769455,0.28865427,1.9972563
2024-04-22 15,24.886362,28.62362,54.67097,68.70016,0.26424688,2.857143e-13
2024-04-22 16,27.447815,30.9804,51.299793,65.54616,0.28720856,2.857143e-13
2024-04-22 17,29.744707,33.07385,52.44574,62.536034,0.30275938,0.08715248
2024-04-22 18,0.20107059,0.20107059,64.30935,0.033279184,0.1,0.2857143
2024-04-22 19,7.778324,9.880611,55.31599,24.932041,0.16347736,1.14791
2024-04-22 20,8.90294,11.103461,53.169144,26.426458,0.17267276,1.2598591
2024-04-22 21,13.58535,16.433472,49.35831,37.019596,0.20958096,1.7173481
2024-04-22 22,15.20961,18.326368,48.45653,40.353115,0.22007786,1.6794975
2024-04-22 23,4.3766694,5.701521,59.05896,13.142976,0.14544411,0.56667775
2024-04-23 00,3.2240758,4.155023,65.01463,7.0840573,0.1262662,0.38862368
我在阅读官网手册时，在第二个时刻，也就是我模式的1日的6点，需要更新下边界，官网中的parame.in如下所示，那他的./fg 原本链接的是0点的wrfinput_d01，我在做的时候在第二个时刻更换为6点的wrfinput_d01，但[size=13.3333px]wrf_input已经给了[size=13.3333px]wrfinput_d01输入文件了，为什么还需要再链接一个相同的呢，或者说我没理解对这个da_file应该链接的文件是什么

da_file            = './fg'

wrf_input          = './wrfinput_d01'

update_lateral_bdy = .false.

update_low_bdy     = .true.

iswater            = 16

Nozigaka46

更新下垫面的时候，da_file 指的应该是第一猜测场first guess，即对应时刻的wrfout

lirui12139

是的，后来我明白了，更改了wrfout作为猜测场，但是我在运行完DA之后去做第二个循环，也就是6点时刻的循环时，将wrfvar_output作为wrfinput重启wrf模式时，他好像还是没有收到同化后的6点时刻的数据，如果您方便的话，可以看我的这个帖子https://bbs.06climate.com/forum.php?mod=viewthread&tid=109933，里面有一些我输出的nc图片内容

最后如果您方便的话，我想加您的qq和您交流以下这个问题，q:1433048582 万分感谢！

任性小王子

lirui12139 发表于 2024-5-10 13:34
是的，后来我明白了，更改了wrfout作为猜测场，但是我在运行完DA之后去做第二个循环，也就是6点时刻的循环 ...

不错的分享，学习啦

一两

lirui12139 发表于 2024-5-10 13:34
是的，后来我明白了，更改了wrfout作为猜测场，但是我在运行完DA之后去做第二个循环，也就是6点时刻的循环 ...

请问这个问题解决了吗

lirui12139

没有解决，因为PM2.5是需要用mosaic机制去转，然后同化气溶胶变量，而不是直接同化颗粒物，比较复杂 被我搁置了