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!
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!