Hi,
I noticed that in mpas_atm_time_integration.F where the Rayleigh damping and horizontal diffusion were applied to the relaxation zones, the calculations of Rayleigh damping coefficient (F1) and horizontal diffusion coefficient (F2) were different from the descriptions (Equations 1 and 2) in Skamarock et al (2018). In the code gamma1=1/(50*dt) and gamma2=1/(10*dt), while in the paper (also WRF model descriptions) gamma1=1/(10*dt) and gamma2=1/(50*dt). Is this a typo in the code?
Thanks,
Yuan
I noticed that in mpas_atm_time_integration.F where the Rayleigh damping and horizontal diffusion were applied to the relaxation zones, the calculations of Rayleigh damping coefficient (F1) and horizontal diffusion coefficient (F2) were different from the descriptions (Equations 1 and 2) in Skamarock et al (2018). In the code gamma1=1/(50*dt) and gamma2=1/(10*dt), while in the paper (also WRF model descriptions) gamma1=1/(10*dt) and gamma2=1/(50*dt). Is this a typo in the code?
Thanks,
Yuan