Ignore:
Timestamp:
Mar 23, 2018 4:30:46 PM (4 years ago)
Author:
suehring
Message:

Remove default surfaces from radiation model and add check for this; revise checks for surface_fraction

File:
1 edited

Legend:

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

    r2925 r2930  
    2525! -----------------
    2626! $Id$
     27! Revise checks for surface_fraction.
     28!
     29! 2925 2018-03-23 14:54:11Z suehring
    2730! Check for further inconsistent settings of surface_fractions.
    2831! Some messages slightly rephrased and error numbers renamed.
     
    27852788       INTEGER(iwp) ::  i      !< loop index along x-direction
    27862789       INTEGER(iwp) ::  j      !< loop index along y-direction
     2790       INTEGER(iwp) ::  n_surf !< number of different surface types at given location
    27872791
    27882792       LOGICAL      ::  check_passed !< flag indicating if a check passed
     
    29912995             ENDIF
    29922996!
    2993 !--          Check for consistency of surface fraction.
    2994 !--          Sum of surface fractions must not exceed one.
    2995              IF ( ANY ( surface_fraction_f%frac(:,j,i) ==                      &
    2996                         surface_fraction_f%fill ) )  THEN
    2997                 message_string = 'If more than one natural surface type is ' //&
    2998                               'given at a location, surface_fraction ' //      &
    2999                               'must be provided.'
    3000                 CALL message( 'netcdf_data_input_mod', 'NDI027',               &
     2997!--          Check for consistency of surface fraction. If more than one type
     2998!--          is set, surface fraction need to be given and the sum must not
     2999!--          be larger than 1.
     3000             n_surf = 0
     3001             IF ( vegetation_type_f%var(j,i) /= vegetation_type_f%fill )       &
     3002                n_surf = n_surf + 1
     3003             IF ( water_type_f%var(j,i)      /= water_type_f%fill )            &
     3004                n_surf = n_surf + 1
     3005             IF ( pavement_type_f%var(j,i)   /= pavement_type_f%fill )         &
     3006                n_surf = n_surf + 1
     3007             
     3008             IF ( n_surf > 1 )  THEN
     3009                IF ( ANY ( surface_fraction_f%frac(:,j,i) ==                   &
     3010                     surface_fraction_f%fill ) )  THEN
     3011                   message_string = 'If more than one surface type is ' //     &
     3012                                 'given at a location, surface_fraction ' //   &
     3013                                 'must be provided.'
     3014                   CALL message( 'netcdf_data_input_mod', 'NDI027',            &
    30013015                                  2, 2, 0, 6, 0 )
    3002              ENDIF
    3003              IF ( SUM ( surface_fraction_f%frac(:,j,i) ) > 1.0_wp )  THEN
    3004                 message_string = 'surface_fraction must not exceed 1'
    3005                 CALL message( 'netcdf_data_input_mod', 'NDI028',               &
    3006                                2, 2, 0, 6, 0 )
     3016                ENDIF
     3017                IF ( SUM ( surface_fraction_f%frac(:,j,i) ) > 1.0_wp )  THEN
     3018                   message_string = 'surface_fraction must not exceed 1'
     3019                   CALL message( 'netcdf_data_input_mod', 'NDI028',            &
     3020                                  2, 2, 0, 6, 0 )
     3021                ENDIF
    30073022             ENDIF
    30083023!
    30093024!--          Check for further mismatches, e.g. vegetation_type is set but
    30103025!--          surface vegetation fraction is zero.
    3011              IF ( ( vegetation_type_f%var(j,i) /= vegetation_type_f%fill  .AND.  &
    3012                     surface_fraction_f%frac(0,j,i) == 0.0_wp )  .OR.             &
    3013                   ( pavement_type_f%var(j,i) /= pavement_type_f%fill     .AND.   &
    3014                     surface_fraction_f%frac(1,j,i) == 0.0_wp )  .OR.             &
    3015                   ( water_type_f%var(j,i) /= water_type_f%fill           .AND.   &
    3016                     surface_fraction_f%frac(2,j,i) == 0.0_wp ) )  THEN
    3017                 WRITE( message_string, * ) 'Mismatch in setting of '     //      &
    3018                                'surface_fraction. Vegetation-, pavement-, or '// &
    3019                                'water surface is given at (i,j) = ( ', i, j,     &
    3020                                ' ), but surface fraction is 0 for the given type.'
    3021                 CALL message( 'netcdf_data_input_mod', 'NDI029',                 &
     3026             IF ( ( vegetation_type_f%var(j,i) /= vegetation_type_f%fill  .AND.&
     3027                 ( surface_fraction_f%frac(0,j,i) == 0.0_wp .OR.               &
     3028                   surface_fraction_f%frac(0,j,i) == surface_fraction_f%fill ) &
     3029                  )  .OR.                                                      &
     3030                  ( pavement_type_f%var(j,i) /= pavement_type_f%fill     .AND. &
     3031                 ( surface_fraction_f%frac(1,j,i) == 0.0_wp .OR.               &
     3032                   surface_fraction_f%frac(1,j,i) == surface_fraction_f%fill ) &
     3033                  )  .OR.                                                      &
     3034                  ( water_type_f%var(j,i) /= water_type_f%fill           .AND. &
     3035                 ( surface_fraction_f%frac(2,j,i) == 0.0_wp .OR.               &
     3036                   surface_fraction_f%frac(2,j,i) == surface_fraction_f%fill ) &
     3037                  ) )  THEN
     3038                WRITE( message_string, * ) 'Mismatch in setting of '     //    &
     3039                             'surface_fraction. Vegetation-, pavement-, or '// &
     3040                             'water surface is given at (i,j) = ( ', i, j,     &
     3041                             ' ), but surface fraction is 0 for the given type.'
     3042                CALL message( 'netcdf_data_input_mod', 'NDI029',               &
    30223043                               2, 2, 0, 6, 0 )
    30233044             ENDIF
     
    30253046!--          Check for further mismatches, e.g. vegetation_type is not set       
    30263047!--          surface vegetation fraction is non-zero.
    3027              IF ( ( vegetation_type_f%var(j,i) == vegetation_type_f%fill  .AND.  &
    3028                     surface_fraction_f%frac(0,j,i) /= 0.0_wp )  .OR.             &
    3029                   ( pavement_type_f%var(j,i) == pavement_type_f%fill     .AND.   &
    3030                     surface_fraction_f%frac(1,j,i) /= 0.0_wp )  .OR.             &
    3031                   ( water_type_f%var(j,i) == water_type_f%fill           .AND.   &
    3032                     surface_fraction_f%frac(2,j,i) /= 0.0_wp ) )  THEN
    3033                 WRITE( message_string, * ) 'Mismatch in setting of '     //      &
    3034                                'surface_fraction. Vegetation-, pavement-, or '// &
    3035                                'water surface is not given at (i,j) = ( ', i, j, &
    3036                                ' ), but surface fraction is not 0 for the ' //   &
    3037                                'given type.'
    3038                 CALL message( 'netcdf_data_input_mod', 'NDI030',                 &
     3048             IF ( ( vegetation_type_f%var(j,i) == vegetation_type_f%fill  .AND.&
     3049                 ( surface_fraction_f%frac(0,j,i) /= 0.0_wp .AND.              &
     3050                   surface_fraction_f%frac(0,j,i) /= surface_fraction_f%fill ) &
     3051                  )  .OR.                                                      &
     3052                  ( pavement_type_f%var(j,i) == pavement_type_f%fill     .AND. &
     3053                 ( surface_fraction_f%frac(1,j,i) /= 0.0_wp .AND.              &
     3054                   surface_fraction_f%frac(1,j,i) /= surface_fraction_f%fill ) &
     3055                  )  .OR.                                                      &
     3056                  ( water_type_f%var(j,i) == water_type_f%fill           .AND. &
     3057                 ( surface_fraction_f%frac(2,j,i) /= 0.0_wp .AND.              &
     3058                   surface_fraction_f%frac(2,j,i) /= surface_fraction_f%fill ) &
     3059                  ) )  THEN
     3060                WRITE( message_string, * ) 'Mismatch in setting of '     //    &
     3061                             'surface_fraction. Vegetation-, pavement-, or '// &
     3062                             'water surface is not given at (i,j) = ( ', i, j, &
     3063                             ' ), but surface fraction is not 0 for the ' //   &
     3064                             'given type.'
     3065                CALL message( 'netcdf_data_input_mod', 'NDI030',               &
    30393066                               2, 2, 0, 6, 0 )
    30403067             ENDIF
Note: See TracChangeset for help on using the changeset viewer.