]> SALOME platform Git repositories - tools/medcoupling.git/blob - doc/MEDMEM/FIELDcreate.py
Salome HOME
Merge from BR_V5_DEV 16Feb09
[tools/medcoupling.git] / doc / MEDMEM / FIELDcreate.py
1 #  Copyright (C) 2007-2008  CEA/DEN, EDF R&D, OPEN CASCADE
2 #
3 #  Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 #  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 #
6 #  This library is free software; you can redistribute it and/or
7 #  modify it under the terms of the GNU Lesser General Public
8 #  License as published by the Free Software Foundation; either
9 #  version 2.1 of the License.
10 #
11 #  This library is distributed in the hope that it will be useful,
12 #  but WITHOUT ANY WARRANTY; without even the implied warranty of
13 #  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 #  Lesser General Public License for more details.
15 #
16 #  You should have received a copy of the GNU Lesser General Public
17 #  License along with this library; if not, write to the Free Software
18 #  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 #
20 #  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 #
22 ######################################################################
23 # This Python script should be executed when the shared library is   #
24 # generated using SWIG 1.3 (or higher) due to the fact that older    #
25 # version could not handle the wrapping of several class constructor #
26 ######################################################################
27 #
28 from libMEDMEM_Swig import *
29
30 MedFile = "pointe.med"
31 meshName = "maa1"
32
33 myMesh = MESH(MED_DRIVER,MedFile,meshName)
34
35 mySupport = SUPPORT(myMesh,"Support on all CELLs",MED_CELL)
36
37 numberOfComponents = 3
38 myField = FIELDDOUBLE(mySupport,numberOfComponents)
39 fieldName = "fieldcelldouble"
40 myField.setName(fieldName)
41
42 for i in range(numberOfComponents):
43     if (i == 0):
44         name = "Vx"
45         desc = "vitesse selon x"
46     elif (i == 1):
47         name = "Vy"
48         desc = "vitesse selon y"
49     else:
50         name = "Vz"
51         desc = "vitesse selon z"
52     unit = "m. s-1"
53     ip1 = i+1
54     myField.setComponentName(ip1,name)
55     myField.setComponentDescription(ip1,desc)
56     myField.setMEDComponentUnit(ip1,unit)
57
58 iterationNumber = 10
59 myField.setIterationNumber(iterationNumber)
60
61 orderNumber = 1
62 myField.setOrderNumber(orderNumber)
63
64 time = 3.435678
65 myField.setTime(time)
66
67 numberOfValue = mySupport.getNumberOfElements(MED_ALL_ELEMENTS)
68
69 for i in range(numberOfValue):
70     ip1 = i+1
71     for j in range(numberOfComponents):
72         jp1 = j+1
73         value = (ip1+jp1)*0.1
74         myField.setValueIJ(ip1,jp1,value)
75
76 id = myField.addDriver(MED_DRIVER)