Changeset 2746
- Timestamp:
- Jan 15, 2018 12:06:04 PM (6 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 12 edited
-
check_parameters.f90 (modified) (2 diffs)
-
data_output_3d.f90 (modified) (3 diffs)
-
header.f90 (modified) (2 diffs)
-
init_3d_model.f90 (modified) (2 diffs)
-
modules.f90 (modified) (2 diffs)
-
netcdf_data_input_mod.f90 (modified) (5 diffs)
-
netcdf_interface_mod.f90 (modified) (3 diffs)
-
plant_canopy_model_mod.f90 (modified) (4 diffs)
-
prognostic_equations.f90 (modified) (3 diffs)
-
radiation_model_mod.f90 (modified) (4 diffs)
-
turbulence_closure_mod.f90 (modified) (3 diffs)
-
urban_surface_mod.f90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/check_parameters.f90
r2743 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2743 2018-01-12 16:03:39Z suehring 27 30 ! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2. 28 31 ! … … 608 611 USE pegrid 609 612 USE plant_canopy_model_mod, & 610 ONLY: pcm_check_data_output, pcm_check_parameters , plant_canopy613 ONLY: pcm_check_data_output, pcm_check_parameters 611 614 612 615 USE pmc_interface, & -
palm/trunk/SOURCE/data_output_3d.f90
r2718 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 190 193 ONLY: air_chemistry, cloud_physics, do3d, do3d_no, do3d_time_count, & 191 194 io_blocks, io_group, land_surface, message_string, & 192 ntdim_3d, nz_do3d, 195 ntdim_3d, nz_do3d, plant_canopy, & 193 196 psolver, simulated_time, time_since_reference_point, & 194 197 urban_surface, varnamelength … … 220 223 221 224 USE plant_canopy_model_mod, & 222 ONLY: pcm_data_output_3d , plant_canopy225 ONLY: pcm_data_output_3d 223 226 224 227 USE radiation_model_mod, & -
palm/trunk/SOURCE/header.f90
r2718 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 405 408 406 409 USE plant_canopy_model_mod, & 407 ONLY: pcm_header , plant_canopy410 ONLY: pcm_header 408 411 409 412 USE pmc_handle_communicator, & -
palm/trunk/SOURCE/init_3d_model.f90
r2718 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 455 458 456 459 USE plant_canopy_model_mod, & 457 ONLY: pcm_init , plant_canopy460 ONLY: pcm_init 458 461 459 462 USE radiation_model_mod, & -
palm/trunk/SOURCE/modules.f90
r2742 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! +plant_canopy 28 ! 29 ! 2742 2018-01-12 14:59:47Z suehring 27 30 ! +tsurf_av 28 31 ! … … 1268 1271 LOGICAL :: outflow_s = .FALSE. !< south domain boundary has non-cyclic outflow? 1269 1272 LOGICAL :: passive_scalar = .FALSE. !< namelist parameter 1273 LOGICAL :: plant_canopy = .FALSE. !< switch for use of plant canopy model 1270 1274 LOGICAL :: precipitation = .FALSE. !< namelist parameter 1271 1275 LOGICAL :: random_heatflux = .FALSE. !< namelist parameter -
palm/trunk/SOURCE/netcdf_data_input_mod.f90
r2718 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Read plant-canopy variables independently on land-surface model usage 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 497 500 USE control_parameters, & 498 501 ONLY: bc_lr_cyc, bc_ns_cyc, land_surface, message_string, & 499 urban_surface502 plant_canopy, urban_surface 500 503 501 504 USE indices, & 502 505 ONLY: nbgp, nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg 506 503 507 504 508 IMPLICIT NONE … … 525 529 ! 526 530 !-- If not static input file is available, skip this routine 527 IF ( .NOT. input_pids_static ) RETURN 528 ! 529 !-- Moreover, skip routine if no land-surface or urban-surface module are 530 !-- applied as no variable is used anyway. 531 IF ( .NOT. input_pids_static ) RETURN 532 ! 533 !-- Read plant canopy variables. 534 IF ( plant_canopy ) THEN 535 DO ii = 0, io_blocks-1 536 IF ( ii == io_group ) THEN 537 #if defined ( __netcdf ) 538 ! 539 !-- Open file in read-only mode 540 CALL open_read_file( TRIM( input_file_static ) // & 541 TRIM( coupling_char ) , id_surf ) 542 ! 543 !-- At first, inquire all variable names. 544 !-- This will be used to check whether an optional input variable 545 !-- exist or not. 546 CALL inquire_num_variables( id_surf, num_vars ) 547 548 ALLOCATE( var_names(1:num_vars) ) 549 CALL inquire_variable_names( id_surf, var_names ) 550 551 552 ! 553 !-- Read leaf area density - resolved vegetation 554 IF ( check_existence( var_names, 'leaf_area_density' ) ) THEN 555 leaf_area_density_f%from_file = .TRUE. 556 CALL get_attribute( id_surf, char_fill, & 557 leaf_area_density_f%fill, & 558 .FALSE., 'leaf_area_density' ) 559 ! 560 !-- Inquire number of vertical vegetation layer 561 CALL get_dimension_length( id_surf, leaf_area_density_f%nz, & 562 'zlad' ) 563 ! 564 !-- Allocate variable for leaf-area density 565 ALLOCATE( leaf_area_density_f%var( & 566 0:leaf_area_density_f%nz-1, & 567 nys:nyn,nxl:nxr) ) 568 569 DO i = nxl, nxr 570 DO j = nys, nyn 571 CALL get_variable( id_surf, 'leaf_area_density', & 572 i, j, & 573 leaf_area_density_f%var(:,j,i) ) 574 ENDDO 575 ENDDO 576 ELSE 577 leaf_area_density_f%from_file = .FALSE. 578 ENDIF 579 580 ! 581 !-- Read basal area density - resolved vegetation 582 IF ( check_existence( var_names, 'basal_area_density' ) ) THEN 583 basal_area_density_f%from_file = .TRUE. 584 CALL get_attribute( id_surf, char_fill, & 585 basal_area_density_f%fill, & 586 .FALSE., 'basal_area_density' ) 587 ! 588 !-- Inquire number of vertical vegetation layer 589 CALL get_dimension_length( id_surf, & 590 basal_area_density_f%nz, & 591 'zlad' ) 592 ! 593 !-- Allocate variable 594 ALLOCATE( basal_area_density_f%var( & 595 0:basal_area_density_f%nz-1, & 596 nys:nyn,nxl:nxr) ) 597 598 DO i = nxl, nxr 599 DO j = nys, nyn 600 CALL get_variable( id_surf, 'basal_area_density', & 601 i, j, & 602 basal_area_density_f%var(:,j,i) ) 603 ENDDO 604 ENDDO 605 ELSE 606 basal_area_density_f%from_file = .FALSE. 607 ENDIF 608 609 ! 610 !-- Read root area density - resolved vegetation 611 IF ( check_existence( var_names, 'root_area_density_lad' ) ) THEN 612 root_area_density_lad_f%from_file = .TRUE. 613 CALL get_attribute( id_surf, char_fill, & 614 root_area_density_lad_f%fill, & 615 .FALSE., 'root_area_density_lad' ) 616 ! 617 !-- Inquire number of vertical soil layers 618 CALL get_dimension_length( id_surf, & 619 root_area_density_lad_f%nz, & 620 'zsoil' ) 621 ! 622 !-- Allocate variable 623 ALLOCATE( root_area_density_lad_f%var & 624 (0:root_area_density_lad_f%nz-1,& 625 nys:nyn,nxl:nxr) ) 626 627 DO i = nxl, nxr 628 DO j = nys, nyn 629 CALL get_variable( id_surf, 'root_area_density_lad', & 630 i, j, & 631 root_area_density_lad_f%var(:,j,i) ) 632 ENDDO 633 ENDDO 634 ELSE 635 root_area_density_lad_f%from_file = .FALSE. 636 ENDIF 637 ! 638 !-- Finally, close input file 639 CALL close_input_file( id_surf ) 640 #endif 641 ENDIF 642 #if defined( __parallel ) 643 CALL MPI_BARRIER( comm2d, ierr ) 644 #endif 645 ENDDO 646 ! 647 !-- Deallocate variable list. Will be re-allocated in case further 648 !-- variables are read from file. 649 IF ( ALLOCATED( var_names ) ) DEALLOCATE( var_names ) 650 651 ENDIF 652 ! 653 !-- Skip the following if no land-surface or urban-surface module are 654 !-- applied. This case, no one of the following variables is used anyway. 531 655 IF ( .NOT. land_surface .OR. .NOT. urban_surface ) RETURN 532 656 ! … … 544 668 545 669 ! 546 !-- At first, inquire all variable names.670 !-- Inquire all variable names. 547 671 !-- This will be used to check whether an optional input variable exist 548 672 !-- or not. … … 953 1077 water_pars_f%from_file = .FALSE. 954 1078 ENDIF 955 956 !957 !-- Read leaf area density - resolved vegetation958 IF ( check_existence( var_names, 'leaf_area_density' ) ) THEN959 leaf_area_density_f%from_file = .TRUE.960 CALL get_attribute( id_surf, char_fill, &961 leaf_area_density_f%fill, &962 .FALSE., 'leaf_area_density' )963 !964 !-- Inquire number of vertical vegetation layer965 CALL get_dimension_length( id_surf, leaf_area_density_f%nz, &966 'zlad' )967 !968 !-- Allocate variable for leaf-area density969 ALLOCATE( leaf_area_density_f%var(0:leaf_area_density_f%nz-1, &970 nys:nyn,nxl:nxr) )971 972 DO i = nxl, nxr973 DO j = nys, nyn974 CALL get_variable( id_surf, 'leaf_area_density', i, j, &975 leaf_area_density_f%var(:,j,i) )976 ENDDO977 ENDDO978 ELSE979 leaf_area_density_f%from_file = .FALSE.980 ENDIF981 982 !983 !-- Read basal area density - resolved vegetation984 IF ( check_existence( var_names, 'basal_area_density' ) ) THEN985 basal_area_density_f%from_file = .TRUE.986 CALL get_attribute( id_surf, char_fill, &987 basal_area_density_f%fill, &988 .FALSE., 'basal_area_density' )989 !990 !-- Inquire number of vertical vegetation layer991 CALL get_dimension_length( id_surf, basal_area_density_f%nz, &992 'zlad' )993 !994 !-- Allocate variable995 ALLOCATE( basal_area_density_f%var(0:basal_area_density_f%nz-1,&996 nys:nyn,nxl:nxr) )997 998 DO i = nxl, nxr999 DO j = nys, nyn1000 CALL get_variable( id_surf, 'basal_area_density', i, j, &1001 basal_area_density_f%var(:,j,i) )1002 ENDDO1003 ENDDO1004 ELSE1005 basal_area_density_f%from_file = .FALSE.1006 ENDIF1007 1008 !1009 !-- Read root area density - resolved vegetation1010 IF ( check_existence( var_names, 'root_area_density_lad' ) ) THEN1011 root_area_density_lad_f%from_file = .TRUE.1012 CALL get_attribute( id_surf, char_fill, &1013 root_area_density_lad_f%fill, &1014 .FALSE., 'root_area_density_lad' )1015 !1016 !-- Inquire number of vertical soil layers1017 CALL get_dimension_length( id_surf, root_area_density_lad_f%nz,&1018 'zsoil' )1019 !1020 !-- Allocate variable1021 ALLOCATE( root_area_density_lad_f%var &1022 (0:root_area_density_lad_f%nz-1, &1023 nys:nyn,nxl:nxr) )1024 1025 DO i = nxl, nxr1026 DO j = nys, nyn1027 CALL get_variable( id_surf, 'root_area_density_lad', i, j,&1028 root_area_density_lad_f%var(:,j,i) )1029 ENDDO1030 ENDDO1031 ELSE1032 root_area_density_lad_f%from_file = .FALSE.1033 ENDIF1034 1035 1079 ! 1036 1080 !-- Read root area density - parametrized vegetation -
palm/trunk/SOURCE/netcdf_interface_mod.f90
r2718 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 496 499 mask_size_l, mask_i, mask_i_global, mask_j, mask_j_global, & 497 500 mask_k_global, message_string, mid, ntdim_2d_xy, & 498 ntdim_2d_xz, ntdim_2d_yz, ntdim_3d, nz_do3d, p rt_time_count,&499 run_description_header, section, simulated_time,&501 ntdim_2d_xz, ntdim_2d_yz, ntdim_3d, nz_do3d, plant_canopy, & 502 prt_time_count, run_description_header, section, simulated_time,& 500 503 simulated_time_at_begin, skip_time_data_output_av, & 501 504 skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz, & … … 520 523 521 524 USE plant_canopy_model_mod, & 522 ONLY: pcm_define_netcdf_grid , plant_canopy525 ONLY: pcm_define_netcdf_grid 523 526 524 527 USE profil_parameter, & -
palm/trunk/SOURCE/plant_canopy_model_mod.f90
r2718 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 175 178 176 179 LOGICAL :: calc_beta_lad_profile = .FALSE. !< switch for calc. of lad from beta func. 177 LOGICAL :: plant_canopy = .FALSE. !< switch for use of canopy model178 180 LOGICAL :: usm_lad_rma = .TRUE. !< use MPI RMA to access LAD for raytracing (instead of global array) 179 181 … … 218 220 !-- Public variables and constants 219 221 PUBLIC pc_heating_rate, canopy_mode, cthf, dt_plant_canopy, lad, lad_s, & 220 pch_index, p lant_canopy, prototype_lad, usm_lad_rma222 pch_index, prototype_lad, usm_lad_rma 221 223 222 224 … … 931 933 SUBROUTINE pcm_parin 932 934 935 USE control_parameters, & 936 ONLY: plant_canopy 933 937 934 938 IMPLICIT NONE -
palm/trunk/SOURCE/prognostic_equations.f90
r2719 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2719 2018-01-02 09:02:06Z maronga 27 30 ! Bugfix for last change. 28 31 ! … … 282 285 inflow_l, intermediate_timestep_count, & 283 286 intermediate_timestep_count_max, large_scale_forcing, & 284 large_scale_subsidence, microphysics_morrison, microphysics_seifert, & 285 microphysics_sat_adjust, neutral, nudging, ocean, outflow_l, & 286 outflow_s, passive_scalar, prho_reference, prho_reference, & 287 large_scale_subsidence, microphysics_morrison, & 288 microphysics_seifert, microphysics_sat_adjust, neutral, nudging,& 289 ocean, outflow_l, outflow_s, passive_scalar, plant_canopy, & 290 prho_reference, prho_reference, & 287 291 prho_reference, pt_reference, pt_reference, pt_reference, & 288 292 scalar_advec, scalar_advec, simulated_time, sloping_surface, & … … 362 366 363 367 USE plant_canopy_model_mod, & 364 ONLY: cthf, p lant_canopy, pcm_tendency368 ONLY: cthf, pcm_tendency 365 369 366 370 USE radiation_model_mod, & -
palm/trunk/SOURCE/radiation_model_mod.f90
r2724 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2724 2018-01-05 12:12:38Z maronga 27 30 ! Set default of average_radiation to .FALSE. 28 31 ! … … 235 238 initializing_actions, io_blocks, io_group, & 236 239 latitude, longitude, large_scale_forcing, lsf_surf, & 237 message_string, microphysics_morrison, p t_surface,&240 message_string, microphysics_morrison, plant_canopy, pt_surface,& 238 241 rho_surface, surface_pressure, time_since_reference_point 239 242 … … 268 271 269 272 USE plant_canopy_model_mod, & 270 ONLY: p lant_canopy, pc_heating_rate, lad_s, usm_lad_rma273 ONLY: pc_heating_rate, lad_s, usm_lad_rma 271 274 272 275 USE pegrid … … 4423 4426 4424 4427 USE plant_canopy_model_mod, & 4425 ONLY: plant_canopy, pch_index, & 4426 pc_heating_rate, lad_s, prototype_lad, usm_lad_rma 4428 ONLY: pch_index, pc_heating_rate, lad_s, prototype_lad, usm_lad_rma 4427 4429 4428 4430 IMPLICIT NONE -
palm/trunk/SOURCE/turbulence_closure_mod.f90
r2718 r2746 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Move flag plant canopy to modules 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 73 76 initializing_actions, intermediate_timestep_count, & 74 77 intermediate_timestep_count_max, kappa, km_constant, les_mw, & 75 ocean, p randtl_number, prho_reference, pt_reference, rans_mode,&76 rans_tke_e, rans_tke_l, simulated_time, timestep_scheme,&77 t urbulence_closure, turbulent_inflow, use_upstream_for_tke,&78 vpt_reference, ws_scheme_sca78 ocean, plant_canopy, prandtl_number, prho_reference, & 79 pt_reference, rans_mode, rans_tke_e, rans_tke_l, simulated_time,& 80 timestep_scheme, turbulence_closure, turbulent_inflow, & 81 use_upstream_for_tke, vpt_reference, ws_scheme_sca 79 82 80 83 USE advec_ws, & … … 103 106 104 107 USE plant_canopy_model_mod, & 105 ONLY: pcm_tendency , plant_canopy108 ONLY: pcm_tendency 106 109 107 110 USE statistics, & -
palm/trunk/SOURCE/urban_surface_mod.f90
r2737 r2746 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Move flag plant canopy to modules 29 ! 30 ! 2737 2018-01-11 14:58:11Z kanani 28 31 ! Removed unused variables t_surf_whole... 29 32 ! … … 231 234 g, pt_surface, large_scale_forcing, lsf_surf, spinup, & 232 235 spinup_pt_mean, spinup_time, time_do3d, dt_do3d, & 233 average_count_3d, varnamelength, urban_surface, kappa 236 average_count_3d, varnamelength, urban_surface, kappa, & 237 plant_canopy 234 238 235 239 USE cpulog, & … … 253 257 254 258 USE plant_canopy_model_mod, & 255 ONLY: pc_heating_rate, plant_canopy,usm_lad_rma259 ONLY: pc_heating_rate, usm_lad_rma 256 260 257 261 USE radiation_model_mod, &
Note: See TracChangeset
for help on using the changeset viewer.
