adjust_ts.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	adjust_ts

Function/Subroutine Documentation

adjust_ts()

subroutine adjust_ts

(

)

Definition at line 46 of file adjust_ts.f90.

!==============================================================================|
   USE all_vars
   USE mod_utils
   use mod_par
   USE bcs
   IMPLICIT NONE
   REAL(SP) :: TAVE,TAVE1,TAVE2
   INTEGER :: I,K,JJ,I1,J,J1,J2,NUM_TAVE,NUM_TAVE1,NUM_TAVE2
!==============================================================================|

   IF (dbg_set(dbg_sbr)) WRITE(ipt,*) "Start: adjust_ts"
   
   IF(numqbc > 0)THEN   

!     IF(RIVER_INFLOW_LOCATION == 'node')THEN
!       DO K=1,KBM1
!         DO I=1,NUMQBC
!           JJ=INODEQ(I)
!           TAVE = 0.0_SP
!           NUM_TAVE = 0
!           DO J=2,NTSN(JJ)-1
!             I1=NBSN(JJ,J)
!        IF(NUMQBC == 1)THEN
!          NUM_TAVE = NUM_TAVE + 1
!          TAVE = TAVE + T1(I1,K)
!        ELSE
!          IF(I == 1)THEN
!            IF(I1 /= INODEQ(I+1))THEN 
!            NUM_TAVE = NUM_TAVE + 1
!            TAVE = TAVE + T1(I1,K)
!        END IF
!          ELSE IF(I == NUMQBC)THEN
!            IF(I1 /= INODEQ(I-1))THEN   
!            NUM_TAVE = NUM_TAVE + 1
!            TAVE = TAVE + T1(I1,K)
!                 END IF
!          ELSE IF(I1 /= INODEQ(I-1) .AND. I1 /= INODEQ(I+1))THEN
!            NUM_TAVE = NUM_TAVE + 1
!                 TAVE = TAVE + T1(I1,K)
!          END IF    
!        END IF
!           END DO
!           T1(JJ,K) = TAVE/FLOAT(NUM_TAVE)
!         END DO
!       END DO
!     ELSE IF(RIVER_INFLOW_LOCATION == 'edge')THEN
!       DO K=1,KBM1
!         DO I=1,NUMQBC
!           J1=N_ICELLQ(I,1)
!           J2=N_ICELLQ(I,2)
!      TAVE1 = 0.0_SP
!      TAVE2 = 0.0_SP
!      NUM_TAVE1 = 0
!      NUM_TAVE2 = 0

!           DO J=2,NTSN(J1)-1
!             I1=NBSN(J1,J)
!        IF(NUMQBC == 1)THEN
!          IF(I1 /= J2)THEN
!            NUM_TAVE1 = NUM_TAVE1 + 1
!                 TAVE1 = TAVE1 + T1(I1,K)
!          END IF
!        ELSE IF(I == 1)THEN
!          IF(I1 /= J2 .AND. I1 /= N_ICELLQ(I+1,1) .AND. &
!            I1 /= N_ICELLQ(I+1,2))THEN
!            NUM_TAVE1 = NUM_TAVE1 + 1
!                 TAVE1 = TAVE1 + T1(I1,K)
!          END IF
!        ELSE IF(I == NUMQBC)THEN
!          IF(I1 /= J2 .AND. I1 /= N_ICELLQ(I-1,1) .AND. &
!            I1 /= N_ICELLQ(I-1,2))THEN
!            NUM_TAVE1 = NUM_TAVE1 + 1
!                 TAVE1 = TAVE1 + T1(I1,K)
!          END IF
!             ELSE IF(I1 /= J2 .AND. &
!           I1 /= N_ICELLQ(I-1,1) .AND. I1 /= N_ICELLQ(I-1,2) .AND.  &
!       I1 /= N_ICELLQ(I+1,1) .AND. I1 /= N_ICELLQ(I+1,2))THEN
!          NUM_TAVE1 = NUM_TAVE1 + 1
!               TAVE1 = TAVE1 + T1(I1,K)
!        END IF
!           END DO
!           T1(J1,K) = TAVE1/FLOAT(NUM_TAVE1)

!           DO J=2,NTSN(J2)-1
!             I1=NBSN(J2,J)
!        IF(NUMQBC == 1)THEN
!          IF(I1 /= J1)THEN
!            NUM_TAVE2 = NUM_TAVE2 + 1
!                 TAVE2 = TAVE2 + T1(I1,K)
!          END IF
!        ELSE IF(I == 1)THEN
!          IF(I1 /= J1 .AND. I1 /= N_ICELLQ(I+1,1) .AND. &
!            I1 /= N_ICELLQ(I+1,2))THEN
!            NUM_TAVE2 = NUM_TAVE2 + 1
!                 TAVE2 = TAVE2 + T1(I1,K)
!          END IF
!        ELSE IF(I == NUMQBC)THEN  
!          IF(I1 /= J1 .AND. I1 /= N_ICELLQ(I-1,1) .AND. &
!            I1 /= N_ICELLQ(I-1,2))THEN
!            NUM_TAVE2 = NUM_TAVE2 + 1
!                 TAVE2 = TAVE2 + T1(I1,K)
!          END IF
!        ELSE IF(I1 /= J1 .AND. &
!           I1 /= N_ICELLQ(I-1,1) .AND. I1 /= N_ICELLQ(I-1,2) .AND.  &
!       I1 /= N_ICELLQ(I+1,1) .AND. I1 /= N_ICELLQ(I+1,2))THEN
!          NUM_TAVE2 = NUM_TAVE2 + 1
!               TAVE2 = TAVE2 + T1(I1,K)
!        END IF
!           END DO
!           T1(J2,K) = TAVE2/FLOAT(NUM_TAVE2)
!      
!         END DO
!       END DO   
!     END IF
   
     DO i=1,numqbc
       IF(river_inflow_location == 'node')THEN
         j = inodeq(i)
         DO k=1,kbm1
           t1(j,k) = max(t1(j,k),tdis(i))
           s1(j,k) = max(s1(j,k),sdis(i))
        END DO
       ELSE IF(river_inflow_location == 'edge')THEN
         j1 = n_icellq(i,1)
         j2 = n_icellq(i,2)
         DO k=1,kbm1
           t1(j1,k) = max(t1(j1,k),tdis(i))
           t1(j2,k) = max(t1(j2,k),tdis(i))
           s1(j1,k) = max(s1(j1,k),sdis(i))
           s1(j2,k) = max(s1(j2,k),sdis(i))
         END DO
       END IF
     END DO

     CALL n2e3d(t1,t)
     CALL n2e3d(s1,s)

   END IF

  IF (dbg_set(dbg_sbr)) WRITE(ipt,*) "End: adjust_ts"
 

   RETURN

Here is the call graph for this function:

Here is the caller graph for this function: