bilan_sisyphe.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\bilan_sisyphe.f
00002 !
00074                      SUBROUTINE BILAN_SISYPHE
00075 !                    ************************
00076 !
00077      &(E,ESOMT,MESH,MSK,MASKEL,T1,T2,S,IELMU,VCUMU,DT,NPTFR,
00078      & INFO,ZFCL_C,ZFCL_S,ZFCL_MS,
00079      & QSCLXC,QSCLYC,NSICLA,VOLTOT,DZF_GF,MASS_GF,LGRAFED,
00080      & NUMLIQ,NFRLIQ,FLBCLA,VF,LT,NIT,NPOIN,VOLU2D,CSF_SABLE,MASDEP,
00081      & MASDEPT,CHARR,SUSP,SLIDE)
00082 !
00083 !***********************************************************************
00084 ! SISYPHE   V6P2                                   21/07/2011
00085 !***********************************************************************
00086 !
00087 !
00088 !
00089 !
00090 !
00091 !
00092 !
00093 !
00094 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00095 !| CSF_SABLE      |-->| 1-POROSITY
00096 !| DT             |-->| TIME STEP
00097 !| DZF_GF         |---| A SUPPRIMER
00098 !| E              |-->| BED EVOLUTION AT A GIVEN TIME STEP
00099 !| ESOMT          |-->| CUMULATED BED EVOLUTION
00100 !| FLBCLA         |-->| BLOCK OF FLUXES AT BOUNDARY FOR EACH CLASS
00101 !| IELMU          |-->| NUMER OF ELEMENTS
00102 !| INFO           |-->| IF YES : INFORMATION IS PRINTED
00103 !| LGRAFED        |---| A SUPPRIMER
00104 !| LT             |-->| CURRENT TIME STEP
00105 !| MASDEP         |-->| VOLUME DEPOSITED ON THE BOTTOM FOR EACH CLASS
00106 !|                |   | FROM THE BEGINNING
00107 !| MASDEPT        |-->| VOLUME DEPOSITED ON THE BOTTOM FOR EACH CLASS
00108 !|                |   | FOR THIS TIME STEP
00109 !| MASKEL         |-->| MASKING OF ELEMENTS
00110 !| MASS_GF        |---|A SUPPRIMER
00111 !| MESH           |<->| MESH STRUCTURE
00112 !| MSK            |-->| IF YES, THERE IS MASKED ELEMENTS
00113 !| NFRLIQ         |-->| NUMBER OF LIQUID BOUNDARIES
00114 !| NIT            |-->| NUMBER OF TIME STEPS
00115 !| NPOIN          |-->| NUMBER OF POINTS
00116 !| NPTFR          |-->| NUMBER OF BOUNDARY NODES
00117 !| NSICLA         |-->| NUMBER OF SIZE CLASSES FOR BED MATERIALS
00118 !| NUMLIQ         |-->| LIQUID BOUNDARY NUMBER OF BOUNDARY POINTS
00119 !| QSCLXC         |<->| TRANSPORT RATE FOR EACH CLASS X-DIRECTION
00120 !| QSCLYC         |<->| TRANSPORT RATE FOR EACH CLASS Y-DIRECTION
00121 !| S              |-->| VOID STRUCTURE
00122 !| T1             |<->| WORK BIEF_OBJ STRUCTURE
00123 !| T2             |<->| WORK BIEF_OBJ STRUCTURE
00124 !| VCUMU          |<->| VOLUME OF SEDIMENT ENTERING THE DOMAIN
00125 !| VF             |-->| IF YES : FINITE VOLUMES IF NO : FINITE ELEMENTS
00126 !| VOLTOT         |-->| VOLUME TOTAL PER CLASS OF SEDIMENT
00127 !| VOLU2D         |-->| INTEGRAL OF TEST FUNCTIONS (NOT ASSEMBLED IN //)
00128 !| ZFCL_C         |<->| BEDLOAD EVOLUTION FOR EACH SEDIMENT CLASS
00129 !| ZFCL_S         |<->| SUSPENDED LOAD EVOLUTION FOR EACH SEDIMENT CLASS
00130 !| ZFCL_MS        |<->| SLIDE EVOLUTION FOR EACH SEDIMENT CLASS
00131 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00132 !
00133       USE BIEF
00134       USE DECLARATIONS_SISYPHE, ONLY : NSICLM,MAXFRO
00135 !
00136       IMPLICIT NONE
00137       INTEGER LNG,LU
00138       COMMON/INFO/LNG,LU
00139 !
00140 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00141 !
00142       INTEGER, INTENT(IN)          :: NPTFR,NFRLIQ,IELMU,NSICLA,LT,NIT
00143       INTEGER, INTENT(IN)          :: NPOIN,NUMLIQ(NPTFR)
00144       DOUBLE PRECISION, INTENT(IN) :: DT
00145       LOGICAL, INTENT(IN)          :: MSK,INFO,VF,SUSP,SLIDE,CHARR
00146 !
00147       LOGICAL, INTENT(IN)          :: LGRAFED
00148       DOUBLE PRECISION, INTENT(INOUT) :: MASS_GF,VCUMU
00149       DOUBLE PRECISION, INTENT(IN)    :: CSF_SABLE,VOLTOT(NSICLA)
00150       DOUBLE PRECISION, INTENT(IN)    :: MASDEP(NSICLA)
00151       DOUBLE PRECISION, INTENT(INOUT) :: MASDEPT(NSICLA)
00152 !
00153 !-----------------------------------------------------------------------
00154 !
00155 !     VECTOR STRUCTURES
00156 !
00157       TYPE(BIEF_OBJ), INTENT(IN)    :: MASKEL,S,ZFCL_C,QSCLXC,QSCLYC
00158       TYPE(BIEF_OBJ), INTENT(IN)    :: E,ESOMT,DZF_GF,VOLU2D,ZFCL_S
00159       TYPE(BIEF_OBJ), INTENT(IN)    :: ZFCL_MS
00160       TYPE(BIEF_OBJ), INTENT(INOUT) :: T1,T2,FLBCLA
00161 !
00162 !-----------------------------------------------------------------------
00163 !
00164 !     MESH STRUCTURES
00165 !
00166       TYPE(BIEF_MESH) :: MESH
00167 !
00168 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00169 !
00170       INTEGER I,IFRLIQ,IPTFR,ICLA
00171       DOUBLE PRECISION RMASSE,RCUMU,RMASCLA(NSICLM)
00172       DOUBLE PRECISION VCUMUCLA(NSICLM),FLUXT,FLUXTCLA,VOLDEP
00173 !     DOUBLE PRECISION MASST
00174       DOUBLE PRECISION FLT_BOUND(MAXFRO),VOLDEPC
00175 !
00176       DOUBLE PRECISION P_DSUM
00177       EXTERNAL         P_DSUM
00178 !
00179 !-----------------------------------------------------------------------
00180 !
00181 !     COMPUTES THE EVOLUTION (E)
00182 !
00183       RMASSE=0.D0
00184       DO I=1,NPOIN
00185         RMASSE=RMASSE+E%R(I)*VOLU2D%R(I)
00186       ENDDO
00187       IF(NCSIZE.GT.1) RMASSE = P_DSUM(RMASSE)
00188 !
00189 !=======================================================================
00190 !
00191 !     COMPUTES THE INTEGRAL OF EVOLUTION AT THE END (ESOMT)
00192 !
00193       IF(LT.EQ.NIT) THEN
00194         RCUMU=0.D0
00195         DO I=1,NPOIN
00196           RCUMU=RCUMU+ESOMT%R(I)*VOLU2D%R(I)
00197         ENDDO
00198         IF(NCSIZE.GT.1) RCUMU = P_DSUM(RCUMU)
00199       ENDIF
00200 !
00201 !=======================================================================
00202 !
00203 !     COMPUTES THE FLUXES AT THE BOUNDARIES
00204 !
00205       IF(CHARR) THEN
00206         CALL OS('X=Y     ',X=T1,Y=FLBCLA%ADR(1)%P)
00207         IF(NSICLA.GT.1) THEN
00208           DO I=2,NSICLA
00209             CALL OS('X=X+Y   ',X=T1,Y=FLBCLA%ADR(I)%P)
00210           ENDDO
00211         ENDIF
00212       ELSE
00213         CALL CPSTVC(FLBCLA%ADR(1)%P,T1)
00214         CALL OS('X=0     ',X=T1)
00215       ENDIF
00216 !
00217       FLUXT=0.D0
00218 !
00219       IF(NFRLIQ.GT.0) THEN
00220         DO IFRLIQ=1,NFRLIQ
00221           FLT_BOUND(IFRLIQ)=0.D0
00222         ENDDO
00223         IF(NPTFR.GT.0) THEN
00224           DO IPTFR=1,NPTFR
00225             IFRLIQ=NUMLIQ(IPTFR)
00226             IF(IFRLIQ.GT.0) THEN
00227               FLT_BOUND(IFRLIQ)=FLT_BOUND(IFRLIQ)+T1%R(IPTFR)
00228             ENDIF
00229           ENDDO
00230         ENDIF
00231         IF(NCSIZE.GT.1) THEN
00232           DO IFRLIQ=1,NFRLIQ
00233             FLT_BOUND(IFRLIQ)=P_DSUM(FLT_BOUND(IFRLIQ))
00234           ENDDO
00235         ENDIF
00236         DO IFRLIQ=1,NFRLIQ
00237           FLUXT=FLUXT+FLT_BOUND(IFRLIQ)
00238         ENDDO
00239       ENDIF
00240 !
00241       VCUMU = VCUMU - FLUXT*DT/CSF_SABLE
00242 !
00243 !     BALANCE IN EXTENDED GRANULOMETRY
00244 !
00245       IF(NSICLA.GT.1) THEN
00246 !
00247         DO ICLA=1,NSICLA
00248 !
00249 !       COMPUTES THE EVOLUTION PER CLASS
00250 !
00251         RMASCLA(ICLA)=0.D0
00252         IF(SUSP.AND.SLIDE.AND.CHARR) THEN
00253           DO I=1,NPOIN
00254             RMASCLA(ICLA)=RMASCLA(ICLA)
00255      &                   +( ZFCL_C%ADR(ICLA)%P%R(I)
00256      &                     +ZFCL_S%ADR(ICLA)%P%R(I)
00257      &                     +ZFCL_MS%ADR(ICLA)%P%R(I) )*VOLU2D%R(I)
00258           ENDDO
00259         ELSEIF(SLIDE.AND.CHARR) THEN
00260           DO I=1,NPOIN
00261             RMASCLA(ICLA)=RMASCLA(ICLA)
00262      &                   +( ZFCL_C%ADR(ICLA)%P%R(I)
00263      &                     +ZFCL_MS%ADR(ICLA)%P%R(I) )*VOLU2D%R(I)
00264           ENDDO
00265         ELSEIF(SUSP.AND.CHARR) THEN
00266           DO I=1,NPOIN
00267             RMASCLA(ICLA)=RMASCLA(ICLA)
00268      &                   +( ZFCL_C%ADR(ICLA)%P%R(I)
00269      &                     +ZFCL_S%ADR(ICLA)%P%R(I) )*VOLU2D%R(I)
00270           ENDDO
00271         ELSEIF(SUSP.AND.SLIDE) THEN
00272           DO I=1,NPOIN
00273             RMASCLA(ICLA)=RMASCLA(ICLA)
00274      &                   +( ZFCL_S%ADR(ICLA)%P%R(I)
00275      &                     +ZFCL_MS%ADR(ICLA)%P%R(I) )*VOLU2D%R(I)
00276           ENDDO
00277         ELSEIF(SUSP) THEN
00278           DO I=1,NPOIN
00279             RMASCLA(ICLA)=RMASCLA(ICLA)
00280      &                   +ZFCL_S%ADR(ICLA)%P%R(I)*VOLU2D%R(I)
00281           ENDDO
00282         ELSEIF(SLIDE) THEN
00283           DO I=1,NPOIN
00284             RMASCLA(ICLA)=RMASCLA(ICLA)
00285      &                   +ZFCL_MS%ADR(ICLA)%P%R(I)*VOLU2D%R(I)
00286           ENDDO
00287         ELSEIF(CHARR) THEN
00288           DO I=1,NPOIN
00289             RMASCLA(ICLA)=RMASCLA(ICLA)
00290      &                   +ZFCL_C%ADR(ICLA)%P%R(I)*VOLU2D%R(I)
00291           ENDDO
00292         ENDIF
00293         IF(NCSIZE.GT.1) RMASCLA(ICLA) = P_DSUM(RMASCLA(ICLA))
00294 !
00295 !       COMPUTES THE FREE FLUXES BY CLASS
00296 !
00297         FLUXTCLA=0.D0
00298         IF(NFRLIQ.GT.0.AND.CHARR) THEN
00299           IF(NPTFR.GT.0) THEN
00300             DO IPTFR=1,NPTFR
00301               IFRLIQ=NUMLIQ(IPTFR)
00302               IF(IFRLIQ.GT.0) THEN
00303                 FLUXTCLA=FLUXTCLA+FLBCLA%ADR(ICLA)%P%R(IPTFR)
00304               ENDIF
00305             ENDDO
00306           ENDIF
00307           IF(NCSIZE.GT.1) FLUXTCLA=P_DSUM(FLUXTCLA)
00308         ENDIF
00309 !
00310         VCUMUCLA(ICLA) = - FLUXTCLA*DT/CSF_SABLE
00311 !
00312         ENDDO
00313 !
00314       ENDIF
00315 !
00316 !=======================================================================
00317 !
00318 !     GRAIN-FEEDING
00319 !     IF (LGRAFED) THEN
00320 !        CALL VECTOR(T1,'=','MASVEC          ',IELMU,
00321 !    &               1.D0,DZF_GF,S,S,S,S,S,MESH,MSK,MASKEL)
00322 !        MASST = BIEF_SUM(T1)
00323 !        IF(NCSIZE.GT.1) MASST = P_DSUM(MASST)
00324 !        MASS_GF = MASS_GF + MASST
00325 !     ENDIF
00326 !     IF(DREDGESIM) ...   ?????
00327 !
00328       VOLDEPC=0.D0
00329       IF(SUSP) THEN
00330         DO I=1,NSICLA
00331           VOLDEPC=VOLDEPC+MASDEPT(I)
00332         ENDDO
00333         VOLDEPC=VOLDEPC/CSF_SABLE
00334       ENDIF
00335 !
00336 !     WRITES OUT THE BALANCE
00337 !
00338       IF(INFO) THEN
00339 !
00340 !       GLOBAL BALANCE
00341 !
00342         WRITE(LU,*)
00343         IF(LNG.EQ.1) THEN
00344           WRITE(LU,1000)
00345           WRITE(LU,1010) RMASSE
00346           IF(NFRLIQ.GT.0) THEN
00347             DO IFRLIQ=1,NFRLIQ
00348               WRITE(LU,1110) IFRLIQ,-FLT_BOUND(IFRLIQ)/CSF_SABLE
00349             ENDDO
00350             WRITE(LU,1111) -FLUXT/CSF_SABLE
00351             IF(SUSP) WRITE(LU,1112) VOLDEPC
00352           ENDIF
00353           WRITE(LU,1033) RMASSE+DT*FLUXT/CSF_SABLE-VOLDEPC
00354         ELSEIF(LNG.EQ.2) THEN
00355           WRITE(LU,2000)
00356           WRITE(LU,2010) RMASSE
00357           IF(NFRLIQ.GT.0) THEN
00358             DO IFRLIQ=1,NFRLIQ
00359               WRITE(LU,2110) IFRLIQ,-FLT_BOUND(IFRLIQ)/CSF_SABLE
00360             ENDDO
00361             WRITE(LU,2111) -FLUXT/CSF_SABLE
00362             IF(SUSP) WRITE(LU,2112) VOLDEPC
00363           ENDIF
00364           WRITE(LU,2033) RMASSE+DT*FLUXT/CSF_SABLE-VOLDEPC
00365         ENDIF
00366 !
00367 !       BALANCE PER CLASS
00368 !
00369         IF(NSICLA.GT.1) THEN
00370           DO I=1,NSICLA
00371             WRITE(LU,*)
00372             IF(LNG.EQ.1) THEN
00373               WRITE(LU,*) 'BILAN POUR LA CLASSE DE SEDIMENT :',I
00374               WRITE(LU,*) 'VOLUME TOTAL DE LA CLASSE :',VOLTOT(I)
00375               WRITE(LU,3011) RMASCLA(I)
00376               WRITE(LU,3032) VCUMUCLA(I)
00377               IF(SUSP) THEN
00378                 WRITE(LU,3033) MASDEPT(I)/CSF_SABLE
00379                 WRITE(LU,1033) RMASCLA(I)-VCUMUCLA(I)
00380      &                                   -MASDEPT(I)/CSF_SABLE
00381               ELSE
00382                 WRITE(LU,1033) RMASCLA(I)-VCUMUCLA(I)
00383               ENDIF
00384             ELSEIF(LNG.EQ.2) THEN
00385               WRITE(LU,*) 'MASS BALANCE FOR SEDIMENT CLASS :',I
00386               WRITE(LU,*) 'TOTAL VOLUME:',VOLTOT(I)
00387               WRITE(LU,3010) RMASCLA(I)
00388               WRITE(LU,3031) VCUMUCLA(I)
00389               IF(SUSP) THEN
00390                 WRITE(LU,3034) MASDEPT(I)/CSF_SABLE
00391                 WRITE(LU,1033) RMASCLA(I)-VCUMUCLA(I)
00392      &                                   -MASDEPT(I)/CSF_SABLE
00393               ELSE
00394                 WRITE(LU,2033) RMASCLA(I)-VCUMUCLA(I)
00395               ENDIF
00396             ENDIF
00397           ENDDO
00398         ENDIF
00399 !
00400 !       FINAL GLOBAL BALANCE
00401 !
00402         IF(LT.EQ.NIT) THEN
00403           VOLDEP=0.D0
00404           DO I=1,NSICLA
00405             VOLDEP=VOLDEP+MASDEP(I)
00406           ENDDO
00407           VOLDEP=VOLDEP/CSF_SABLE
00408           WRITE(LU,*)
00409           IF(LNG.EQ.1) THEN
00410             WRITE(LU,1030) RCUMU
00411             WRITE(LU,1031) VCUMU
00412             IF(SUSP) WRITE(LU,1032) VOLDEP
00413             WRITE(LU,1033) RCUMU-VCUMU-VOLDEP
00414           ELSEIF(LNG.EQ.2) THEN
00415             WRITE(LU,2030) RCUMU
00416             WRITE(LU,2031) VCUMU
00417             IF(SUSP) WRITE(LU,2032) VOLDEP
00418             WRITE(LU,2033) RCUMU-VCUMU-VOLDEP
00419           ENDIF
00420         ENDIF
00421 !
00422 !       IF(LGRAFED) THEN
00423 !         IF(LNG.EQ.1) THEN
00424 !           WRITE(LU, 4000) MASST
00425 !           WRITE(LU, 4010) MASS_GF
00426 !         ENDIF
00427 !         IF(LNG.EQ.2) THEN
00428 !           WRITE(LU, 4001) MASST
00429 !           WRITE(LU, 4011) MASS_GF
00430 !         ENDIF
00431 !       ENDIF
00432 !
00433       ENDIF
00434 !
00435 1000  FORMAT(1X,'BILAN DE MASSE (EN VOLUME, VIDES INCLUS) : ')
00436 1010  FORMAT(1X,'SOMME DES EVOLUTIONS : ',G16.7,' M3')
00437 !1020  FORMAT(1X,'FLUX IMPOSE          : ', G16.7,' M3/S'
00438 !     &         ,'  ( M3/S  >0 = ENTRANT )')
00439 !1021  FORMAT(1X,'FLUX LIBRE           : ', G16.7,' M3/S'
00440 !     &         ,'  ( M3/S  >0 = ENTRANT )')
00441 1030  FORMAT(1X,'SOMME DES EVOLUTIONS CUMULEES : ',G16.7)
00442 1031  FORMAT(1X,'VOLUME ENTRE AUX FRONTIERES   : ',G16.7,' M3'
00443      &         ,'  ( SI <0 VOLUME SORTI )')
00444 1032  FORMAT(1X,'VOLUME DEPOSE SUR LE FOND     : ',G16.7,' M3'
00445      &         ,'  ( SI <0 VOLUME ERODE )')
00446 1033  FORMAT(1X,'VOLUME PERDU                  : ',G16.7,' M3'
00447      &         ,'  ( SI <0 VOLUME SORTI )')
00448 1110  FORMAT(1X,'FRONTIERE ',1I3,' FLUX EN CHARRIAGE = ',G16.7,
00449      &          ' M3/S >0=ENTRANT, AVEC VIDES')
00450 1111  FORMAT(1X,'TOTAL         FLUX EN CHARRIAGE = ',G16.7,
00451      &          ' M3/S >0=ENTRANT, AVEC VIDES')
00452 1112  FORMAT(1X,'DEPOT DE SUSPENSION SUR LE FOND = ',G16.7,
00453      &          ' ( M3/S )')
00454 2000  FORMAT(1X,'MASS-BALANCE (IN VOLUME, INCLUDING VOID): ')
00455 2010  FORMAT(1X,'SUM OF THE EVOLUTIONS : ',G16.7,' M3')
00456 !2020  FORMAT(1X,'PRESCRIBED FLOW       : ',G16.7,' M3/S'
00457 !     &         ,'  M3/S >0=ENTERING, WITH VOIDS')
00458 !2021  FORMAT(1X,'FREE FLOW             : ',G16.7,' M3/S'
00459 !     &         ,'  M3/S >0=ENTERING, WITH VOIDS')
00460 2030  FORMAT(1X,'SUM OF THE CUMULATED EVOLUTIONS : ',G16.7)
00461 2031  FORMAT(1X,'VOLUME THAT ENTERED THE DOMAIN  : ',G16.7,' M3'
00462      &         ,'  ( IF <0 EXIT )')
00463 2032  FORMAT(1X,'VOLUME DEPOSITED ON THE BOTTOM  : ',G16.7,' M3'
00464      &         ,'  ( IF <0 ERODED )')
00465 2033  FORMAT(1X,'LOST VOLUME                     : ',G16.7,' M3'
00466      &         ,'  ( IF <0 EXIT )')
00467 2110  FORMAT(1X,'BOUNDARY ',1I3,' BEDLOAD FLUX = ',G16.7,
00468      &          '  ( M3/S  >0 = ENTERING )')
00469 2111  FORMAT(1X,'TOTAL        BEDLOAD FLUX = ',G16.7,
00470      &          '  ( M3/S  >0 = ENTERING )')
00471 2112  FORMAT(1X,'DEPOSIT ON BOTTOM         = ',G16.7,
00472      &          '  ( M3/S )')
00473 3010  FORMAT(1X,'SUM OF THE EVOLUTIONS FOR THIS CLASS: ',G16.7)
00474 3031  FORMAT(1X,'VOLUME THAT ENTERED THE DOMAIN FOR THIS CLASS: '
00475      &       ,G16.7,' M3')
00476 3034  FORMAT(1X,'VOLUME DEPOSITED ON BOTTOM FOR THIS CLASS:     '
00477      &       ,G16.7,' M3')
00478 3011  FORMAT(1X,'SOMME DES EVOLUTIONS POUR CETTE CLASSE : ',G16.7)
00479 3032  FORMAT(1X,'VOLUME ENTRE DANS LE DOMAINE POUR CETTE CLASSE : '
00480      &       ,G16.7,' M3')
00481 3033  FORMAT(1X,'VOLUME DEPOSE SUR LE FOND POUR CETTE CLASSE    : '
00482      &       ,G16.7,' M3')
00483 !4000  FORMAT(1X,'GRAIN-FEEDING A CET INSTANT       : ',G16.7)
00484 !4010  FORMAT(1X,'GRAIN-FEEDING JUSQU''A MAINTENANT : ',G16.7)
00485 !4001  FORMAT(1X,'GRAIN-FEEDING THIS MOMENT : ',G16.7)
00486 !4011  FORMAT(1X,'GRAIN-FEEDING UNTIL NOW   : ',G16.7)
00487 !
00488 !-----------------------------------------------------------------------
00489 !
00490       RETURN
00491       END

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