running mp_thompson using the hail aware option

[ldfowler58]

I implemented module_mp_thompson.F from WRF version 4.6.1 in MPAS and tried to run mp_thompson with the option is_hail_aware = .false. (In WRF, mp_physics = 8). My run crashes on the very first timestep at line 2537. When using the option mp_physics = 8, dimNRHG is equal to NRHG1 (=1). This implies that the third dimension of the array tcg_racg is also equal to 1.

line 2537 is prr_rcg(k) = tcg_racg(idx_g1,idx_g,idx_bg(k),idx_r1,idx_r) where idx_bg(k) = idx_bg1 which is equal to 5 and not 1. This is the reason why I get the following error:
forrtl: severe (408): fort: (2): Subscript #3 of the array TCG_RACG has value 5 which is greater than the upper bound of 1

when running the Thompson cloud microphysics scheme in MPAS. Has someone made a similar post in the past.

Thanks,
Laura

 

[Ming Chen]

Laura,

I didn't see any reports regarding this issue in WRF. It seems that this is specific for MPAS.

Note that in WRF, is_hail_aware is hard wired to be .false. and it is .true. under the condition:

if (PRESENT(ng)) then
is_hail_aware = .TRUE.
dimNRHG = NRHG
end if

In the official release of MPASV8.2, I couldn't fine this option, --- please let me know if I miss something.