34 integer ncompo, nnodes
36 integer ntria3, nquad4
38 character*64 fname, finame, lfname
40 character*16 cpname, cpunit
50 parameter(fname =
"./UsesCase_MEDfield_1.med")
51 parameter(lfname=
"./UsesCase_MEDmesh_1.med")
52 parameter(mname =
"2D unstructured mesh")
53 parameter(finame =
"TEMPERATURE_FIELD")
54 parameter(cpname =
"TEMPERATURE")
55 parameter(cpunit =
"C")
56 parameter(dtunit =
" ")
57 parameter(nnodes = 15, ncompo = 1 )
58 parameter(ntria3 = 8, nquad4 = 4)
61 data verval / 0., 100., 200., 300., 400.,
62 & 500., 600., 700., 800., 900,
63 & 1000., 1100, 1200., 1300., 1500. /
64 data tria3v / 1000., 2000., 3000., 4000.,
65 & 5000., 6000., 7000., 8000. /
66 data quad4v / 10000., 20000., 30000., 4000. /
70 call mfiope(fid,fname,med_acc_creat,cret)
71 if (cret .ne. 0 )
then
72 print *,
'ERROR : file creation'
78 call mlnliw(fid,mname,lfname,cret)
79 if (cret .ne. 0 )
then
80 print *,
'ERROR : create mesh link ...'
90 if (cret .ne. 0 )
then
91 print *,
'ERROR : create field ...'
97 call mfdrvw(fid,finame,med_no_dt,med_no_it,dt,med_node,
98 & med_none,med_full_interlace,med_all_constituent,
100 if (cret .ne. 0 )
then
101 print *,
'ERROR : write field values on vertices'
108 call mfdrvw(fid,finame,med_no_dt,med_no_it,dt,med_cell,
109 & med_tria3,med_full_interlace,med_all_constituent,
110 & ntria3,tria3v,cret)
111 if (cret .ne. 0 )
then
112 print *,
'ERROR : write field values on MED_TRIA3'
118 call mfdrvw(fid,finame,med_no_dt,med_no_it,dt,med_cell,
119 & med_quad4,med_full_interlace,med_all_constituent,
120 & nquad4,quad4v,cret)
121 if (cret .ne. 0 )
then
122 print *,
'ERROR : write field values on MED_QUAD4'
129 if (cret .ne. 0 )
then
130 print *,
'ERROR : close file'