Changeset 2932 for palm/trunk/SOURCE/radiation_model_mod.f90
- Timestamp:
- Mar 26, 2018 9:39:22 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/radiation_model_mod.f90
r2930 r2932 28 28 ! ----------------- 29 29 ! $Id$ 30 ! renamed radiation_par to radiation_parameters 31 ! 32 ! 2930 2018-03-23 16:30:46Z suehring 30 33 ! Remove default surfaces from radiation model, does not make much sense to 31 34 ! apply radiation model without energy-balance solvers; Further, add check for … … 2571 2574 ! Description: 2572 2575 ! ------------ 2573 !> Parin for &radiation_par for radiation model2576 !> Parin for &radiation_parameters for radiation model 2574 2577 !------------------------------------------------------------------------------! 2575 2578 SUBROUTINE radiation_parin … … 2592 2595 average_radiation, & 2593 2596 surf_reflections, svfnorm_report_thresh 2594 2597 2598 NAMELIST /radiation_parameters/ albedo, albedo_type, albedo_lw_dir, & 2599 albedo_lw_dif, albedo_sw_dir, albedo_sw_dif, & 2600 constant_albedo, dt_radiation, emissivity, & 2601 lw_radiation, net_radiation, & 2602 radiation_scheme, skip_time_do_radiation, & 2603 sw_radiation, unscheduled_radiation_calls, & 2604 split_diffusion_radiation, & 2605 max_raytracing_dist, min_irrf_value, & 2606 nrefsteps, mrt_factors, rma_lad_raytrace, & 2607 dist_max_svf, & 2608 average_radiation, & 2609 surf_reflections, svfnorm_report_thresh 2610 2595 2611 line = ' ' 2596 2612 2597 2613 ! 2598 !-- Try to find radiation model package2614 !-- Try to find radiation model namelist 2599 2615 REWIND ( 11 ) 2600 2616 line = ' ' 2601 DO WHILE ( INDEX( line, '&radiation_par ' ) == 0 )2617 DO WHILE ( INDEX( line, '&radiation_parameters' ) == 0 ) 2602 2618 READ ( 11, '(A)', END=10 ) line 2603 2619 ENDDO … … 2606 2622 ! 2607 2623 !-- Read user-defined namelist 2608 READ ( 11, radiation_par )2624 READ ( 11, radiation_parameters ) 2609 2625 2610 2626 ! 2611 2627 !-- Set flag that indicates that the radiation model is switched on 2612 2628 radiation = .TRUE. 2629 2630 GOTO 12 2631 ! 2632 !-- Try to find old namelist 2633 10 REWIND ( 11 ) 2634 line = ' ' 2635 DO WHILE ( INDEX( line, '&radiation_par' ) == 0 ) 2636 READ ( 11, '(A)', END=12 ) line 2637 ENDDO 2638 BACKSPACE ( 11 ) 2639 2640 ! 2641 !-- Read user-defined namelist 2642 READ ( 11, radiation_par ) 2643 2644 message_string = 'namelist radiation_par is deprecated and will be ' // & 2645 'removed in near future. Please &use namelist ' // & 2646 'radiation_parameters instead' 2647 CALL message( 'radiation_parin', 'PA0487', 0, 1, 0, 6, 0 ) 2648 2649 2650 ! 2651 !-- Set flag that indicates that the radiation model is switched on 2652 radiation = .TRUE. 2653 2654 12 CONTINUE 2655 2656 2613 2657 2614 2658 !-- Set radiation_interactions flag according to urban_ and land_surface flag 2615 IF ( urban_surface .OR. land_surface ) radiation_interactions = .TRUE. 2616 2617 10 CONTINUE 2659 IF ( urban_surface .OR. land_surface ) radiation_interactions = .TRUE. 2618 2660 2619 2620 2661 END SUBROUTINE radiation_parin 2621 2662 … … 5451 5492 CASE DEFAULT 5452 5493 WRITE(message_string, *) 'ERROR: the surface type ',td , ' is not supported for calculating SVF' 5453 CALL message( 'radiation_calc_svf', 'PA0 XXX', 1, 2, 0, 6, 0 )5494 CALL message( 'radiation_calc_svf', 'PA0488', 1, 2, 0, 6, 0 ) 5454 5495 END SELECT 5455 5496
Note: See TracChangeset
for help on using the changeset viewer.