sensor_energy_bilan.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "parit_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sensor_energy_bilan (jft, jlt, ei, ek, off, ipart, itask, sensors)

Function/Subroutine Documentation

sensor_energy_bilan()

subroutine sensor_energy_bilan	(	integer, intent(in)	jft,
		integer, intent(in)	jlt,
		intent(in)	ei,
		intent(in)	ek,
		intent(in)	off,
		integer, dimension(*), intent(in)	ipart,
		integer, intent(in)	itask,
		type (sensors_), intent(inout)	sensors )

Definition at line 50 of file sensor_energy_bilan.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE sensor_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "parit_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(in) :: JFT,JLT
      INTEGER, INTENT(in) :: ITASK
      INTEGER, DIMENSION(*), INTENT(in) ::  IPART
C     REAL
      my_real, DIMENSION(*), INTENT(in) :: ei,ek
      my_real, DIMENSION(*), INTENT(in) :: off
      type (sensors_),INTENT(INOUT) :: SENSORS
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,MX,II,J,IC,NEL,IJK
      INTEGER :: K,IPART_SENS,LOCAL_INDEX
      my_real, DIMENSION(MVSIZ,2) :: fstparit   
      INTEGER :: NGR_SENSOR
C-----------------------------------------------
      ngr_sensor=sensors%NGR_SENSOR(itask+1)
      IF(iparit>0) THEN
         IF(sensors%SENS_GROUP(ngr_sensor)%NUM_PART>0) THEN
            IF(sensors%SENS_GROUP(ngr_sensor)%NUM_PART==1) THEN
                DO i=jft,jlt
                    IF( off(i)/=zero ) THEN
                        fstparit(i,1) = ei(i)
                        fstparit(i,2) = ek(i)
                    ELSE
                        fstparit(i,1) = zero
                        fstparit(i,2) = zero
                    ENDIF
                ENDDO
                nel = jlt-jft+1
                k = sensors%SENS_GROUP(ngr_sensor)%PART(1,3)    !   id of the sensor
                ijk = itask+1
                CALL sum_6_float(1,nel,fstparit(1,1),sensors%SENSOR_STRUCT(k)%FBSAV6_SENS(1,1,ijk),2)
                CALL sum_6_float(1,nel,fstparit(1,2),sensors%SENSOR_STRUCT(k)%FBSAV6_SENS(2,1,ijk),2)
            ELSEIF(sensors%SENS_GROUP(ngr_sensor)%NUM_PART>1) THEN
 
               DO j=1,sensors%SENS_GROUP(ngr_sensor)%NUM_PART
                    ipart_sens = sensors%SENS_GROUP(ngr_sensor)%PART(j,1)   !   id of the part linked to the sensor
                    local_index = 0     
                    DO i=jft,jlt
                        mx=ipart(i)
                        IF( (ipart_sens==mx).AND.(off(i)/=zero) ) THEN
                            local_index = local_index + 1
                            fstparit(local_index,1) = ei(i)
                            fstparit(local_index,2) = ek(i)   
                        ENDIF
                    ENDDO
 
                    k = sensors%SENS_GROUP(ngr_sensor)%PART(j,3)   !   id of the sensor
                    ijk = itask+1
                    CALL sum_6_float(1,local_index,fstparit(1,1),sensors%SENSOR_STRUCT(k)%FBSAV6_SENS(1,1,ijk),2)
                    CALL sum_6_float(1,local_index,fstparit(1,2),sensors%SENSOR_STRUCT(k)%FBSAV6_SENS(2,1,ijk),2)
                ENDDO
            ENDIF
         ENDIF
      ENDIF          
C---
      RETURN