condin_adj.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\condin_adj.f
00002 !
00089                      SUBROUTINE CONDIN_ADJ
00090 !                    *********************
00091 !
00092      &(ALIRE,NRES,TROUVE)
00093 !
00094 !***********************************************************************
00095 ! TELEMAC2D   V6P1                                   21/08/2010
00096 !***********************************************************************
00097 !
00098 !
00099 !
00100 !
00101 !
00102 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00103 !| ALIRE          |-->| ARRAY FOR ALL VARIABLES 0:VARIABLE DISCARDED
00104 !|                |   |                         1:VARIABLE TO BE READ
00105 !| NRES           |-->| LOGICAL UNIT OF RESULTS FILE
00106 !| TROUVE         |<--| ARRAY FOR ALL VARIABLES 0:VARIABLE NOT FOUND
00107 !|                |   |                         1:VARIABLE FOUND
00108 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00109 !
00110       USE BIEF
00111       USE DECLARATIONS_TELEMAC2D
00112 !
00113       IMPLICIT NONE
00114       INTEGER LNG,LU
00115       COMMON/INFO/LNG,LU
00116 !
00117 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00118 !
00119       INTEGER, INTENT(IN)    :: ALIRE(*),NRES
00120       INTEGER, INTENT(INOUT) :: TROUVE(*)
00121 !
00122 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00123 !
00124       INTEGER I,ITER
00125       DOUBLE PRECISION HIST(1),AT1
00126 !
00127 !-----------------------------------------------------------------------
00128 !
00129 !     CONDIN FOR ADJOINT PROBLEM: T=T(N+1) : P=0,Q=0,R=0
00130 !
00131       CALL OS( 'X=C     ',PP,PP,PP,0.D0)
00132       CALL OS( 'X=C     ',QQ,QQ,QQ,0.D0)
00133       CALL OS( 'X=C     ',RR,RR,RR,0.D0)
00134 !
00135 !     JUST IN CASE CV1,.. IS WRITTEN IN THE RESULT FILE
00136 !
00137       CALL OS( 'X=C     ',CV1,CV1,CV1,0.D0)
00138       CALL OS( 'X=C     ',CV2,CV2,CV2,0.D0)
00139       CALL OS( 'X=C     ',CV3,CV3,CV3,0.D0)
00140 !
00141 !     READS THE LAST TIME IN THE TELEMAC RESULT FILE (NRES)
00142 !     INITIALISES U,V AND H
00143 !
00144       REWIND NRES
00145 !
00146       CALL BIEF_SUITE(VARSOR,VARCL,ITER,NRES,'SERAFIN ',
00147      &           HIST,0,NPOIN,AT,TEXTE,VARCLA,
00148      &           NVARCL,TROUVE,ALIRE,LISTIN,.TRUE.,MAXVAR)
00149 !
00150 !     GIVES MEASUREMENTS HD,UD AND VD AT THE LAST TIME STEP
00151 !     (ITER AND AT GIVEN BY THE PREVIOUS CALL TO SUITE)
00152 !
00153       CALL MESURES(ITER,AT)
00154 !     INITIALISES HH, UU, VV
00155 !
00156       CALL OS( 'X=Y     ' , HH   , H , H , 0.D0 )
00157       CALL OS( 'X=Y     ' , UU   , U , U , 0.D0 )
00158       CALL OS( 'X=Y     ' , VV   , V , V , 0.D0 )
00159       CALL OS( 'X=C     ' , HIT1 , HIT1 , HIT1 , 0.D0 )
00160       CALL OS( 'X=C     ' , UIT1 , UIT1 , UIT1 , 0.D0 )
00161       CALL OS( 'X=C     ' , VIT1 , VIT1 , VIT1 , 0.D0 )
00162 !
00163 !     READS TELEMAC2D RESULTS (RESULT FILE - UNIT NRES)
00164 !     THIS IS TO HAVE UN AT THE LAST TIME STEP INTO U.
00165 !
00166 !     BEWARE : ASSUMES THAT NVARRES HAS ALREADY BEEN COMPUTED
00167 !
00168       DO I=1,2*(NVARRES+1)
00169         BACKSPACE NRES
00170       ENDDO
00171       CALL LITENR(VARSOR,VARCL,NRES,'STD',HIST,0,NPOIN,AT1,TEXTE,
00172      &           TEXRES,NVARRES,VARCLA,0,TROUVE,ALIRE,W,.FALSE.,MAXVAR)
00173 !
00174 !-----------------------------------------------------------------------
00175 !
00176       AT = AT + DT
00177 !
00178 !-----------------------------------------------------------------------
00179 !
00180       RETURN
00181       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0