f/test24.f
C* This file is part of MED.
C*
C* COPYRIGHT (C) 1999 - 2021 EDF R&D, CEA/DEN
C* MED is free software: you can redistribute it and/or modify
C* it under the terms of the GNU Lesser General Public License as published by
C* the Free Software Foundation, either version 3 of the License, or
C* (at your option) any later version.
C*
C* MED is distributed in the hope that it will be useful,
C* but WITHOUT ANY WARRANTY; without even the implied warranty of
C* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
C* GNU Lesser General Public License for more details.
C*
C* You should have received a copy of the GNU Lesser General Public License
C* along with MED. If not, see <http://www.gnu.org/licenses/>.
C*
C *******************************************************************************
C * - Nom du fichier : test24.f
C *
C * - Description : lecture de mailles MED_POLYGONE dans le maillage MED
C * du fichier test23.med
C *
C ******************************************************************************
program test23
C
implicit none
include 'med.hf'
C
integer*8 fid
integer cret,mdim,nmaa,npoly,i,j,k,taille
integer edim,nstep,stype,atype, chgt, tsf
character*64 maa
character*200 desc
integer ni, n, isize;
parameter(ni=4, n=3)
integer index(ni),ind1,ind2
character*16 nom(n)
integer num(n),fam(n)
integer con(16)
integer type
character*16 nomcoo(2)
character*16 unicoo(2)
character(16) :: dtunit
C
C Ouverture du fichier test23.med en lecture seule
call mfiope(fid,'test23.med',med_acc_rdonly, cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur ouverture du fichier'
call efexit(-1)
endif
print *,'Ouverture du fichier test23.med'
C
C Lecture du nombre de maillages
call mmhnmh(fid,nmaa,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur lecture nombre de maillage'
call efexit(-1)
endif
print *,'Nombre de maillages : ',nmaa
C
C Lecture de toutes les mailles MED_POLYGONE
C dans chaque maillage
do 10 i=1,nmaa
C
C Info sur chaque maillage
call mmhmii(fid,i,maa,edim,mdim,type,desc,
& dtunit,stype,nstep,atype,
& nomcoo,unicoo,cret)
if (cret .ne. 0 ) then
print *,'Erreur lecture infos maillage'
call efexit(-1)
endif
print *,cret
print *,'Maillage : ',maa
print *,'Dimension : ',mdim
C
C Combien de mailles polygones
call mmhnme(fid,maa,med_no_dt,med_no_it,med_cell,med_polygon,
& med_index_node,med_nodal,chgt,tsf,isize,cret)
npoly = isize - 1;
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur lecture du nombre de polygone'
call efexit(-1)
endif
print *,'Nombre de mailles MED_POLYGONE : ',npoly
C
C Taille des connectivites
call mmhnme(fid,maa,med_no_dt,med_no_it,med_cell,med_polygon,
& med_connectivity,med_nodal,chgt,tsf,taille,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur lecture infos polygones'
call efexit(-1)
endif
print *,'Taille de la connectivite : ',taille
C
C Lecture de la connectivite
call mmhpgr(fid,maa,med_no_dt,med_no_it,med_cell,
& med_nodal,index,con,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur lecture des connectivites polygones'
call efexit(-1)
endif
print *,'Lecture de la connectivite des polygones'
C
C Lecture des noms
call mmhear(fid,maa,med_no_dt,med_no_it,
& med_cell,med_polygon,nom,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur lecture des noms des polygones'
call efexit(-1)
endif
print *,'Lecture des noms'
C
C Lecture des numeros
call mmhfnr(fid,maa,med_no_dt,med_no_it,med_cell,med_polygon,
& num,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur lecture des numeros des polygones'
call efexit(-1)
endif
print *,'Lecture des numeros'
C
C Lecture des numeros de familles
call mmhfnr(fid,maa,med_no_dt,med_no_it,med_cell,med_polygon,
& fam,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur lecture des numeros de famille des
& polygones'
call efexit(-1)
endif
print *,'Lecture des numeros de famille'
C
C Affichage des resultats
print *,'Affichage des resultats'
do 20 j=1,npoly
C
print *,'>> Maille polygone ',j
print *,'---- Connectivite ---- : '
ind1 = index(j)
ind2 = index(j+1)
do 30 k=ind1,ind2-1
print *,con(k)
30 continue
c print *,'---- Nom ---- : ',nom(j)
print *,'---- Numero ----: ',num(j)
print *,'---- Numero de famille ---- : ',fam(j)
C
20 continue
C
10 continue
C
C Fermeture du fichier
call mficlo(fid,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur fermeture du fichier'
call efexit(-1)
endif
print *,'Fermeture du fichier'
C
end
mmhnme
subroutine mmhnme(fid, name, numdt, numit, entype, geotype, datype, cmode, chgt, tsf, n, cret)
Cette routine permet de lire le nombre d'entités dans un maillage pour une étape de calcul donnée.
Definition: medmesh.f:551
mmhfnr
subroutine mmhfnr(fid, name, numdt, numit, entype, geotype, num, cret)
Cette routine permet la lecture des numéros de famille d'un type d'entité d'un maillage.
Definition: medmesh.f:487
mficlo
subroutine mficlo(fid, cret)
Fermeture d'un fichier MED.
Definition: medfile.f:82
mmhnmh
subroutine mmhnmh(fid, n, cret)
Cette routine permet de lire le nombre de maillages dans un fichier.
Definition: medmesh.f:41
mfiope
subroutine mfiope(fid, name, access, cret)
Ouverture d'un fichier MED.
Definition: medfile.f:42
mmhmii
subroutine mmhmii(fid, it, name, sdim, mdim, mtype, desc, dtunit, stype, nstep, atype, aname, aunit, cret)
Cette routine permet de lire les informations relatives à un maillage dans un fichier.
Definition: medmesh.f:110
test23
program test23
Definition: test23.f:24
mmhear
subroutine mmhear(fid, mname, numdt, numit, entype, geotype, ename, cret)
Cette routine permet de lire les noms d'un type d'entité d'un maillage.
Definition: medmesh.f:529
mmhpgr
subroutine mmhpgr(fid, name, numdt, numit, entype, cmode, index, con, cret)
Cette routine permet la lecture des connectivités de polygones.
Definition: medmesh.f:912