find_dt_target.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "com06_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	find_dt_target (ms, stifn, target_dt_tab, per_adm_tab, dt, tmp, dtsca, totmas, nval, nnod)

Function/Subroutine Documentation

find_dt_target()

subroutine find_dt_target	(		ms,
			stifn,
			target_dt_tab,
			per_adm_tab,
			dt,
			tmp,
			dtsca,
			totmas,
		integer	nval,
		integer	nnod )

Definition at line 29 of file find_dt_target.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   A n a l y s e   M o d u l e
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER NVAL,NNOD
        my_real
     .     ms(*),stifn(*),totmas,target_dt_tab(*),per_adm_tab(*),dtsca,dt(*),tmp(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "com06_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER I,N,COMPT,K
        INTEGER :: CPT
        my_real sumk,sumk_old,summ,summ_old,per_adm,target_dt,seuil,chunk
C=======================================================================
C
C--------------------------------------------------------------------------------------
C     EXTRACTED FROM ADD_MASS_STAT in starter - used in both starter and engine
C--------------------------------------------------------------------------------------
C
        sumk = zero
        summ = zero
        sumk_old = zero
        summ_old = zero
        compt = 1
        seuil = per_adm_tab(1)
C
        DO i=1,nnod
          n=nint(tmp(i))
        ENDDO
C
        DO i=1,nnod
          IF (i > 1) THEN
            IF (dt(i) > dt(i-1)) THEN
              sumk_old = sumk
              summ_old = summ
            ENDIF
          ENDIF
          n=nint(tmp(i))
          per_adm = (dt(i)*sumk_old - summ_old)/(max(em20,totmas))
          IF (i>1) THEN
            DO WHILE ((per_adm > seuil).AND.(compt<=nval))
              target_dt_tab(compt) = dtsca*sqrt(two*(totmas*seuil+summ_old)/max(em20,sumk_old))
              compt = compt+1
              IF(compt<=nval) seuil = per_adm_tab(compt)
            ENDDO
            IF (compt > nval) EXIT
          ENDIF
          sumk = sumk + stifn(n)
          summ = summ + ms(n)
C-- If threshold are not crossed target dt values must be computed for the last node
          IF (i==nnod) THEN
            DO k=compt,nval
              target_dt_tab(k) = dtsca*sqrt(two*(totmas*seuil+summ)/max(em20,sumk))
              IF(k+1 <= nval) seuil = per_adm_tab(k+1)
            ENDDO
          ENDIF
        ENDDO
C
 
        RETURN