thmonv.F File Reference

#include "implicit_f.inc"
#include "task_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	thmonv (j1, j2, ithbuf, l1, l2, wa, fsav, fsavvent, ivolu, iform)

Function/Subroutine Documentation

thmonv()

subroutine thmonv	(	integer	j1,
		integer	j2,
		integer, dimension(*)	ithbuf,
		integer	l1,
		integer	l2,
			wa,
			fsav,
			fsavvent,
		integer, dimension(nimv,*)	ivolu,
		integer	iform )

Definition at line 31 of file thmonv.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "task_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
      INTEGER ITHBUF(*),IVOLU(NIMV,*),IFORM
      INTEGER J1,J2,L1,L2,K,IAD2
      my_real
     .        wa(*),fsav(nthvki,*),fsavvent(5,*)
C
      INTEGER I,J,L,II,IADVENT,IV,KV,NVENT,ID_VENT,K1
      CHARACTER TITR1*40,ID_TITR1*2
C
C   gather deja effectue, seul p0 ecrit 
      IF (ispmd==0) THEN
        ii = 0
        DO j=j1,j2
          DO l=l1,l2
            ii=ii+1
            wa(ii) = zero
          ENDDO
        ENDDO
c
        ii = 0
        id_vent = 0
        iad2=j1+3*((j2-j1)+1)
        DO j=j1,j2
              i=ithbuf(j)
              CALL fretitl2(titr1,ithbuf(iad2),40)
              id_titr1(1:2) = titr1(21:22)
              READ(id_titr1,fmt='(I2)',err=100) id_vent
100           CONTINUE
c
              DO l=l1,l2
                k=ithbuf(l) 
                k1 = (k-101)/5+1
                ii=ii+1
                IF(k<=nthvki)THEN
                  wa(ii)=fsav(k,i)
                ELSEIF(k>=100)THEN
C                 variables additionnelles AOUTi,BOUTi,UOUTi,MOUTi,HOUTi
                  kv=mod(k-101,5)+1
                  iv=id_vent
                  iadvent=ivolu(16,i)
                  nvent=ivolu(11,i)
                  IF (iv<=nvent) THEN 
                    wa(ii) =fsavvent(kv,iadvent+iv)
                  END IF
                END IF
              ENDDO
              iad2 = iad2 + 40
        ENDDO
        IF(ii>0)CALL wrtdes(wa,wa,ii,iform,1)
      ENDIF
C
      RETURN