Ignore:
Timestamp:
Mar 6, 2016 6:36:17 PM (8 years ago)
Author:
raasch
Message:

NetCDF routines modularized; new parameter netcdf_deflate; further changes in the pmc

File:
1 edited

Legend:

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

    r1692 r1783  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! name change of netcdf routines and module + related changes,
     22! switch back of netcdf data format moved from time integration routine to here
    2223!
    2324! Former revisions:
     
    9697   
    9798    USE control_parameters,                                                    &
    98         ONLY:  cloud_physics, domask, domask_no, domask_time_count,            &
    99                mask_i, mask_j, mask_k, mask_size, mask_size_l,                 &
    100                mask_start_l, max_masks, message_string, mid,                   &
    101                netcdf_data_format, nz_do3d, simulated_time
     99        ONLY:  cloud_physics, domask, domask_no, domask_time_count, mask_i,    &
     100               mask_j, mask_k, mask_size, mask_size_l, mask_start_l,           &
     101               max_masks, message_string, mid, nz_do3d, simulated_time
    102102    USE cpulog,                                                                &
    103103        ONLY:  cpu_log, log_point
    104104
    105 
    106    
    107105    USE indices,                                                               &
    108106        ONLY:  nbgp, nxl, nxr, nyn, nys, nzb, nzt
     
    110108    USE kinds
    111109   
    112     USE netcdf
    113    
    114     USE netcdf_control
     110    USE NETCDF
     111   
     112    USE netcdf_interface,                                                      &
     113        ONLY:  id_set_mask, id_var_domask, id_var_time_mask, nc_stat,          &
     114               netcdf_data_format, netcdf_handle_error
    115115   
    116116    USE particle_attributes,                                                   &
     
    136136    INTEGER(iwp) ::  k        !<
    137137    INTEGER(iwp) ::  n        !<
     138    INTEGER(iwp) ::  netcdf_data_format_save !<
    138139    INTEGER(iwp) ::  psi      !<
    139140    INTEGER(iwp) ::  sender   !<
     
    158159
    159160    CALL cpu_log (log_point(49),'data_output_mask','start')
     161
     162!
     163!-- Parallel netcdf output is not tested so far for masked data, hence
     164!-- netcdf_data_format is switched back to non-paralell output.
     165    netcdf_data_format_save = netcdf_data_format
     166    IF ( netcdf_data_format == 5 ) netcdf_data_format = 3
     167    IF ( netcdf_data_format == 6 ) netcdf_data_format = 4
    160168
    161169!
     
    183191                               start = (/ domask_time_count(mid,av) /),       &
    184192                               count = (/ 1 /) )
    185        CALL handle_netcdf_error( 'data_output_mask', 460 )
     193       CALL netcdf_handle_error( 'data_output_mask', 460 )
    186194    ENDIF
    187195
     
    564572               count = (/ mask_size_l(mid,1), mask_size_l(mid,2),  &
    565573                          mask_size_l(mid,3), 1 /) )
    566           CALL handle_netcdf_error( 'data_output_mask', 461 )
     574          CALL netcdf_handle_error( 'data_output_mask', 461 )
    567575       ELSE
    568576#endif
     
    611619                  count = (/ mask_size(mid,1), mask_size(mid,2), &
    612620                             mask_size(mid,3), 1 /) )
    613              CALL handle_netcdf_error( 'data_output_mask', 462 )
     621             CALL netcdf_handle_error( 'data_output_mask', 462 )
    614622
    615623          ELSE
     
    656664            count = (/ mask_size_l(mid,1), mask_size_l(mid,2), &
    657665                       mask_size_l(mid,3), 1 /) )
    658        CALL handle_netcdf_error( 'data_output_mask', 463 )
     666       CALL netcdf_handle_error( 'data_output_mask', 463 )
    659667#endif
    660668
     
    672680#endif
    673681
     682!
     683!-- Switch back to original format given by user (see beginning of this routine)
     684    netcdf_data_format = netcdf_data_format_save
    674685
    675686    CALL cpu_log( log_point(49), 'data_output_mask', 'stop' )
Note: See TracChangeset for help on using the changeset viewer.