Ignore:
Timestamp:
Sep 16, 2008 3:29:03 PM (13 years ago)
Author:
raasch
Message:

further adjustments for SGI and other small changes

File:
1 edited

Legend:

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

    r90 r197  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Time coordinate t=0 stored on netcdf-file only if an output is required for
     7! this time for at least one of the profiles
    78!
    89! Former revisions:
     
    3940
    4041    INTEGER ::  i, id, ilc, ils, j, k, sr
     42    LOGICAL ::  output_for_t0
    4143    REAL    ::  uxma, uxmi
    4244
     
    101103#if defined( __netcdf )
    102104!
    103 !--          Store initial time (t=0) to time axis         
    104              nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr, (/ 0.0 /), &
    105                                      start = (/ 1 /), count = (/ 1 /) )
    106              IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 329 )
     105!--          Store initial time (t=0) to time axis, but only if an output
     106!--          is required for at least one of the profiles
     107             output_for_t0 = .FALSE.
     108             DO  i = 1, dopr_n
     109                IF ( dopr_initial_index(i) /= 0 )  THEN
     110                   nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr,  &
     111                                           (/ 0.0 /), start = (/ 1 /), &
     112                                           count = (/ 1 /) )
     113                   IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 329 )
     114                   output_for_t0 = .TRUE.
     115                   EXIT
     116                ENDIF
     117             ENDDO
    107118
    108119!
     
    331342          ENDDO   ! Loop over dopr_n for initial profiles
    332343
    333           IF ( netcdf_output )  dopr_time_count = dopr_time_count + 1
     344          IF ( netcdf_output  .AND.  output_for_t0 )  THEN
     345             dopr_time_count = dopr_time_count + 1
     346          ENDIF
    334347
    335348       ENDIF   ! Initial profiles
Note: See TracChangeset for help on using the changeset viewer.