Bug on WRF v3.9.1.1 initialization for real data cases

Any issues with the actual running of the WRF.

Bug on WRF v3.9.1.1 initialization for real data cases

Postby duoying » Wed Jul 11, 2018 10:33 am

For WRF v3.9.1.1, there’s a bug in dyn_em/module_initialize_real.F. After the Do loop at line 1681, it’s necessary to add the following line to avoid the dry pressure grid%pd_gc at middle grids from decreasing continually.

IF ( skip_middle_points_t ( ids , ide , jds , jde , i , j , em_width , hold_ups ) ) CYCLE

1680 DO j = jts, MIN(jte,jde-1)
1681 DO i = its, MIN(ite,ide-1)
1682 IF ( skip_middle_points_t ( ids , ide , jds , jde , i , j , em_width , hold_ups ) ) CYCLE ! add this line
1683 grid%psfc_gc(i,j) = grid%pd_gc(i,1,j)
1684 grid%pd_gc(i,1,j) = grid%pslv_gc(i,j) - ( grid%p_gc(i,1,j) - grid%pd_gc(i,1,j) )
1685 grid%ht_gc(i,j) = grid%ght_gc(i,1,j)
1686 grid%ght_gc(i,1,j) = 0.
1687 END DO
1688 END DO

The dry pressure of middling points adjacent to boundary are used in vertical interpolation for wind. In our case of 1-year simulation, the continually decreasing surface dry pressure finally crashes the initialization at May.

Another way to fix it, comment out Line 1716:

1712 ! Put things back to normal.
1713
1714 DO j = jts, MIN(jte,jde-1)
1715 DO i = its, MIN(ite,ide-1)
1716 ! IF ( skip_middle_points_t ( ids , ide , jds , jde , i , j , em_width , hold_ups ) ) CYCLE ! disable this line
1717 grid%pd_gc(i,1,j) = grid%psfc_gc(i,j)
1718 grid%ght_gc(i,1,j) = grid%ht_gc(i,j)
1719 END DO
1720 END DO

Not sure if this bug exists in recently released WRF v4.0.
duoying
 
Posts: 1
Joined: Wed Jul 11, 2018 9:55 am

Return to Runtime Problems

Who is online

Users browsing this forum: Google [Bot] and 5 guests

cron