fct_t.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	fct_t

Function/Subroutine Documentation

fct_t()

subroutine fct_t

(

)

Definition at line 45 of file fct_t.f90.

  !#  if defined (1)

  !==============================================================================|
  USE all_vars
  USE mod_utils
  USE bcs
  USE mod_obcs
  IMPLICIT NONE
  REAL(SP):: TMAX,TMIN
  INTEGER :: I,J,K
  !==============================================================================|


  IF(heating_type == 'body') RETURN

  IF(dbg_set(dbg_sbr)) WRITE(ipt,*)"Start: fct_t"

  nodes: DO i=1,m
     ! SKIP OPEN BOUNDARY NODES
     IF(iobcn > 0)THEN
        DO j=1,iobcn
           IF(i == i_obc_n(j)) cycle nodes
        END DO
     END IF

     ! SKIP RIVER INFLOW POINTS
     IF(numqbc > 0)THEN
        DO j=1,numqbc
           IF(river_inflow_location == 'node')THEN
              IF(i == inodeq(j)) cycle nodes
           END IF
           IF(river_inflow_location == 'edge')THEN
              IF(i == n_icellq(j,1) .OR. i == n_icellq(j,2)) cycle nodes
           END IF
        END DO
     END IF

     ! SKIP GROUND WATER INFLOW POINTS
     IF(bfwdis(i) .GT. 0.0_sp .and. groundwater_temp_on) cycle nodes

     DO k=1,kbm1
        tmax = maxval(t1(nbsn(i,1:ntsn(i)),k))
        tmin = minval(t1(nbsn(i,1:ntsn(i)),k))

        IF(k == 1)THEN
           tmax = max(tmax,(t1(i,k)*dz(i,k+1)+t1(i,k+1)*dz(i,k))/  &
                (dz(i,k)+dz(i,k+1)))
           tmin = min(tmin,(t1(i,k)*dz(i,k+1)+t1(i,k+1)*dz(i,k))/  &
                (dz(i,k)+dz(i,k+1)))
        ELSE IF(k == kbm1)THEN
           tmax = max(tmax,(t1(i,k)*dz(i,k-1)+t1(i,k-1)*dz(i,k))/  &
                (dz(i,k)+dz(i,k-1)))
           tmin = min(tmin,(t1(i,k)*dz(i,k-1)+t1(i,k-1)*dz(i,k))/  &
                (dz(i,k)+dz(i,k-1)))
        ELSE
           tmax = max(tmax,(t1(i,k)*dz(i,k-1)+t1(i,k-1)*dz(i,k))/  &
                (dz(i,k)+dz(i,k-1)),                             &
                (t1(i,k)*dz(i,k+1)+t1(i,k+1)*dz(i,k))/           &
                (dz(i,k)+dz(i,k+1)))
           tmin = min(tmin,(t1(i,k)*dz(i,k-1)+t1(i,k-1)*dz(i,k))/  &
                (dz(i,k)+dz(i,k-1)),                             &
                (t1(i,k)*dz(i,k+1)+t1(i,k+1)*dz(i,k))/           &
                (dz(i,k)+dz(i,k+1)))
        END IF

        IF(tmin-tf1(i,k) > 0.0_sp)tf1(i,k) = tmin
        IF(tf1(i,k)-tmax > 0.0_sp)tf1(i,k) = tmax

     END DO

  END DO nodes

  IF(dbg_set(dbg_sbr)) WRITE(ipt,*)"End: fct_t"
  !#  endif

Here is the call graph for this function:

Here is the caller graph for this function: