Hello everyone,
When I was working on the cycle mode of WRFDA, I conducted data assimilation at 0:00, 6:00, 12:00, and 18:00 respectively. Moving to 6:00, I updated the lower boundary of WRFOUT according to the manual, then assimilated DA, updated the lateral boundary, and ran the data for the next 6 hours; But there are some issues with my current simulation results. I feel like the mode restarted when I was looping, instead of continuing the previous hour;
For example, the following is the data value of a certain site, which seems to restart the mode at times 0, 6, 12, and 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
So, I have two unclear points.
Firstly, when I first ran from 0:00 on the 1st to 6:00 on the 1st, the data at 6:00 was not actually calculated. Therefore, should I set the end time to 7:00 on the 1st, and then use the integrated 6:00 data to update the boundaries and assimilate the data, and then go through the loop?
Secondly, while reading the official website manual, at the second moment, which is 6 o'clock on the 1st of my mode, I need to update the lower boundary. The parame.in namelist on the official website is shown below The original link for fg was wrfinput-d01 at 0:00, but I changed it to wrfinput_d01 at 6:00 in the second moment. However, wrf_input has already provided the input file for wrfinput-d01. Why do I need to link another file that is the same? Or, I don't understand what file should be linked to this da_file
da_file = './fg'
wrf_input = './wrfinput_d01'
update_lateral_bdy = .false.
update_low_bdy = .true.
iswater = 16
When I was working on the cycle mode of WRFDA, I conducted data assimilation at 0:00, 6:00, 12:00, and 18:00 respectively. Moving to 6:00, I updated the lower boundary of WRFOUT according to the manual, then assimilated DA, updated the lateral boundary, and ran the data for the next 6 hours; But there are some issues with my current simulation results. I feel like the mode restarted when I was looping, instead of continuing the previous hour;
For example, the following is the data value of a certain site, which seems to restart the mode at times 0, 6, 12, and 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
So, I have two unclear points.
Firstly, when I first ran from 0:00 on the 1st to 6:00 on the 1st, the data at 6:00 was not actually calculated. Therefore, should I set the end time to 7:00 on the 1st, and then use the integrated 6:00 data to update the boundaries and assimilate the data, and then go through the loop?
Secondly, while reading the official website manual, at the second moment, which is 6 o'clock on the 1st of my mode, I need to update the lower boundary. The parame.in namelist on the official website is shown below The original link for fg was wrfinput-d01 at 0:00, but I changed it to wrfinput_d01 at 6:00 in the second moment. However, wrf_input has already provided the input file for wrfinput-d01. Why do I need to link another file that is the same? Or, I don't understand what file should be linked to this da_file
da_file = './fg'
wrf_input = './wrfinput_d01'
update_lateral_bdy = .false.
update_low_bdy = .true.
iswater = 16