1 : C DECO SOURCE AM 08/02/14 21:20:05 6050
2 : SUBROUTINE DECO
3 : IMPLICIT INTEGER(I-N)
4 : IMPLICIT REAL*8(A-H,O-Z)
5 : *--------------------------------------------------------------------*
6 : * *
7 : * Opérateur DEnsite de COurant *
8 : * ____________________________ *
9 : * *
10 : * Calcul de la densite de courant de Foucault *
11 : * *
12 : * Syntaxe n›1: DNS1 = DECO MOD1 DEP1 ( CAR1 ) ; *
13 : * *
14 : * Entr{es: *
15 : * *
16 : * MOD1 objet MMODEL *
17 : * DEP1 CHPOINT (FONCTION DE COURANT) *
18 : * CAR1 MCHAML de sous-type CARACTERISTIQUES *
19 : * *
20 : * Sortie: *
21 : * *
22 : * DNS1 MCHAML de sous-type GRADIENT *
23 : * *
24 : * *
25 : * Auteurs, date de cr{ation: *
26 : * *
27 : * *
28 : * Y. STEPHAN le 22/09/97 (copie de GRAD) *
29 : *--------------------------------------------------------------------*
30 : -INC CCOPTIO
31 : -INC CCHAMP
32 : *
33 : * NINFOS est le nombre d'informations contenues dans INFELE
34 : * ( VOIR LE S-P ELQUOI )
35 : *
36 : PARAMETER ( NINFOS=15 )
37 : *
38 : * MC 19/01/98 : Fonction indisponible en défo planes généralisées
39 : IF(IFOUR.EQ.-3) THEN
39 : CALL ERREUR(710)
40 : RETURN
41 : ENDIF
42 : *
43 : IPCHL1=0
44 : IPCHE1=0
45 : IRET1=0
46 : IRET2=0
47 : IRET3=0
48 : ICONV=0
49 : *
50 : *
51 : * LECTURE D'UN MMODEL
52 : *
52 : CALL LIROBJ('MMODEL',IPMODL,1,IRET)
53 : IF (IERR.NE.0) RETURN
54 : *
55 : IF (IRET.NE.0) THEN
56 : *
57 : * 1-}RE SYNTAXE
58 : * _____________
59 : *
60 : *
61 : * LECTURE D'UN CHPOINT
62 : *
62 : CALL LIROBJ('CHPOINT',IPCHP1,0,IRET1)
63 : IF(IERR.NE.0) RETURN
64 : IF (IRET1.NE.0) THEN
65 : *
66 : * CHPOINT ---> CHAMELEM AUX NOEUDS
67 : *
67 : CALL CHAME1(0,IPMODL,IPCHP1,' ',IPCHE2,1)
68 : IF (IERR.NE.0) RETURN
69 : ENDIF
70 : *
71 : * LECTURE D'UN MCHAML (CARACTERISTIQUE)
72 : *
72 : CALL LIROBJ('MCHAML',IPCHE1,0,IRET3)
73 : IF (IERR.NE.0) RETURN
74 : *
75 : ENDIF
76 : *
77 : * GESTION DES OP{RANDES
78 : * IPCHL1 = 0 (MCHAML calcule)
79 : *
79 : CALL DECO1(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)
80 : *
81 : IF (IRET.EQ.1) THEN
82 : IF (ICONV.EQ.0) THEN
82 : IF (IRET1.NE.0) CALL DTCHAM(IPCHE2)
82 : CALL ECROBJ('MCHAML',IPCHL1)
83 : ELSE
83 : CALL POPEL3(IPAFF1,1,IPCHL1,IPCHEL2)
83 : IF (IPCHE1.NE.0) CALL DTCHAM(IPCHE1)
83 : CALL DTCHAM(IPCHL1)
83 : CALL DTMODL(IPMODL)
83 : CALL ECROBJ('CHAMELEM',IPCHEL2)
84 : ENDIF
85 : ENDIF
86 : RETURN
87 : END
88 :
89 :
90 :
91 :
92 :
93 :
94 :
95 :
© Cast3M 2003 - All rights reserved.
Disclaimer