Esope source files

   1 : C MEDRCH    SOURCE    CHAT      11/10/26    22:15:25     7176           
   2 :       subroutine medrch(ifid,mlmot1,mlmot2,nbnoin,nommaa,nbpta,icret)
   3 :       implicit real*8(a-h,o-z)
   4 :       implicit integer (i-n)
   5 : *
   6 : * lecture des champs s'il y a lieu
   7 : *      
   8 : -INC CCOPTIO
   9 : -INC SMELEME
  10 : -INC SMCHPOI
  11 : -INC SMCHAML
  12 : -INC SMTABLE
  13 : -INC SMLMOTS
  14 :       character*4 cha4a,cha4b,cha4c
  15 :       character*8 cha8a,ityc,cha8b,cha8c
  16 :       character*32 nocham,cha32a,cha32b,cha32c
  17 :       segment ichamp
  18 :          character*16 nomcom(nc+1),iunit(nc+1)
  19 :       endsegment 
  20 :       segment iprofi        
  21 :          character*32 nompro(m)
  22 :          integer melem(m)
  23 :       endsegment
  24 :       character*32 nommaa      
  25 :       parameter (MED_MAILLE=0,MED_FACE=1,MED_ARETE=2,MED_NOEUD=3)
  26 :       parameter (MED_NO_PFLMOD=0, MED_GLOBAL=1, MED_COMPACT=2)
  27 :                   
  28 : *      write(6,*) ' medrch : appel efncha'
  29 :       ind=0
  29 :       call efncha4(ifid,ind,nbch,icret)
  30 :       if(icret.ne.0) then
  31 : *        write(6,*) ' plouf 10'
  31 :         call efexit4(-1)
  32 :         return
  33 :       endif
  34 : *      write(6,*) ' sortie efncha nombre de champs ', nbch
  35 : *
  36 : * quelques initialisations 
  37 : *
  38 : *
  39 :       nat=1
  40 :       nsoupo=1
  41 :       nbnn=1
  42 :       nbsous=0
  43 :       nbref=0
  44 :       n1ptel=0
  45 :       n1el=0
  46 :       n2ptel=0
  47 :       n2el=0
  48 :       m=20
  49 :       nprof=0
  50 :       segini iprofi 
  51 :       do ia=1,nbch
  52 : * combien de composantes
  53 : *        write(6,*) ' lecture du nb de composantes du champ ' , ia
  53 :         call efncha4(ifid,ia,nc,icret)
  54 : *        write(6,*) ' nb de composante du champ ' , nc 
  55 :         if(icret.ne.0) then
  56 : *          write(6,*) ' plouf 11'     
  56 :           call efexit4(-1)  
  57 :           return
  58 :         endif     
  59 :         segini ichamp
  60 : *        write(6,*) ' avant appelefchai'  
  60 : ******        call efchai(ifid,ia,nocham,itypch,nomcom,iunit,nc,icret) 
  60 :         call efchai4(ifid,ia,nocham,itypch,nomcom,iunit,nc,icret)
  61 : *        write(6,*) ' info issu de efchai '
  62 : *        write(6,*) ' nom du champ ', nocham , ' nb composantes ' ,nc
  63 : *        write(6,*) 'nom des composantes ' ,(nomcom(iou),iou=1,nc) 
  64 : *        write(6,*) ' type du champ ' ,itypch       
  65 :         if(icret.ne.0) then
  66 : *          write(6,*) ' plouf 12'
  66 :           call efexit4(-1)
  67 :           return
  68 :         endif   
  69 : *
  70 : *   est-ce un champ par point?si oui on le lit
  71 : * 
  72 : *        write(6,*) ' avant appel a medchp '      
  72 :         call medchp(ifid,nocham,itypch,nc,nbnoin,nommaa,iprofi,nprof,
  73 :      $         ichamp,nbpta,ich,icret)
  74 :         if(icret.ne.0) then
  75 : *          write(6,*) ' plouf 121'
  75 :           call efexit4(-1)
  76 :           return
  77 :         endif        
  78 : *        if(ich.ne.0) then
  79 : *           write(6,*) ' on vient de lire un chpoint de pointeur', ich
  80 : *        else 
  81 : *           write(6,*) ' le champ n est pas un chpoint'
  82 : *        endif        
  82 : *        call medcen(ifid,nocham,itypch,nc,nbnoin,nommaa,iprofi,nprof,
  83 : *     $         ichamp,nbpta,ich,icret)
  84 :       enddo
  85 :       return
  86 :       end
  87 :          
  88 :  
  89 :