Esope source files

   1 : C ISOVAL    SOURCE    PV        09/03/12    21:25:36     6325           
   2 :       SUBROUTINE ISOVAL
   3 :       IMPLICIT REAL*8 (A-H,O-Z)
   4 :       IMPLICIT INTEGER (I-N)
   5 : C***********************************************************************
   6 : C NOM         : ISOVAL
   7 : C DESCRIPTION : Construit le maillage de l'isovaleur d'un champ par 
   8 : C               éléments
   9 : C
  10 : C
  11 : C
  12 : C LANGAGE     : ESOPE
  13 : C AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  14 : C               mél : gounand@semt2.smts.cea.fr
  15 : C***********************************************************************
  16 : C APPELES          : 
  17 : C APPELES (E/S)    : 
  18 : C APPELES (BLAS)   :
  19 : C APPELES (CALCUL) :
  20 : C APPELE PAR       :
  21 : C***********************************************************************
  22 : C SYNTAXE GIBIANE    :
  23 : C ENTREES            :
  24 : C ENTREES/SORTIES    :
  25 : C SORTIES            :
  26 : C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  27 : C***********************************************************************
  28 : C VERSION    : v1, 15/12/2008, version initiale
  29 : C HISTORIQUE : v1, 15/12/2008, création
  30 : C HISTORIQUE :
  31 : C HISTORIQUE :
  32 : C***********************************************************************
  33 : C Prière de PRENDRE LE TEMPS de compléter les commentaires
  34 : C en cas de modification de ce sous-programme afin de faciliter
  35 : C la maintenance !
  36 : C***********************************************************************
  37 : -INC CCOPTIO
  38 : -INC CCREEL
  39 : -INC SMCHAML
  40 :       CHARACTER*8 MCHA
  41 :       INTEGER IMPR,IRET
  42 : *     
  43 : * Executable statements
  44 : *
  45 :       IMPR=0
  46 :       IF (IMPR.GT.2) WRITE(IOIMP,*) 'Entrée dans isoval.eso' 
  47 : *
  48 :       MCHA='MCHAML  '
  48 :       CALL LIROBJ(MCHA,MCHELM,1,IRETOU)
  49 : *  37 2
  50 : * On ne trouve pas d'objet de type %m1:8
  51 :       IF (IRETOU.NE.1) THEN
  52 :          MOTERR(1:8)=MCHA
  52 :          CALL ERREUR(37)
  53 :          RETURN
  54 :       ENDIF
  55 : *
  55 :       CALL LIRREE(XISO,1,IRETOU)
  56 : *  37 2
  57 : * On ne trouve pas d'objet de type %m1:8
  58 :       IF (IRETOU.NE.1) THEN
  59 :          MOTERR(1:8)='FLOTTANT'
  59 :          CALL ERREUR(37)
  60 :          RETURN
  61 :       ENDIF
  62 : *
  63 : * Recherche du maximum
  64 : *    
  65 :       IPLMOT=0
  66 :       IPLACE=0
  67 :       KPLUS=1
  68 :       LABSO=1
  68 :       CALL MAXICH(MCHELM,IPLMOT,MCHA,IPLACE,XMAX,KPLUS,LABSO)
  69 :       XTOL=MAX(XMAX*XZPREC,XPETIT)
  70 : *
  71 : * Calcul de l'isovaleur
  72 : *
  72 :       CALL ISOVA1(MCHELM,XISO,XTOL,MELEME)
  73 : *
  74 : * On renvoie le résultat
  75 : *
  75 :       CALL ECROBJ('MAILLAGE',MELEME)
  76 :       RETURN
  77 : *
  78 : * End of subroutine ISOVAL
  79 : *
  80 :       END
  81 : 
  82 :  
  83 :