wrihyd.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\wrihyd.f
00002 !
00058                      SUBROUTINE WRIHYD
00059 !                    *****************
00060 !
00061      &(TITRE , ITSTRT , ITSTOP , ITSTEP , NPOIN2 , MBND   ,
00062      & NSEG  , NOLAY  , NOMGEO , NOMLIM ,
00063      & F     , NSTEPA , NOMSOU , NOMMAB , NOMCOU ,
00064      & NOMINI, NOMVEB , NOMMAF , NOMSAL , NOMTEM , NOMVEL , NOMVIS ,
00065      & NHYD  , SALI_DEL,TEMP_DEL,VELO_DEL,DIFF_DEL,MARDAT,MARTIM)
00066 !
00067 !***********************************************************************
00068 ! TELEMAC2D   V6P2                                   21/08/2010
00069 !***********************************************************************
00070 !
00071 !
00072 !
00073 !
00074 !
00075 !
00076 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00077 !| DIFF_DEL       |-->| IF YES, WRITES DIFFUSION FILE FOR DELWAQ
00078 !| F              |-->| ARRAY TO STORE FRACTION OF DEPTH PER LAYER
00079 !| ITSTEP         |-->| TIME STEP
00080 !| ITSTOP         |-->| STOP TIME
00081 !| ITSTRT         |-->| START TIME
00082 !| MARDAT         |-->| DATE (YEAR, MONTH,DAY)
00083 !| MARTIM         |-->| TIME (HOUR, MINUTE,SECOND)
00084 !| MBND           |-->| SEQUENTIAL COUNTER OPEN BOUNDARIES
00085 !| NHYD           |-->| DELWAQ STEERING FILE CANAL
00086 !| NOLAY          |-->| NUMBER OF PLANES
00087 !| NOMCOU         |-->| FLUX FILE
00088 !| NOMGEO         |-->| RESULT FILE OF THE SIMULATION
00089 !| NOMINI         |-->| HORIZONTAL SURFACE FILE
00090 !| NOMLIM         |-->| BOUNDARY FILE OF THE SIMULATION
00091 !| NOMMAB         |-->| AREA FILE
00092 !| NOMMAF         |-->| NODE DISTANCE FILE
00093 !| NOMSAL         |-->| SALINITY FOR DELWAQ FILE
00094 !| NOMSOU         |-->| VOLUME FILE
00095 !| NOMTEM         |-->| TEMPERATURE FOR DELWAQ FILE
00096 !| NOMVEB         |-->| NODE EXCHANGE FILE
00097 !| NOMVEL         |-->| VELOCITY FILE
00098 !| NOMVIS         |-->| DIFFUSION FILE
00099 !| NPOIN2         |-->| NUMBER OF 2D POINTS IN THE MESH
00100 !| NSEG           |-->| NUMBER OF 2D SEGMENTS IN THE MESH
00101 !| NSTEPA         |-->| NUMBER OF TIME-STEPS FOR TIME AGGREGATION
00102 !| SALI_DEL       |-->| IF YES, THERE IS SALINITY
00103 !| TEMP_DEL       |-->| IF YES, THERE IS TEMPERATURE
00104 !| TITRE          |-->| TITLE OF STUDY
00105 !| VELO_DEL       |-->| IF YES, WRITES VELOCITY FILE FOR DELWAQ
00106 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00107 !
00108       IMPLICIT NONE
00109       INTEGER LNG,LU
00110       COMMON/INFO/LNG,LU
00111 !
00112 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00113 !
00114       INTEGER,          INTENT(IN) :: NHYD,ITSTRT,ITSTOP,ITSTEP,NPOIN2
00115       INTEGER,          INTENT(IN) :: NSEG,NOLAY,NSTEPA,MBND
00116       INTEGER,          INTENT(IN) :: MARDAT(3),MARTIM(3)
00117       CHARACTER(LEN=72),  INTENT(IN) :: TITRE
00118       CHARACTER(LEN=144), INTENT(IN) :: NOMSOU,NOMMAB,NOMCOU,NOMSAL
00119       CHARACTER(LEN=144), INTENT(IN) :: NOMINI,NOMVEB,NOMMAF,NOMVEL
00120       CHARACTER(LEN=144), INTENT(IN) :: NOMGEO,NOMLIM,NOMTEM,NOMVIS
00121       DOUBLE PRECISION, INTENT(IN) :: F(NPOIN2,NOLAY)
00122       LOGICAL,          INTENT(IN) :: SALI_DEL,TEMP_DEL
00123       LOGICAL,          INTENT(IN) :: VELO_DEL,DIFF_DEL
00124 !
00125 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00126 !
00127       INTEGER ILAY,IWAQ,I
00128       INTEGER IYEAR,IMONTH,IDAY,IHOUR,IMIN,ISEC
00129       DOUBLE PRECISION REFER_DAY,JULIAN_DAY
00130 !
00131       DOUBLE PRECISION JULTIM
00132       EXTERNAL         JULTIM
00133 !
00134 !-----------------------------------------------------------------------
00135 !
00136       IF( MARDAT(1) .EQ. 1900 ) THEN
00137         IF(LNG.EQ.1) THEN
00138           WRITE(LU,*) ' Avertissement grave :'
00139           WRITE(LU,*) ' VOUS UTILISEZ L ANNEE DE REFERENCE PAR DEFAUT'
00140           WRITE(LU,*) ' DE 1900 !!!'
00141           WRITE(LU,*) ' PUISQUE LES TEMPS DE DELWAQ SONT COMPTES A'
00142           WRITE(LU,*) ' PARTIR DE CETTE ANNEE ET S ETALENT SUR PLUS DE'
00143           WRITE(LU,*) ' 68 ANS, VOUS NE POURREZ PAS EXECUTER'
00144           WRITE(LU,*) ' APRES 1968 !!!'
00145           WRITE(LU,*) ' SOLUTION :'
00146           WRITE(LU,*) ' EDITEZ MANUELLEMENT LE FICHIER ~~~.hyd POUR'
00147           WRITE(LU,*) ' MODIFIER TOUTES LES ANNEES MENTIONNEES DE'
00148           WRITE(LU,*) ' 50 OU 100 ANS PLUS TARD !!!'
00149           WRITE(LU,*) ' VOUS N AVEZ PAS BESOIN DE REFAIRE TOURNER'
00150           WRITE(LU,*) ' L HYDRODYNAMIQUE !!!'
00151         ENDIF
00152         IF(LNG.EQ.2) THEN
00153           WRITE(LU,*) ' Severe warning:'
00154           WRITE(LU,*) ' YOU USE THE DEFAULT REFERENCE YEAR OF 1900!!!'
00155           WRITE(LU,*) ' SINCE DELWAQ TIMES ARE COUNTED FROM THIS YEAR,'
00156           WRITE(LU,*) ' AND SPAN AT MOST SOME 68 YEARS, YOU MAY NOT BE'
00157           WRITE(LU,*) ' ABLE TO RUN BEYOND 1968!!!'
00158           WRITE(LU,*) ' Remedy:'
00159           WRITE(LU,*) ' MANUALLY EDIT THE ~~~.hyd FILE TO ADJUST ALL'
00160           WRITE(LU,*) ' MENTIONED YEARS TO 50 OR 100 YEARS LATER!!!'
00161           WRITE(LU,*) ' YOU NEED not RERUN THE HYDRODYNAMICS!!!'
00162         ENDIF
00163       ENDIF
00164 !
00165       WRITE ( NHYD, '(A)' )
00166      &    "task      full-coupling                              "
00167       WRITE ( NHYD, '(A)' )
00168      &    "                                                     "
00169       WRITE ( NHYD, '(A)' )
00170      &    "#                                                    "
00171       WRITE ( NHYD, '(A)' )
00172      &    "# telemac data                                       "
00173       WRITE ( NHYD, '(A)' )
00174      &    "#                                                    "
00175       WRITE ( NHYD, '(A)' )
00176      &    "                                                     "
00177       WRITE ( NHYD, '(A)' )
00178      &    "geometry  finite-elements                            "
00179       WRITE ( NHYD, '(A)' )
00180      &    "                                                     "
00181       WRITE ( NHYD, '(A)' )
00182      &    "horizontal-aggregation       no                      "
00183       WRITE ( NHYD, '(A)' )
00184      &    "minimum-vert-diffusion-used  no                      "
00185       WRITE ( NHYD, '(A)' )
00186      &    "vertical-diffusion           calculated              "
00187       WRITE ( NHYD, '(A)' )
00188      &    "description                                          "
00189 !  CHARACTER(LEN=72),  INTENT(IN) :: TITRE
00190 !  PREVIOUS TESTS FOR LEN GREATER THAN 80 ARE USELESS
00191       IWAQ = LEN_TRIM(TITRE)
00192       IF ( IWAQ .EQ. 0 ) THEN
00193         WRITE ( NHYD, '(A)' )
00194      &    "   '                                    '            "
00195         WRITE ( NHYD, '(A)' )
00196      &    "   '                                    '            "
00197         WRITE ( NHYD, '(A)' )
00198      &    "   '                                    '            "
00199       ELSEIF ( IWAQ .LE. 40 ) THEN
00200         WRITE ( NHYD, '(A,A,A)' ) "   '",TITRE(1:IWAQ),"'"
00201         WRITE ( NHYD, '(A)' )
00202      &    "   '                                    '            "
00203         WRITE ( NHYD, '(A)' )
00204      &    "   '                                    '            "
00205       ELSE
00206         WRITE ( NHYD, '(A,A,A)' ) "   '",TITRE(1:40),"'"
00207         WRITE ( NHYD, '(A,A,A)' ) "   '",TITRE(41:IWAQ),"'"
00208         WRITE ( NHYD, '(A)' )
00209      &    "   '                                    '            "
00210       ENDIF
00211       WRITE ( NHYD, '(A)' )
00212      &    "end-description                                      "
00213       WRITE ( NHYD, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00214      &    "reference-time           '",MARDAT(1),MARDAT(2),MARDAT(3),
00215      &                                 MARTIM(1),MARTIM(2),MARTIM(3),"'"
00216       REFER_DAY  = JULTIM(MARDAT(1),MARDAT(2),MARDAT(3),
00217      &                    MARTIM(1),MARTIM(2),MARTIM(3),0.D0)
00218       JULIAN_DAY = REFER_DAY + DBLE(ITSTRT)/(86400.D0*36525.D0)
00219       CALL GREGTIM( JULIAN_DAY, IYEAR, IMONTH, IDAY, IHOUR, IMIN, ISEC )
00220       WRITE ( NHYD, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00221      &    "hydrodynamic-start-time  '",IYEAR,IMONTH,IDAY,
00222      &                                 IHOUR,IMIN  ,ISEC, "'"
00223       JULIAN_DAY = REFER_DAY + DBLE(ITSTOP)/(86400.D0*36525.D0)
00224       CALL GREGTIM( JULIAN_DAY, IYEAR, IMONTH, IDAY, IHOUR, IMIN, ISEC )
00225       WRITE ( NHYD, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00226      &    "hydrodynamic-stop-time   '",IYEAR,IMONTH,IDAY,
00227      &                                 IHOUR,IMIN  ,ISEC, "'"
00228 !      IDAY  = NSTEPA/86400
00229 !      IHOUR = (NSTEPA-IDAY*86400)/3600
00230 !      IMIN  = (NSTEPA-IDAY*86400-IHOUR*3600)/60
00231 !      ISEC  =  NSTEPA-IDAY*86400-IHOUR*3600-IMIN*60
00232 !      WRITE ( NHYD, '(A,I2.2,I2.2,I2.2,A)' )
00233 !     &    "hydrodynamic-timestep    '00000000",IHOUR,IMIN,ISEC,"'"
00234       WRITE ( NHYD, '(A,I14,A)' )
00235      &    "hydrodynamic-timestep    '",NSTEPA,"'"
00236       WRITE ( NHYD, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00237      &    "conversion-ref-time      '",MARDAT(1),MARDAT(2),MARDAT(3),
00238      &                                 MARTIM(1),MARTIM(2),MARTIM(3),"'"
00239       JULIAN_DAY = REFER_DAY + DBLE(ITSTRT)/(86400.D0*36525.D0)
00240       CALL GREGTIM( JULIAN_DAY, IYEAR, IMONTH, IDAY, IHOUR, IMIN, ISEC )
00241       WRITE ( NHYD, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00242      &    "conversion-start-time    '",IYEAR,IMONTH,IDAY,
00243      &                                 IHOUR,IMIN  ,ISEC, "'"
00244       JULIAN_DAY = REFER_DAY + DBLE(ITSTOP)/(86400.D0*36525.D0)
00245       CALL GREGTIM( JULIAN_DAY, IYEAR, IMONTH, IDAY, IHOUR, IMIN, ISEC )
00246       WRITE ( NHYD, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00247      &    "conversion-stop-time     '",IYEAR,IMONTH,IDAY,
00248      &                                 IHOUR,IMIN  ,ISEC, "'"
00249 !      IDAY  = NSTEPA/86400
00250 !      IHOUR = (NSTEPA-IDAY*86400)/3600
00251 !      IMIN  = (NSTEPA-IDAY*86400-IHOUR*3600)/60
00252 !      ISEC  =  NSTEPA-IDAY*86400-IHOUR*3600-IMIN*60
00253 !      WRITE ( NHYD, '(A,I2.2,I2.2,I2.2,A)' )
00254 !     &    "conversion-timestep      '00000000",IHOUR,IMIN,ISEC,"'"
00255       WRITE ( NHYD, '(A,I14,A)' )
00256      &    "conversion-timestep      '",NSTEPA,"'"
00257       WRITE ( NHYD, '(A,I6)'  )
00258      &    "grid-cells-first-direction ",NPOIN2
00259       WRITE ( NHYD, '(A,I6,A)')
00260      &    "grid-cells-second-direction",NSEG+MBND," # nr of exchanges!"
00261       WRITE ( NHYD, '(A,I6)' )
00262      &    "number-hydrodynamic-layers ",NOLAY
00263       WRITE ( NHYD, '(A,I6)' )
00264      &    "number-water-quality-layers",NOLAY
00265       IWAQ = LEN_TRIM(NOMGEO)
00266       I    = IWAQ
00267       DO WHILE(I.GE.1)
00268         IF((NOMGEO(I:I).NE.'/').AND.(NOMGEO(I:I).NE.
00269 
00270 
00271 
00272 
00273 
00274 
00275 '\')) THEN          I = I-1        ELSE          EXIT        ENDIF      ENDDO      IF(IWAQ.NE.0) THEN        WRITE ( NHYD, '(A,A,A)
00276 ' )     &    "hydrodynamic-file        '",NOMGEO(I+1:IWAQ),"
00277 
00278 '"      ELSE        WRITE ( NHYD, '(A)' ) "hydrodynamic-file        '
00279 
00280 '"      ENDIF      WRITE ( NHYD, '(A)
00281 
00282 
00283 ' )     &    "aggregation-file         none                        "      IF(IWAQ.NE.0) THEN        WRITE ( NHYD, '(A,A,A)
00284 ' )     &    "grid-indices-file        '",NOMGEO(I+1:IWAQ),"
00285 
00286 '"      ELSE        WRITE ( NHYD, '(A)' ) "grid-indices-file        '
00287 
00288 
00289 
00290 
00291 '"      ENDIF      IWAQ = LEN_TRIM(NOMLIM)      I    = IWAQ      DO WHILE(I.GE.1)        IF((NOMLIM(I:I).NE.'/').AND.(NOMLIM(I:I).NE.')) THEN
00292           I = I-1
00293         ELSE
00294           EXIT
00295         ENDIF
00296       ENDDO
00297       IF(IWAQ.NE.0) THEN
00298         WRITE ( NHYD, '(A,A,A)' )
00299      &    "grid-coordinates-file    '",NOMLIM(I+1:IWAQ),"'"
00300       ELSE
00301         WRITE ( NHYD, '(A)' ) "grid-coordinates-file    ''"
00302       ENDIF
00303       IWAQ = LEN_TRIM(NOMSOU)
00304       I    = IWAQ
00305       DO WHILE(I.GE.1)
00306         IF((NOMSOU(I:I).NE.'/').AND.(NOMSOU(I:I).NE.
00307 
00308 
00309 
00310 
00311 
00312 
00313 '\')) THEN          I = I-1        ELSE          EXIT        ENDIF      ENDDO      IF(IWAQ.NE.0) THEN        WRITE ( NHYD, '(A,A,A)
00314 ' )     &    "volumes-file             '",NOMSOU(I+1:IWAQ),"
00315 
00316 '"      ELSE        WRITE ( NHYD, '(A)' ) "volumes-file             '
00317 
00318 
00319 
00320 
00321 '"      ENDIF      IWAQ = LEN_TRIM(NOMMAB)      I    = IWAQ      DO WHILE(I.GE.1)        IF((NOMMAB(I:I).NE.'/').AND.(NOMMAB(I:I).NE.')) THEN
00322           I = I-1
00323         ELSE
00324           EXIT
00325         ENDIF
00326       ENDDO
00327       IF(IWAQ.NE.0) THEN
00328         WRITE ( NHYD, '(A,A,A)' )
00329      &    "areas-file               '",NOMMAB(I+1:IWAQ),"'"
00330       ELSE
00331         WRITE ( NHYD, '(A)' ) "areas-file               ''"
00332       ENDIF
00333       IWAQ = LEN_TRIM(NOMCOU)
00334       I    = IWAQ
00335       DO WHILE(I.GE.1)
00336         IF((NOMCOU(I:I).NE.'/').AND.(NOMCOU(I:I).NE.
00337 
00338 
00339 
00340 
00341 
00342 
00343 '\')) THEN          I = I-1        ELSE          EXIT        ENDIF      ENDDO      IF(IWAQ.NE.0) THEN        WRITE ( NHYD, '(A,A,A)
00344 ' )     &    "flows-file               '",NOMCOU(I+1:IWAQ),"
00345 
00346 '"      ELSE        WRITE ( NHYD, '(A)' ) "flows-file               '
00347 
00348 
00349 
00350 
00351 '"      ENDIF      IWAQ = LEN_TRIM(NOMVEB)      I    = IWAQ      DO WHILE(I.GE.1)        IF((NOMVEB(I:I).NE.'/').AND.(NOMVEB(I:I).NE.')) THEN
00352           I = I-1
00353         ELSE
00354           EXIT
00355         ENDIF
00356       ENDDO
00357       IF(IWAQ.NE.0) THEN
00358         WRITE ( NHYD, '(A,A,A)' )
00359      &    "pointers-file            '",NOMVEB(I+1:IWAQ),"'"
00360       ELSE
00361         WRITE ( NHYD, '(A)' ) "pointers-file            ''"
00362       ENDIF
00363       IWAQ = LEN_TRIM(NOMMAF)
00364       I    = IWAQ
00365       DO WHILE(I.GE.1)
00366         IF((NOMMAF(I:I).NE.'/').AND.(NOMMAF(I:I).NE.
00367 
00368 
00369 
00370 
00371 
00372 
00373 '\')) THEN          I = I-1        ELSE          EXIT        ENDIF      ENDDO      IF(IWAQ.NE.0) THEN        WRITE ( NHYD, '(A,A,A)
00374 ' )     &    "lengths-file             '",NOMMAF(I+1:IWAQ),"
00375 
00376 '"      ELSE        WRITE ( NHYD, '(A)' ) "lengths-file             '
00377 
00378 
00379 
00380 
00381 
00382 '"      ENDIF      IF(SALI_DEL) THEN        IWAQ = LEN_TRIM(NOMSAL)        I    = IWAQ        DO WHILE(I.GE.1)          IF((NOMSAL(I:I).NE.'/').AND.(NOMSAL(I:I).NE.')) THEN
00383             I = I-1
00384           ELSE
00385             EXIT
00386           ENDIF
00387         ENDDO
00388         IF(IWAQ.NE.0) THEN
00389           WRITE ( NHYD, '(A,A,A)' )
00390      &    "salinity-file            '",NOMSAL(I+1:IWAQ),"'"
00391         ELSE
00392           WRITE ( NHYD, '(A)' ) "salinity-file            ''"
00393         ENDIF
00394       ELSE
00395         WRITE ( NHYD, '(A)' )
00396      &    "salinity-file            none                        "
00397       ENDIF
00398       IF(TEMP_DEL) THEN
00399         IWAQ = LEN_TRIM(NOMTEM)
00400         I    = IWAQ
00401         DO WHILE(I.GE.1)
00402           IF((NOMTEM(I:I).NE.'/').AND.(NOMTEM(I:I).NE.
00403 
00404 
00405 
00406 
00407 
00408 
00409 '\')) THEN            I = I-1          ELSE            EXIT          ENDIF        ENDDO        IF(IWAQ.NE.0) THEN          WRITE ( NHYD, '(A,A,A)
00410 ' )     &    "temperature-file         '",NOMTEM(I+1:IWAQ),"
00411 
00412 '"        ELSE          WRITE ( NHYD, '(A)' ) "temperature-file         '
00413 
00414 
00415 '"        ENDIF      ELSE        WRITE ( NHYD, '(A)
00416 
00417 
00418 
00419 
00420 
00421 
00422 ' )     &    "temperature-file         none                        "      ENDIF      IF(DIFF_DEL) THEN        IWAQ = LEN_TRIM(NOMVIS)        I    = IWAQ        DO WHILE(I.GE.1)          IF((NOMVIS(I:I).NE.'/').AND.(NOMVIS(I:I).NE.')) THEN
00423             I = I-1
00424           ELSE
00425             EXIT
00426           ENDIF
00427         ENDDO
00428         IF(IWAQ.NE.0) THEN
00429           WRITE ( NHYD, '(A,A,A)' )
00430      &    "vert-diffusion-file      '",NOMVIS(I+1:IWAQ),"'"
00431         ELSE
00432           WRITE ( NHYD, '(A)' ) "vert-diffusion-file      ''"
00433         ENDIF
00434       ELSE
00435         WRITE ( NHYD, '(A)' )
00436      &    "vert-diffusion-file      none                        "
00437       ENDIF
00438       IF(VELO_DEL) THEN
00439         IWAQ = LEN_TRIM(NOMVEL)
00440         I    = IWAQ
00441         DO WHILE(I.GE.1)
00442           IF((NOMVEL(I:I).NE.'/').AND.(NOMVEL(I:I).NE.
00443 
00444 
00445 
00446 
00447 
00448 
00449 '\')) THEN            I = I-1          ELSE            EXIT          ENDIF        ENDDO        IF(IWAQ.NE.0) THEN          WRITE ( NHYD, '(A,A,A)
00450 ' )     &    "velocity-file            '",NOMVEL(I+1:IWAQ),"
00451 
00452 '"        ELSE          WRITE ( NHYD, '(A)' ) "velocity-file            '
00453 
00454 
00455 '"        ENDIF      ELSE        WRITE ( NHYD, '(A)
00456 
00457 
00458 
00459 
00460 
00461 ' )     &    "velocity-file            none                        "      ENDIF      IWAQ = LEN_TRIM(NOMINI)      I    = IWAQ      DO WHILE(I.GE.1)        IF((NOMINI(I:I).NE.'/').AND.(NOMINI(I:I).NE.')) THEN
00462           I = I-1
00463         ELSE
00464           EXIT
00465         ENDIF
00466       ENDDO
00467       IF(IWAQ.NE.0) THEN
00468         WRITE ( NHYD, '(A,A,A)' )
00469      &    "surfaces-file            '",NOMINI(I+1:IWAQ),"'"
00470       ELSE
00471         WRITE ( NHYD, '(A)' ) "surfaces-file            ''"
00472       ENDIF
00473       WRITE ( NHYD, '(A)' )
00474      &    "total-grid-file          none                        "
00475       WRITE ( NHYD, '(A)' )
00476      &    "discharges-file          none                        "
00477       WRITE ( NHYD, '(A)' )
00478      &    "chezy-coefficients-file  none                        "
00479       WRITE ( NHYD, '(A)' )
00480      &    "shear-stresses-file      none                        "
00481       WRITE ( NHYD, '(A)' )
00482      &    "walking-discharges-file  none                        "
00483       IF ( NOLAY .GT. 1 ) THEN
00484         WRITE ( NHYD, '(A)' )
00485      &       "minimum-vert-diffusion                            "
00486         WRITE ( NHYD, '(A)' )
00487      &       "   upper-layer       0.0000E+00                   "
00488         WRITE ( NHYD, '(A)' )
00489      &       "   lower-layer       0.0000E+00                   "
00490         WRITE ( NHYD, '(A)' )
00491      &       "   interface-depth   0.0000E+00                   "
00492         WRITE ( NHYD, '(A)' )
00493      &       "end-minimum-vert-diffusion                        "
00494       ENDIF
00495       WRITE ( NHYD, '(A)' )
00496      &    "constant-dispersion                                  "
00497       WRITE ( NHYD, '(A)' )
00498      &    "   first-direction    0.0000                         "
00499       WRITE ( NHYD, '(A)' )
00500      &    "   second-direction   0.0000                         "
00501       WRITE ( NHYD, '(A)' )
00502      &    "   third-direction    0.0000                         "
00503       WRITE ( NHYD, '(A)' )
00504      &    "end-constant-dispersion                              "
00505       WRITE ( NHYD, '(A)' )
00506      &    "hydrodynamic-layers                               "
00507 !      DO ILAY=1,NOLAY
00508 !     FROM TOP TO BOTTOM IN DELWAQ
00509       DO ILAY=NOLAY,1,-1
00510         WRITE ( NHYD, '(F10.4)' ) F(1,ILAY)
00511       ENDDO
00512       WRITE ( NHYD, '(A)' )
00513      &    "end-hydrodynamic-layers                           "
00514       WRITE ( NHYD, '(A)' )
00515      &    "water-quality-layers                              "
00516       DO ILAY=1,NOLAY
00517         WRITE ( NHYD, '(F10.4)' ) 1.0
00518       ENDDO
00519       WRITE ( NHYD, '(A)' )
00520      &    "end-water-quality-layers                          "
00521       WRITE ( NHYD, '(A)' )
00522      &    "discharges                                           "
00523       WRITE ( NHYD, '(A)' )
00524      &    "end-discharges                                       "
00525 !
00526 !-----------------------------------------------------------------------
00527 !
00528       RETURN
00529       END

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