Changeset 1308 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Mar 13, 2014 2:58:42 PM (10 years ago)

Author:

fricke

Message:

Adjustments for parallel NetCDF output for lccrayh/lccrayb (Cray XC30 systems)

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (4 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! +netcdf_data_format_save
+! Calculate fixed number of output time levels for parallel netcdf output.
+! For masked data, parallel netcdf output is not tested so far, hence
+! netcdf_data_format is switched back to non-paralell output.
 !
 ! Former revisions:
@@ -337 +340 @@
     CHARACTER (LEN=100) ::  action
 
-    INTEGER ::  i, ilen, iremote = 0, j, k, kk, position, prec
+    INTEGER ::  i, ilen, iremote = 0, j, k, kk, netcdf_data_format_save, &
+                position, prec
     LOGICAL ::  found, ldum
     REAL    ::  gradient, remote = 0.0, simulation_time_since_reference
@@ -3265 +3269 @@
 !
 !--    Generate masks for masked data output
+!--    Parallel netcdf output is not tested so far for masked data, hence 
+!--    netcdf_data_format is switched back to non-paralell output.
+       netcdf_data_format_save = netcdf_data_format
+       IF ( netcdf_data_format > 4 )  THEN
+          IF ( netcdf_data_format == 5 ) netcdf_data_format = 3
+          IF ( netcdf_data_format == 6 ) netcdf_data_format = 4
+          message_string = 'netCDF file formats '//                            &
+                           '5 (parallel netCDF 4) and ' //                     &
+                           '6 (parallel netCDF 4 Classic model) '//            &
+                           '&are currently not supported (not yet tested) ' // &
+                           'for masked data.&Using respective non-parallel' // &  
+                           ' output for masked data.'
+          CALL message( 'check_parameters', 'PA0383', 0, 0, 0, 6, 0 )
+       ENDIF
        CALL init_masks
+       netcdf_data_format = netcdf_data_format_save
     ENDIF
 
@@ -3294 +3313 @@
     ENDIF
 #endif
+
+!
+!-- Calculate fixed number of output time levels for parallel netcdf output.
+!-- The time dimension has to be limited for paralell output, otherwise the
+!-- performance drops significantly.
+    IF ( netcdf_output )  THEN
+       IF ( netcdf_data_format > 4 )  THEN
+          ntdim_3d(0) = INT( ( end_time - skip_time_do3d ) / dt_do3d )
+          IF ( do3d_at_begin ) ntdim_3d(0) = ntdim_3d(0) + 1
+          ntdim_3d(1) = INT( ( end_time - skip_time_data_output_av ) &
+                              / dt_data_output_av )
+          ntdim_2d_xy(0) = INT( ( end_time - skip_time_do2d_xy ) / dt_do2d_xy )
+          ntdim_2d_xz(0) = INT( ( end_time - skip_time_do2d_xz ) / dt_do2d_xz )
+          ntdim_2d_yz(0) = INT( ( end_time - skip_time_do2d_yz ) / dt_do2d_yz )
+          IF ( do2d_at_begin )  THEN
+             ntdim_2d_xy(0) = ntdim_2d_xy(0) + 1
+             ntdim_2d_xz(0) = ntdim_2d_xz(0) + 1
+             ntdim_2d_yz(0) = ntdim_2d_yz(0) + 1
+          ENDIF
+          ntdim_2d_xy(1) = ntdim_3d(1)
+          ntdim_2d_xz(1) = ntdim_3d(1)
+          ntdim_2d_yz(1) = ntdim_3d(1)
+       ENDIF
+    ENDIF
 
 #if ! defined( __check )