Moist idealized baroclinic wave test case produces NaNs

J

jagu0521

New member
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!
 
Could you attach your full namelist.init_atmosphere and namelist.atmosphere files? Which version of MPAS-A are you running?
 
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.
 

Attachments

Thanks for looking over the namelists! Unfortunately, I am still getting the same issue with the above namelist.atmosphere. The following is an excerpt of the log.atmosphere.0000.out that was produced:

Begin timestep 0000-01-01_00:00:00
--- time to update background surface albedo, greeness fraction.
--- time to calculate additional physics_diagnostics = F
coupled RK3 dynamics-transport integration

global min, max w -0.161751E-02 0.377960E-02
global min, max u -5.45398 5.45399
global min, max scalar 1 0.00000 0.278022E-01
global min, max scalar 2 0.00000 0.00000
global min, max scalar 3 0.00000 0.00000
Timing for integration step: 0.558665 s

Begin timestep 0000-01-01_00:07:30
--- time to calculate additional physics_diagnostics = F
coupled RK3 dynamics-transport integration

global min, max w -0.174202E-02 0.183917E-03
global min, max u -5.45385 5.45386
global min, max scalar 1 0.00000 0.278022E-01
global min, max scalar 2 0.00000 0.00000
global min, max scalar 3 0.00000 0.00000
Timing for integration step: 0.556557 s

Begin timestep 0000-01-01_00:15:00
--- time to calculate additional physics_diagnostics = F
coupled RK3 dynamics-transport integration

global min, max w NaN NaN
global min, max u NaN NaN
global min, max scalar 1 0.00000 0.278022E-01
global min, max scalar 2 0.00000 0.00000
global min, max scalar 3 0.00000 0.00000
Timing for integration step: 0.556421 s

Begin timestep 0000-01-01_00:22:30
--- time to calculate additional physics_diagnostics = F
coupled RK3 dynamics-transport integration

global min, max w NaN NaN
global min, max u NaN NaN
global min, max scalar 1 0.00000 0.00000
global min, max scalar 2 0.00000 0.00000
global min, max scalar 3 0.00000 0.00000
Timing for integration step: 0.558373 s

I have run the dry baroclinic wave test case with no issue before and it is only when changing moisture to true in the subroutine, that I am getting NaNs for both u and w. Is there something else in the subroutine or namelists that I should be modifying?
 
In addition to the changes Michael suggested, please also modify your options as listed below, then try again.

Code: 
config_dt = 450.0
config_scalar_advection  = false
config_positive_definite = false
config_coef_3rd_order = 1.0
config_monotonic = false

Let us know whether your case works or not. Thanks.
 
Thank you for the update. I am sorry to know that this case still doesn't work for you. I will repeat this case and get back to you once it is done.
 
You must log in or register to reply here.
Back
Top