Ignore:
Timestamp:
Jul 27, 2007 9:09:17 AM (17 years ago)
Author:
raasch
Message:

preliminary version for coupled runs

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/check_open.f90

    r83 r102  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Different filenames are used in case of a coupled simulation,
     7! coupling_char added to all relevant filenames
    78!
    89! Former revisions:
     
    146147       CASE ( 11 )
    147148
    148           OPEN ( 11, FILE='PARIN', FORM='FORMATTED', STATUS='OLD' )
     149          OPEN ( 11, FILE='PARIN'//coupling_char, FORM='FORMATTED', &
     150                     STATUS='OLD' )
    149151
    150152       CASE ( 13 )
    151153
    152154          IF ( myid_char == '' )  THEN
    153              OPEN ( 13, FILE='BININ'//myid_char, FORM='UNFORMATTED', &
    154                         STATUS='OLD' )
    155           ELSE
    156              OPEN ( 13, FILE='BININ/'//myid_char, FORM='UNFORMATTED', &
    157                         STATUS='OLD' )
     155             OPEN ( 13, FILE='BININ'//coupling_char//myid_char, &
     156                        FORM='UNFORMATTED', STATUS='OLD' )
     157          ELSE
     158             OPEN ( 13, FILE='BININ'//coupling_char//'/'//myid_char, &
     159                        FORM='UNFORMATTED', STATUS='OLD' )
    158160          ENDIF
    159161
     
    161163
    162164          IF ( myid_char == '' )  THEN
    163              OPEN ( 14, FILE='BINOUT'//myid_char, FORM='UNFORMATTED', &
    164                         POSITION='APPEND' )
     165             OPEN ( 14, FILE='BINOUT'//coupling_char//myid_char, &
     166                        FORM='UNFORMATTED', POSITION='APPEND' )
    165167          ELSE
    166168             IF ( myid == 0  .AND. .NOT. openfile(file_id)%opened_before )  THEN
    167                 CALL local_system( 'mkdir  BINOUT' )
     169                CALL local_system( 'mkdir  BINOUT' // coupling_char )
    168170             ENDIF
    169171#if defined( __parallel )
     
    173175             CALL MPI_BARRIER( comm2d, ierr )
    174176#endif
    175              OPEN ( 14, FILE='BINOUT/'//myid_char_14, FORM='UNFORMATTED', &
    176                         POSITION='APPEND' )
     177             OPEN ( 14, FILE='BINOUT'//coupling_char//'/'//myid_char_14, &
     178                        FORM='UNFORMATTED', POSITION='APPEND' )
    177179          ENDIF
    178180
    179181       CASE ( 15 )
    180182
    181           OPEN ( 15, FILE='RUN_CONTROL', FORM='FORMATTED' )
     183          OPEN ( 15, FILE='RUN_CONTROL'//coupling_char, FORM='FORMATTED' )
    182184
    183185       CASE ( 16 )
    184186
    185           OPEN ( 16, FILE='LIST_PROFIL', FORM='FORMATTED' )
     187          OPEN ( 16, FILE='LIST_PROFIL'//coupling_char, FORM='FORMATTED' )
    186188
    187189       CASE ( 17 )
    188190
    189           OPEN ( 17, FILE='LIST_PROFIL_1D', FORM='FORMATTED' )
     191          OPEN ( 17, FILE='LIST_PROFIL_1D'//coupling_char, FORM='FORMATTED' )
    190192
    191193       CASE ( 18 )
    192194
    193           OPEN ( 18, FILE='CPU_MEASURES', FORM='FORMATTED' )
     195          OPEN ( 18, FILE='CPU_MEASURES'//coupling_char, FORM='FORMATTED' )
    194196
    195197       CASE ( 19 )
    196198
    197           OPEN ( 19, FILE='HEADER', FORM='FORMATTED' )
     199          OPEN ( 19, FILE='HEADER'//coupling_char, FORM='FORMATTED' )
    198200
    199201       CASE ( 20 )
    200202
    201203          IF ( myid == 0  .AND. .NOT. openfile(file_id)%opened_before )  THEN
    202              CALL local_system( 'mkdir  DATA_LOG' )
     204             CALL local_system( 'mkdir  DATA_LOG' // coupling_char )
    203205          ENDIF
    204206          IF ( myid_char == '' )  THEN
    205              OPEN ( 20, FILE='DATA_LOG/_0000', FORM='UNFORMATTED', &
    206                         POSITION='APPEND' )
     207             OPEN ( 20, FILE='DATA_LOG'//coupling_char//'/_0000', &
     208                        FORM='UNFORMATTED', POSITION='APPEND' )
    207209          ELSE
    208210#if defined( __parallel )
     
    212214             CALL MPI_BARRIER( comm2d, ierr )
    213215#endif
    214              OPEN ( 20, FILE='DATA_LOG/'//myid_char, FORM='UNFORMATTED', &
    215                         POSITION='APPEND' )
     216             OPEN ( 20, FILE='DATA_LOG'//coupling_char//'/'//myid_char, &
     217                        FORM='UNFORMATTED', POSITION='APPEND' )
    216218          ENDIF
    217219
     
    219221
    220222          IF ( data_output_2d_on_each_pe )  THEN
    221              OPEN ( 21, FILE='PLOT2D_XY'//myid_char, FORM='UNFORMATTED', &
    222                         POSITION='APPEND' )
    223           ELSE
    224              OPEN ( 21, FILE='PLOT2D_XY', FORM='UNFORMATTED', &
    225                         POSITION='APPEND' )
     223             OPEN ( 21, FILE='PLOT2D_XY'//coupling_char//myid_char, &
     224                        FORM='UNFORMATTED', POSITION='APPEND' )
     225          ELSE
     226             OPEN ( 21, FILE='PLOT2D_XY'//coupling_char, &
     227                        FORM='UNFORMATTED', POSITION='APPEND' )
    226228          ENDIF
    227229
     
    250252!--          Create output file for local parameters
    251253             IF ( iso2d_output )  THEN
    252                 OPEN ( 27, FILE='PLOT2D_XY_LOCAL', FORM='FORMATTED', &
    253                            DELIM='APOSTROPHE' )
     254                OPEN ( 27, FILE='PLOT2D_XY_LOCAL'//coupling_char, &
     255                           FORM='FORMATTED', DELIM='APOSTROPHE' )
    254256                openfile(27)%opened = .TRUE.
    255257             ENDIF
     
    260262
    261263          IF ( data_output_2d_on_each_pe )  THEN
    262              OPEN ( 22, FILE='PLOT2D_XZ'//myid_char, FORM='UNFORMATTED', &
    263                         POSITION='APPEND' )
    264           ELSE
    265              OPEN ( 22, FILE='PLOT2D_XZ', FORM='UNFORMATTED', &
     264             OPEN ( 22, FILE='PLOT2D_XZ'//coupling_char//myid_char, &
     265                        FORM='UNFORMATTED', POSITION='APPEND' )
     266          ELSE
     267             OPEN ( 22, FILE='PLOT2D_XZ'//coupling_char, FORM='UNFORMATTED', &
    266268                        POSITION='APPEND' )
    267269          ENDIF
     
    289291!
    290292!--          Create output file for local parameters
    291              OPEN ( 28, FILE='PLOT2D_XZ_LOCAL', FORM='FORMATTED', &
    292                         DELIM='APOSTROPHE' )
     293             OPEN ( 28, FILE='PLOT2D_XZ_LOCAL'//coupling_char, &
     294                        FORM='FORMATTED', DELIM='APOSTROPHE' )
    293295             openfile(28)%opened = .TRUE.
    294296
     
    298300
    299301          IF ( data_output_2d_on_each_pe )  THEN
    300              OPEN ( 23, FILE='PLOT2D_YZ'//myid_char, FORM='UNFORMATTED', &
    301                         POSITION='APPEND' )
    302           ELSE
    303              OPEN ( 23, FILE='PLOT2D_YZ', FORM='UNFORMATTED', &
     302             OPEN ( 23, FILE='PLOT2D_YZ'//coupling_char//myid_char, &
     303                        FORM='UNFORMATTED', POSITION='APPEND' )
     304          ELSE
     305             OPEN ( 23, FILE='PLOT2D_YZ'//coupling_char, FORM='UNFORMATTED', &
    304306                        POSITION='APPEND' )
    305307          ENDIF
     
    327329!
    328330!--          Create output file for local parameters
    329              OPEN ( 29, FILE='PLOT2D_YZ_LOCAL', FORM='FORMATTED', &
    330                         DELIM='APOSTROPHE' )
     331             OPEN ( 29, FILE='PLOT2D_YZ_LOCAL'//coupling_char, &
     332                        FORM='FORMATTED', DELIM='APOSTROPHE' )
    331333             openfile(29)%opened = .TRUE.
    332334
     
    335337       CASE ( 30 )
    336338
    337           OPEN ( 30, FILE='PLOT3D_DATA'//myid_char, FORM='UNFORMATTED' )
     339          OPEN ( 30, FILE='PLOT3D_DATA'//coupling_char//myid_char, &
     340                     FORM='UNFORMATTED' )
    338341!
    339342!--       Write coordinate file for AVS
     
    508511             WRITE ( suffix, '(''_'',I1)' )  file_id - 40
    509512          ENDIF
    510           OPEN ( file_id, FILE='PLOT1D_DATA'//TRIM( suffix ), FORM='FORMATTED' )
     513          OPEN ( file_id, FILE='PLOT1D_DATA'//coupling_char//TRIM( suffix ), &
     514                          FORM='FORMATTED' )
    511515!
    512516!--       Write contents comments at the top of the file
     
    521525             WRITE ( suffix, '(''_'',I1)' )  file_id - 50
    522526          ENDIF
    523           OPEN ( file_id, FILE='PLOTTS_DATA'//TRIM( suffix ), FORM='FORMATTED',&
    524                           RECL=496 )
     527          OPEN ( file_id, FILE='PLOTTS_DATA'//coupling_char//TRIM( suffix ), &
     528                          FORM='FORMATTED', RECL=496 )
    525529!
    526530!--       Write PROFIL parameter file for output of time series
     
    539543!
    540544!--       Write RAHMEN parameter
    541           OPEN ( 90, FILE='PLOTTS_PAR'//TRIM( suffix ), FORM='FORMATTED', &
    542                      DELIM='APOSTROPHE' )
     545          OPEN ( 90, FILE='PLOTTS_PAR'//coupling_char//TRIM( suffix ), &
     546                     FORM='FORMATTED', DELIM='APOSTROPHE' )
    543547          WRITE ( 90, RAHMEN )
    544548!
     
    575579
    576580          IF ( myid_char == '' )  THEN
    577              OPEN ( 80, FILE='PARTICLE_INFOS'//myid_char, FORM='FORMATTED', &
    578                         POSITION='APPEND' )
     581             OPEN ( 80, FILE='PARTICLE_INFOS'//coupling_char//myid_char, &
     582                        FORM='FORMATTED', POSITION='APPEND' )
    579583          ELSE
    580584             IF ( myid == 0  .AND.  .NOT. openfile(80)%opened_before )  THEN
    581                 CALL local_system( 'mkdir  PARTICLE_INFOS' )
     585                CALL local_system( 'mkdir  PARTICLE_INFOS' // coupling_char )
    582586             ENDIF
    583587#if defined( __parallel )
     
    592596             ENDIF
    593597#endif
    594              OPEN ( 80, FILE='PARTICLE_INFOS/'//myid_char, FORM='FORMATTED', &
    595                         POSITION='APPEND' )
     598             OPEN ( 80, FILE='PARTICLE_INFOS'//coupling_char//'/'//myid_char, &
     599                        FORM='FORMATTED', POSITION='APPEND' )
    596600          ENDIF
    597601
     
    602606       CASE ( 81 )
    603607
    604              OPEN ( 81, FILE='PLOTSP_X_PAR', FORM='FORMATTED', &
     608             OPEN ( 81, FILE='PLOTSP_X_PAR'//coupling_char, FORM='FORMATTED', &
    605609                        DELIM='APOSTROPHE', RECL=1500, POSITION='APPEND' )
    606610
    607611       CASE ( 82 )
    608612
    609              OPEN ( 82, FILE='PLOTSP_X_DATA', FORM='FORMATTED', &
     613             OPEN ( 82, FILE='PLOTSP_X_DATA'//coupling_char, FORM='FORMATTED', &
    610614                        POSITION = 'APPEND' )
    611615
    612616       CASE ( 83 )
    613617
    614              OPEN ( 83, FILE='PLOTSP_Y_PAR', FORM='FORMATTED', &
     618             OPEN ( 83, FILE='PLOTSP_Y_PAR'//coupling_char, FORM='FORMATTED', &
    615619                        DELIM='APOSTROPHE', RECL=1500, POSITION='APPEND' )
    616620
    617621       CASE ( 84 )
    618622
    619              OPEN ( 84, FILE='PLOTSP_Y_DATA', FORM='FORMATTED', &
     623             OPEN ( 84, FILE='PLOTSP_Y_DATA'//coupling_char, FORM='FORMATTED', &
    620624                        POSITION='APPEND' )
    621625
     
    623627
    624628          IF ( myid_char == '' )  THEN
    625              OPEN ( 85, FILE='PARTICLE_DATA'//myid_char, FORM='UNFORMATTED', &
    626                         POSITION='APPEND' )
     629             OPEN ( 85, FILE='PARTICLE_DATA'//coupling_char//myid_char, &
     630                        FORM='UNFORMATTED', POSITION='APPEND' )
    627631          ELSE
    628632             IF ( myid == 0  .AND.  .NOT. openfile(85)%opened_before )  THEN
    629                 CALL local_system( 'mkdir  PARTICLE_DATA' )
     633                CALL local_system( 'mkdir  PARTICLE_DATA' // coupling_char )
    630634             ENDIF
    631635#if defined( __parallel )
     
    635639             CALL MPI_BARRIER( comm2d, ierr )
    636640#endif
    637              OPEN ( 85, FILE='PARTICLE_DATA/'//myid_char, FORM='UNFORMATTED', &
    638                         POSITION='APPEND' )
     641             OPEN ( 85, FILE='PARTICLE_DATA'//coupling_char//'/'//myid_char, &
     642                        FORM='UNFORMATTED', POSITION='APPEND' )
    639643          ENDIF
    640644
     
    656660!--       Set filename depending on unit number
    657661          IF ( file_id == 101 )  THEN
    658              filename = 'DATA_2D_XY_NETCDF'
     662             filename = 'DATA_2D_XY_NETCDF' // coupling_char
    659663             av = 0
    660664          ELSE
    661              filename = 'DATA_2D_XY_AV_NETCDF'
     665             filename = 'DATA_2D_XY_AV_NETCDF' // coupling_char
    662666             av = 1
    663667          ENDIF
     
    716720!--       Set filename depending on unit number
    717721          IF ( file_id == 102 )  THEN
    718              filename = 'DATA_2D_XZ_NETCDF'
     722             filename = 'DATA_2D_XZ_NETCDF' // coupling_char
    719723             av = 0
    720724          ELSE
    721              filename = 'DATA_2D_XZ_AV_NETCDF'
     725             filename = 'DATA_2D_XZ_AV_NETCDF' // coupling_char
    722726             av = 1
    723727          ENDIF
     
    776780!--       Set filename depending on unit number
    777781          IF ( file_id == 103 )  THEN
    778              filename = 'DATA_2D_YZ_NETCDF'
     782             filename = 'DATA_2D_YZ_NETCDF' // coupling_char
    779783             av = 0
    780784          ELSE
    781              filename = 'DATA_2D_YZ_AV_NETCDF'
     785             filename = 'DATA_2D_YZ_AV_NETCDF' // coupling_char
    782786             av = 1
    783787          ENDIF
     
    834838       CASE ( 104 )
    835839!
     840!--       Set filename
     841          filename = 'DATA_1D_PR_NETCDF' // coupling_char
     842
     843!
    836844!--       Inquire, if there is a NetCDF file from a previuos run. This should
    837845!--       be opened for extension, if its variables match the actual run.
    838           INQUIRE( FILE='DATA_1D_PR_NETCDF', EXIST=netcdf_extend )
     846          INQUIRE( FILE=filename, EXIST=netcdf_extend )
    839847
    840848          IF ( netcdf_extend )  THEN
    841849!
    842850!--          Open an existing NetCDF file for output
    843              nc_stat = NF90_OPEN( 'DATA_1D_PR_NETCDF', NF90_WRITE, id_set_pr )
     851             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_pr )
    844852             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 29 )
    845853!
     
    854862                nc_stat = NF90_CLOSE( id_set_pr )
    855863                IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 30 )
    856                 CALL local_system( 'rm DATA_1D_PR_NETCDF' )
     864                CALL local_system( 'rm ' // TRIM( filename ) )
    857865             ENDIF
    858866
     
    864872             IF ( netcdf_64bit )  THEN
    865873#if defined( __netcdf_64bit )
    866                 nc_stat = NF90_CREATE( 'DATA_1D_PR_NETCDF',                    &
     874                nc_stat = NF90_CREATE( filename,                               &
    867875                                       OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
    868876                                       id_set_pr )
     
    870878                IF ( myid == 0 )  PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
    871879                                               'offset allowed on this machine'
    872                 nc_stat = NF90_CREATE( 'DATA_1D_PR_NETCDF', NF90_NOCLOBBER, &
    873                                        id_set_pr )
     880                nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_pr )
    874881#endif
    875882             ELSE
    876                 nc_stat = NF90_CREATE( 'DATA_1D_PR_NETCDF', NF90_NOCLOBBER, &
    877                                        id_set_pr )
     883                nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_pr )
    878884             ENDIF
    879885             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 31 )
     
    886892       CASE ( 105 )
    887893!
     894!--       Set filename
     895          filename = 'DATA_1D_TS_NETCDF' // coupling_char
     896
     897!
    888898!--       Inquire, if there is a NetCDF file from a previuos run. This should
    889899!--       be opened for extension, if its variables match the actual run.
    890           INQUIRE( FILE='DATA_1D_TS_NETCDF', EXIST=netcdf_extend )
     900          INQUIRE( FILE=filename, EXIST=netcdf_extend )
    891901
    892902          IF ( netcdf_extend )  THEN
    893903!
    894904!--          Open an existing NetCDF file for output
    895              nc_stat = NF90_OPEN( 'DATA_1D_TS_NETCDF', NF90_WRITE, id_set_ts )
     905             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_ts )
    896906             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 32 )
    897907!
     
    906916                nc_stat = NF90_CLOSE( id_set_ts )
    907917                IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 33 )
    908                 CALL local_system( 'rm DATA_1D_TS_NETCDF' )
     918                CALL local_system( 'rm ' // TRIM( filename ) )
    909919             ENDIF
    910920
     
    916926             IF ( netcdf_64bit )  THEN
    917927#if defined( __netcdf_64bit )
    918                 nc_stat = NF90_CREATE( 'DATA_1D_TS_NETCDF',                    &
     928                nc_stat = NF90_CREATE( filename,                               &
    919929                                       OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
    920930                                       id_set_ts )
     
    922932                IF ( myid == 0 )  PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
    923933                                               'offset allowed on this machine'
    924                 nc_stat = NF90_CREATE( 'DATA_1D_TS_NETCDF', NF90_NOCLOBBER, &
    925                                        id_set_ts )
     934                nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_ts )
    926935#endif
    927936             ELSE
    928                 nc_stat = NF90_CREATE( 'DATA_1D_TS_NETCDF', NF90_NOCLOBBER, &
    929                                        id_set_ts )
     937                nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_ts )
    930938             ENDIF
    931939             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 34 )
     
    941949!--       Set filename depending on unit number
    942950          IF ( file_id == 106 )  THEN
    943              filename = 'DATA_3D_NETCDF'
     951             filename = 'DATA_3D_NETCDF' // coupling_char
    944952             av = 0
    945953          ELSE
    946              filename = 'DATA_3D_AV_NETCDF'
     954             filename = 'DATA_3D_AV_NETCDF' // coupling_char
    947955             av = 1
    948956          ENDIF
     
    10001008       CASE ( 107 )
    10011009!
     1010!--       Set filename
     1011          filename = 'DATA_1D_SP_NETCDF' // coupling_char
     1012
     1013!
    10021014!--       Inquire, if there is a NetCDF file from a previuos run. This should
    10031015!--       be opened for extension, if its variables match the actual run.
    1004           INQUIRE( FILE='DATA_1D_SP_NETCDF', EXIST=netcdf_extend )
     1016          INQUIRE( FILE=filename, EXIST=netcdf_extend )
    10051017
    10061018          IF ( netcdf_extend )  THEN
    10071019!
    10081020!--          Open an existing NetCDF file for output
    1009              nc_stat = NF90_OPEN( 'DATA_1D_SP_NETCDF', NF90_WRITE, id_set_sp )
     1021             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_sp )
    10101022             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 38 )
    10111023!
     
    10201032                nc_stat = NF90_CLOSE( id_set_sp )
    10211033                IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 39 )
    1022                 CALL local_system( 'rm DATA_1D_SP_NETCDF' )
     1034                CALL local_system( 'rm ' // TRIM( filename ) )
    10231035             ENDIF
    10241036
     
    10301042             IF ( netcdf_64bit )  THEN
    10311043#if defined( __netcdf_64bit )
    1032                 nc_stat = NF90_CREATE( 'DATA_1D_SP_NETCDF',                    &
     1044                nc_stat = NF90_CREATE( filename,                               &
    10331045                                       OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
    10341046                                       id_set_sp )
     
    10361048                IF ( myid == 0 )  PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
    10371049                                               'offset allowed on this machine'
    1038                 nc_stat = NF90_CREATE( 'DATA_1D_SP_NETCDF', NF90_NOCLOBBER, &
    1039                                        id_set_sp )
     1050                nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_sp )
    10401051#endif
    10411052             ELSE
    1042                 nc_stat = NF90_CREATE( 'DATA_1D_SP_NETCDF', NF90_NOCLOBBER, &
    1043                                        id_set_sp )
     1053                nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_sp )
    10441054             ENDIF
    10451055             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 40 )
     
    10541064
    10551065          IF ( myid_char == '' )  THEN
    1056              filename = 'DATA_PRT_NETCDF'
    1057           ELSE
    1058              filename = 'DATA_PRT_NETCDF/' // myid_char
     1066             filename = 'DATA_PRT_NETCDF' // coupling_char
     1067          ELSE
     1068             filename = 'DATA_PRT_NETCDF' // coupling_char // '/' // myid_char
    10591069          ENDIF
    10601070!
     
    10911101                IF ( myid == 0  .AND. .NOT. openfile(file_id)%opened_before ) &
    10921102                THEN    ! needs modification in case of non-extendable sets
    1093                    CALL local_system( 'mkdir  DATA_PRT_NETCDF/' )
     1103                   CALL local_system( 'mkdir  DATA_PRT_NETCDF' // &
     1104                                       coupling_char // '/' )
    10941105                ENDIF
    10951106#if defined( __parallel )
     
    11261137       CASE ( 109 )
    11271138!
     1139!--       Set filename
     1140          filename = 'DATA_1D_PTS_NETCDF' // coupling_char
     1141
     1142!
    11281143!--       Inquire, if there is a NetCDF file from a previuos run. This should
    11291144!--       be opened for extension, if its variables match the actual run.
    1130           INQUIRE( FILE='DATA_1D_PTS_NETCDF', EXIST=netcdf_extend )
     1145          INQUIRE( FILE=filename, EXIST=netcdf_extend )
    11311146
    11321147          IF ( netcdf_extend )  THEN
    11331148!
    11341149!--          Open an existing NetCDF file for output
    1135              nc_stat = NF90_OPEN( 'DATA_1D_PTS_NETCDF', NF90_WRITE, id_set_pts )
     1150             nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_pts )
    11361151             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 393 )
    11371152!
     
    11461161                nc_stat = NF90_CLOSE( id_set_pts )
    11471162                IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 394 )
    1148                 CALL local_system( 'rm DATA_1D_PTS_NETCDF' )
     1163                CALL local_system( 'rm ' // TRIM( filename ) )
    11491164             ENDIF
    11501165
     
    11561171             IF ( netcdf_64bit )  THEN
    11571172#if defined( __netcdf_64bit )
    1158                 nc_stat = NF90_CREATE( 'DATA_1D_PTS_NETCDF',                   &
     1173                nc_stat = NF90_CREATE( filename,                               &
    11591174                                       OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
    11601175                                       id_set_pts )
     
    11621177                IF ( myid == 0 )  PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
    11631178                                               'offset allowed on this machine'
    1164                 nc_stat = NF90_CREATE( 'DATA_1D_PTS_NETCDF', NF90_NOCLOBBER, &
    1165                                        id_set_pts )
     1179                nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_pts )
    11661180#endif
    11671181             ELSE
    1168                 nc_stat = NF90_CREATE( 'DATA_1D_PTS_NETCDF', NF90_NOCLOBBER, &
    1169                                        id_set_pts )
     1182                nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_pts )
    11701183             ENDIF
    11711184             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 395 )
Note: See TracChangeset for help on using the changeset viewer.