test/VisuPrs/Animation/D2.py

   1 # Copyright (C) 2010-2016  CEA/DEN, EDF R&D
   2 #
   3 # This library is free software; you can redistribute it and/or
   4 # modify it under the terms of the GNU Lesser General Public
   5 # License as published by the Free Software Foundation; either
   6 # version 2.1 of the License, or (at your option) any later version.
   7 #
   8 # This library is distributed in the hope that it will be useful,
   9 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 # Lesser General Public License for more details.
  12 #
  13 # You should have received a copy of the GNU Lesser General Public
  14 # License along with this library; if not, write to the Free Software
  15 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 #
  17 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 #
  19
  20 #This case corresponds to: /visu/animation/D2 case
  21 #%Create animation for Gauss Points for 'pression' field of the the given MED file and dumps picture files in TIFF format %
  22
  23 import sys
  24 import os
  25 from paravistest import *
  26 from presentations import *
  27 from pvsimple import *
  28
  29 # Directory for saving snapshots
  30 picturedir = get_picture_dir("Animation/D2")
  31
  32 theFileName = datadir +  "TimeStamps.med"
  33 print(" --------------------------------- ")
  34 print("file ", theFileName)
  35 print(" --------------------------------- ")
  36
  37 OpenDataFile(theFileName)
  38 aProxy = GetActiveSource()
  39 if aProxy is None:
  40         raise RuntimeError("Error: can't import file.")
  41 else: print("OK")
  42
  43 print("Creating a Viewer.........................", end=' ')
  44 aView = GetRenderView()
  45 time_value = get_time(aProxy, 0)
  46 aView.ViewTime = time_value
  47 UpdatePipeline(time=time_value, proxy=aProxy)
  48 reset_view(aView)
  49 Render(aView)
  50
  51 if aView is None : print("Error")
  52 else : print("OK")
  53
  54 # Gauss Points  creation
  55 prs= GaussPointsOnField(aProxy,EntityType.CELL,'pression' , 2,  multiplier=0.01)
  56 prs.Visibility=1
  57 aView.ResetCamera()
  58 print("Creating an Animation.....................", end=' ')
  59 my_format = "tiff"
  60 print("Current format to save snapshots: ",my_format)
  61 # Add path separator to the end of picture path if necessery
  62 if not picturedir.endswith(os.sep):
  63     picturedir += os.sep
  64
  65 # Select only the current field:
  66 aProxy.AllArrays = ['TS0/dom/ComSup0/pression@@][@@P0']
  67 time_value = get_time(aProxy, 0)
  68 aView.ViewTime = time_value
  69 UpdatePipeline(time=time_value, proxy=aProxy)
  70
  71 # Animation creation and saving into set of files into picturedir
  72 scene = AnimateReader(aProxy,aView,picturedir+"D2_dom."+my_format)
  73 nb_frames = len(scene.TimeKeeper.TimestepValues)
  74
  75 pics = os.listdir(picturedir)
  76 if len(pics) != nb_frames:
  77    print("FAILED!!! Number of made pictures is equal to ", len(pics), " instead of ", nb_frames)
  78
  79 for pic in pics:
  80     os.remove(picturedir+pic)
  81
  82 # Prepare animation  performance
  83 scene.PlayMode = 1 #  set RealTime mode for animation performance
  84 # set period
  85 scene.Duration = 30 # correspond to set the speed of animation in VISU
  86 scene.GoToFirst()
  87 print("Animation.................................", end=' ')
  88 scene.Play()
  89 scene.GoToFirst()