gewitterblitz
New member
Hi,
I am trying to implement a 'theta-tendency nudging' method to force convective updrafts in my idealized WRF simulation. I want the tendency term to be applied throughout the simulation duration or maybe use a start and end time later on. As per my current understanding, doing so through 'module_initialize_ideal.F' would not be possible since I want the time-varying forcing. Therefore, I decided to add the relevant code directly in the 'solve_em.F' file, at line 1046 (WRF V4.6), right before Step 3, the acoustic steps.
I tried implementing the t_tend nudging using the following pseudo-code:
`t_tend(i,k,j) = t_tend(i,k,j) + some 3d gaussian function`, where i, k, and j range from (grid%i_start(ij),grid%i_end(ij)), (grid%j_start(ij),grid%j_end(ij)), and (1,k_end), respectively.
I added these lines of code after calls to the `first_rk_step_part1`, `first_rk_step_part2`, `rk_tendency`, and `rk_addtend_dry` subroutines. Since all these subroutines have t_tend defined as an INOUT variable, I thought I could simply modify t_tend directly after these subroutines are called within the solve_em subroutine.
I set `delt=0` for the quarter_ss ideal_case (ideal case = 2) and compile the code without any error. However, I do not notice any forcing being applied applied over the domain in the simulation output.
Could someone please help me identify the issues with my current implementation? Can the `t_tend` term not be modified within solve_em.F directly? Alternative suggestions to implement these changes are most welcome.
Thank you for your time and consideration!
I am trying to implement a 'theta-tendency nudging' method to force convective updrafts in my idealized WRF simulation. I want the tendency term to be applied throughout the simulation duration or maybe use a start and end time later on. As per my current understanding, doing so through 'module_initialize_ideal.F' would not be possible since I want the time-varying forcing. Therefore, I decided to add the relevant code directly in the 'solve_em.F' file, at line 1046 (WRF V4.6), right before Step 3, the acoustic steps.
I tried implementing the t_tend nudging using the following pseudo-code:
`t_tend(i,k,j) = t_tend(i,k,j) + some 3d gaussian function`, where i, k, and j range from (grid%i_start(ij),grid%i_end(ij)), (grid%j_start(ij),grid%j_end(ij)), and (1,k_end), respectively.
I added these lines of code after calls to the `first_rk_step_part1`, `first_rk_step_part2`, `rk_tendency`, and `rk_addtend_dry` subroutines. Since all these subroutines have t_tend defined as an INOUT variable, I thought I could simply modify t_tend directly after these subroutines are called within the solve_em subroutine.
I set `delt=0` for the quarter_ss ideal_case (ideal case = 2) and compile the code without any error. However, I do not notice any forcing being applied applied over the domain in the simulation output.
Could someone please help me identify the issues with my current implementation? Can the `t_tend` term not be modified within solve_em.F directly? Alternative suggestions to implement these changes are most welcome.
Thank you for your time and consideration!