Mengjuan Liu
New member
It seems the staggering of deformations, i.e. defor11, defor12, defor23 ... in the registry is different from how they calculated in module_diffusion_em.F. I wonder what's the intention of that? For example, defor33 is located at the center, same as defor11, defor22, but it is defined in the registry to be vertically staggered.
In Regsitry.EM_COMMON
state real defor11 ikj misc 1 - r "defor11" "DEFORMATION 11" "s-1"
state real defor22 ikj misc 1 - r "defor22" "DEFORMATION 22" "s-1"
state real defor12 ikj misc 1 - r "defor12" "DEFORMATION 12" "s-1"
state real defor33 ikj misc 1 z r "defor33" "DEFORMATION 33" "s-1"
state real defor13 ikj misc 1 z r "defor13" "DEFORMATION 13" "s-1"
state real defor23 ikj misc 1 z r "defor23" "DEFORMATION 23" "s-1"
In the subroutine cal_deform_and_div in dyn_em/module_diffusion_em.F, for example defor33 is calculated as
DO j = j_start, j_end
DO k = kts, ktf
DO i = i_start, i_end
tmp1(i,k,j) = ( w(i,k+1,j) - w(i,k,j) ) * rdzw(i,k,j)
END DO
END DO
END DO
DO j = j_start, j_end
DO k = kts, ktf
DO i = i_start, i_end
defor33(i,k,j) = 2.0 * tmp1(i,k,j)
END DO
END DO
END DO
Thanks!
In Regsitry.EM_COMMON
state real defor11 ikj misc 1 - r "defor11" "DEFORMATION 11" "s-1"
state real defor22 ikj misc 1 - r "defor22" "DEFORMATION 22" "s-1"
state real defor12 ikj misc 1 - r "defor12" "DEFORMATION 12" "s-1"
state real defor33 ikj misc 1 z r "defor33" "DEFORMATION 33" "s-1"
state real defor13 ikj misc 1 z r "defor13" "DEFORMATION 13" "s-1"
state real defor23 ikj misc 1 z r "defor23" "DEFORMATION 23" "s-1"
In the subroutine cal_deform_and_div in dyn_em/module_diffusion_em.F, for example defor33 is calculated as
DO j = j_start, j_end
DO k = kts, ktf
DO i = i_start, i_end
tmp1(i,k,j) = ( w(i,k+1,j) - w(i,k,j) ) * rdzw(i,k,j)
END DO
END DO
END DO
DO j = j_start, j_end
DO k = kts, ktf
DO i = i_start, i_end
defor33(i,k,j) = 2.0 * tmp1(i,k,j)
END DO
END DO
END DO
Thanks!