Domain decomposition of latitude and longitude array

[JasonWok]

Update:
This problem has been resolved.
This post will not be deleted in case of anyone finds it useful in the future.

---

I have one innermost domain with e_we = e_sn = 52.
And I used 4 processes to compute the WRF:
mpiexec_mpt -spawn -np 4 ./wrf.exe

In phys/module_sf_bep_bem.F, link, how can I get actual decomposed domain (latitude and longitude information) for each process without the grey/padding grids?

For example I have the following outputs for rsl.out.0001. And the definition of xlatall and xlongall are stated here.:

its, ite          27          51 jts, jte           1          26                                                                                                                                                                                                         
 ims, ime          17          57 jms, jme          -4          36                                                                                                                                                                                                         
 shape of xlatall          41          41 shape of xlongall          41          41