Changeset 3378 for palm/trunk/SOURCE/urban_surface_mod.f90
- Timestamp:
- Oct 19, 2018 12:34:59 PM (6 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/urban_surface_mod.f90 (modified) (12 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r3371 r3378 28 28 ! ----------------- 29 29 ! $Id$ 30 ! merge from radiation branch (r3362) into trunk 31 ! (moh.hefny): 32 ! - check the requested output variables if they are correct 33 ! - added unscheduled_radiation_calls switch to control force_radiation_call 34 ! - minor formate changes 35 ! 36 ! 3371 2018-10-18 13:40:12Z knoop 30 37 ! Set flag indicating that albedo at urban surfaces is already initialized 31 38 ! … … 399 406 idcsf, ndcsf, kdcsf, pct, & 400 407 startland, endland, startwall, endwall, skyvf, skyvft, nzub, & 401 nzut, n zpt, npcbl, pcbl408 nzut, npcbl, pcbl, unscheduled_radiation_calls 402 409 403 410 USE statistics, & … … 2399 2406 IMPLICIT NONE 2400 2407 2401 CHARACTER (len=*),INTENT(IN) :: variable !: 2402 CHARACTER (len=*),INTENT(OUT) :: unit !: 2403 2404 CHARACTER (len=varnamelength) :: var 2405 2408 CHARACTER(LEN=*),INTENT(IN) :: variable !: 2409 CHARACTER(LEN=*),INTENT(OUT) :: unit !: 2410 2411 INTEGER(iwp) :: i,j !< index 2412 CHARACTER(LEN=varnamelength) :: var 2413 INTEGER(iwp), PARAMETER :: nl1 = 32 !< number of directional usm variables 2414 CHARACTER(LEN=varnamelength), DIMENSION(nl1) :: varlist1 = & !< list of directional usm variables 2415 (/'usm_rad_net','usm_rad_insw','usm_rad_inlw','usm_rad_inswdir', & 2416 'usm_rad_inswdif','usm_rad_inswref','usm_rad_inlwdif','usm_wshf', & 2417 'usm_rad_inlwref','usm_rad_outsw','usm_rad_outlw','usm_rad_hf', & 2418 'usm_rad_ressw','usm_rad_reslw','usm_wghf','usm_wghf_window', & 2419 'usm_wghf_green','usm_iwghf','usm_iwghf_window','usm_surfz', & 2420 'usm_surfwintrans','usm_surfcat','usm_surfalb','usm_surfemis', & 2421 'usm_t_surf','usm_t_surf_window','usm_t_surf_green','usm_t_green', & 2422 'usm_t_surf_10cm','usm_t_wall','usm_t_window','usm_t_green'/) 2423 2424 INTEGER(iwp), PARAMETER :: nl2 = 9 !< number of other variables 2425 CHARACTER(LEN=varnamelength), DIMENSION(nl2) :: varlist2 = & !< list of other usm variables 2426 (/'usm_skyvf','usm_skyvft','usm_svf','usm_dif','usm_rad_pc_inlw', & 2427 'usm_rad_pc_insw','usm_rad_pc_inswdir','usm_rad_pc_inswdif', & 2428 'usm_rad_pc_inswref'/) 2429 2430 INTEGER(iwp), PARAMETER :: nd = 5 !< number of directions 2431 CHARACTER(LEN=6), DIMENSION(nd), PARAMETER :: dirname = & !< direction names 2432 (/'_roof ','_south','_north','_west ','_east '/) 2433 LOGICAL :: lfound !< flag if the variable is found 2434 2435 2436 lfound = .FALSE. 2406 2437 var = TRIM(variable) 2438 2439 !-- check if variable exists 2440 ! directional variables 2441 DO i = 1, nl1 2442 DO j = 1, nd 2443 IF ( TRIM(var) == TRIM(varlist1(i))//TRIM(dirname(j)) ) THEN 2444 lfound = .TRUE. 2445 EXIT 2446 ENDIF 2447 IF ( lfound ) EXIT 2448 ENDDO 2449 ENDDO 2450 IF ( lfound ) GOTO 10 2451 ! other variables 2452 DO i = 1, nl2 2453 IF ( TRIM(var) == TRIM(varlist2(i)) ) THEN 2454 lfound = .TRUE. 2455 EXIT 2456 ENDIF 2457 ENDDO 2458 IF ( .NOT. lfound ) THEN 2459 unit = 'illegal' 2460 RETURN 2461 ENDIF 2462 10 CONTINUE 2463 2464 !-- check if variable exists 2407 2465 IF ( var(1:12) == 'usm_rad_net_' .OR. var(1:13) == 'usm_rad_insw_' .OR. & 2408 2466 var(1:13) == 'usm_rad_inlw_' .OR. var(1:16) == 'usm_rad_inswdir_' .OR. & … … 2593 2651 2594 2652 CASE ( 'usm_surfz' ) 2595 !-- array of lw radiation falling to local surface after i-th reflection2653 !-- array of surface height (z) 2596 2654 IF ( idsint == iup_u ) THEN 2597 2655 DO m = 1, surf_usm_h%ns … … 5332 5390 ENDIF 5333 5391 ENDDO 5334 !!!!!!!!!!!!!HACK!!!!!!!!!!!!!!!!!!!5335 ! t_window_v_p(l)%t = t_wall_v_p(l)%t5336 ! surf_usm_v(l)%tt_window_m = surf_usm_v(l)%tt_wall_m5337 ! t_green_v_p(l)%t = t_wall_v_p(l)%t5338 ! surf_usm_v(l)%tt_green_m = surf_usm_v(l)%tt_wall_m5339 !!!!!!!!!!!!!HACK!!!!!!!!!!!!!!!!!!!5340 5392 ENDDO 5341 5393 … … 6914 6966 INTEGER(iwp) :: wtc 6915 6967 REAL(wp), DIMENSION(n_surface_params) :: wtp 6916 6917 6968 LOGICAL :: ascii_file = .FALSE. 6918 6919 6969 INTEGER(iwp), DIMENSION(0:17, nysg:nyng, nxlg:nxrg) :: usm_par 6920 6970 REAL(wp), DIMENSION(1:14, nysg:nyng, nxlg:nxrg) :: usm_val … … 7742 7792 !-- in case of fast changes in the skin temperature, it is required to 7743 7793 !-- update the radiative fluxes in order to keep the solution stable 7744 IF ( ( ABS( t_surf_h_p(m) - t_surf_h(m) ) > 1.0_wp ) .OR.&7794 IF ( ( ( ABS( t_surf_h_p(m) - t_surf_h(m) ) > 1.0_wp ) .OR. & 7745 7795 ( ABS( t_surf_green_h_p(m) - t_surf_green_h(m) ) > 1.0_wp ) .OR. & 7746 ( ABS( t_surf_window_h_p(m) - t_surf_window_h(m) ) > 1.0_wp ) ) THEN 7796 ( ABS( t_surf_window_h_p(m) - t_surf_window_h(m) ) > 1.0_wp ) ) & 7797 .AND. unscheduled_radiation_calls ) THEN 7747 7798 force_radiation_call_l = .TRUE. 7748 7799 ENDIF … … 7897 7948 7898 7949 !-- add RK3 term 7899 t_surf_v_p(l)%t(m) = t_surf_v_p(l)%t(m) + dt_3d * tsc(3) * &7950 t_surf_v_p(l)%t(m) = t_surf_v_p(l)%t(m) + dt_3d * tsc(3) * & 7900 7951 surf_usm_v(l)%tt_surface_m(m) 7901 7952 t_surf_window_v_p(l)%t(m) = t_surf_window_v_p(l)%t(m) + dt_3d * tsc(3) * & … … 7907 7958 !-- store also vpt_surface, which is, due to the lack of moisture on roofs simply 7908 7959 !-- assumed to be the surface temperature. 7909 surf_usm_v(l)%pt_surface(m) = ( surf_usm_v(l)%frac(ind_veg_wall,m) * t_surf_v_p(l)%t(m) &7910 + surf_usm_v(l)%frac(ind_wat_win,m) * t_surf_window_v_p(l)%t(m) &7960 surf_usm_v(l)%pt_surface(m) = ( surf_usm_v(l)%frac(ind_veg_wall,m) * t_surf_v_p(l)%t(m) & 7961 + surf_usm_v(l)%frac(ind_wat_win,m) * t_surf_window_v_p(l)%t(m) & 7911 7962 + surf_usm_v(l)%frac(ind_pav_green,m) * t_surf_green_v_p(l)%t(m) ) & 7912 7963 / exner(k) … … 7920 7971 stend_window = ( t_surf_window_v_p(l)%t(m) - t_surf_window_v(l)%t(m) - dt_3d * tsc(3) *& 7921 7972 surf_usm_v(l)%tt_surface_window_m(m) ) / ( dt_3d * tsc(2) ) 7922 stend_green = ( t_surf_green_v_p(l)%t(m) - t_surf_green_v(l)%t(m) - dt_3d * tsc(3) * &7973 stend_green = ( t_surf_green_v_p(l)%t(m) - t_surf_green_v(l)%t(m) - dt_3d * tsc(3) * & 7923 7974 surf_usm_v(l)%tt_surface_green_m(m) ) / ( dt_3d * tsc(2) ) 7924 7975 … … 7929 7980 surf_usm_v(l)%tt_surface_window_m(m) = stend_window 7930 7981 surf_usm_v(l)%tt_surface_green_m(m) = stend_green 7931 ELSEIF ( intermediate_timestep_count < &7982 ELSEIF ( intermediate_timestep_count < & 7932 7983 intermediate_timestep_count_max ) THEN 7933 surf_usm_v(l)%tt_surface_m(m) = -9.5625_wp * stend + &7984 surf_usm_v(l)%tt_surface_m(m) = -9.5625_wp * stend + & 7934 7985 5.3125_wp * surf_usm_v(l)%tt_surface_m(m) 7935 surf_usm_v(l)%tt_surface_green_m(m) = -9.5625_wp * stend_green + 7986 surf_usm_v(l)%tt_surface_green_m(m) = -9.5625_wp * stend_green + & 7936 7987 5.3125_wp * surf_usm_v(l)%tt_surface_green_m(m) 7937 surf_usm_v(l)%tt_surface_window_m(m) = -9.5625_wp * stend_window + 7988 surf_usm_v(l)%tt_surface_window_m(m) = -9.5625_wp * stend_window + & 7938 7989 5.3125_wp * surf_usm_v(l)%tt_surface_window_m(m) 7939 7990 ENDIF … … 7943 7994 !-- update the radiative fluxes in order to keep the solution stable 7944 7995 7945 IF ( ( ABS( t_surf_v_p(l)%t(m) - t_surf_v(l)%t(m) ) > 1.0_wp ) .OR.&7996 IF ( ( ( ABS( t_surf_v_p(l)%t(m) - t_surf_v(l)%t(m) ) > 1.0_wp ) .OR. & 7946 7997 ( ABS( t_surf_green_v_p(l)%t(m) - t_surf_green_v(l)%t(m) ) > 1.0_wp ) .OR. & 7947 ( ABS( t_surf_window_v_p(l)%t(m) - t_surf_window_v(l)%t(m) ) > 1.0_wp ) ) THEN 7998 ( ABS( t_surf_window_v_p(l)%t(m) - t_surf_window_v(l)%t(m) ) > 1.0_wp ) ) & 7999 .AND. unscheduled_radiation_calls ) THEN 7948 8000 force_radiation_call_l = .TRUE. 7949 8001 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
