diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Dragon/src/MCTPSP.f | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Dragon/src/MCTPSP.f')
| -rw-r--r-- | Dragon/src/MCTPSP.f | 64 |
1 files changed, 64 insertions, 0 deletions
diff --git a/Dragon/src/MCTPSP.f b/Dragon/src/MCTPSP.f new file mode 100644 index 0000000..2a8cb66 --- /dev/null +++ b/Dragon/src/MCTPSP.f @@ -0,0 +1,64 @@ +*DECK MCTPSP + SUBROUTINE MCTPSP(IPTRK,POS,IREG,IEV) +*----------------------------------------------------------------------- +* +*Purpose: +* Store position and region index in TRACKING table for PSP display. +* +*Copyright: +* Copyright (C) 2008 Ecole Polytechnique de Montreal +* This library is free software; you can redistribute it and/or +* modify it under the terms of the GNU Lesser General Public +* License as published by the Free Software Foundation; either +* version 2.1 of the License, or (at your option) any later version +* +*Author(s): R. Le Tellier +* +*Parameters: input +* IPTRK pointer to the TRACKING data structure. +* POS point global coordinates. +* IREG region/surface index. +* IEV event index. +* +*----------------------------------------------------------------------- +* + USE GANLIB + IMPLICIT NONE +*---- +* SUBROUTINE ARGUMENTS +*---- + TYPE(C_PTR) IPTRK + INTEGER IREG,IEV + DOUBLE PRECISION POS(3) +*---- +* LOCAL VARIABLES +*---- + INTEGER ILONG,ITYLCM,NPOINT +*---- +* ALLOCATABLE ARRAYS +*---- + INTEGER, ALLOCATABLE, DIMENSION(:) :: REGI,EVENT + DOUBLE PRECISION, ALLOCATABLE, DIMENSION(:,:) :: COORD +* + CALL LCMSIX(IPTRK,'MCpoints',1) + CALL LCMLEN(IPTRK,'REGI',ILONG,ITYLCM) + NPOINT=ILONG+1 + ALLOCATE(COORD(3,NPOINT),REGI(NPOINT),EVENT(NPOINT)) + IF (ILONG.GT.0) THEN + CALL LCMGET(IPTRK,'COORD',COORD) + CALL LCMGET(IPTRK,'REGI',REGI) + CALL LCMGET(IPTRK,'EVENT',EVENT) + ENDIF + COORD(1,NPOINT)=POS(1) + COORD(2,NPOINT)=POS(2) + COORD(3,NPOINT)=POS(3) + REGI(NPOINT)=IREG + EVENT(NPOINT)=IEV + CALL LCMPUT(IPTRK,'COORD',3*NPOINT,4,COORD) + CALL LCMPUT(IPTRK,'REGI',NPOINT,1,REGI) + CALL LCMPUT(IPTRK,'EVENT',NPOINT,1,EVENT) + DEALLOCATE(EVENT,REGI,COORD) + CALL LCMSIX(IPTRK,' ',2) +* + RETURN + END |