From 50307a6a53ef1ca66dffb87e60426142faa312c7 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Sat, 19 Jun 2021 23:40:54 -0400 Subject: [PATCH] Calculate the original angular momentum magnitude --- src/user/user_udio_reader.f90 | 1 + 1 file changed, 1 insertion(+) diff --git a/src/user/user_udio_reader.f90 b/src/user/user_udio_reader.f90 index ea10222e1..add51bc2d 100644 --- 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