f/test27.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 : test27.f
C *
C * - Description : creation de maillages structures (grille cartesienne |
C * grille standard ) dans le fichier test27.med
C *
C *****************************************************************************
program test27
C
implicit none
include 'med.hf'
C
C
integer*8 fid
integer cret
C ** la dimension du maillage **
integer mdim,sdim
C ** nom du maillage de longueur maxi MED_TAILLE_NOM **
character*64 maa
C ** le nombre de noeuds **
integer nnoe
C ** table des coordonnees **
real*8 coo(8)
character*16 nomcoo(2), unicoo(2)
character*200 desc
integer strgri(2)
C ** grille cartesienne **
integer axe,nind
real*8 indice(4)
C
C
data coo /0.0,0.0,1.0,0.0,0.0,1.0,1.0,1.0/
data nomcoo /"x","y"/, unicoo /"cm","cm"/
C
C Creation du fichier test27.med
call mfiope(fid,'test27.med',med_acc_rdwr, cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur creation du fichier'
call efexit(-1)
endif
print *,'Creation du fichier test27.med'
C
C Creation d'un maillage MED_NON_STRUCTURE
mdim = 2
sdim = 2
maa = 'maillage vide'
desc = 'un maillage vide'
call mmhcre(fid,maa,mdim,sdim,med_unstructured_mesh,
& desc,"",med_sort_dtit,med_cartesian,
& nomcoo,unicoo,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur creation du maillage'
call efexit(-1)
endif
C
C Creation d'une grille cartesienne
mdim = 2
maa = 'grille cartesienne'
desc = 'un exemple de grille cartesienne'
call mmhcre(fid,maa,mdim,sdim,med_structured_mesh,
& desc,"",med_sort_dtit,med_cartesian,
& nomcoo,unicoo,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur creation du maillage'
call efexit(-1)
endif
print *,'Creation d un maillage MED_STRUCTURE'
C
C On specifie la nature du maillage structure
call mmhgtw(fid,maa,med_cartesian_grid,cret)
print *,cret
print *,'On definit la nature de la grille :
& MED_GRILLE_CARTESIENNE'
if (cret .ne. 0 ) then
print *,'Erreur ecriture de la nature de la grille'
call efexit(-1)
endif
C
C On definit les indices de la grille selon chaque dimension
indice(1) = 1.1d0
indice(2) = 1.2d0
indice(3) = 1.3d0
indice(4) = 1.4d0
nind = 4
axe = 1
call mmhgcw(fid,maa,med_no_dt,med_no_it,med_undef_dt,
& axe,nind,indice,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur ecriture des indices'
call efexit(-1)
endif
print *,'Ecriture des indices des coordonnees selon axe X'
C
indice(1) = 2.1d0
indice(2) = 2.2d0
indice(3) = 2.3d0
indice(4) = 2.4d0
nind = 4
axe = 2
call mmhgcw(fid,maa,med_no_dt,med_no_it,med_undef_dt,
& axe,nind,indice,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur ecriture des indices'
call efexit(-1)
endif
print *,'Ecriture des indices des coordonnees selon axe Y'
C
C Creation d'une grille MED_CURVILINEAR_GRID de dimension 2
maa = 'grille curviligne'
mdim = 2
desc = 'un exemple de grille curviligne'
call mmhcre(fid,maa,mdim,sdim,med_structured_mesh,
& desc,"",med_sort_dtit,med_cartesian,
& nomcoo,unicoo,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur creation de maillage'
call efexit(-1)
endif
print *,'Nouveau maillage MED_STRUCTURE'
C
call mmhgtw(fid,maa,med_curvilinear_grid,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur ecriture de la nature de la grille'
call efexit(-1)
endif
print *,'On definit la nature du maillage : MED_GRILLE_STANDARD'
C
C On ecrit les coordonnes de la grille
nnoe = 4
call mmhcow(fid,maa,med_no_dt,med_no_it,med_undef_dt,
& med_full_interlace,nnoe,coo,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur ecriture des coordonnees des noeuds'
call efexit(-1)
endif
print *,'Ecriture des coordonnees de la grille'
C
C On definit la structure des coordonnees de la grille
strgri(1) = 2
strgri(2) = 2
call mmhgsw(fid,maa,med_no_dt,med_no_it,med_undef_dt,
& strgri,cret)
print *,cret
if (cret .ne. 0 ) then
print *,'Erreur ecriture de la structure'
call efexit(-1)
endif
print *,'Ecriture de la structure de la grille : / 2,2 /'
C
C On ferme le 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
mficlo
subroutine mficlo(fid, cret)
Fermeture d'un fichier MED.
Definition: medfile.f:82
mmhcow
subroutine mmhcow(fid, name, numdt, numit, dt, swm, n, coo, cret)
Cette routine permet d'écrire dans un maillage le tableau des coordonnées des noeuds,...
Definition: medmesh.f:299
test27
program test27
Definition: test27.f:25
mmhcre
subroutine mmhcre(fid, name, sdim, mdim, mtype, desc, dtunit, stype, atype, aname, aunit, cret)
Cette routine permet de créer un maillage dans un fichier.
Definition: medmesh.f:20
mfiope
subroutine mfiope(fid, name, access, cret)
Ouverture d'un fichier MED.
Definition: medfile.f:42
mmhgsw
subroutine mmhgsw(fid, name, numdt, numit, dt, st, cret)
Cette routine définit la structure (nombre de points sur chaque axe du repère) d'un maillage structur...
Definition: medmesh.f:259
mmhgcw
subroutine mmhgcw(fid, name, numdt, numit, dt, axis, size, index, cret)
Cette routine permet l'écriture des coordonnées des noeuds d'un maillage structuré selon un axe du re...
Definition: medmesh.f:383
mmhgtw
subroutine mmhgtw(fid, name, gtype, cret)
Cette routine permet de définir le type d'un maillage structuré (MED_STRUCTURED_MESH).
Definition: medmesh.f:223