@@ -25,7 +25,7 @@ module SatellitePhenologyMod
2525 !
2626 ! !PUBLIC MEMBER FUNCTIONS:
2727 public :: GetSatellitePhenologyInputs ! put the data into the correct format
28- public :: SetSPModeCanopyStructs ! CLMSP Ecosystem dynamics: phenology, vegetation
28+ public :: SetSPModeCanopyStructs ! CLM(BGC)-SP phenology and vegetation
2929 public :: SatellitePhenologyInit ! Dynamically allocate memory
3030 public :: interpMonthlyVeg ! interpolate monthly vegetation data
3131 public :: readAnnualVegetation ! Read in annual vegetation (needed for Dry-deposition)
@@ -114,10 +114,10 @@ subroutine GetSatellitePhenologyInputs(bounds, num_filter, filter, canopystate_i
114114 !- ----------------------------------------------------------------------
115115
116116 associate( &
117- tlai_driver = > canopystate_inst% tlai_input_patch , & ! Output: [real(r8) (:) ] one-sided leaf area index, no burying by snow
118- tsai_driver = > canopystate_inst% tsai_input_patch , & ! Output: [real(r8) (:) ] one-sided stem area index, no burying by snow
119- htop_driver = > canopystate_inst% htop_input_patch , & ! Output: [real(r8) (:) ] canopy top (m)
120- hbot_driver = > canopystate_inst% hbot_input_patch & ! Output: [real(r8) (:) ] canopy bottom (m)
117+ tlai_tinterp = > canopystate_inst% tlai_input_patch , & ! Output: [real(r8) (:) ] one-sided leaf area index, no burying by snow
118+ tsai_tinterp = > canopystate_inst% tsai_input_patch , & ! Output: [real(r8) (:) ] one-sided stem area index, no burying by snow
119+ htop_tinterp = > canopystate_inst% htop_input_patch , & ! Output: [real(r8) (:) ] canopy top (m)
120+ hbot_tinterp = > canopystate_inst% hbot_input_patch & ! Output: [real(r8) (:) ] canopy bottom (m)
121121 )
122122
123123 if (use_lai_streams) then
@@ -145,12 +145,12 @@ subroutine GetSatellitePhenologyInputs(bounds, num_filter, filter, canopystate_i
145145 ! bottom height HBOT <- mhvb1 and mhvb2
146146
147147 if (.not. use_lai_streams) then
148- tlai_driver (p) = timwt(1 )* mlai2t(p,1 ) + timwt(2 )* mlai2t(p,2 )
148+ tlai_tinterp (p) = timwt(1 )* mlai2t(p,1 ) + timwt(2 )* mlai2t(p,2 )
149149 endif
150150
151- tsai_driver (p) = timwt(1 )* msai2t(p,1 ) + timwt(2 )* msai2t(p,2 )
152- htop_driver (p) = timwt(1 )* mhvt2t(p,1 ) + timwt(2 )* mhvt2t(p,2 )
153- hbot_driver (p) = timwt(1 )* mhvb2t(p,1 ) + timwt(2 )* mhvb2t(p,2 )
151+ tsai_tinterp (p) = timwt(1 )* msai2t(p,1 ) + timwt(2 )* msai2t(p,2 )
152+ htop_tinterp (p) = timwt(1 )* mhvt2t(p,1 ) + timwt(2 )* mhvt2t(p,2 )
153+ hbot_tinterp (p) = timwt(1 )* mhvb2t(p,1 ) + timwt(2 )* mhvb2t(p,2 )
154154
155155 end do
156156 end associate
@@ -159,7 +159,7 @@ end subroutine GetSatellitePhenologyInputs
159159
160160 ! ==============================================================================
161161
162- subroutine SetSPModeCanopyStructs (bounds , num_filter , filter , &
162+ subroutine SetSatellitePhenologyBGCCanopyStructs (bounds , num_filter , filter , &
163163 waterdiagnosticbulk_inst , canopystate_inst )
164164 !
165165 ! !DESCRIPTION:
@@ -187,7 +187,8 @@ subroutine SetSPModeCanopyStructs(bounds, num_filter, filter, &
187187 real (r8 ) :: fb ! fraction of canopy layer covered by snow
188188
189189 if (use_fates_sp) then
190- return ! we should not be here
190+ write (iulog,* ) ' SetSatellitePhenologyBGCCanopyStructs' , ' should not be calling this method when use_fates_sp == .true.'
191+ call endrun(msg= errMsg(sourcefile, __LINE__))
191192 end if
192193
193194 associate( &
@@ -250,7 +251,7 @@ subroutine SetSPModeCanopyStructs(bounds, num_filter, filter, &
250251
251252 end associate
252253
253- end subroutine SetSPModeCanopyStructs
254+ end subroutine SetSatellitePhenologyBGCCanopyStructs
254255
255256 ! ==============================================================================
256257 subroutine interpMonthlyVeg (bounds , canopystate_inst )
0 commit comments