Running 60km-20km mesh

Wei Huang

unread,

Aug 30, 2018, 12:12:22 PM8/30/18

to MPAS-Atmosphere Help

With Michael's help, I have been able to generate the static.nc file, and then generated the x4.535554.init.nc file.

But when I ran it, the model blow up quickly:

[whuang@cy001 run.535554]$ grep global log.atmosphere.0000.out | more

global min, max w -29.0344649209343 34.8392375112011

global min, max u -3312.45592704106 3456.19273924604

global min, max w -3107788.19600904 9762270.87691499

global min, max u -1904324993.19938 1304398041.48502

global min, max w -6150.06917814954 7744.30310326945

global min, max u -15047645.7754188 16731533.2048295

global min, max w -0.817013509721307 0.773054297597560

global min, max u -9605.25682877766 6311.81543508287

global min, max w -0.306602117267682 0.450222308125158

global min, max u -92.7939317368138 92.7571427539108

global min, max w -0.300626324228148 0.487439273389991

global min, max u -90.3329891640772 90.5421181666273

global min, max w -0.273457618782564 0.399769697706082

global min, max u -89.6490892159091 87.8639388836426

global min, max w -0.285363029305200 0.382125987627251

global min, max u -89.6638228363405 86.5007533220126

global min, max w -0.155157847525878 0.313438963651474

global min, max u -59.0884139742292 43.7549736556630

global min, max w NaN NaN

global min, max u NaN NaN

global min, max w NaN NaN

global min, max u NaN NaN

global min, max w NaN NaN

global min, max u NaN NaN

global min, max w NaN NaN

Here is my namelist.atmosphere:

[whuang@cy001 run.535554]$ cat namelist.atmosphere

&nhyd_model

config_time_integration_order = 2

config_dt = 3.0

config_start_time = '2017-10-25_00:00:00'

config_run_duration = '1_00:00:00'

config_split_dynamics_transport = true

config_number_of_sub_steps = 2

config_dynamics_split_steps = 3

config_h_mom_eddy_visc2 = 0.0

config_h_mom_eddy_visc4 = 0.0

config_v_mom_eddy_visc2 = 0.0

config_h_theta_eddy_visc2 = 0.0

config_h_theta_eddy_visc4 = 0.0

config_v_theta_eddy_visc2 = 0.0

config_horiz_mixing = '2d_smagorinsky'

config_len_disp = 120000.0

config_visc4_2dsmag = 0.05

config_w_adv_order = 3

config_theta_adv_order = 3

config_scalar_adv_order = 3

config_u_vadv_order = 3

config_w_vadv_order = 3

config_theta_vadv_order = 3

config_scalar_vadv_order = 3

config_scalar_advection = true

config_positive_definite = false

config_monotonic = true

config_coef_3rd_order = 0.25

config_epssm = 0.1

config_smdiv = 0.1

config_apvm_upwinding = 0.0

/

&damping

config_zd = 22000.0

config_xnutr = 0.2

/

&io

config_pio_num_iotasks = 0

config_pio_stride = 1

/

&decomposition

config_block_decomp_file_prefix = '../data/x4.535554.graph.info.part.'

/

&restart

config_do_restart = false

/

&printout

config_print_global_minmax_vel = true

config_print_detailed_minmax_vel = false

/

&IAU

config_IAU_option = 'off'

config_IAU_window_length_s = 21600.

/

&physics

config_sst_update = false

config_sstdiurn_update = false

config_deepsoiltemp_update = false

config_radtlw_interval = '00:30:00'

config_radtsw_interval = '00:30:00'

config_bucket_update = 'none'

config_physics_suite = 'mesoscale_reference'

/

&soundings

config_sounding_interval = 'none'

/

and my namelist.init_atmophere:

[whuang@cy001 run.535554]$ cat namelist.init_atmosphere

&nhyd_model

config_init_case = 7

config_start_time = '2017-10-25_00:00:00'

config_stop_time = '2017-10-25_00:00:00'

config_theta_adv_order = 3

config_coef_3rd_order = 0.25

/

&dimensions

config_nvertlevels = 41

config_nsoillevels = 4

config_nfglevels = 38

config_nfgsoillevels = 4

/

&data_sources

config_geog_data_path = '/lustre/whuang/data4wrf/pre_geog/'

config_met_prefix = 'FILE'

config_sfc_prefix = 'SST'

config_fg_interval = 86400

config_landuse_data = 'USGS'

config_topo_data = 'GTOPO30'

config_use_spechumd = false

/

&vertical_grid

config_ztop = 30000.0

config_nsmterrain = 1

config_smooth_surfaces = true

config_dzmin = 0.3

config_nsm = 30

config_tc_vertical_grid = true

/

&interpolation_control

config_extrap_airtemp = 'linear'

/

&preproc_stages

config_static_interp = false

config_native_gwd_static = false

config_vertical_grid = true

config_met_interp = true

config_input_sst = false

config_frac_seaice = true

/

&io

config_pio_num_iotasks = 0

config_pio_stride = 1

/

&decomposition

config_block_decomp_file_prefix = '../data/x4.535554.graph.info.part.'

/

I checked initial u,w min/max value:

ncl 0> f = addfile("x4.535554.init.nc", "r")

ncl 1> u = f->u

ncl 2> w = f->w

ncl 4> print("min u:" + min(u) + ", max u:" + max(u))

(0) min u:-92.82956594606307, max u:92.71529482341505

ncl 5> print("min w:" + min(w) + ", max w:" + max(w))

(0) min w:-0.7512803869532644, max w:0.6458138433241239

Which seems reasonable.

Could you please help me figure out what is wrong here?

Thanks,

Wei

Bill Skamarock

unread,

Aug 30, 2018, 12:55:35 PM8/30/18

to Wei Huang, MPAS-Atmosphere Help

Hi Wei,

A couple of things.

(1) You need to set the namelist.atmosphere parameter "config_len_disp" to the smallest nominal horizontal cell-to-cell distance (as stated in the User's Guide, page 29). The default namelist.atmosphere sets this parameter for a 120 km mesh, (i.e. 120,000 meters). You need to set it to 15,000 meters.

(2) You should be able to use a 90 second tilmestep with this mesh.

If, with these 2 changes, you still see unphysical velocities after the first time step (as indicated in your output) then that suggests there is a problem with the initial data.

I do not think that ncl commands min and max will see NaNs, you might want to try the NCL function isnan_ieee to see if there are any in the initial data.

Regards,

Bill

--
You received this message because you are subscribed to the Google Groups "MPAS-Atmosphere Help" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mpas-atmosphere-help+unsub...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

--

Dr. William C. Skamarock, Senior Scientist
National Center for Atmospheric Research
Mesoscale and Microscale Meteorology Division
P.O. Box 3000
Boulder, Colorado 80307-3000
Phone: (303) 497-8161
Fax: (303) 497-8171
emails...@ucar.edu
http://www2.mmm.ucar.edu/people/skamarock

Wei Huang

unread,

Aug 30, 2018, 1:29:13 PM8/30/18

to MPAS-Atmosphere Help

Bill,

With those two changes, it runs to end smoothly.

It seems the key is "config_len_disp", which I did not understand. You explanation helped my know it more.

Certainly the time_step change makes the integration more efficiently.

I know NCL won't print NaN with min/max function.

The value I printed is from the initial condition, which I try to say that the initial u/w is a valid range,

which means my problem is in the model setup. Here it is, as you have pointed it out.