Calculate the original angular momentum magnitude

MintonGroup · Jun 20, 2021 · 50307a6 · 50307a6
1 parent d47485a
commit 50307a6
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/user/user_udio_reader.f90 b/src/user/user_udio_reader.f90
@@ -196,6 +196,7 @@ module subroutine user_udio_reader(param, unit, iotype, v_list, iostat, iomsg)
                   param_value = user_get_token(line, ifirst, ilast, iostat) 
                   read(param_value, *) param%Ltot_orig(i)
                end do
+                  param%Lmag_orig = norm2(param%Ltot_orig(:))
             case("LORBIT_ORIG")
                read(param_value, *) param%Lorbit_orig(1)
                do i = 2, NDIM