Moist idealized baroclinic wave test case produces NaNs

[jagu0521]

Hi everyone!

I am using the JW idealized baroclinic wave test case and I am having trouble getting the moist case to run. I have edited the subroutine in mpas_init_atm_cases.F by uncommenting logical, parameter :: moisture = .true. and commenting out logical, parameter :: moisture = .false. without changing anything else in the subroutine. init_atmosphere runs with no errors and x1.40962.init.nc shows that moisture has been initialized properly. In namelist.atmosphere, I have decreased the config_dt to 120.0, and set config_scalar_advection, config_positive_definite, and config_monotonic all to true. When I run atmosphere.model, I get NaNs for u and w almost immediately:
(excerpt from log.atmosphere.0000.out)

min/max of meshScalingDel2 = 1.00000 1.00000
min/max of meshScalingDel4 = 1.00000 1.00000
Timing for diagnostic computation: 0.211806 s
Timing for stream output: 0.596264 s

Begin timestep 0000-01-01_00:00:00
coupled RK3 dynamics-transport integration

global min, max w -0.735288E-02 0.451069E-04
global min, max u -5.45494 5.45494
Timing for integration step: 0.439597 s

Begin timestep 0000-01-01_00:00:05
coupled RK3 dynamics-transport integration

global min, max w -0.951185E-02 0.681508E-04
global min, max u -5.45492 5.45492
Timing for integration step: 0.438246 s

Begin timestep 0000-01-01_00:00:10
coupled RK3 dynamics-transport integration

global min, max w NaN NaN
global min, max u NaN NaN
Timing for integration step: 0.438428 s

Why is this happening? Should I be changing or adding other flags in namelist.atmosphere?

Any advice would be appreciated, thank you!

 

[mgduda]

Could you attach your full namelist.init_atmosphere and namelist.atmosphere files? Which version of MPAS-A are you running?

 

[jagu0521]

I am using MPAS-Atmosphere Version 8.3.1 and both of the namelist files that I am using are attached.

 

[mgduda]

Thanks very much for attaching your namelists. It may be that the values for config_v_theta_eddy_visc2 and config_v_mom_eddy_visc2 that you've specified are too large. Could you try running with the attached namelist? Also, if you do plan to eventually use the Kessler microphysics, you will likely need to patch in the fix from MPAS-Dev/MPAS-Model PR #1452, though other microphysics schemes should work without any issue.