OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cbatempel.F File Reference
#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine cbatempel (jft, jlt, ng, ixc, temp, tempel)

Function/Subroutine Documentation

◆ cbatempel()

subroutine cbatempel ( integer jft,
integer jlt,
integer ng,
integer, dimension(nixc,*) ixc,
temp,
tempel )

Definition at line 28 of file cbatempel.F.

30C-----------------------------------------------
31C I m p l i c i t T y p e s
32C-----------------------------------------------
33#include "implicit_f.inc"
34C-----------------------------------------------
35C D u m m y A r g u m e n t s
36C-----------------------------------------------
37 INTEGER JFT, JLT,NG,IXC(NIXC,*)
38 my_real temp(*),tempel(*)
39C-----------------------------------------------
40C L o c a l V a r i a b l e s
41C-----------------------------------------------
42 INTEGER I,KSI ,ETA ,SF1,SF2 ,SF3 ,SF4
43 my_real
44 . vpg(2,4),pg1,pg
45 parameter(pg=.577350269189626)
46 parameter(pg1=-.577350269189626)
47C--------------------------
48C INITIALISATION
49C-------------------------
50 DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
51C-----------------------------------------------
52C --- Temperature in Gauss point NG
53 ksi = vpg(1,ng)
54 eta = vpg(2,ng)
55
56 sf1 = (1-ksi)*(1-eta)
57 sf2 = (1+ksi)*(1-eta)
58 sf3 = (1+ksi)*(1+eta)
59 sf4 = (1-ksi)*(1+eta)
60
61 DO i=jft,jlt
62 tempel(i) = fourth * (sf1 *temp(ixc(2,i)) + sf2*temp(ixc(3,i)) +
63 . sf3 *temp(ixc(4,i)) + sf4*temp(ixc(5,i)) )
64 ENDDO
65C
66 RETURN
#define my_real
Definition cppsort.cpp:32