fvtemp.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	fvtemp (ityp, efac_in, cpa_in, cpb_in, cpc_in, cpd_in, cpe_in, cpf_in, rmw_in, temp0_in, temp_in)

Function/Subroutine Documentation

fvtemp()

subroutine fvtemp	(	integer	ityp,
		intent(in)	efac_in,
		intent(in)	cpa_in,
		intent(in)	cpb_in,
		intent(in)	cpc_in,
		intent(in)	cpd_in,
		intent(in)	cpe_in,
		intent(in)	cpf_in,
		intent(in)	rmw_in,
		intent(in)	temp0_in,
		intent(inout)	temp_in )

Definition at line 32 of file fvtemp.F.

       USE cast_mod, ONLY : double_to_my_real
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ITYP
      my_real, intent(in) ::
     .        efac_in, cpa_in, cpb_in, cpc_in, cpd_in, cpe_in, cpf_in, rmw_in, temp0_in
      my_real, intent(inout) ::  temp_in
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER K
      DOUBLE PRECISION :: CVA, HH, TEMP2, EFAC
      DOUBLE PRECISION :: CPA, CPB, CPC, CPD, CPE, CPF, RMW, TEMP0, TEMP
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      cpa=cpa_in
      cpb=cpb_in
      cpc=cpc_in
      cpd=cpd_in
      cpe=cpe_in
      cpf=cpf_in
      rmw=rmw_in
      temp0=temp0_in
      temp=temp_in
      efac = efac_in
 
      IF(ityp == 6) THEN
         cva=cpa-rmw
         IF(cpc == zero) THEN
            IF(cpb == zero) THEN
               temp=efac/cva
            ELSE IF( cpb .NE. zero) THEN
               temp=(sqrt(max(zero,cva*cva+two*cpb*efac))-cva)/cpb
            ELSE 
               temp = temp0
            ENDIF
         ELSE
            temp=temp0
            DO k=1,5
               hh=cpa+half*cpb*temp+third*cpc*temp*temp
               temp=efac/(hh-rmw)
            ENDDO
         ENDIF
      ELSEIF(ityp == 8) THEN
         temp=temp0
         DO k=1,5
            temp2=temp*temp
            hh=cpa+half*cpb*temp+third*cpc*temp2+
     .         fourth*cpd*temp*temp2-cpe/temp2+one_fifth*cpf*temp2*temp2   
            IF( (hh-rmw) > zero) THEN
              temp=efac/(hh-rmw)
            ELSE
              temp_in=zero
              RETURN
            ENDIF
         ENDDO
      ENDIF
      ! TEMP__IN**4 will be used, so we need to limit the value TEMP_IN 
      ! to avoid overflow
      ! in sp sqrt(sqrt(HUGE(TEMP_IN))) is 4.2E+9
      ! in dp sqrt(sqrt(HUGE(TEMP_IN))) is 1.2E+77
      temp_in = double_to_my_real(temp,max_value = sqrt(sqrt(huge(temp_in))))
      RETURN