Skip to content

Commit 6590929

Browse files
author
Loren
committed
minimal changes (Diagnostics_,(Base.F90,Custom.F90)) for custom magnetic outputs
1 parent 6790d65 commit 6590929

File tree

2 files changed

+1072
-90
lines changed

2 files changed

+1072
-90
lines changed

‎src/Diagnostics/Diagnostics_Base.F90‎

Lines changed: 113 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -81,11 +81,119 @@ Module Diagnostics_Base
8181
!//////////////////////////////////////////////////////////
8282
! Custom Outputs: range from ...
8383
Integer, Parameter :: custom_offset = dcheck_off+100 !2200
84-
Integer, Parameter :: cross_helicity = custom_offset + 1 ! v dot B
85-
Integer, Parameter :: turb_cross_helicity = custom_offset+2
86-
Integer, Parameter :: ell0_vr = custom_offset+3
87-
Integer, Parameter :: ell0_tvar = custom_offset+4
88-
Integer, Parameter :: ell0_dpdr = custom_offset+5
84+
Integer, Parameter :: set_offset = 15
85+
Integer, Parameter :: set_offset2 = 12
86+
Integer, Parameter :: set_offset3 = 9
87+
88+
! total
89+
Integer, Parameter :: induct_work_r = custom_offset + 1
90+
Integer, Parameter :: induct_work_t = custom_offset + 2
91+
Integer, Parameter :: induct_work_p = custom_offset + 3
92+
93+
Integer, Parameter :: ishear_work_r = custom_offset + 4
94+
Integer, Parameter :: ishear_work_t = custom_offset + 5
95+
Integer, Parameter :: ishear_work_p = custom_offset + 6
96+
97+
Integer, Parameter :: iadvec_work_r = custom_offset + 7
98+
Integer, Parameter :: iadvec_work_t = custom_offset + 8
99+
Integer, Parameter :: iadvec_work_p = custom_offset + 9
100+
101+
Integer, Parameter :: icomp_work_r = custom_offset + 10
102+
Integer, Parameter :: icomp_work_t = custom_offset + 11
103+
Integer, Parameter :: icomp_work_p = custom_offset + 12
104+
105+
Integer, Parameter :: idiff_work_r = custom_offset + 13
106+
Integer, Parameter :: idiff_work_t = custom_offset + 14
107+
Integer, Parameter :: idiff_work_p = custom_offset + 15
108+
109+
!pmp
110+
Integer, Parameter :: induct_work_pmp_r = custom_offset + 1 + set_offset
111+
Integer, Parameter :: induct_work_pmp_t = custom_offset + 2 + set_offset
112+
Integer, Parameter :: induct_work_pmp_p = custom_offset + 3 + set_offset
113+
114+
Integer, Parameter :: ishear_work_pmp_r = custom_offset + 4 + set_offset
115+
Integer, Parameter :: ishear_work_pmp_t = custom_offset + 5 + set_offset
116+
Integer, Parameter :: ishear_work_pmp_p = custom_offset + 6 + set_offset
117+
118+
Integer, Parameter :: iadvec_work_pmp_r = custom_offset + 7 + set_offset
119+
Integer, Parameter :: iadvec_work_pmp_t = custom_offset + 8 + set_offset
120+
Integer, Parameter :: iadvec_work_pmp_p = custom_offset + 9 + set_offset
121+
122+
Integer, Parameter :: icomp_work_pmp_r = custom_offset + 10 + set_offset
123+
Integer, Parameter :: icomp_work_pmp_t = custom_offset + 11 + set_offset
124+
Integer, Parameter :: icomp_work_pmp_p = custom_offset + 12 + set_offset
125+
126+
!ppm
127+
Integer, Parameter :: induct_work_ppm_r = custom_offset + 1 + set_offset + set_offset2
128+
Integer, Parameter :: induct_work_ppm_t = custom_offset + 2 + set_offset + set_offset2
129+
Integer, Parameter :: induct_work_ppm_p = custom_offset + 3 + set_offset + set_offset2
130+
131+
Integer, Parameter :: ishear_work_ppm_r = custom_offset + 4 + set_offset + set_offset2
132+
Integer, Parameter :: ishear_work_ppm_t = custom_offset + 5 + set_offset + set_offset2
133+
Integer, Parameter :: ishear_work_ppm_p = custom_offset + 6 + set_offset + set_offset2
134+
135+
Integer, Parameter :: iadvec_work_ppm_r = custom_offset + 7 + set_offset + set_offset2
136+
Integer, Parameter :: iadvec_work_ppm_t = custom_offset + 8 + set_offset + set_offset2
137+
Integer, Parameter :: iadvec_work_ppm_p = custom_offset + 9 + set_offset + set_offset2
138+
139+
Integer, Parameter :: icomp_work_ppm_r = custom_offset + 10 + set_offset + set_offset2
140+
Integer, Parameter :: icomp_work_ppm_t = custom_offset + 11 + set_offset + set_offset2
141+
Integer, Parameter :: icomp_work_ppm_p = custom_offset + 12 + set_offset + set_offset2
142+
143+
!mmm
144+
Integer, Parameter :: induct_work_mmm_r = custom_offset + 1 + set_offset + 2*set_offset2
145+
Integer, Parameter :: induct_work_mmm_t = custom_offset + 2 + set_offset + 2*set_offset2
146+
Integer, Parameter :: induct_work_mmm_p = custom_offset + 3 + set_offset + 2*set_offset2
147+
148+
Integer, Parameter :: ishear_work_mmm_r = custom_offset + 4 + set_offset + 2*set_offset2
149+
Integer, Parameter :: ishear_work_mmm_t = custom_offset + 5 + set_offset + 2*set_offset2
150+
Integer, Parameter :: ishear_work_mmm_p = custom_offset + 6 + set_offset + 2*set_offset2
151+
152+
Integer, Parameter :: iadvec_work_mmm_r = custom_offset + 7 + set_offset + 2*set_offset2
153+
Integer, Parameter :: iadvec_work_mmm_t = custom_offset + 8 + set_offset + 2*set_offset2
154+
Integer, Parameter :: iadvec_work_mmm_p = custom_offset + 9 + set_offset + 2*set_offset2
155+
156+
Integer, Parameter :: icomp_work_mmm_r = custom_offset + 10 + set_offset + 2*set_offset2
157+
Integer, Parameter :: icomp_work_mmm_t = custom_offset + 11 + set_offset + 2*set_offset2
158+
Integer, Parameter :: icomp_work_mmm_p = custom_offset + 12 + set_offset + 2*set_offset2
159+
160+
Integer, Parameter :: idiff_work_mm_r = custom_offset + 13 + set_offset + 2*set_offset2
161+
Integer, Parameter :: idiff_work_mm_t = custom_offset + 14 + set_offset + 2*set_offset2
162+
Integer, Parameter :: idiff_work_mm_p = custom_offset + 15 + set_offset + 2*set_offset2
163+
164+
!mpp
165+
Integer, Parameter :: induct_work_mpp_r = custom_offset + 1 + 2*set_offset + 2*set_offset2
166+
Integer, Parameter :: induct_work_mpp_t = custom_offset + 2 + 2*set_offset + 2*set_offset2
167+
Integer, Parameter :: induct_work_mpp_p = custom_offset + 3 + 2*set_offset + 2*set_offset2
168+
169+
Integer, Parameter :: ishear_work_mpp_r = custom_offset + 4 + 2*set_offset + 2*set_offset2
170+
Integer, Parameter :: ishear_work_mpp_t = custom_offset + 5 + 2*set_offset + 2*set_offset2
171+
Integer, Parameter :: ishear_work_mpp_p = custom_offset + 6 + 2*set_offset + 2*set_offset2
172+
173+
Integer, Parameter :: iadvec_work_mpp_r = custom_offset + 7 + 2*set_offset + 2*set_offset2
174+
Integer, Parameter :: iadvec_work_mpp_t = custom_offset + 8 + 2*set_offset + 2*set_offset2
175+
Integer, Parameter :: iadvec_work_mpp_p = custom_offset + 9 + 2*set_offset + 2*set_offset2
176+
177+
!ppp
178+
Integer, Parameter :: induct_work_ppp_r = custom_offset + 1 + 2*set_offset + 2*set_offset2 + set_offset3
179+
Integer, Parameter :: induct_work_ppp_t = custom_offset + 2 + 2*set_offset + 2*set_offset2 + set_offset3
180+
Integer, Parameter :: induct_work_ppp_p = custom_offset + 3 + 2*set_offset + 2*set_offset2 + set_offset3
181+
182+
Integer, Parameter :: ishear_work_ppp_r = custom_offset + 4 + 2*set_offset + 2*set_offset2 + set_offset3
183+
Integer, Parameter :: ishear_work_ppp_t = custom_offset + 5 + 2*set_offset + 2*set_offset2 + set_offset3
184+
Integer, Parameter :: ishear_work_ppp_p = custom_offset + 6 + 2*set_offset + 2*set_offset2 + set_offset3
185+
186+
Integer, Parameter :: iadvec_work_ppp_r = custom_offset + 7 + 2*set_offset + 2*set_offset2 + set_offset3
187+
Integer, Parameter :: iadvec_work_ppp_t = custom_offset + 8 + 2*set_offset + 2*set_offset2 + set_offset3
188+
Integer, Parameter :: iadvec_work_ppp_p = custom_offset + 9 + 2*set_offset + 2*set_offset2 + set_offset3
189+
190+
Integer, Parameter :: icomp_work_ppp_r = custom_offset + 10 + 2*set_offset + 2*set_offset2 + set_offset3
191+
Integer, Parameter :: icomp_work_ppp_t = custom_offset + 11 + 2*set_offset + 2*set_offset2 + set_offset3
192+
Integer, Parameter :: icomp_work_ppp_p = custom_offset + 12 + 2*set_offset + 2*set_offset2 + set_offset3
193+
194+
Integer, Parameter :: idiff_work_pp_r = custom_offset + 13 + 2*set_offset + 2*set_offset2 + set_offset3
195+
Integer, Parameter :: idiff_work_pp_t = custom_offset + 14 + 2*set_offset + 2*set_offset2 + set_offset3
196+
Integer, Parameter :: idiff_work_pp_p = custom_offset + 15 + 2*set_offset + 2*set_offset2 + set_offset3
89197

90198
include "turbKE_codes.F"
91199
include "axial_field_codes.F"

0 commit comments

Comments
 (0)