Ignore:
Timestamp:
May 9, 2018 8:42:38 AM (6 years ago)
Author:
maronga
Message:

series of bugfixes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r2977 r3014  
    2828! -----------------
    2929! $Id$
     30! Added pc_transpiration_rate
     31!
     32! 2977 2018-04-17 10:27:57Z kanani
    3033! Implement changes from branch radiation (r2948-2971) with minor modifications.
    3134! (moh.hefny):
     
    264267!>       by zero, e.g. in case fraq(0,m) + fraq(1,m) = 0?!
    265268!> @todo Use unit 90 for OPEN/CLOSE of input files (FK)
     269!> @todo Move plant canopy stuff into plant canopy code
    266270!------------------------------------------------------------------------------!
    267271 MODULE urban_surface_mod
     
    309313   
    310314    USE plant_canopy_model_mod,                                                &
    311         ONLY:  pc_heating_rate
     315        ONLY:  pc_heating_rate, pc_transpiration_rate
    312316   
    313317    USE radiation_model_mod,                                                   &
     
    23612365        INTEGER(iwp)                                           :: ids,idsint,idsidx,isurf,isvf,isurfs,isurflt
    23622366        INTEGER(iwp)                                           :: is,js,ks,i,j,k,iwl,istat, l, m
    2363         INTEGER(iwp)                                           ::  k_topo    !< topography top index
     2367        INTEGER(iwp)                                           :: k_topo    !< topography top index
    23642368
    23652369        dirstart = (/ startland, startwall, startwall, startwall, startwall /)
     
    45394543!--             in case of cthf /= 0 => we need to allocate it for our use here
    45404544                ALLOCATE( pc_heating_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     4545
     4546                pc_heating_rate = 0.0_wp
     4547
     4548            ENDIF
     4549
     4550            IF ( .NOT.  ALLOCATED( pc_transpiration_rate) )  THEN
     4551!--             then pc_heating_rate is allocated in init_plant_canopy
     4552!--             in case of cthf /= 0 => we need to allocate it for our use here
     4553                ALLOCATE( pc_transpiration_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     4554
     4555                pc_transpiration_rate = 0.0_wp
     4556
     4557
    45414558            ENDIF
    45424559        ENDIF
Note: See TracChangeset for help on using the changeset viewer.