bil3d.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\bil3d.f
00002 !
00131                      SUBROUTINE BIL3D
00132 !                    ****************
00133 !
00134      &(LT,IKLBORL,NPTFR,NETAG,NELEM)
00135 !
00136 !***********************************************************************
00137 ! TELEMAC3D   V7P0                                   21/08/2010
00138 !***********************************************************************
00139 !
00140 !
00141 !
00142 !
00143 !
00144 !
00145 !
00146 !
00147 !
00148 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00149 !| IKLBORL        |-->| CONNECTIVITY TABLE OF LATERAL BOUNDARIES
00150 !| LT             |-->| CURRENT TIME STEP NUMBER
00151 !| NELEM          |-->| NUMBER OF ELEMENTS IN THE 2D MESH
00152 !| NETAG          |-->| NUMBER OF LAYERS (I.E. NPLAN-1)
00153 !| NPTFR          |-->| NUMBER OF 2D BOUNDARY POINTS
00154 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00155 !
00156       USE BIEF
00157       USE INTERFACE_TELEMAC3D, EX_BIL3D => BIL3D
00158       USE DECLARATIONS_TELEMAC
00159       USE DECLARATIONS_TELEMAC3D
00160 !
00161       IMPLICIT NONE
00162       INTEGER LNG,LU
00163       COMMON/INFO/LNG,LU
00164 !
00165 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00166 !
00167       INTEGER, INTENT(IN) :: LT,NPTFR,NETAG,NELEM
00168       INTEGER, INTENT(IN) :: IKLBORL(NELEBX,4)
00169 !
00170 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00171 !
00172       DOUBLE PRECISION FLUTOT,SUR12,A1,A2,FLURAIN,D,Z_1,Z_2,REALRAIN
00173 !
00174       INTEGER I,L1,L2,L3,L4,N1,N2,N3,N4,IVBIL,ILIQ,IELEB
00175       INTEGER IPTFR,IETAGE,ITRAC
00176 !                            25=MAXTRA+5
00177       DOUBLE PRECISION FLUDI(25),FLUS1(25),FLUXTOTAL
00178       INTRINSIC SQRT
00179 !
00180       DOUBLE PRECISION P_DSUM
00181       EXTERNAL         P_DSUM
00182 !
00183 !=======================================================================
00184 !
00185 !     COMPUTES THE FLUXES AT THE LIQUID BOUNDARIES
00186 !
00187       IF(NFRLIQ.GT.0) THEN
00188         DO I=1,NFRLIQ
00189           FLUX_BOUNDARIES(I)=0.D0
00190         ENDDO
00191         IF(NPTFR.GT.0) THEN
00192           DO I=1,NPTFR
00193             ILIQ=NUMLIQ%I(I)
00194             IF(ILIQ.GT.0) THEN
00195               FLUX_BOUNDARIES(ILIQ)=FLUX_BOUNDARIES(ILIQ)+FLBOR%R(I)
00196             ENDIF
00197           ENDDO
00198         ENDIF
00199         IF(NCSIZE.GT.1) THEN
00200           DO I=1,NFRLIQ
00201             FLUX_BOUNDARIES(I)=P_DSUM(FLUX_BOUNDARIES(I))
00202           ENDDO
00203         ENDIF
00204       ENDIF
00205 !
00206 !=======================================================================
00207 !
00208 !   COMPUTES THE ADVECTION FLUXES ON THE VARIOUS OPEN BOUNDARIES
00209 !
00210 !=======================================================================
00211 !
00212 !     TOTAL FLUX, STARTING WITH THE LIQUID BOUNDARIES
00213 !
00214       FLUXTOTAL=0.D0
00215       IF(NFRLIQ.GT.0) THEN
00216         DO I=1,NFRLIQ
00217           FLUXTOTAL=FLUXTOTAL+FLUX_BOUNDARIES(I)
00218         ENDDO
00219       ENDIF
00220 !
00221 !     RAIN AND EVAPORATION
00222 !
00223       IF(RAIN) THEN
00224         FLURAIN=0.D0
00225         DO I=1,NPOIN2
00226           FLURAIN=FLURAIN+PLUIE%R(I)
00227         ENDDO
00228         IF(NCSIZE.GT.1) FLURAIN = P_DSUM(FLURAIN)
00229         FLUXTOTAL=FLUXTOTAL-FLURAIN
00230       ENDIF
00231 !
00232       IF(NSCE.GE.1) THEN
00233         DO I=1,NSCE
00234           FLUXTOTAL = FLUXTOTAL - QSCE2(I)
00235         ENDDO
00236 !       IF(NTRAC.GT.0) THEN
00237 !         DO I=1,NSCE
00238 !         DO IVBIL=2,1+NTRAC
00239 !         NOW DONE IN CVDF3D AND MURD3D
00240 !                                                      TASCE(MAXSCE,NTRAC)
00241 !             FLUX%R(IVBIL) = FLUX%R(IVBIL) - QSCE2(I)*TASCE2(I,IVBIL-1)
00242 !         ENDDO
00243 !         ENDDO
00244 !       ENDIF
00245       ENDIF
00246 !
00247 !=======================================================================
00248 !
00249 !  COMPUTES THE FLUXES BY DIFFUSION AND SOURCES OF TRACERS
00250 !
00251 !=======================================================================
00252 !
00253       SUR12 = 1.D0/12.D0
00254 !
00255 !   ====================================
00256 !   FLUX BY SOURCES OF TRACERS
00257 !   ====================================
00258 !
00259       IF(NTRAC.NE.0) THEN
00260         DO ITRAC=1,NTRAC
00261           FLUS1(5+ITRAC) = 0.D0
00262           IF(S1TA%ADR(ITRAC)%P%TYPR.NE.'0') THEN
00263             DO I=1,NPOIN3
00264 !             JMH NOTE : IS NOT VOLU WITH SUPG
00265 !                        BUT A VALUE WEIGHTED WITH VOLUN
00266 !                        INVESTIGATE WHEN THE CASE ARISES...
00267               FLUS1(5+ITRAC)=FLUS1(5+ITRAC)
00268      &                        +S1TA%ADR(ITRAC)%P%R(I)*VOLU%R(I)*
00269      &                           TA%ADR(ITRAC)%P%R(I)
00270             ENDDO
00271             IF(NCSIZE.GT.1) THEN
00272               FLUS1(5+ITRAC) = P_DSUM(FLUS1(5+ITRAC))
00273             ENDIF
00274           ENDIF
00275         ENDDO
00276       ENDIF
00277 !
00278 !   ===============================
00279 !   FLUX BY DIFFUSION OF TRACERS
00280 !   ===============================
00281 !
00282 !   IN THE CASE OF A SEDIMENT, THIS FLUX INCLUDES THE SETTLING VELOCITY
00283 !   SEE HOW ATABOF AND BTABOF ARE DONE IN FLUSED
00284 !
00285       IF(NTRAC.NE.0) THEN
00286 !
00287         DO ITRAC=1,NTRAC
00288 !
00289           FLUDI(5+ITRAC) = 0.D0
00290 !
00291 !         BOTTOM AND FREE SURFACE
00292 !
00293           IF(ATABOF%ADR(ITRAC)%P%TYPR.NE.'0') THEN
00294 !           WITH MASS-LUMPING LIKE IN DIFF3D
00295 !           CORRECTION CV+JMH 28/10/2013  SEE ALSO DIFF3D
00296             IF(SIGMAG.OR.OPTBAN.EQ.1) THEN
00297               DO I=1,NPOIN2
00298 !               1) FLUX THROUGH THE BOTTOM NOT TAKEN INTO ACCOUNT FOR TIDAL FLATS
00299 !               2) TRACER AT ACTUAL BOTTOM PLANE TAKEN INTO ACCOUNT, THEN IPBOT
00300 !                  IN ADDRESS OF TA.
00301                 IF(IPBOT%I(I).NE.NPLAN-1) THEN
00302                   FLUDI(5+ITRAC) = FLUDI(5+ITRAC)
00303      &            + ATABOF%ADR(ITRAC)%P%R(I)*VOLU2D%R(I)
00304      &            *TA%ADR(ITRAC)%P%R(IPBOT%I(I)*NPOIN2+I)
00305                 ENDIF
00306               ENDDO
00307              ELSE
00308               DO I=1,NPOIN2
00309                 FLUDI(5+ITRAC) = FLUDI(5+ITRAC)
00310      &          + ATABOF%ADR(ITRAC)%P%R(I)*VOLU2D%R(I)
00311      &                             *TA%ADR(ITRAC)%P%R(I)
00312               ENDDO
00313              ENDIF
00314           ENDIF
00315 !
00316           IF(ATABOS%ADR(ITRAC)%P%TYPR.NE.'0') THEN
00317 !           WITH MASS-LUMPING LIKE IN DIFF3D
00318             DO I=1,NPOIN2
00319               FLUDI(5+ITRAC) = FLUDI(5+ITRAC)
00320      &        + ATABOS%ADR(ITRAC)%P%R(I)*VOLU2D%R(I)
00321      &        * TA%ADR(ITRAC)%P%R(I+NETAG*NPOIN2)
00322             ENDDO
00323           ENDIF
00324 !
00325           IF(BTABOF%ADR(ITRAC)%P%TYPR.NE.'0') THEN
00326             DO I=1,NPOIN2
00327               FLUDI(5+ITRAC) = FLUDI(5+ITRAC)
00328      &        + VOLU2D%R(I)*BTABOF%ADR(ITRAC)%P%R(I)
00329             ENDDO
00330           ENDIF
00331 !
00332           IF(BTABOS%ADR(ITRAC)%P%TYPR.NE.'0') THEN
00333             DO I=1,NPOIN2
00334               FLUDI(5+ITRAC) = FLUDI(5+ITRAC)
00335      &        +VOLU2D%R(I)*BTABOS%ADR(ITRAC)%P%R(I)
00336             ENDDO
00337           ENDIF
00338 !
00339 !         VALUE OF TRACER IN RAIN
00340 !
00341           IF(RAIN) THEN
00342             REALRAIN=0.D0
00343             DO I=1,NPOIN2
00344               REALRAIN=REALRAIN+MAX(PLUIE%R(I),0.D0)
00345             ENDDO
00346             FLUDI(5+ITRAC)=FLUDI(5+ITRAC)+REALRAIN*TRAIN(ITRAC)
00347           ENDIF
00348 !
00349 !         LATERAL BOUNDARIES
00350 !
00351           DO IETAGE=1,NETAG
00352 !
00353             IF(ATABOL%ADR(ITRAC)%P%TYPR.NE.'0') THEN
00354               DO IPTFR=1,NPTFR
00355 !
00356                 IELEB=(IETAGE-1)*NPTFR+IPTFR
00357                 L1 = IKLBORL(IELEB,1)
00358                 L2 = IKLBORL(IELEB,2)
00359                 L3 = IKLBORL(IELEB,3)
00360                 L4 = IKLBORL(IELEB,4)
00361                 N1 = NBOR3%I(L1)
00362                 N2 = NBOR3%I(L2)
00363                 N3 = NBOR3%I(L3)
00364                 N4 = NBOR3%I(L4)
00365 !
00366                 FLUDI(5+ITRAC) = FLUDI(5+ITRAC) + SUR12 *
00367      &          ( ATABOL%ADR(ITRAC)%P%R(L1)*TA%ADR(ITRAC)%P%R(N1)
00368      &          + ATABOL%ADR(ITRAC)%P%R(L2)*TA%ADR(ITRAC)%P%R(N2)
00369      &          + ATABOL%ADR(ITRAC)%P%R(L3)*TA%ADR(ITRAC)%P%R(N3)
00370      &          + ATABOL%ADR(ITRAC)%P%R(L4)*TA%ADR(ITRAC)%P%R(N4) )
00371 !
00372               ENDDO
00373             ENDIF
00374 !
00375             IF(BTABOL%ADR(ITRAC)%P%TYPR.NE.'0') THEN
00376               DO IPTFR=1,NPTFR
00377 !
00378                 IELEB=(IETAGE-1)*NPTFR+IPTFR
00379                 L1 = IKLBORL(IELEB,1)
00380                 L2 = IKLBORL(IELEB,2)
00381                 L3 = IKLBORL(IELEB,3)
00382                 L4 = IKLBORL(IELEB,4)
00383                 N1 = NBOR3%I(L1)
00384                 N2 = NBOR3%I(L2)
00385                 N3 = NBOR3%I(L3)
00386                 N4 = NBOR3%I(L4)
00387                 D = SQRT((X(N2)-X(N1))**2 + (Y(N2)-Y(N1))**2)
00388                 Z_1 = ZPROP%R(N4) - ZPROP%R(N1)
00389                 Z_2 = ZPROP%R(N3) - ZPROP%R(N2)
00390                 A1 = D * (Z_1+Z_1+Z_2)
00391                 A2 = D * (Z_2+Z_2+Z_1)
00392 !
00393                 FLUDI(5+ITRAC) = FLUDI(5+ITRAC) + SUR12 *
00394      &          ( BTABOL%ADR(ITRAC)%P%R(L1)*A1
00395      &           +BTABOL%ADR(ITRAC)%P%R(L2)*A2
00396      &           +BTABOL%ADR(ITRAC)%P%R(L3)*A2
00397      &           +BTABOL%ADR(ITRAC)%P%R(L4)*A1)
00398 !
00399               ENDDO
00400             ENDIF
00401 !
00402           ENDDO
00403 !
00404           IF(NCSIZE.GT.1) FLUDI(5+ITRAC) = P_DSUM(FLUDI(5+ITRAC))
00405 !
00406         ENDDO
00407 !
00408       ENDIF
00409 !
00410 !=======================================================================
00411 !
00412 ! COMPUTES CUMULATIVE FLUXES AND WRITES OUT
00413 !
00414 !=======================================================================
00415 !
00416       FLUXTOTCUM = FLUXTOTCUM + FLUXTOTAL
00417 !
00418       IF(INFOGR) THEN
00419         IF(LNG.EQ.1) THEN
00420           WRITE(LU,601) MASSEN_WATER,MASSE_WATER,DT*FLUXTOTAL,
00421      &                  MASSEN_WATER-MASSE_WATER-DT*FLUXTOTAL
00422         ENDIF
00423         IF(LNG.EQ.2) THEN
00424           WRITE(LU,602) MASSEN_WATER,MASSE_WATER,DT*FLUXTOTAL,
00425      &                  MASSEN_WATER-MASSE_WATER-DT*FLUXTOTAL
00426         ENDIF
00427         IF(RAIN) THEN
00428           IF(LNG.EQ.1) WRITE(LU,603) DT*FLURAIN
00429           IF(LNG.EQ.2) WRITE(LU,604) DT*FLURAIN
00430         ENDIF
00431         IF(NFRLIQ.GT.0) THEN
00432           DO I=1,NFRLIQ
00433             IF(LNG.EQ.1) WRITE(LU,3020) I,-FLUX_BOUNDARIES(I)
00434             IF(LNG.EQ.2) WRITE(LU,4020) I,-FLUX_BOUNDARIES(I)
00435           ENDDO
00436         ENDIF
00437       ENDIF
00438 !
00439 !-----------------------------------------------------------------------
00440 !
00441       IF(NTRAC.GT.0) THEN
00442 !
00443         DO IVBIL=6,5+NTRAC
00444 !
00445           FLUTOT = FLUX%R(IVBIL) - FLUDI(IVBIL) + FLUS1(IVBIL)
00446           FLUCUM%R(IVBIL) = FLUCUM%R(IVBIL) + FLUTOT
00447 !
00448           IF(INFOGR) THEN
00449             IF(SEDI.AND.(IVBIL.EQ.NTRAC+5)) THEN
00450               IF(LNG.EQ.1) WRITE(LU,611) FLUX%R(IVBIL),
00451      &        -FLUDI(IVBIL),MASSEN%R(IVBIL),MASSE%R(IVBIL),DT*FLUTOT,
00452      &                      MASSEN%R(IVBIL)-MASSE%R(IVBIL)-DT*FLUTOT
00453               IF(LNG.EQ.2) WRITE(LU,612) FLUX%R(IVBIL),
00454      &        -FLUDI(IVBIL),MASSEN%R(IVBIL),MASSE%R(IVBIL),DT*FLUTOT,
00455      &                      MASSEN%R(IVBIL)-MASSE%R(IVBIL)-DT*FLUTOT
00456             ELSE
00457               IF(LNG.EQ.1) THEN
00458                 WRITE(LU,621) IVBIL-5,FLUX%R(IVBIL),
00459      &          -FLUDI(IVBIL),MASSEN%R(IVBIL),MASSE%R(IVBIL),DT*FLUTOT,
00460      &                        MASSEN%R(IVBIL)-MASSE%R(IVBIL)-DT*FLUTOT
00461               ENDIF
00462               IF(LNG.EQ.2) THEN
00463                 WRITE(LU,622) IVBIL-5,FLUX%R(IVBIL),
00464      &          -FLUDI(IVBIL),MASSEN%R(IVBIL),MASSE%R(IVBIL),DT*FLUTOT,
00465      &                        MASSEN%R(IVBIL)-MASSE%R(IVBIL)-DT*FLUTOT
00466               ENDIF
00467             ENDIF
00468           ENDIF
00469 !
00470         ENDDO
00471 !
00472       ENDIF
00473 !
00474 !-----------------------------------------------------------------------
00475 !
00476 601   FORMAT(/,'  EAU',
00477      &       /,'MASSE AU PAS DE TEMPS PRECEDENT               : ',G16.7,
00478      &       /,'MASSE AU PAS DE TEMPS EN COURS                : ',G16.7,
00479      &       /,'MASSE SORTIE PAR LES LIMITES PENDANT CE TEMPS : ',G16.7,
00480      &       /,'ERREUR SUR LA MASSE AU COURS DU PAS DE TEMPS  : ',G16.7)
00481 !
00482 602   FORMAT(/,'  WATER',
00483      &       /,'MASS AT THE PREVIOUS TIME STEP                : ',G16.7,
00484      &       /,'MASS AT THE PRESENT TIME STEP                 : ',G16.7,
00485      &       /,'MASS LEAVING THE DOMAIN DURING THIS TIME STEP : ',G16.7,
00486      &       /,'ERROR ON THE MASS DURING THIS TIME STEP       : ',G16.7)
00487 !
00488 603   FORMAT(  'BILAN PLUIE-EVAPORATION                       : ',G16.7)
00489 604   FORMAT(  'BALANCE RAIN-EVAPORATION                      : ',G16.7)
00490 !
00491 611   FORMAT(/,'  SEDIMENT EN SUSPENSION ',
00492      &       /,'FLUX CONVECTIF A TRAVERS LES BORDS            : ',G16.7,
00493      &       /,'FLUX DIFFUSIF + DEPOT                         : ',G16.7,
00494      &       /,'MASSE AU PAS DE TEMPS PRECEDENT               : ',G16.7,
00495      &       /,'MASSE AU PAS DE TEMPS EN COURS                : ',G16.7,
00496      &       /,'MASSE SORTIE PAR LES LIMITES PENDANT CE TEMPS : ',G16.7,
00497      &       /,'ERREUR SUR LA MASSE AU COURS DU PAS DE TEMPS  : ',G16.7)
00498 !
00499 612   FORMAT(/,'  SEDIMENT IN SUSPENSION ',
00500      &       /,'ADVECTIVE FLUX THROUGH THE BOUNDARIES         : ',G16.7,
00501      &       /,'DIFFUSIVE FLUX + DEPOSITION                   : ',G16.7,
00502      &       /,'MASS AT THE PREVIOUS TIME STEP                : ',G16.7,
00503      &       /,'MASS AT THE PRESENT TIME STEP                 : ',G16.7,
00504      &       /,'MASS LEAVING THE DOMAIN DURING THIS TIME STEP : ',G16.7,
00505      &       /,'ERROR ON THE MASS DURING THIS TIME STEP       : ',G16.7)
00506 !
00507 621   FORMAT(/,'  TRACEUR  ',I2,
00508      &       /,'FLUX CONVECTIF A TRAVERS BORDS OU SOURCES     : ',G16.7,
00509      &       /,'FLUX DIFFUSIF A TRAVERS LES BORDS             : ',G16.7,
00510      &       /,'MASSE AU PAS DE TEMPS PRECEDENT               : ',G16.7,
00511      &       /,'MASSE AU PAS DE TEMPS EN COURS                : ',G16.7,
00512      &       /,'MASSE SORTIE (FRONTIERES OU SOURCE)           : ',G16.7,
00513      &       /,'ERREUR SUR LA MASSE AU COURS DU PAS DE TEMPS  : ',G16.7)
00514 !
00515 622   FORMAT(/,'  TRACER ',I2,
00516      &       /,'ADVECTIVE FLUX THROUGH BOUNDARIES OR SOURCES  : ',G16.7,
00517      &       /,'DIFFUSIVE FLUX THROUGH THE BOUNDARIES         : ',G16.7,
00518      &       /,'MASS AT THE PREVIOUS TIME STEP                : ',G16.7,
00519      &       /,'MASS AT THE PRESENT TIME STEP                 : ',G16.7,
00520      &       /,'MASS EXITING (BOUNDARIES OR SOURCE)           : ',G16.7,
00521      &       /,'ERROR ON THE MASS DURING THIS TIME STEP       : ',G16.7)
00522 3020  FORMAT('FLUX FRONTIERE ',I4,' : ', G16.7 ,' M3/S',
00523      &          '  ( >0 : ENTRANT  <0 : SORTANT )')
00524 4020  FORMAT('FLUX BOUNDARY ',I4,': ', G16.7 ,' M3/S',
00525      &          '  ( >0 : ENTERING  <0 : EXITING )')
00526 !
00527 !=======================================================================
00528 !
00529 ! PRINTOUTS SPECIFIC TO THE LAST TIMESTEP
00530 !
00531 !=======================================================================
00532 !
00533       IF(LT.EQ.NIT) THEN
00534 !
00535         WRITE(LU,*)
00536         CALL MITTIT(11,AT,LT)
00537         WRITE (LU,'(A4,F16.4)') 'T = ',AT
00538 !
00539         IF(LNG.EQ.1) WRITE(LU,701) MASINI_WATER,MASSE_WATER,
00540      &     DT*FLUXTOTCUM, MASINI_WATER-MASSE_WATER-DT*FLUXTOTCUM
00541         IF(LNG.EQ.2) WRITE(LU,702) MASINI_WATER,MASSE_WATER,
00542      &     DT*FLUXTOTCUM, MASINI_WATER-MASSE_WATER-DT*FLUXTOTCUM
00543 !
00544 !-----------------------------------------------------------------------
00545 !
00546         IF (NTRAC.GT.0) THEN
00547 !
00548           DO IVBIL=6,5+NTRAC
00549 !
00550             IF(SEDI.AND.(IVBIL.EQ.NTRAC+5)) THEN
00551               IF(LNG.EQ.1) WRITE(LU,711)
00552      &              MASINI%R(IVBIL),MASSE%R(IVBIL),DT*FLUCUM%R(IVBIL),
00553      &              MASINI%R(IVBIL)-MASSE%R(IVBIL)-DT*FLUCUM%R(IVBIL)
00554               IF(LNG.EQ.2) WRITE(LU,712)
00555      &              MASINI%R(IVBIL),MASSE%R(IVBIL),DT*FLUCUM%R(IVBIL),
00556      &              MASINI%R(IVBIL)-MASSE%R(IVBIL)-DT*FLUCUM%R(IVBIL)
00557             ELSE
00558               IF(LNG.EQ.1) WRITE(LU,721) IVBIL-5,
00559      &              MASINI%R(IVBIL),MASSE%R(IVBIL),DT*FLUCUM%R(IVBIL),
00560      &              MASINI%R(IVBIL)-MASSE%R(IVBIL)-DT*FLUCUM%R(IVBIL)
00561               IF(LNG.EQ.2) WRITE(LU,722) IVBIL-5,
00562      &              MASINI%R(IVBIL),MASSE%R(IVBIL),DT*FLUCUM%R(IVBIL),
00563      &              MASINI%R(IVBIL)-MASSE%R(IVBIL)-DT*FLUCUM%R(IVBIL)
00564             ENDIF
00565 !
00566           ENDDO
00567 !
00568         ENDIF
00569 !
00570       ENDIF
00571 !
00572 !-----------------------------------------------------------------------
00573 !
00574 701   FORMAT(//,'--- EAU ---',
00575      &       /,'MASSE INITIALE (DEBUT DE CE CALCUL) : ',G16.7,
00576      &       /,'MASSE FINALE                        : ',G16.7,
00577      &       /,'MASSE SORTIE DU DOMAINE (OU SOURCE) : ',G16.7,
00578      &       /,'MASSE PERDUE                        : ',G16.7)
00579 !
00580 702   FORMAT(//,'--- WATER ---',
00581      &       /,'INITIAL MASS                        : ',G16.7,
00582      &       /,'FINAL MASS                          : ',G16.7,
00583      &       /,'MASS LEAVING THE DOMAIN (OR SOURCE) : ',G16.7,
00584      &       /,'MASS LOSS                           : ',G16.7)
00585 !
00586 711   FORMAT(//,'--- SEDIMENT ---',
00587      &       /,'MASSE INITIALE (DEBUT DE CE CALCUL) : ',G16.7,
00588      &       /,'MASSE FINALE                        : ',G16.7,
00589      &       /,'MASSE SORTIE                        : ',G16.7,
00590      &       /,'MASSE PERDUE                        : ',G16.7)
00591 !
00592 712   FORMAT(//,'--- SEDIMENT ---',
00593      &       /,'INITIAL MASS                        : ',G16.7,
00594      &       /,'FINAL MASS                          : ',G16.7,
00595      &       /,'MASS LEAVING THE DOMAIN             : ',G16.7,
00596      &       /,'MASS LOSS                           : ',G16.7)
00597 !
00598 721   FORMAT(//,'--- TRACEUR ',I2,' ---',
00599      &       /,'MASSE INITIALE (DEBUT DE CE CALCUL) : ',G16.7,
00600      &       /,'MASSE FINALE                        : ',G16.7,
00601      &       /,'MASSE SORTIE (FRONTIERES OU SOURCE) : ',G16.7,
00602      &       /,'MASSE PERDUE                        : ',G16.7)
00603 !
00604 722   FORMAT(//,'--- TRACER',I2,' ---',
00605      &       /,'INITIAL MASS                        : ',G16.7,
00606      &       /,'FINAL MASS                          : ',G16.7,
00607      &       /,'MASS EXITING (BOUNDARIES OR SOURCE) : ',G16.7,
00608      &       /,'MASS LOSS                           : ',G16.7)
00609 !
00610 !-----------------------------------------------------------------------
00611 !
00612       RETURN
00613       END

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