gredelhyd_autop.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\gretel\gredelhyd_autop.f
00002 !
00058                      PROGRAM GREDELHYD_AUTOP
00059 !                    ***********************
00060 !
00061 !
00062 !***********************************************************************
00063 ! PARALLEL   V6P2                                   21/08/2010
00064 !***********************************************************************
00065 !
00066 !
00067 !
00068 !
00069 !
00070 !
00071 !
00072 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00073 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00074 !
00075       IMPLICIT NONE
00076       INTEGER LNG,LU
00077       COMMON/INFO/LNG,LU
00078       INTEGER LI
00079 !
00080       CHARACTER(LEN=30) GEO
00081 !
00082       INTEGER ERR,NELEM,ECKEN,NDUM,I,J,NBV1,NBV2,PARAM(10)
00083       INTEGER NPLAN,NPOIN2,NPROC,I_S,I_SP,I_LEN,IDUM,NPTFR,NSEG2,MBND
00084       INTEGER IYEAR,IMONTH,IDAY,IHOUR,IMIN,ISEC
00085 !
00086       INTEGER, DIMENSION(:)  , ALLOCATABLE :: LIHBOR             ! LIHBOR(NPTFR)
00087       INTEGER, DIMENSION(:)  , ALLOCATABLE :: NBOR               ! NBOR(*)
00088       INTEGER, DIMENSION(:)  , ALLOCATABLE :: NBOR0,LIHBOR0      ! NBOR0(NPTFR),LIHBOR0(NPTFR)
00089 !
00090       REAL RDUM
00091       REAL,    DIMENSION(:)  , ALLOCATABLE :: F
00092 !
00093       DOUBLE PRECISION REFER_DAY,JULIAN_DAY
00094       DOUBLE PRECISION JULTIM
00095       EXTERNAL         JULTIM
00096 !
00097       LOGICAL IS
00098 !
00099       CHARACTER*30 RES
00100       CHARACTER*50 RESPAR
00101       CHARACTER*11 EXTENS
00102       CHARACTER*30 CONLIM
00103       EXTERNAL    EXTENS
00104 !
00105       INTEGER ITSTRT,ITSTOP,NSTEPA
00106       INTEGER MARDAT(3),MARTIM(3)
00107       CHARACTER*144 TITRE
00108       CHARACTER*144 NOMGEO,NOMLIM
00109       CHARACTER*144 NOMSOU,NOMMAB,NOMCOU,NOMSAL,NOMTEM
00110       CHARACTER*144 NOMINI,NOMVEB,NOMMAF,NOMVEL,NOMVIS
00111       LOGICAL SALI_DEL,TEMP_DEL
00112       LOGICAL VELO_DEL,DIFF_DEL
00113 !
00114       LI=5
00115       LU=6
00116       LNG=2
00117 !HW
00118 !JAJ INTRODUCE YOURSELF WITH THE RELEASE DATE
00119 !
00120       WRITE(LU,*) 'I AM GREDELHYD, COUSIN OF GRETEL FROM BAW HAMBURG'
00121       WRITE(LU,*)
00122 !
00123       WRITE (LU, ADVANCE='NO',
00124      &    FMT='(/,'' GLOBAL GEOMETRY FILE: '')')
00125 !      REWIND(LI)
00126       READ(LI,*) GEO
00127       WRITE(LU,*) GEO
00128 !
00129 ! READS FILENAMES AND THE NUMBER OF PROCESSORS / PARTITIONS
00130 !
00131       WRITE (LU, ADVANCE='NO', FMT='(/,'' RESULT FILE: '')')
00132       READ(LI,*) RES
00133       WRITE(LU,*) RES
00134 !
00135       WRITE (LU,ADVANCE='NO',FMT='(/,'' NUMBER OF PROCESSORS: '')')
00136       READ (LI,*) NPROC
00137       WRITE(LU,*) NPROC
00138       INQUIRE (FILE=GEO,EXIST=IS)
00139       IF (.NOT.IS) THEN
00140         WRITE (LU,*) 'FILE DOES NOT EXIST: ', GEO
00141         CALL PLANTE(1)
00142         STOP
00143       END IF
00144 !
00145       I_S  = LEN (RES)
00146       I_SP = I_S + 1
00147       DO I=1,I_S
00148         IF(RES(I_SP-I:I_SP-I) .NE. ' ') EXIT
00149       ENDDO
00150       I_LEN=I_SP - I
00151 !
00152 !     GEOMETRY FILE, READ UNTIL 10 PARAMETERS:
00153 !
00154       OPEN(2,FILE=GEO,FORM='UNFORMATTED',STATUS='OLD',ERR=990)
00155       READ(2,ERR=990)
00156       READ(2,ERR=990) NBV1,NBV2
00157       DO I=1,NBV1+NBV2
00158         READ(2,ERR=990)
00159       ENDDO ! I
00160       GO TO 992
00161 990   WRITE(LU,*) 'ERROR WHEN OPENING OR READING FILE: ',GEO
00162       CALL PLANTE(1)
00163       STOP
00164 992   CONTINUE
00165 !     READS THE 10 PARAMETERS AND THE DATE
00166       READ(2) (PARAM(I),I=1,10)
00167       IF(PARAM(10).EQ.1) READ(2) (PARAM(I),I=1,6)
00168 !
00169 !     RESULTS FILE:
00170 !
00171       OPEN(3,FILE=RES,FORM='FORMATTED',ERR=991)
00172       GO TO 993
00173 991   WRITE(LU,*) 'ERROR WHEN OPENING FILE: ',RES
00174       CALL PLANTE(1)
00175       STOP
00176 993   CONTINUE
00177 !
00178 !     1) READS THE BEGINNING OF THE FIRST RESULTS FILE
00179 !
00180 !CC      RESPAR=RES // EXTENS(2**IDIMS-1,0)
00181 !
00182       RESPAR=RES(1:I_LEN) // EXTENS(NPROC-1,0)
00183 !
00184       INQUIRE (FILE=RESPAR,EXIST=IS)
00185       IF (.NOT.IS) THEN
00186         WRITE (LU,*) 'FILE DOES NOT EXIST: ', RESPAR
00187         WRITE (LU,*) 'CHECK THE NUMBER OF PROCESSORS'
00188         WRITE (LU,*) 'AND THE RESULT FILE CORE NAME'
00189         CALL PLANTE(1)
00190         STOP
00191       END IF
00192 !
00193       OPEN(4,FILE=RESPAR,FORM='FORMATTED',ERR=994)
00194       GO TO 995
00195 994   WRITE(LU,*) 'ERROR WHEN OPENING FILE: ',RESPAR
00196       CALL PLANTE(1)
00197       STOP
00198 995   CONTINUE
00199 !
00200       READ(4,'(I6)')NPLAN
00201       CLOSE(4)
00202 !
00203       ALLOCATE(F(NPLAN),STAT=ERR)
00204       CALL CHECK_ALLOCATE(ERR, 'F')
00205 !
00206 !  5 : 4 PARAMETERS
00207 !
00208       READ(2) NELEM,NPOIN2,ECKEN,NDUM
00209       WRITE(LU,*) '4 PARAMETERS IN GEOMETRY FILE'
00210       WRITE(LU,*) 'NELEM=',  NELEM
00211       WRITE(LU,*) 'NPOIN2=', NPOIN2
00212       WRITE(LU,*) 'ECKEN=',  ECKEN
00213       WRITE(LU,*) 'NDUM=',   NDUM
00214 !
00215 !----------------------------------------------------------------------
00216 !
00217       IF(NPLAN.LE.1) THEN
00218         CONLIM = "T2DCLI"
00219       ELSE
00220         CONLIM = "T3DCLI"
00221       ENDIF
00222 !
00223       OPEN(4,FILE=CONLIM,FORM='FORMATTED',ERR=996)
00224       GO TO 997
00225  996  WRITE(LU,*) 'ERROR WHEN OPENING FILE: ',CONLIM
00226       CALL PLANTE(1)
00227       STOP
00228  997  CONTINUE
00229 !
00230       ALLOCATE(LIHBOR0(NPOIN2),STAT=ERR)
00231       CALL CHECK_ALLOCATE(ERR, 'LIHBOR')
00232       ALLOCATE(NBOR0(NPOIN2),STAT=ERR)
00233       CALL CHECK_ALLOCATE(ERR, 'NBOR')
00234       DO I=1,NPOIN2
00235         READ(4,*,END=989) LIHBOR0(I),IDUM,IDUM,RDUM,RDUM,RDUM,RDUM,
00236      &                    IDUM,RDUM,RDUM,RDUM,NBOR0(I),IDUM
00237       ENDDO
00238 !
00239       CLOSE(4)
00240  989  NPTFR=I-1
00241 !
00242       ALLOCATE(LIHBOR(NPTFR),STAT=ERR)
00243       CALL CHECK_ALLOCATE(ERR, 'LIHBOR')
00244       ALLOCATE(NBOR(NPTFR),STAT=ERR)
00245       CALL CHECK_ALLOCATE(ERR, 'NBOR')
00246 !
00247       MBND=0
00248 !
00249       DO I=1,NPTFR
00250         NBOR(I)   = NBOR0(I)
00251         LIHBOR(I) = LIHBOR0(I)
00252         IF (LIHBOR(I).NE.2) THEN
00253           MBND = MBND + 1
00254         ENDIF
00255       ENDDO
00256 !
00257 !     WITH PRISMS, DIFFERENT FROM 2D VALUES, OTHERWISE
00258 !
00259       NSEG2 = (3*NELEM+NPTFR)/2
00260 !
00261 !
00262       OPEN(4,FILE=RESPAR,FORM='FORMATTED',ERR=984)
00263       GO TO 985
00264 984   WRITE(LU,*) 'ERROR WHEN OPENING FILE: ',RESPAR
00265       CALL PLANTE(1)
00266       STOP
00267 985   CONTINUE
00268 !
00269       READ(4,'(I6)')NPLAN
00270       READ(4,'(I3)')J
00271       READ(4,'(A)')TITRE(1:J)
00272       READ(4,'(I4)')MARDAT(1)
00273       READ(4,'(I2)')MARDAT(2)
00274       READ(4,'(I2)')MARDAT(3)
00275       READ(4,'(I2)')MARTIM(1)
00276       READ(4,'(I2)')MARTIM(2)
00277       READ(4,'(I2)')MARTIM(3)
00278       READ(4,'(I14)')ITSTRT
00279       READ(4,'(I14)')ITSTOP
00280       READ(4,'(I14)')NSTEPA
00281       READ(4,'(I6)')NPLAN
00282 !
00283       WRITE(3, '(A)' )
00284      &    "task      full-coupling                              "
00285       WRITE(3, '(A)' )
00286      &    "                                                     "
00287       WRITE(3, '(A)' )
00288      &    "#                                                    "
00289       WRITE(3, '(A)' )
00290      &    "# telemac data                                       "
00291       WRITE(3, '(A)' )
00292      &    "#                                                    "
00293       WRITE(3, '(A)' )
00294      &    "                                                     "
00295       WRITE(3, '(A)' )
00296      &    "geometry  finite-elements                            "
00297       WRITE(3, '(A)' )
00298      &    "                                                     "
00299       WRITE(3, '(A)' )
00300      &    "horizontal-aggregation       no                      "
00301       WRITE(3, '(A)' )
00302      &    "minimum-vert-diffusion-used  no                      "
00303       WRITE(3, '(A)' )
00304      &    "vertical-diffusion           calculated              "
00305       WRITE(3, '(A)' )
00306      &    "description                                          "
00307 !      J = LEN_TRIM(TITRE)
00308       IF ( J .GT. 40 ) THEN
00309         WRITE (3, '(A,A,A)' ) "   '",TITRE(1:40),"'"
00310         IF ( J .GT. 80 ) THEN
00311           WRITE (3, '(A,A,A)' ) "   '",TITRE(41:80),"'"
00312           IF ( J .GT. 120 ) THEN
00313             WRITE (3, '(A,A,A)' ) "   '",TITRE(81:120),"'"
00314           ELSE
00315             WRITE (3, '(A,A,A)' ) "   '",TITRE(81:J),"'"
00316           ENDIF
00317         ELSE
00318            WRITE (3, '(A,A,A)' ) "   '",TITRE(41:J),"'"
00319            WRITE (3, '(A)' )
00320      &   "   '                                    '            "
00321         ENDIF
00322       ELSE
00323         WRITE (3, '(A,A,A)' ) "   '",TITRE(1:J),"'"
00324         WRITE (3, '(A)' )
00325      &   "   '                                    '            "
00326         WRITE (3, '(A)' )
00327      &   "   '                                    '            "
00328       ENDIF
00329 !      WRITE(3, '(A,A,A)' )
00330 !     &    "   '",TITRE(1:J),"'"
00331 !      WRITE(3, '(A)' )
00332 !     &    "   '                                    '            "
00333 !      WRITE(3, '(A)' )
00334 !     &    "   '                                    '            "
00335       WRITE(3, '(A)' )
00336      &    "end-description                                      "
00337       WRITE(3, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00338      &"reference-time           '",MARDAT(1),MARDAT(2),MARDAT(3),
00339      &                             MARTIM(1),MARTIM(2),MARTIM(3),"'"
00340       REFER_DAY  = JULTIM(MARDAT(1),MARDAT(2),MARDAT(3),
00341      &                              MARTIM(1),MARTIM(2),MARTIM(3),0.D0)
00342       JULIAN_DAY = REFER_DAY + DBLE(ITSTRT)/(86400.D0*36525.D0)
00343       CALL GREGTIM( JULIAN_DAY, IYEAR, IMONTH, IDAY,
00344      &                                    IHOUR, IMIN,   ISEC )
00345       WRITE(3, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00346      &    "hydrodynamic-start-time  '",IYEAR,IMONTH,IDAY,
00347      &                                 IHOUR,IMIN  ,ISEC, "'"
00348       JULIAN_DAY = REFER_DAY + DBLE(ITSTOP)/(86400.D0*36525.D0)
00349       CALL GREGTIM( JULIAN_DAY, IYEAR, IMONTH, IDAY,
00350      &                                    IHOUR, IMIN,   ISEC )
00351       WRITE(3, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00352      &    "hydrodynamic-stop-time   '",IYEAR,IMONTH,IDAY,
00353      &                                 IHOUR,IMIN  ,ISEC, "'"
00354       WRITE(3, '(A,I14,A)' )
00355      &    "hydrodynamic-timestep    '",NSTEPA,"'"
00356       WRITE(3, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00357      &    "conversion-ref-time      '",MARDAT(1),MARDAT(2),MARDAT(3),
00358      &                                 MARTIM(1),MARTIM(2),MARTIM(3),"'"
00359       JULIAN_DAY = REFER_DAY + DBLE(ITSTRT)/(86400.D0*36525.D0)
00360       CALL GREGTIM( JULIAN_DAY, IYEAR, IMONTH, IDAY,
00361      &                                    IHOUR, IMIN,   ISEC )
00362       WRITE(3, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00363      &    "conversion-start-time    '",IYEAR,IMONTH,IDAY,
00364      &                                 IHOUR,IMIN  ,ISEC, "'"
00365       JULIAN_DAY = REFER_DAY + DBLE(ITSTOP)/(86400.D0*36525.D0)
00366       CALL GREGTIM( JULIAN_DAY, IYEAR, IMONTH, IDAY,
00367      &                                    IHOUR, IMIN,   ISEC )
00368       WRITE(3, '(A,I4.4,I2.2,I2.2,I2.2,I2.2,I2.2,A)' )
00369      &    "conversion-stop-time     '",IYEAR,IMONTH,IDAY,
00370      &                                 IHOUR,IMIN  ,ISEC, "'"
00371       WRITE(3, '(A,I14,A)' )
00372      &    "conversion-timestep      '",NSTEPA,"'"
00373       WRITE(3, '(A,I6)'  )
00374      &    "grid-cells-first-direction ",NPOIN2
00375       WRITE(3, '(A,I6,A)')
00376      &    "grid-cells-second-direction",NSEG2+MBND," # nr of exchanges!"
00377       WRITE(3, '(A,I6)' )
00378      &    "number-hydrodynamic-layers ",NPLAN
00379       WRITE(3, '(A,I6)' )
00380      &    "number-water-quality-layers",NPLAN
00381       READ(4,'(I3)')J
00382       READ(4,'(A)')NOMGEO(1:J)
00383       WRITE(3, '(A,A,A)' )
00384      &    "hydrodynamic-file        '",NOMGEO(1:J),"'"
00385       WRITE(3, '(A)' )
00386      &    "aggregation-file         none                        "
00387       WRITE(3, '(A,A,A)' )
00388      &    "grid-indices-file        '",NOMGEO(1:J),"'"
00389       READ(4,'(I3)')J
00390       READ(4,'(A)')NOMLIM(1:J)
00391       WRITE(3, '(A,A,A)' )
00392      &    "grid-coordinates-file    '",NOMLIM(1:J),"'"
00393       READ(4,'(I3)')J
00394       READ(4,'(A)') NOMSOU(1:J)
00395       I = J
00396       DO WHILE((NOMSOU(I:I).NE.'/').AND.(NOMSOU(I:I).NE.
00397 
00398 
00399 
00400 '\')     &                             .AND.(I.GE.1))        I = I-1      ENDDO      WRITE(3, '(A,A,A)
00401 ' )     &    "volumes-file             '",NOMSOU(I+1:J),"
00402 '"      READ(4,'(I3)
00403 ')J      READ(4,'(A)
00404 
00405 ') NOMMAB(1:J)      I = J      DO WHILE((NOMMAB(I:I).NE.'/').AND.(NOMMAB(I:I).NE.')
00406      &                             .AND.(I.GE.1))
00407         I = I-1
00408       ENDDO
00409       WRITE(3, '(A,A,A)' )
00410      &    "areas-file               '",NOMMAB(I+1:J),"'"
00411       READ(4,'(I3)')J
00412       READ(4,'(A)') NOMCOU(1:J)
00413       I = J
00414       DO WHILE((NOMCOU(I:I).NE.'/').AND.(NOMCOU(I:I).NE.
00415 
00416 
00417 
00418 '\')     &                             .AND.(I.GE.1))        I = I-1      ENDDO      WRITE(3, '(A,A,A)
00419 ' )     &    "flows-file               '",NOMCOU(I+1:J),"
00420 '"      READ(4,'(I3)
00421 ')J      READ(4,'(A)
00422 
00423 ') NOMVEB(1:J)      I = J      DO WHILE((NOMVEB(I:I).NE.'/').AND.(NOMVEB(I:I).NE.')
00424      &                             .AND.(I.GE.1))
00425         I = I-1
00426       ENDDO
00427       WRITE(3, '(A,A,A)' )
00428      &    "pointers-file            '",NOMVEB(I+1:J),"'"
00429       READ(4,'(I3)')J
00430       READ(4,'(A)')NOMMAF(1:J)
00431       I = J
00432       DO WHILE((NOMMAF(I:I).NE.'/').AND.(NOMMAF(I:I).NE.
00433 
00434 
00435 
00436 '\')     &                             .AND.(I.GE.1))        I = I-1      ENDDO      WRITE(3, '(A,A,A)
00437 ' )     &    "lengths-file             '",NOMMAF(I+1:J),"
00438 '"      READ(4,'(L1)
00439 
00440 ') SALI_DEL      IF(SALI_DEL) THEN        READ(4,'(I3)
00441 ')J        READ(4,'(A)
00442 
00443 ') NOMSAL(1:J)        I = J        DO WHILE((NOMSAL(I:I).NE.'/').AND.(NOMSAL(I:I).NE.')
00444      &                               .AND.(I.GE.1))
00445           I = I-1
00446         ENDDO
00447         WRITE(3, '(A,A,A)' )
00448      &    "salinity-file            '",NOMSAL(I+1:J),"'"
00449       ELSE
00450         WRITE(3, '(A)' )
00451      &    "salinity-file            none                        "
00452       ENDIF
00453       READ(4,'(L1)') TEMP_DEL
00454       IF(TEMP_DEL) THEN
00455         READ(4,'(I3)')J
00456         READ(4,'(A)') NOMTEM(1:J)
00457         I = J
00458         DO WHILE((NOMTEM(I:I).NE.'/').AND.(NOMTEM(I:I).NE.
00459 
00460 
00461 
00462 '\')     &                               .AND.(I.GE.1))          I = I-1        ENDDO        WRITE(3, '(A,A,A)
00463 ' )     &    "temperature-file         '",NOMTEM(I+1:J),"
00464 
00465 '"      ELSE        WRITE(3, '(A)
00466 
00467 
00468 ' )     &    "temperature-file         none                        "      ENDIF      READ(4,'(L1)
00469 
00470 ') DIFF_DEL      IF(DIFF_DEL) THEN        READ(4,'(I3)
00471 ')J        READ(4,'(A)
00472 
00473 ') NOMVIS(1:J)        I = J        DO WHILE((NOMVIS(I:I).NE.'/').AND.(NOMVIS(I:I).NE.')
00474      &                               .AND.(I.GE.1))
00475           I = I-1
00476         ENDDO
00477         WRITE(3, '(A,A,A)' )
00478      &    "vert-diffusion-file      '",NOMVIS(I+1:J),"'"
00479       ELSE
00480         WRITE(3, '(A)' )
00481      &    "vert-diffusion-file      none                        "
00482       ENDIF
00483       READ(4,'(L1)') VELO_DEL
00484       IF(VELO_DEL) THEN
00485         READ(4,'(I3)')J
00486         READ(4,'(A)') NOMVEL(1:J)
00487         I = J
00488         DO WHILE((NOMVEL(I:I).NE.'/').AND.(NOMVEL(I:I).NE.
00489 
00490 
00491 
00492 '\')     &                               .AND.(I.GE.1))          I = I-1        ENDDO        WRITE(3, '(A,A,A)
00493 ' )     &    "velocity-file            '",NOMVEL(I+1:J),"
00494 
00495 '"      ELSE        WRITE(3, '(A)
00496 
00497 
00498 ' )     &    "velocity-file            none                        "      ENDIF      READ(4,'(I3)
00499 ')J      READ(4,'(A)
00500 
00501 ') NOMINI(1:J)      I = J      DO WHILE((NOMINI(I:I).NE.'/').AND.(NOMINI(I:I).NE.')
00502      &                             .AND.(I.GE.1))
00503         I = I-1
00504       ENDDO
00505       WRITE(3, '(A,A,A)' )
00506      &    "surfaces-file            '",NOMINI(I+1:J),"'"
00507 !
00508       WRITE(3, '(A)' )
00509      &    "total-grid-file          none                        "
00510       WRITE(3, '(A)' )
00511      &    "discharges-file          none                        "
00512       WRITE(3, '(A)' )
00513      &    "chezy-coefficients-file  none                        "
00514       WRITE(3, '(A)' )
00515      &    "shear-stresses-file      none                        "
00516       WRITE(3, '(A)' )
00517      &    "walking-discharges-file  none                        "
00518       IF ( NPLAN .GT. 1 ) THEN
00519         WRITE(3, '(A)' )
00520      &      "minimum-vert-diffusion                            "
00521         WRITE(3, '(A)' )
00522      &      "   upper-layer       0.0000E+00                   "
00523         WRITE(3, '(A)' )
00524      &      "   lower-layer       0.0000E+00                   "
00525         WRITE(3, '(A)' )
00526      &      "   interface-depth   0.0000E+00                   "
00527         WRITE(3, '(A)' )
00528      &      "end-minimum-vert-diffusion                        "
00529       ENDIF
00530       WRITE(3, '(A)' )
00531      &    "constant-dispersion                                  "
00532       WRITE(3, '(A)' )
00533      &    "   first-direction    0.0000                         "
00534       WRITE(3, '(A)' )
00535      &    "   second-direction   0.0000                         "
00536       WRITE(3, '(A)' )
00537      &    "   third-direction    0.0000                         "
00538       WRITE(3, '(A)' )
00539      &    "end-constant-dispersion                              "
00540       WRITE(3, '(A)' )
00541      &    "hydrodynamic-layers                               "
00542       DO I=1,NPLAN
00543         READ(4,'(F10.4)')F(I)
00544       ENDDO
00545       DO I=1,NPLAN
00546         WRITE(3, '(F10.4)' ) F(I)
00547       ENDDO
00548       WRITE(3, '(A)' )
00549      &    "end-hydrodynamic-layers                           "
00550       WRITE(3, '(A)' )
00551      &    "water-quality-layers                              "
00552       DO I=1,NPLAN
00553         WRITE(3, '(F10.4)' ) 1.0
00554       ENDDO
00555       WRITE(3, '(A)' )
00556      &    "end-water-quality-layers                          "
00557       WRITE(3, '(A)' )
00558      &    "discharges                                           "
00559       WRITE(3, '(A)' )
00560      &    "end-discharges                                       "
00561 !
00562       WRITE(LU,*) 'END OF PROGRAM '
00563 !
00564       CLOSE(2)
00565       CLOSE(3)
00566       CLOSE(4)
00567 !
00568       STOP 0
00569       END PROGRAM GREDELHYD_AUTOP

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