From 737a90bdfc8eb45a19d4c771c88cf7bc678e6c0b Mon Sep 17 00:00:00 2001
From: Paul-Gerhard Reinhard <mpt218@homework.(none)>
Date: Mon, 9 May 2016 14:29:40 +0200
Subject: [PATCH] isovector external field corrected
---
Code/dynamic.f90 | 1 +
Code/external.f90 | 11 +++++++----
2 files changed, 8 insertions(+), 4 deletions(-)
diff --git a/Code/dynamic.f90 b/Code/dynamic.f90
index f7b6208..eb5e401 100644
--- a/Code/dynamic.f90
+++ b/Code/dynamic.f90
@@ -51,6 +51,7 @@ CONTAINS
END IF
! external boost
IF(texternal) THEN
+ CALL init_external
CALL extboost(text_timedep)
END IF
! Create protocol files
diff --git a/Code/external.f90 b/Code/external.f90
index 3a2c844..1eac1c0 100644
--- a/Code/external.f90
+++ b/Code/external.f90
@@ -14,13 +14,15 @@ MODULE External
CONTAINS
!***********************************************************************
SUBROUTINE getin_external
+ NAMELIST/extern/ amplq0,radext,widext,isoext,ipulse,omega,tau0,taut, &
+ textfield_periodic
+ READ(5,extern)
+ END SUBROUTINE getin_external
+ SUBROUTINE init_external
REAL(db) :: facn,facp,facr,xlim,ylim,zlim
INTEGER :: ix,iy,iz
CHARACTER(14),PARAMETER :: pulsetype(0:2)=(/ 'Instantaneous ', &
'Gaussian ','Cosine squared' /)
- NAMELIST/extern/ amplq0,radext,widext,isoext,ipulse,omega,tau0,taut, &
- textfield_periodic
- READ(5,extern)
IF(ipulse<0.OR.ipulse>2) STOP &
' External field: called with invalid pulse type'
IF(wflag) THEN
@@ -49,6 +51,7 @@ CONTAINS
facn=-1.0D0/(mass_number-charge_number)
facp=1.0D0/charge_number
ENDIF
+ WRITE(*,*) 'EXTERNAL: ',facn,facp
ALLOCATE(extfield(nx,ny,nz,2))
xlim=nx*dx
ylim=ny*dy
@@ -68,7 +71,7 @@ CONTAINS
ENDDO
ENDDO
ENDDO
- END SUBROUTINE getin_external
+ END SUBROUTINE init_external
!***********************************************************************
SUBROUTINE extfld(time)
REAL(db) :: time
--
GitLab