Changeset 4387 for palm/trunk/SOURCE

Timestamp:

Jan 28, 2020 11:44:20 AM (5 years ago)

Author:

banzhafs

Message:

Added subroutine get_variable_string_generic to netcdf_data_input_mod

File:

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

palm/trunk/SOURCE/netcdf_data_input_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+!
 ! 
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Added subroutine get_variable_string_generic ( )
+! and added to interface get_variable to circumvent
+! unknown application-specific restrictions
+! in existing function get_variable_string ( ),
+! which is retained for backward compatibility (ECC)
+!
+! 4370 2020-01-10 14:00:44Z raasch
 ! collective read switched off on NEC Aurora to avoid hang situations
-! 
+!
 ! 4362 2020-01-07 17:15:02Z suehring
 ! Input of plant canopy variables from static driver moved to plant-canopy 
@@ -706 +713 @@
        MODULE PROCEDURE get_variable_5d_real_dynamic   ! 2B removed as z is out of emission_values
        MODULE PROCEDURE get_variable_string
+       MODULE PROCEDURE get_variable_string_generic    ! (ecc) generic string function
+
     END INTERFACE get_variable
 
@@ -4362 +4371 @@
 !> input file.  
 !------------------------------------------------------------------------------!
-    SUBROUTINE get_variable_string( id, variable_name, var_string, names_number)
+
+    SUBROUTINE get_variable_string( id, variable_name, var_string, names_number )
 #if defined( __netcdf )
 
@@ -4428 +4438 @@
 #endif
     END SUBROUTINE get_variable_string
+
+
+!
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Generalized routine for reading strings from a netCDF variable
+!> to replace existing get_variable_string ( )
+!>
+!> Improvements:
+!>   - Expanded string size limit from 25 to 512
+!>   - No longer removes spaces between text magically (this seems to have
+!>     been aimed at a very specific application, but I don't know what)
+!>   - Simplified implementation
+!>
+!> Caveats:
+!>   - Somehow I could not get the subroutine to work with str_array(:,:)
+!>     so I reverted to a hard-coded str_array(:,512), hopefully large enough
+!>     for most general applications.  This also means the character variable
+!>     used for str_array must be of size (:,512)
+!>     (ecc 20200128)    
+!------------------------------------------------------------------------------!
+
+ SUBROUTINE get_variable_string_generic ( id, var_name, str_array, num_str, str_len )
+
+    IMPLICIT NONE
+
+    CHARACTER(LEN=*),                INTENT(IN)    :: var_name       !> netCDF variable name
+    CHARACTER(LEN=512), ALLOCATABLE, INTENT(INOUT) :: str_array(:)   !> output string array
+
+    INTEGER(iwp)              :: buf_len   !> string buffer size
+    INTEGER(iwp)              :: id_var    !> netCDF variable ID
+    INTEGER(iwp)              :: k         !> generic counter
+
+    INTEGER(iwp), INTENT(IN)  :: id        !> netCDF file ID
+    INTEGER(iwp), INTENT(IN)  :: num_str   !> number of string elements in array
+    INTEGER(iwp), INTENT(IN)  :: str_len   !> size of each string element
+
+#if defined( __netcdf )
+
+!
+!-- set buffer length to up to hard-coded string size
+
+    buf_len = MIN( ABS(str_len), 512 )
+
+!
+!-- allocate necessary memories for string variables
+
+    ALLOCATE(str_array(num_str))
+!
+!-- get variable id
+
+    nc_stat = NF90_INQ_VARID( id, TRIM(var_name), id_var )
+!
+!-- extract string variables
+
+    DO k = 1, num_str
+       str_array(k) = ''
+       nc_stat = NF90_GET_VAR( id, id_var, str_array(k),  &
+                      start = (/ 1, k /), count = (/ buf_len, 1 /)  )
+       CALL handle_error ( 'get_variable_string_generic', 702 )
+    ENDDO
+
+#endif
+
+ END SUBROUTINE get_variable_string_generic
+
 
 !------------------------------------------------------------------------------!