Changeset 2698 for palm/trunk/SOURCE/plant_canopy_model_mod.f90
 Timestamp:
 Dec 14, 2017 6:46:24 PM (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/plant_canopy_model_mod.f90
r2696 r2698 20 20 ! Current revisions: 21 21 !  22 ! 22 ! Bugfix in get_topography_top_index 23 23 ! 24 24 ! Former revisions: … … 151 151 152 152 USE surface_mod, & 153 ONLY: get_topography_top_index 153 ONLY: get_topography_top_index_ji 154 154 155 155 … … 400 400 DO j = nys, nyn 401 401 IF ( pch_index_ji(j,i) /= 0 ) THEN 402 k_topo = get_topography_top_index ( j, i, 's' )402 k_topo = get_topography_top_index_ji( j, i, 's' ) 403 403 DO k = k_topo, k_topo + pch_index_ji(j,i) 404 404 local_pf(i,j,k) = pc_heating_rate(kk_topo,j,i) … … 415 415 DO j = nys, nyn 416 416 IF ( pch_index_ji(j,i) /= 0 ) THEN 417 k_topo = get_topography_top_index ( j, i, 's' )417 k_topo = get_topography_top_index_ji( j, i, 's' ) 418 418 DO k = k_topo, k_topo + pch_index_ji(j,i) 419 419 local_pf(i,j,k) = lad_s(kk_topo,j,i) … … 803 803 ! Check whether topography and local vegetation on top exceed 804 804 ! height of the model domain. 805 k = get_topography_top_index ( j, i, 's' )805 k = get_topography_top_index_ji( j, i, 's' ) 806 806 IF ( k + pch_index_ji(j,i) >= nzt + 1 ) THEN 807 807 message_string = 'Local vegetation height on top of ' // & … … 1110 1110 ! 1111 1111 ! Determine topographytop index on ugrid 1112 k_wall = get_topography_top_index ( j, i, 'u' )1112 k_wall = get_topography_top_index_ji( j, i, 'u' ) 1113 1113 DO k = k_wall+1, k_wall + pch_index_ji(j,i) 1114 1114 … … 1174 1174 ! 1175 1175 ! Determine topographytop index on vgrid 1176 k_wall = get_topography_top_index ( j, i, 'v' )1176 k_wall = get_topography_top_index_ji( j, i, 'v' ) 1177 1177 1178 1178 DO k = k_wall+1, k_wall + pch_index_ji(j,i) … … 1239 1239 ! 1240 1240 ! Determine topographytop index on wgrid 1241 k_wall = get_topography_top_index ( j, i, 'w' )1241 k_wall = get_topography_top_index_ji( j, i, 'w' ) 1242 1242 1243 1243 DO k = k_wall+1, k_wall + pch_index_ji(j,i)  1 … … 1291 1291 ! 1292 1292 ! Determine topographytop index on scalargrid 1293 k_wall = get_topography_top_index ( j, i, 's' )1293 k_wall = get_topography_top_index_ji( j, i, 's' ) 1294 1294 1295 1295 DO k = k_wall+1, k_wall + pch_index_ji(j,i) … … 1308 1308 ! 1309 1309 ! Determine topographytop index on scalargrid 1310 k_wall = get_topography_top_index ( j, i, 's' )1310 k_wall = get_topography_top_index_ji( j, i, 's' ) 1311 1311 1312 1312 DO k = k_wall+1, k_wall + pch_index_ji(j,i) … … 1338 1338 ! 1339 1339 ! Determine topographytop index on scalargrid 1340 k_wall = get_topography_top_index ( j, i, 's' )1340 k_wall = get_topography_top_index_ji( j, i, 's' ) 1341 1341 1342 1342 DO k = k_wall+1, k_wall + pch_index_ji(j,i) … … 1367 1367 ! 1368 1368 ! Determine topographytop index on scalargrid 1369 k_wall = get_topography_top_index ( j, i, 's' )1369 k_wall = get_topography_top_index_ji( j, i, 's' ) 1370 1370 1371 1371 DO k = k_wall+1, k_wall + pch_index_ji(j,i) … … 1461 1461 ! 1462 1462 ! Determine topographytop index on ugrid 1463 k_wall = get_topography_top_index ( j, i, 'u' )1463 k_wall = get_topography_top_index_ji( j, i, 'u' ) 1464 1464 DO k = k_wall + 1, k_wall + pch_index_ji(j,i) 1465 1465 … … 1521 1521 ! 1522 1522 ! Determine topographytop index on vgrid 1523 k_wall = get_topography_top_index ( j, i, 'v' )1523 k_wall = get_topography_top_index_ji( j, i, 'v' ) 1524 1524 1525 1525 DO k = k_wall + 1, k_wall + pch_index_ji(j,i) … … 1581 1581 ! 1582 1582 ! Determine topographytop index on wgrid 1583 k_wall = get_topography_top_index ( j, i, 'w' )1583 k_wall = get_topography_top_index_ji( j, i, 'w' ) 1584 1584 1585 1585 DO k = k_wall + 1, k_wall + pch_index_ji(j,i)  1 … … 1628 1628 ! 1629 1629 ! Determine topographytop index on scalar grid 1630 k_wall = get_topography_top_index ( j, i, 's' )1630 k_wall = get_topography_top_index_ji( j, i, 's' ) 1631 1631 1632 1632 DO k = k_wall + 1, k_wall + pch_index_ji(j,i) … … 1641 1641 ! 1642 1642 ! Determine topographytop index on scalar grid 1643 k_wall = get_topography_top_index ( j, i, 's' )1643 k_wall = get_topography_top_index_ji( j, i, 's' ) 1644 1644 1645 1645 DO k = k_wall + 1, k_wall + pch_index_ji(j,i) … … 1667 1667 ! 1668 1668 ! Determine topographytop index on scalar grid 1669 k_wall = get_topography_top_index ( j, i, 's' )1669 k_wall = get_topography_top_index_ji( j, i, 's' ) 1670 1670 1671 1671 DO k = k_wall + 1, k_wall + pch_index_ji(j,i) … … 1693 1693 ! 1694 1694 ! Determine topographytop index on scalar grid 1695 k_wall = get_topography_top_index ( j, i, 's' )1695 k_wall = get_topography_top_index_ji( j, i, 's' ) 1696 1696 1697 1697 DO k = k_wall + 1, k_wall + pch_index_ji(j,i)
Note: See TracChangeset
for help on using the changeset viewer.