- Timestamp:
- Apr 4, 2018 6:01:41 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/netcdf_data_input_mod.f90
r2946 r2947 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Checks for dynamic input revised 28 ! 29 ! 2946 2018-04-04 17:01:23Z suehring 27 30 ! Bugfix for revision 2945, perform checks only if dynamic input file is 28 31 ! available. … … 1964 1967 1965 1968 USE indices, & 1966 ONLY: nx, nxl, nxlu, nxr, ny, nyn, nys, nysv, nzb, nz t1969 ONLY: nx, nxl, nxlu, nxr, ny, nyn, nys, nysv, nzb, nz, nzt 1967 1970 1968 1971 IMPLICIT NONE … … 2024 2027 CALL get_dimension_length( id_dynamic, init_3d%ny, 'y' ) 2025 2028 CALL get_dimension_length( id_dynamic, init_3d%nyv, 'yv' ) 2029 2030 ! 2031 !-- Check for correct horizontal and vertical dimension. Please note, 2032 !-- checks are performed directly here and not called from 2033 !-- check_parameters as some varialbes are still not allocated there. 2034 !-- Moreover, please note, u- and v-grid has 1 grid point less on 2035 !-- Inifor grid. 2036 IF ( init_3d%nx-1 /= nx .OR. init_3d%nxu-1 /= nx - 1 .OR. & 2037 init_3d%ny-1 /= ny .OR. init_3d%nyv-1 /= ny - 1 ) THEN 2038 message_string = 'Number of inifor horizontal grid points '// & 2039 'does not match the number of numeric grid '//& 2040 'points.' 2041 CALL message( 'netcdf_data_input_mod', 'NDI003', 1, 2, 0, 6, 0 ) 2042 ENDIF 2043 2044 IF ( init_3d%nzu-1 /= nz ) THEN 2045 message_string = 'Number of inifor vertical grid points ' // & 2046 'does not match the number of numeric grid '//& 2047 'points.' 2048 CALL message( 'netcdf_data_input_mod', 'NDI003', 1, 2, 0, 6, 0 ) 2049 ENDIF 2026 2050 ! 2027 2051 !-- Read vertical dimensions. Later, these are required for eventual … … 2782 2806 ONLY: initializing_actions, forcing, message_string 2783 2807 2784 USE indices, &2785 ONLY: nx, ny, nz2786 2787 2808 IMPLICIT NONE 2788 2809 … … 2803 2824 TRIM( coupling_char ) 2804 2825 CALL message( 'netcdf_data_input_mod', 'NDI010', 1, 2, 0, 6, 0 ) 2805 ENDIF2806 !2807 !-- Check for correct horizontal and vertical dimension.2808 !-- Please note, u- and v-grid has 1 grid point less on Inifor grid.2809 IF ( input_pids_dynamic ) THEN2810 2811 IF ( init_3d%nx-1 /= nx .OR. init_3d%nxu-1 /= nx - 1 .OR. &2812 init_3d%ny-1 /= ny .OR. init_3d%nyv-1 /= ny - 1 ) THEN2813 message_string = 'Number of inifor horizontal grid points does '//&2814 'not match the number of numeric grid points.'2815 CALL message( 'netcdf_data_input_mod', 'NDI003', 1, 2, 0, 6, 0 )2816 ENDIF2817 2818 IF ( init_3d%nzu /= nz ) THEN2819 message_string = 'Number of inifor vertical grid points does ' // &2820 'not match the number of numeric grid points.'2821 CALL message( 'netcdf_data_input_mod', 'NDI003', 1, 2, 0, 6, 0 )2822 ENDIF2823 2824 2826 ENDIF 2825 2827
Note: See TracChangeset
for help on using the changeset viewer.