From 2daafd3962294008684fabd66d62e3370fe594bc Mon Sep 17 00:00:00 2001
From: Sergey Yurchenko <sergei.yurchenko@gmail.com>
Date: Fri, 17 May 2024 15:42:09 +0100
Subject: [PATCH] Update pot_CH3OH_Bowman.f90

---
 pot_CH3OH_Bowman.f90 | 18 ++++++++++++++----
 1 file changed, 14 insertions(+), 4 deletions(-)

diff --git a/pot_CH3OH_Bowman.f90 b/pot_CH3OH_Bowman.f90
index 8763ed4..e855550 100644
--- a/pot_CH3OH_Bowman.f90
+++ b/pot_CH3OH_Bowman.f90
@@ -62,19 +62,29 @@ function MLpoten(ncoords,natoms,local,xyz,force) result(f)
    real(ark),intent(in)   ::  local(ncoords)
    real(ark),intent(in)   ::  xyz(natoms,3)
    real(ark),intent(in)   ::  force(:)
-   real(ark)              ::  f
+   real(ark)              ::  f,xyz_user(3,natoms)
    integer(ik) :: m,mr,N
    !
    m  = int(force(1),ik)
    mr = int(force(2),ik)
    !
    N = size(force(:))
+   ! 
+   ! convert TROVE xyz to uzer's XYZ
+   !
+   xyz_user(:,1) = xyz(3,:)
+   xyz_user(:,2) = xyz(4,:)
+   xyz_user(:,3) = xyz(5,:)
+   xyz_user(:,4) = xyz(6,:)
+   xyz_user(:,5) = xyz(1,:)
+   xyz_user(:,6) = xyz(2,:)
+   !
+   xyz_user = xyz_user/bohr
    !
-   call potshell(force(3:N),m,mr,xyz,f)
+   call potshell(force(3:N),m,mr,xyz_user,f)
    !
  end function MLpoten
- 
- 
+
  
  !
  subroutine potshell(coef,m,mr,xn,V)