About 4D array in Registry.EM

[studymore]

New 4D arrays are needed in my research.Now one questions come out. First, I want to let you know how I do as follows:
1, Added the 4D array MyNewArray in Registry.EM.

state real - ikjftb MyNewArray l - - -
state real test1 ikjftb MyNewArray l - i01rhusdf=(bdy_interp:dt) "test1" "My test1" "no unit"
state real test2 ikjftb MyNewArray l - i01rhusdf=(bdy_interp:dt) "test2" "My test2" "no unit"
2, Added test1 and test2 to microphysics scheme "linshceme" in package in Registry.EM

package kesslerscheme mp_physics==1 - moist:qv,qc,qr;MyNewArray:test1,test2

In linshceme, array test1 and test2 are needed for calculation. Then I did some modification in solve_em.F, module_mirophysics_driver.F and module_mp_lin.F. In namelist.input, max_dom = 2. When I ran wrf.exe, error came out as follows:

......
Timing for Writing wrfout_d01_2004-08-12_00:00:00 for domain 1: 0.14800 elapsed seconds.
Timing for processing lateral boundary for domain 1: 0.02800 elapsed seconds.
WRF NUMBER OF TILES = 1
Timing for Writing wrfout_d02_2004-08-12_00:00:00 for domain 2: 0.26500 elapsed seconds.
WRF NUMBER OF TILES = 1
Segmentation fault

Someone can help me ?
Thanks a lot!!!!!!!!!!!

[studymore]

If I delete "MyNewArray:test1,test2" in the followings:

package kesslerscheme mp_physics==1 - moist:qv,qc,qr;MyNewArray:test1,test2

wrf.exe could run successfully and calculation in module_mp_lin.F is correct. However, I could not find test1 and test2 in wrfout_d01... and wrfout_d02... files.

why?

[studymore]

if max_dom is set to 1, then wrf.exe could run successfully too. Problems may be caused by nest.

[alexgc]

so did you find a fix for this?

[jimmyc]

But I don't see a 4d array being constructed. I see you specified a 4d array and only have two arrays instead of 4....