Esope source files

   1 : C TRJDTS    SOURCE    CHAT      05/01/13    03:49:35     5004           
   2 :       SUBROUTINE TRJDTS(NSAUV,MLREE6,TMIN,KSAUV,DTSTOC,DTCUMU)
   3 : ***************************************************************************
   4 : ***  Permet de selectionner le pas de tps de sauvegarde ds la
   5 : *** liste 'MLREE6' considéré lors de la prochaine sauvegarde de resultats.
   6 : ***
   7 : *** APPELES 1 = aucun
   8 : *** APPELES 2 = aucun
   9 : ***
  10 : *** E = 'NSAUV' taille de la liste des tps de sauvegarde
  11 : ***     'MLREE6' liste des tps de sauvegarde
  12 : ***     'TMIN' instant de depart du lacher de la particule
  13 : ***
  14 : *** E/S = 'KSAUV' indice liste des tps de sauvegarde considéré 
  15 : ***       'DTSTOC' pas de tps de sauvegarde considéré
  16 : ***       'DTCUMU' cumul des pas de tps entre deux sauvegardes
  17 : ****************************************************************************** 
  18 : 
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8 (A-H,O-Z)
  21 : -INC SMLREEL
  22 : -INC CCREEL
  23 :       POINTEUR MLREE6.MLREEL
  24 :       IF (NSAUV.LT.0) THEN
  25 : ***      sauvegarde des qu'un saut de particule est effectué
  26 :          DTSTOC=0.D0
  27 :       ELSEIF (NSAUV.EQ.0) THEN
  28 : ***      sauvegarde à intervalle de tps régulier (1er tps de 'MLREE6')
  29 :          DTSTOC=MLREE6.PROG(1)
  30 :       ELSEIF (KSAUV.LT.NSAUV) THEN
  31 : ***      sauvegarde selon intervalles des tps de la liste 'MLREE6'
  32 :          KSAUV=KSAUV+1
  33 :          IF (KSAUV.EQ.1) THEN
  34 :             TPS1=TMIN
  35 :          ELSE
  36 :             TPS1=MLREE6.PROG(KSAUV-1)
  37 :          ENDIF
  38 :          TPS2=MLREE6.PROG(KSAUV)
  39 :          DTSTOC=TPS2-TPS1
  40 :       ELSE
  41 : ***      les pts ne sont plus sauvegardés lorsque 'MLREE6' entierement lue
  42 :          DTSTOC=XGRAND
  43 :       ENDIF
  44 : ***   remise à 0 du cumul des pas de tps apres la sauvegarde
  45 :       DTCUMU=0.D0
  46 : 
  47 :       RETURN
  48 :       END
  49 :  
  50 :