Changeset 4441 for palm/trunk/SOURCE/land_surface_model_mod.f90
- Timestamp:
- Mar 4, 2020 7:20:35 PM (4 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (32 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/land_surface_model_mod.f90
r4429 r4441 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! Change order of dimension in surface arrays %frac, %emissivity and %albedo 23 ! to allow for better vectorization in the radiation interactions. 23 24 ! 24 25 ! Former revisions: … … 2502 2503 !-- fraction. 2503 2504 !-- Set default values at each surface element. 2504 ALLOCATE ( surf_lsm_h%albedo_type( 0:2,1:surf_lsm_h%ns) )2505 ALLOCATE ( surf_lsm_h%emissivity( 0:2,1:surf_lsm_h%ns) )2505 ALLOCATE ( surf_lsm_h%albedo_type(1:surf_lsm_h%ns,0:2) ) 2506 ALLOCATE ( surf_lsm_h%emissivity(1:surf_lsm_h%ns,0:2) ) 2506 2507 ! 2507 2508 !-- Initialize albedo type according to its default type, in order to set values 2508 2509 !-- independent on default albedo_type in radiation model. 2509 surf_lsm_h%albedo_type( ind_veg_wall,:) = &2510 surf_lsm_h%albedo_type(:,ind_veg_wall) = & 2510 2511 INT( vegetation_pars(ind_v_at,vegetation_type) ) 2511 surf_lsm_h%albedo_type( ind_wat_win,:) = &2512 surf_lsm_h%albedo_type(:,ind_wat_win) = & 2512 2513 INT( water_pars(ind_w_at,water_type) ) 2513 surf_lsm_h%albedo_type( ind_pav_green,:) = &2514 surf_lsm_h%albedo_type(:,ind_pav_green) = & 2514 2515 INT( pavement_pars(ind_p_at,pavement_type) ) 2515 2516 surf_lsm_h%emissivity = emissivity 2516 2517 DO l = 0, 3 2517 ALLOCATE ( surf_lsm_v(l)%albedo_type( 0:2,1:surf_lsm_v(l)%ns) )2518 ALLOCATE ( surf_lsm_v(l)%emissivity( 0:2,1:surf_lsm_v(l)%ns) )2518 ALLOCATE ( surf_lsm_v(l)%albedo_type(1:surf_lsm_v(l)%ns,0:2) ) 2519 ALLOCATE ( surf_lsm_v(l)%emissivity(1:surf_lsm_v(l)%ns,0:2) ) 2519 2520 ! 2520 2521 !-- Initialize albedo type according to its default type, in order to 2521 2522 !-- set values independent on default albedo_type in radiation model. 2522 surf_lsm_v(l)%albedo_type( ind_veg_wall,:) = &2523 surf_lsm_v(l)%albedo_type(:,ind_veg_wall) = & 2523 2524 INT( vegetation_pars(ind_v_at,vegetation_type) ) 2524 surf_lsm_v(l)%albedo_type( ind_wat_win,:) = &2525 surf_lsm_v(l)%albedo_type(:,ind_wat_win) = & 2525 2526 INT( water_pars(ind_w_at,water_type) ) 2526 surf_lsm_v(l)%albedo_type( ind_pav_green,:) = &2527 surf_lsm_v(l)%albedo_type(:,ind_pav_green) = & 2527 2528 INT( pavement_pars(ind_p_at,pavement_type) ) 2528 2529 surf_lsm_v(l)%emissivity = emissivity … … 2531 2532 !-- Allocate arrays for relative surface fraction. 2532 2533 !-- 0 - vegetation fraction, 2 - water fraction, 1 - pavement fraction 2533 ALLOCATE( surf_lsm_h%frac( 0:2,1:surf_lsm_h%ns) )2534 ALLOCATE( surf_lsm_h%frac(1:surf_lsm_h%ns,0:2) ) 2534 2535 surf_lsm_h%frac = 0.0_wp 2535 2536 DO l = 0, 3 2536 ALLOCATE( surf_lsm_v(l)%frac( 0:2,1:surf_lsm_v(l)%ns) )2537 ALLOCATE( surf_lsm_v(l)%frac(1:surf_lsm_v(l)%ns,0:2) ) 2537 2538 surf_lsm_v(l)%frac = 0.0_wp 2538 2539 ENDDO … … 2588 2589 2589 2590 surf_lsm_h%vegetation_surface = .TRUE. 2590 surf_lsm_h%frac( ind_veg_wall,:) = 1.0_wp2591 surf_lsm_h%frac(:,ind_veg_wall) = 1.0_wp 2591 2592 DO l = 0, 3 2592 2593 surf_lsm_v(l)%vegetation_surface = .TRUE. 2593 surf_lsm_v(l)%frac( ind_veg_wall,:) = 1.0_wp2594 surf_lsm_v(l)%frac(:,ind_veg_wall) = 1.0_wp 2594 2595 ENDDO 2595 2596 … … 2597 2598 2598 2599 surf_lsm_h%water_surface = .TRUE. 2599 surf_lsm_h%frac( ind_wat_win,:) = 1.0_wp2600 surf_lsm_h%frac(:,ind_wat_win) = 1.0_wp 2600 2601 ! 2601 2602 !-- Note, vertical water surface does not really make sense. 2602 2603 DO l = 0, 3 2603 2604 surf_lsm_v(l)%water_surface = .TRUE. 2604 surf_lsm_v(l)%frac( ind_wat_win,:) = 1.0_wp2605 surf_lsm_v(l)%frac(:,ind_wat_win) = 1.0_wp 2605 2606 ENDDO 2606 2607 … … 2608 2609 2609 2610 surf_lsm_h%pavement_surface = .TRUE. 2610 surf_lsm_h%frac( ind_pav_green,:) = 1.0_wp2611 surf_lsm_h%frac(:,ind_pav_green) = 1.0_wp 2611 2612 DO l = 0, 3 2612 2613 surf_lsm_v(l)%pavement_surface = .TRUE. 2613 surf_lsm_v(l)%frac( ind_pav_green,:) = 1.0_wp2614 surf_lsm_v(l)%frac(:,ind_pav_green) = 1.0_wp 2614 2615 ENDDO 2615 2616 … … 2745 2746 IF ( surface_fraction_f%frac(ind_veg_wall,j,i) /= & 2746 2747 surface_fraction_f%fill ) THEN 2747 surf_lsm_h%frac( ind_veg_wall,m) = &2748 surf_lsm_h%frac(m,ind_veg_wall) = & 2748 2749 surface_fraction_f%frac(ind_veg_wall,j,i) 2749 2750 ENDIF 2750 2751 IF ( surface_fraction_f%frac(ind_pav_green,j,i) /= & 2751 2752 surface_fraction_f%fill ) THEN 2752 surf_lsm_h%frac( ind_pav_green,m) = &2753 surf_lsm_h%frac(m,ind_pav_green) = & 2753 2754 surface_fraction_f%frac(ind_pav_green,j,i) 2754 2755 ENDIF 2755 2756 IF ( surface_fraction_f%frac(ind_wat_win,j,i) /= & 2756 2757 surface_fraction_f%fill ) THEN 2757 surf_lsm_h%frac( ind_wat_win,m) = &2758 surf_lsm_h%frac(m,ind_wat_win) = & 2758 2759 surface_fraction_f%frac(ind_wat_win,j,i) 2759 2760 ENDIF … … 2761 2762 !-- Check if sum of relative fractions is zero. This case, give an 2762 2763 !-- error message. 2763 IF ( SUM ( surf_lsm_h%frac( :,m) ) == 0.0_wp ) THEN2764 IF ( SUM ( surf_lsm_h%frac(m,:) ) == 0.0_wp ) THEN 2764 2765 WRITE( message_string, * ) & 2765 2766 'surface fractions at grid point (j,i) = (', & … … 2773 2774 !-- fractions to one. Note, at the moment no tile approach is 2774 2775 !-- implemented, so that relative fractions are either 1 or zero. 2775 IF ( SUM ( surf_lsm_h%frac( :,m) ) > 1.0_wp .OR. &2776 SUM ( surf_lsm_h%frac( :,m) ) < 1.0_wp ) THEN2777 surf_lsm_h%frac( :,m) = surf_lsm_h%frac(:,m) / &2778 SUM ( surf_lsm_h%frac( :,m) )2776 IF ( SUM ( surf_lsm_h%frac(m,:) ) > 1.0_wp .OR. & 2777 SUM ( surf_lsm_h%frac(m,:) ) < 1.0_wp ) THEN 2778 surf_lsm_h%frac(m,:) = surf_lsm_h%frac(m,:) / & 2779 SUM ( surf_lsm_h%frac(m,:) ) 2779 2780 2780 2781 ENDIF … … 2791 2792 IF ( surface_fraction_f%frac(ind_veg_wall,j,i) /= & 2792 2793 surface_fraction_f%fill ) THEN 2793 surf_lsm_v(l)%frac( ind_veg_wall,m) = &2794 surf_lsm_v(l)%frac(m,ind_veg_wall) = & 2794 2795 surface_fraction_f%frac(ind_veg_wall,j,i) 2795 2796 ENDIF 2796 2797 IF ( surface_fraction_f%frac(ind_pav_green,j,i) /= & 2797 2798 surface_fraction_f%fill ) THEN 2798 surf_lsm_v(l)%frac( ind_pav_green,m) = &2799 surf_lsm_v(l)%frac(m,ind_pav_green) = & 2799 2800 surface_fraction_f%frac(ind_pav_green,j,i) 2800 2801 ENDIF 2801 2802 IF ( surface_fraction_f%frac(ind_wat_win,j,i) /= & 2802 2803 surface_fraction_f%fill ) THEN 2803 surf_lsm_v(l)%frac( ind_wat_win,m) = &2804 surf_lsm_v(l)%frac(m,ind_wat_win) = & 2804 2805 surface_fraction_f%frac(ind_wat_win,j,i) 2805 2806 ENDIF … … 2807 2808 !-- Check if sum of relative fractions is zero. This case, give an 2808 2809 !-- error message. 2809 IF ( SUM ( surf_lsm_v(l)%frac( :,m) ) == 0.0_wp ) THEN2810 IF ( SUM ( surf_lsm_v(l)%frac(m,:) ) == 0.0_wp ) THEN 2810 2811 WRITE( message_string, * ) & 2811 2812 'surface fractions at grid point (j,i) = (', & … … 2819 2820 !-- fractions to one. Note, at the moment no tile approach is 2820 2821 !-- implemented, so that relative fractions are either 1 or zero. 2821 IF ( SUM ( surf_lsm_v(l)%frac( :,m) ) > 1.0_wp .OR. &2822 SUM ( surf_lsm_v(l)%frac( :,m) ) < 1.0_wp ) THEN2823 surf_lsm_v(l)%frac( :,m) = surf_lsm_v(l)%frac(:,m) / &2824 SUM ( surf_lsm_v(l)%frac( :,m) )2822 IF ( SUM ( surf_lsm_v(l)%frac(m,:) ) > 1.0_wp .OR. & 2823 SUM ( surf_lsm_v(l)%frac(m,:) ) < 1.0_wp ) THEN 2824 surf_lsm_v(l)%frac(m,:) = surf_lsm_v(l)%frac(m,:) / & 2825 SUM ( surf_lsm_v(l)%frac(m,:) ) 2825 2826 2826 2827 ENDIF … … 2834 2835 2835 2836 IF ( vegetation_type_f%var(j,i) /= vegetation_type_f%fill ) & 2836 surf_lsm_h%frac( ind_veg_wall,m) = 1.0_wp2837 surf_lsm_h%frac(m,ind_veg_wall) = 1.0_wp 2837 2838 IF ( pavement_type_f%var(j,i) /= pavement_type_f%fill ) & 2838 surf_lsm_h%frac( ind_pav_green,m) = 1.0_wp2839 surf_lsm_h%frac(m,ind_pav_green) = 1.0_wp 2839 2840 IF ( water_type_f%var(j,i) /= water_type_f%fill ) & 2840 surf_lsm_h%frac( ind_wat_win,m) = 1.0_wp2841 surf_lsm_h%frac(m,ind_wat_win) = 1.0_wp 2841 2842 ENDDO 2842 2843 DO l = 0, 3 … … 2848 2849 2849 2850 IF ( vegetation_type_f%var(j,i) /= vegetation_type_f%fill ) & 2850 surf_lsm_v(l)%frac( ind_veg_wall,m) = 1.0_wp2851 surf_lsm_v(l)%frac(m,ind_veg_wall) = 1.0_wp 2851 2852 IF ( pavement_type_f%var(j,i) /= pavement_type_f%fill ) & 2852 surf_lsm_v(l)%frac( ind_pav_green,m) = 1.0_wp2853 surf_lsm_v(l)%frac(m,ind_pav_green) = 1.0_wp 2853 2854 IF ( water_type_f%var(j,i) /= water_type_f%fill ) & 2854 surf_lsm_v(l)%frac( ind_wat_win,m) = 1.0_wp2855 surf_lsm_v(l)%frac(m,ind_wat_win) = 1.0_wp 2855 2856 ENDDO 2856 2857 ENDDO … … 3218 3219 surf_lsm_h%f_sw_in(m) = f_shortwave_incoming 3219 3220 surf_lsm_h%c_surface(m) = c_surface 3220 surf_lsm_h%albedo_type( ind_veg_wall,m) = albedo_type3221 surf_lsm_h%emissivity( ind_veg_wall,m) = emissivity3221 surf_lsm_h%albedo_type(m,ind_veg_wall) = albedo_type 3222 surf_lsm_h%emissivity(m,ind_veg_wall) = emissivity 3222 3223 3223 3224 surf_lsm_h%vegetation_type(m) = vegetation_type … … 3247 3248 surf_lsm_v(l)%f_sw_in(m) = f_shortwave_incoming 3248 3249 surf_lsm_v(l)%c_surface(m) = c_surface 3249 surf_lsm_v(l)%albedo_type( ind_veg_wall,m) = albedo_type3250 surf_lsm_v(l)%emissivity( ind_veg_wall,m) = emissivity3250 surf_lsm_v(l)%albedo_type(m,ind_veg_wall) = albedo_type 3251 surf_lsm_v(l)%emissivity(m,ind_veg_wall) = emissivity 3251 3252 3252 3253 surf_lsm_v(l)%vegetation_type(m) = vegetation_type … … 3285 3286 surf_lsm_h%f_sw_in(m) = vegetation_pars(ind_v_f_sw_in,st) 3286 3287 surf_lsm_h%c_surface(m) = vegetation_pars(ind_v_c_surf,st) 3287 surf_lsm_h%albedo_type( ind_veg_wall,m) = INT( vegetation_pars(ind_v_at,st) )3288 surf_lsm_h%emissivity( ind_veg_wall,m) = vegetation_pars(ind_v_emis,st)3288 surf_lsm_h%albedo_type(m,ind_veg_wall) = INT( vegetation_pars(ind_v_at,st) ) 3289 surf_lsm_h%emissivity(m,ind_veg_wall) = vegetation_pars(ind_v_emis,st) 3289 3290 3290 3291 surf_lsm_h%vegetation_type(m) = st … … 3314 3315 surf_lsm_v(l)%f_sw_in(m) = vegetation_pars(ind_v_f_sw_in,st) 3315 3316 surf_lsm_v(l)%c_surface(m) = vegetation_pars(ind_v_c_surf,st) 3316 surf_lsm_v(l)%albedo_type( ind_veg_wall,m) = INT( vegetation_pars(ind_v_at,st) )3317 surf_lsm_v(l)%emissivity( ind_veg_wall,m) = vegetation_pars(ind_v_emis,st)3317 surf_lsm_v(l)%albedo_type(m,ind_veg_wall) = INT( vegetation_pars(ind_v_at,st) ) 3318 surf_lsm_v(l)%emissivity(m,ind_veg_wall) = vegetation_pars(ind_v_emis,st) 3318 3319 3319 3320 surf_lsm_v(l)%vegetation_type(m) = st … … 3393 3394 IF ( vegetation_pars_f%pars_xy(ind_v_at,j,i) /= & 3394 3395 vegetation_pars_f%fill ) & 3395 surf_lsm_h%albedo_type( ind_veg_wall,m) = &3396 surf_lsm_h%albedo_type(m,ind_veg_wall) = & 3396 3397 INT( vegetation_pars_f%pars_xy(ind_v_at,j,i) ) 3397 3398 IF ( vegetation_pars_f%pars_xy(ind_v_emis,j,i) /= & 3398 3399 vegetation_pars_f%fill ) & 3399 surf_lsm_h%emissivity( ind_veg_wall,m) = &3400 surf_lsm_h%emissivity(m,ind_veg_wall) = & 3400 3401 vegetation_pars_f%pars_xy(ind_v_emis,j,i) 3401 3402 ENDIF … … 3469 3470 IF ( vegetation_pars_f%pars_xy(ind_v_at,j,i) /= & 3470 3471 vegetation_pars_f%fill ) & 3471 surf_lsm_v(l)%albedo_type( ind_veg_wall,m) = &3472 surf_lsm_v(l)%albedo_type(m,ind_veg_wall) = & 3472 3473 INT( vegetation_pars_f%pars_xy(ind_v_at,j,i) ) 3473 3474 IF ( vegetation_pars_f%pars_xy(ind_v_emis,j,i) /= & 3474 3475 vegetation_pars_f%fill ) & 3475 surf_lsm_v(l)%emissivity( ind_veg_wall,m) = &3476 surf_lsm_v(l)%emissivity(m,ind_veg_wall) = & 3476 3477 vegetation_pars_f%pars_xy(ind_v_emis,j,i) 3477 3478 ENDIF … … 3523 3524 surf_lsm_h%lambda_surface_u(m) = 1.0E10_wp 3524 3525 surf_lsm_h%c_surface(m) = 0.0_wp 3525 surf_lsm_h%albedo_type( ind_wat_win,m) = albedo_type3526 surf_lsm_h%emissivity( ind_wat_win,m) = emissivity3526 surf_lsm_h%albedo_type(m,ind_wat_win) = albedo_type 3527 surf_lsm_h%emissivity(m,ind_wat_win) = emissivity 3527 3528 3528 3529 surf_lsm_h%water_type(m) = water_type … … 3544 3545 surf_lsm_v(l)%lambda_surface_u(m) = 1.0E10_wp 3545 3546 surf_lsm_v(l)%c_surface(m) = 0.0_wp 3546 surf_lsm_v(l)%albedo_type( ind_wat_win,m) = albedo_type3547 surf_lsm_v(l)%emissivity( ind_wat_win,m) = emissivity3547 surf_lsm_v(l)%albedo_type(m,ind_wat_win) = albedo_type 3548 surf_lsm_v(l)%emissivity(m,ind_wat_win) = emissivity 3548 3549 3549 3550 surf_lsm_v(l)%water_type(m) = water_type … … 3577 3578 surf_lsm_h%lambda_surface_u(m) = water_pars(ind_w_lambda_u,st) 3578 3579 surf_lsm_h%c_surface(m) = 0.0_wp 3579 surf_lsm_h%albedo_type( ind_wat_win,m) = INT( water_pars(ind_w_at,st) )3580 surf_lsm_h%emissivity( ind_wat_win,m) = water_pars(ind_w_emis,st)3580 surf_lsm_h%albedo_type(m,ind_wat_win) = INT( water_pars(ind_w_at,st) ) 3581 surf_lsm_h%emissivity(m,ind_wat_win) = water_pars(ind_w_emis,st) 3581 3582 3582 3583 surf_lsm_h%water_type(m) = st … … 3605 3606 water_pars(ind_w_lambda_u,st) 3606 3607 surf_lsm_v(l)%c_surface(m) = 0.0_wp 3607 surf_lsm_v(l)%albedo_type( ind_wat_win,m) = &3608 surf_lsm_v(l)%albedo_type(m,ind_wat_win) = & 3608 3609 INT( water_pars(ind_w_at,st) ) 3609 surf_lsm_v(l)%emissivity( ind_wat_win,m) = &3610 surf_lsm_v(l)%emissivity(m,ind_wat_win) = & 3610 3611 water_pars(ind_w_emis,st) 3611 3612 … … 3664 3665 IF ( water_pars_f%pars_xy(ind_w_at,j,i) /= & 3665 3666 water_pars_f%fill ) & 3666 surf_lsm_h%albedo_type( ind_wat_win,m) = &3667 surf_lsm_h%albedo_type(m,ind_wat_win) = & 3667 3668 INT( water_pars_f%pars_xy(ind_w_at,j,i) ) 3668 3669 3669 3670 IF ( water_pars_f%pars_xy(ind_w_emis,j,i) /= & 3670 3671 water_pars_f%fill ) & 3671 surf_lsm_h%emissivity( ind_wat_win,m) = &3672 surf_lsm_h%emissivity(m,ind_wat_win) = & 3672 3673 water_pars_f%pars_xy(ind_w_emis,j,i) 3673 3674 ENDIF … … 3724 3725 IF ( water_pars_f%pars_xy(ind_w_at,j,i) /= & 3725 3726 water_pars_f%fill ) & 3726 surf_lsm_v(l)%albedo_type( ind_wat_win,m) = &3727 surf_lsm_v(l)%albedo_type(m,ind_wat_win) = & 3727 3728 INT( water_pars_f%pars_xy(ind_w_at,j,i) ) 3728 3729 3729 3730 IF ( water_pars_f%pars_xy(ind_w_emis,j,i) /= & 3730 3731 water_pars_f%fill ) & 3731 surf_lsm_v(l)%emissivity( ind_wat_win,m) = &3732 surf_lsm_v(l)%emissivity(m,ind_wat_win) = & 3732 3733 water_pars_f%pars_xy(ind_w_emis,j,i) 3733 3734 ENDIF … … 3810 3811 * 0.25_wp 3811 3812 3812 surf_lsm_h%albedo_type( ind_pav_green,m) = albedo_type3813 surf_lsm_h%emissivity( ind_pav_green,m) = emissivity3813 surf_lsm_h%albedo_type(m,ind_pav_green) = albedo_type 3814 surf_lsm_h%emissivity(m,ind_pav_green) = emissivity 3814 3815 3815 3816 surf_lsm_h%pavement_type(m) = pavement_type … … 3848 3849 * 0.25_wp 3849 3850 3850 surf_lsm_v(l)%albedo_type( ind_pav_green,m) = albedo_type3851 surf_lsm_v(l)%emissivity( ind_pav_green,m) = emissivity3851 surf_lsm_v(l)%albedo_type(m,ind_pav_green) = albedo_type 3852 surf_lsm_v(l)%emissivity(m,ind_pav_green) = emissivity 3852 3853 3853 3854 surf_lsm_v(l)%pavement_type(m) = pavement_type … … 3908 3909 * dz_soil(nzb_soil) & 3909 3910 * 0.25_wp 3910 surf_lsm_h%albedo_type( ind_pav_green,m) = INT( pavement_pars(ind_p_at,st) )3911 surf_lsm_h%emissivity( ind_pav_green,m) = pavement_pars(ind_p_emis,st)3911 surf_lsm_h%albedo_type(m,ind_pav_green) = INT( pavement_pars(ind_p_at,st) ) 3912 surf_lsm_h%emissivity(m,ind_pav_green) = pavement_pars(ind_p_emis,st) 3912 3913 3913 3914 surf_lsm_h%pavement_type(m) = st … … 3961 3962 * dz_soil(nzb_soil) & 3962 3963 * 0.25_wp 3963 surf_lsm_v(l)%albedo_type( ind_pav_green,m) = &3964 surf_lsm_v(l)%albedo_type(m,ind_pav_green) = & 3964 3965 INT( pavement_pars(ind_p_at,st) ) 3965 surf_lsm_v(l)%emissivity( ind_pav_green,m) = &3966 surf_lsm_v(l)%emissivity(m,ind_pav_green) = & 3966 3967 pavement_pars(ind_p_emis,st) 3967 3968 … … 4030 4031 IF ( pavement_pars_f%pars_xy(ind_p_at,j,i) /= & 4031 4032 pavement_pars_f%fill ) & 4032 surf_lsm_h%albedo_type( ind_pav_green,m) = &4033 surf_lsm_h%albedo_type(m,ind_pav_green) = & 4033 4034 INT( pavement_pars_f%pars_xy(ind_p_at,j,i) ) 4034 4035 IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /= & 4035 4036 pavement_pars_f%fill ) & 4036 surf_lsm_h%emissivity( ind_pav_green,m) = &4037 surf_lsm_h%emissivity(m,ind_pav_green) = & 4037 4038 pavement_pars_f%pars_xy(ind_p_emis,j,i) 4038 4039 ENDIF … … 4090 4091 IF ( pavement_pars_f%pars_xy(ind_p_at,j,i) /= & 4091 4092 pavement_pars_f%fill ) & 4092 surf_lsm_v(l)%albedo_type( ind_pav_green,m) = &4093 surf_lsm_v(l)%albedo_type(m,ind_pav_green) = & 4093 4094 INT( pavement_pars_f%pars_xy(ind_p_at,j,i) ) 4094 4095 4095 4096 IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /= & 4096 4097 pavement_pars_f%fill ) & 4097 surf_lsm_v(l)%emissivity( ind_pav_green,m) = &4098 surf_lsm_v(l)%emissivity(m,ind_pav_green) = & 4098 4099 pavement_pars_f%pars_xy(ind_p_emis,j,i) 4099 4100 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
