Changeset 3941 for palm/trunk/SOURCE/land_surface_model_mod.f90
- Timestamp:
- Apr 30, 2019 9:48:33 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/land_surface_model_mod.f90
r3933 r3941 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Check that at least one surface type is set at surface element. 28 ! 29 ! 3933 2019-04-25 12:33:20Z kanani 27 30 ! Remove unused subroutine and allocation of pt_2m, this is done in surface_mod 28 31 ! now (surfaces%pt_2m) … … 2277 2280 ENDIF 2278 2281 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 temperature2281 2282 surf%pt_surface(m) = surf_t_surface_p%var_1d(m) / exner(nzb) 2282 2283 … … 2881 2882 IF ( water_type_f%var(j,i) /= water_type_f%fill ) & 2882 2883 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 2883 2897 ENDDO 2884 2898 ! … … 2959 2973 surf_lsm_v(l)%water_surface(m) = .TRUE. 2960 2974 ! 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 2963 2986 ENDDO 2964 2987 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.