summaryrefslogtreecommitdiff
path: root/Trivac/data/Ktests_proc/pkinet1.c2m
diff options
context:
space:
mode:
authorstainer_t <thomas.stainer@oecd-nea.org>2025-09-08 13:48:49 +0200
committerstainer_t <thomas.stainer@oecd-nea.org>2025-09-08 13:48:49 +0200
commit7dfcc480ba1e19bd3232349fc733caef94034292 (patch)
tree03ee104eb8846d5cc1a981d267687a729185d3f3 /Trivac/data/Ktests_proc/pkinet1.c2m
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Trivac/data/Ktests_proc/pkinet1.c2m')
-rwxr-xr-xTrivac/data/Ktests_proc/pkinet1.c2m45
1 files changed, 45 insertions, 0 deletions
diff --git a/Trivac/data/Ktests_proc/pkinet1.c2m b/Trivac/data/Ktests_proc/pkinet1.c2m
new file mode 100755
index 0000000..e4b79ac
--- /dev/null
+++ b/Trivac/data/Ktests_proc/pkinet1.c2m
@@ -0,0 +1,45 @@
+***********************************************
+* Procedure : pkinet1.c2m *
+* Author : D. Sekki (04/08) *
+***********************************************
+PARAMETER KINET2 MACRO TRACK SYSTEM FLUX ::
+ ::: LINKED_LIST KINET2 MACRO TRACK SYSTEM FLUX ; ;
+MODULE INIKIN: KINSOL: BIVACA: TRIVAA: GREP: DELETE: END: ;
+LINKED_LIST SYSTEMP KINET1 ;
+INTEGER itr iedit nDel := 0 0 6 ;
+REAL delta := 0.01 ;
+INTEGER ity maxItr nadi ;
+STRING fTemp pTemp ;
+REAL precf ;
+ :: >>fTemp<< >>pTemp<< >>maxItr<<
+ >>precf<< >>nadi<< ;
+GREP: SYSTEM :: GETVAL 'STATE-VECTOR' 4 >>ity<< ;
+KINET1 := INIKIN: MACRO TRACK SYSTEM FLUX :: EDIT 1
+ NDEL <<nDel>>
+ BETA 0.000266 0.001491 0.001316
+ 0.002849 0.000896 0.000182
+ LAMBDA 0.0127 0.0317 0.1150
+ 0.3110 1.4000 3.8700
+ CHID 1.0 1.0 1.0 1.0 1.0 1.0
+ 0.0 0.0 0.0 0.0 0.0 0.0 ;
+IF ity 1 = THEN
+ SYSTEMP := BIVACA: MACRO TRACK :: EDIT 0 UNIT ;
+ELSEIF ity 2 = ity 3 = + THEN
+ SYSTEMP := TRIVAA: MACRO TRACK :: EDIT 0 UNIT ;
+ELSE
+ ECHO "Unknown type=" ity ;
+ENDIF ;
+WHILE itr maxItr < DO
+ EVALUATE itr := itr 1 + ;
+ IF itr maxItr = THEN
+ EVALUATE iedit := 1 ;
+ ENDIF ;
+ KINET1 := KINSOL: KINET1 MACRO TRACK SYSTEMP ::
+ EDIT <<iedit>> DELTA <<delta>>
+ SCHEME FLUX <<fTemp>> PREC <<pTemp>>
+ EXTE <<precf>> ADI <<nadi>> ;
+ENDWHILE ;
+KINET2 := KINET1 ;
+KINET1 SYSTEMP := DELETE: KINET1 SYSTEMP ;
+END: ;
+QUIT .