Changeset 2317 for palm/trunk/SOURCE/urban_surface_mod.f90
- Timestamp:
- Jul 20, 2017 5:27:19 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r2296 r2317 21 21 ! Current revisions: 22 22 ! ------------------ 23 ! 23 ! Get topography top index via Function call 24 24 ! 25 25 ! Former revisions: … … 603 603 ! 604 604 !-- Find topography top index 605 k_topo = MAXLOC( MERGE( 1, 0, & 606 BTEST( wall_flags_0(:,j,i), 12 ) & 607 ), DIM = 1 & 608 ) - 1 605 k_topo = get_topography_top_index( j, i, 's' ) 606 609 607 DO k = nzt+1, 0, -1 610 608 IF ( lad_s(k,j,i) /= 0.0_wp ) THEN … … 706 704 DO j = ijdb(3,ids), ijdb(4,ids) 707 705 708 k_topo = MAXLOC( MERGE( 1, 0, & 709 BTEST( wall_flags_0(:,j,i), 12 ) & 710 ), DIM = 1 & 711 ) - 1 712 k_topo2 = MAXLOC( MERGE( 1, 0, & 713 BTEST( wall_flags_0(:,j-jdir(ids),i-idir(ids)), 12 ) & 714 ), DIM = 1 & 715 ) - 1 706 k_topo = get_topography_top_index( j, i, 's' ) 707 k_topo2 = get_topography_top_index( j-jdir(ids), i-idir(ids), 's' ) 716 708 717 709 k = nzut - MAX( k_topo, k_topo2 ) … … 732 724 ! 733 725 !-- Find topography top index 734 k_topo = MAXLOC( MERGE( 1, 0, & 735 BTEST( wall_flags_0(:,j,i), 12 ) & 736 ), DIM = 1 & 737 ) - 1 726 k_topo = get_topography_top_index( j, i, 's' ) 727 738 728 DO k = k_topo + 1, pct(j,i) 739 729 ipcgb = ipcgb + 1 … … 813 803 DO i = ijdb(1,ids), ijdb(2,ids) 814 804 DO j = ijdb(3,ids), ijdb(4,ids) 815 k_topo = MAXLOC( MERGE( 1, 0, & 816 BTEST( wall_flags_0(:,j,i), 12 ) & 817 ), DIM = 1 & 818 ) - 1 819 k_topo2 = MAXLOC( MERGE( 1, 0, & 820 BTEST( wall_flags_0(:,j-jdir(ids),i-idir(ids)), 12 ) & 821 ), DIM = 1 & 822 ) - 1 805 k_topo = get_topography_top_index( j, i, 's' ) 806 k_topo2 = get_topography_top_index( j-jdir(ids), i-idir(ids), 's' ) 807 823 808 DO k = MAX(k_topo,k_topo2)+1, nzut 824 809 isurf = isurf + 1 … … 1929 1914 DO i = nxl, nxr 1930 1915 DO j = nys, nyn 1931 k = MAXLOC( & 1932 MERGE( 1, 0, & 1933 BTEST( wall_flags_0(:,j,i), 12 ) & 1934 ), DIM = 1 & 1935 ) - 1 1916 k = get_topography_top_index( j, i, 's' ) 1936 1917 1937 1918 usm_lad(k:nzut, j, i) = lad_s(0:nzut-k, j, i) … … 3806 3787 ! 3807 3788 !-- Following expression equals former kk = k - nzb_s_inner(j,i) 3808 kk = k - ( MAXLOC( & 3809 MERGE( 1, 0, & 3810 BTEST( wall_flags_0(:,j,i), 12 ) & 3811 ), DIM = 1 & 3812 ) - 1 & 3813 ) !- lad arrays are defined flat 3789 kk = k - ( get_topography_top_index( j, i, 's' ) ) !- lad arrays are defined flat 3814 3790 pc_heating_rate(kk, j, i) = (pcbinsw(ipcgb) + pcbinlw(ipcgb)) & 3815 3791 * pchf_prep(k) * pt(k, j, i) !-- = dT/dt
Note: See TracChangeset
for help on using the changeset viewer.