Hi guys,
I am interested in tracer_opt = 2, but I can not find any files to learn its mechanism.
So I have some questions about it.
1) The tr17_5 and tr17_6 are stratospheric tracers, right? What is the difference between them?
2) The domains in namelist.input is
max_dom = 2,
s_we = 1, 1,
e_we = 141, 340,
s_sn = 1, 1,
e_sn = 94, 205,
i_parent_start = 1, 14,
j_parent_start = 1, 13,
parent_grid_ratio = 1, 3.
I want to release tracers only in the entire domain 2. And I set this
IF (config_flags%tracer_opt .eq. 2) THEN
DO j = jds+13, jde + 94, 1
DO i = ids+14, ide+ 113, 1
IF ( its .LE. i .and. ite .GE. i .and. jts .LE. j .and. jte .GE. j ) THEN
tracer(i, 1, j, P_tr17_5) = 1.
tracer(i, 1, j, P_tr17_6) = 1.
END IF
END DO
END DO
END IF
in solve_em.F and module_initialize_real.F files.
The i_parent_start of domain 2 is 14, so I set ids+ 14. Then the ide+ 113 is form e_we(d01, 141) - i_parent_start(d02, 14) - 14(right boundary of domain 2). And j is the same as it.
I'm not sure if the range of i and j is correct. Perhaps setting ide + 113 to ide + 127(i_parent_end) would be better?
3) I have already run the result, and I display the result of 100 hPa as follow:
Why are there few tracers north of 30°N? Perhaps due to the stratospheric tracers do not display tracers in the stratosphere, but only in the troposphere?
Or that's only caused by the settings in solve_em.F and module_initialize_real.F files.
Looking forward to your reply!
Thanks!
I am interested in tracer_opt = 2, but I can not find any files to learn its mechanism.
So I have some questions about it.
1) The tr17_5 and tr17_6 are stratospheric tracers, right? What is the difference between them?
2) The domains in namelist.input is
max_dom = 2,
s_we = 1, 1,
e_we = 141, 340,
s_sn = 1, 1,
e_sn = 94, 205,
i_parent_start = 1, 14,
j_parent_start = 1, 13,
parent_grid_ratio = 1, 3.
I want to release tracers only in the entire domain 2. And I set this
IF (config_flags%tracer_opt .eq. 2) THEN
DO j = jds+13, jde + 94, 1
DO i = ids+14, ide+ 113, 1
IF ( its .LE. i .and. ite .GE. i .and. jts .LE. j .and. jte .GE. j ) THEN
tracer(i, 1, j, P_tr17_5) = 1.
tracer(i, 1, j, P_tr17_6) = 1.
END IF
END DO
END DO
END IF
in solve_em.F and module_initialize_real.F files.
The i_parent_start of domain 2 is 14, so I set ids+ 14. Then the ide+ 113 is form e_we(d01, 141) - i_parent_start(d02, 14) - 14(right boundary of domain 2). And j is the same as it.
I'm not sure if the range of i and j is correct. Perhaps setting ide + 113 to ide + 127(i_parent_end) would be better?
3) I have already run the result, and I display the result of 100 hPa as follow:
Why are there few tracers north of 30°N? Perhaps due to the stratospheric tracers do not display tracers in the stratosphere, but only in the troposphere?
Or that's only caused by the settings in solve_em.F and module_initialize_real.F files.
Looking forward to your reply!
Thanks!