Changeset 3030


Ignore:
Timestamp:
May 23, 2018 2:37:00 PM (6 years ago)
Author:
raasch
Message:

variable if renamed ivar

File:
1 edited

Legend:

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

    r2718 r3030  
    2525! -----------------
    2626! $Id$
     27! variable if renamed ivar
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    159162    INTEGER(iwp) ::  ngp      !<
    160163    INTEGER(iwp) ::  i        !<
    161     INTEGER(iwp) ::  if       !<
     164    INTEGER(iwp) ::  ivar     !<
    162165    INTEGER(iwp) ::  j        !<
    163166    INTEGER(iwp) ::  k        !<
     
    223226!
    224227!-- Loop over all variables to be written.
    225     if = 1
    226 
    227     DO  WHILE ( domask(mid,av,if)(1:1) /= ' ' )
     228    ivar = 1
     229
     230    DO  WHILE ( domask(mid,av,ivar)(1:1) /= ' ' )
    228231!
    229232!--    Reallocate local_pf on PE 0 since its shape changes during MPI exchange
    230        IF ( netcdf_data_format < 5   .AND.  myid == 0  .AND.  if > 1 )  THEN
     233       IF ( netcdf_data_format < 5   .AND.  myid == 0  .AND.  ivar > 1 )  THEN
    231234          DEALLOCATE( local_pf )
    232235          ALLOCATE( local_pf(mask_size_l(mid,1),mask_size_l(mid,2), &
     
    240243!--    Store the variable chosen.
    241244       resorted = .FALSE.
    242        SELECT CASE ( TRIM( domask(mid,av,if) ) )
     245       SELECT CASE ( TRIM( domask(mid,av,ivar) ) )
    243246
    244247          CASE ( 'e' )
     
    511514!--          Radiation quantity
    512515             IF ( radiation )  THEN
    513                 CALL radiation_data_output_mask(av, domask(mid,av,if), found,  &
     516                CALL radiation_data_output_mask(av, domask(mid,av,ivar), found,&
    514517                                                local_pf )
    515518             ENDIF
     
    518521!--          User defined quantity
    519522             IF ( .NOT. found )  THEN
    520                 CALL user_data_output_mask(av, domask(mid,av,if), found,       &
     523                CALL user_data_output_mask(av, domask(mid,av,ivar), found,     &
    521524                                           local_pf )
    522525             ENDIF
     
    526529             IF ( .NOT. found )  THEN
    527530                WRITE ( message_string, * ) 'no output available for: ', &
    528                                             TRIM( domask(mid,av,if) )
     531                                            TRIM( domask(mid,av,ivar) )
    529532                CALL message( 'data_output_mask', 'PA0327', 0, 0, 0, 6, 0 )
    530533             ENDIF
     
    559562!
    560563!--       (1) a. Parallel I/O using netCDF 4 (not yet tested)
    561           nc_stat = NF90_PUT_VAR( id_set_mask(mid,av),  &
    562                id_var_domask(mid,av,if),  &
    563                local_pf,  &
    564                start = (/ mask_start_l(mid,1), mask_start_l(mid,2),  &
    565                           mask_start_l(mid,3), domask_time_count(mid,av) /),  &
    566                count = (/ mask_size_l(mid,1), mask_size_l(mid,2),  &
     564          nc_stat = NF90_PUT_VAR( id_set_mask(mid,av),                         &
     565               id_var_domask(mid,av,ivar), local_pf,                           &
     566               start = (/ mask_start_l(mid,1), mask_start_l(mid,2),            &
     567                          mask_start_l(mid,3), domask_time_count(mid,av) /),   &
     568               count = (/ mask_size_l(mid,1), mask_size_l(mid,2),              &
    567569                          mask_size_l(mid,3), 1 /) )
    568570          CALL netcdf_handle_error( 'data_output_mask', 461 )
     
    608610             ENDDO
    609611
    610              nc_stat = NF90_PUT_VAR( id_set_mask(mid,av),  &
    611                   id_var_domask(mid,av,if), total_pf, &
    612                   start = (/ 1, 1, 1, domask_time_count(mid,av) /), &
    613                   count = (/ mask_size(mid,1), mask_size(mid,2), &
     612             nc_stat = NF90_PUT_VAR( id_set_mask(mid,av),                      &
     613                  id_var_domask(mid,av,ivar), total_pf,                        &
     614                  start = (/ 1, 1, 1, domask_time_count(mid,av) /),            &
     615                  count = (/ mask_size(mid,1), mask_size(mid,2),               &
    614616                             mask_size(mid,3), 1 /) )
    615617             CALL netcdf_handle_error( 'data_output_mask', 462 )
     
    652654!--    (2) For serial execution of PALM, the single processor (PE0) holds all
    653655!--    data and writes them directly to file.
    654        nc_stat = NF90_PUT_VAR( id_set_mask(mid,av),  &
    655             id_var_domask(mid,av,if),       &
    656             local_pf, &
    657             start = (/ 1, 1, 1, domask_time_count(mid,av) /), &
    658             count = (/ mask_size_l(mid,1), mask_size_l(mid,2), &
    659                        mask_size_l(mid,3), 1 /) )
     656       nc_stat = NF90_PUT_VAR( id_set_mask(mid,av),                            &
     657                               id_var_domask(mid,av,ivar), local_pf,           &
     658                             start = (/ 1, 1, 1, domask_time_count(mid,av) /), &
     659                             count = (/ mask_size_l(mid,1), mask_size_l(mid,2),&
     660                               mask_size_l(mid,3), 1 /) )
    660661       CALL netcdf_handle_error( 'data_output_mask', 463 )
    661662#endif
    662663
    663        if = if + 1
     664       ivar = ivar + 1
    664665
    665666    ENDDO
Note: See TracChangeset for help on using the changeset viewer.