Ignore:
Timestamp:
Apr 30, 2019 9:48:33 AM (6 years ago)
Author:
suehring
Message:

Move checks for correct dimensions in static input file; improve checks concerning buildings; check whether at least one surface type is set at a natural-type surface element

File:
1 edited

Legend:

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

    r3933 r3941  
    2525! -----------------
    2626! $Id$
     27! Check that at least one surface type is set at surface element.
     28!
     29! 3933 2019-04-25 12:33:20Z kanani
    2730! Remove unused subroutine and allocation of pt_2m, this is done in surface_mod
    2831! now (surfaces%pt_2m)
     
    22772280       ENDIF
    22782281
    2279 
    2280 !        pt(k+k_off,j+j_off,i+i_off) = surf_t_surface_p%var_1d(m) / exner(nzb)  !is actually no air temperature
    22812282       surf%pt_surface(m)          = surf_t_surface_p%var_1d(m) / exner(nzb)
    22822283
     
    28812882                IF ( water_type_f%var(j,i)      /= water_type_f%fill )         &
    28822883                   surf_lsm_h%water_surface(m) = .TRUE.
     2884!
     2885!--             Check if at least one type is set.
     2886                IF ( .NOT. surf_lsm_h%vegetation_surface(m)  .AND.             &
     2887                     .NOT. surf_lsm_h%pavement_surface(m)    .AND.             &
     2888                     .NOT. surf_lsm_h%water_surface(m) )  THEN
     2889                   WRITE( message_string, * ) 'Horizontal surface element ' // &
     2890                                       ' at i, j = ',  i, j,                   &
     2891                                       ' is neither a vegetation, ' //         &
     2892                                       'pavement, nor a water surface.'
     2893                   CALL message( 'land_surface_model_mod', 'PA0619',          &
     2894                                  2, 2, myid, 6, 0 )
     2895                ENDIF
     2896               
    28832897             ENDDO
    28842898!
     
    29592973                      surf_lsm_v(l)%water_surface(m) = .TRUE.
    29602974!
    2961 !--                Check if surface element has the flag %building_covered_all,
    2962 !--                this case, set vegetation_type appropriate
     2975!--                Check if at least one type is set.
     2976                   IF ( .NOT. surf_lsm_v(l)%vegetation_surface(m)  .AND.       &
     2977                        .NOT. surf_lsm_v(l)%pavement_surface(m)    .AND.       &
     2978                        .NOT. surf_lsm_v(l)%water_surface(m) )  THEN
     2979                      WRITE( message_string, * ) 'Vertical surface element ' //&
     2980                                       ' at i, j = ',  i, j,                   &
     2981                                       ' is neither a vegetation, ' //         &
     2982                                       'pavement, nor a water surface.'
     2983                      CALL message( 'land_surface_model_mod', 'PA0619',        &
     2984                                     2, 2, myid, 6, 0 )
     2985                   ENDIF
    29632986                ENDDO
    29642987             ENDDO
Note: See TracChangeset for help on using the changeset viewer.