diff --git a/src/science/casa-cnp/biogeochem_casa.F90 b/src/science/casa-cnp/biogeochem_casa.F90
index ab0f4ebb1..ef4f860b6 100644
--- a/src/science/casa-cnp/biogeochem_casa.F90
+++ b/src/science/casa-cnp/biogeochem_casa.F90
@@ -52,7 +52,7 @@ SUBROUTINE biogeochem(ktau,dels,idoY,LALLOC,veg,soil,casabiome,casapool,casaflux
 ! zero annual sums
 IF (idoy==1) CALL casa_cnpflux(casaflux,casapool,casabal,.TRUE.)
 
-IF (cable_user%PHENOLOGY_SWITCH.EQ.'MODIS' .OR. cable_runtime%esm15 ) THEN
+IF (cable_user%PHENOLOGY_SWITCH.EQ.'MODIS' ) THEN
        CALL phenology(idoy,veg,phen)
 ENDIF
 CALL avgsoil(veg,soil,casamet)
diff --git a/src/science/casa-cnp/casa_feedback.F90 b/src/science/casa-cnp/casa_feedback.F90
index c916b8401..0636ba843 100644
--- a/src/science/casa-cnp/casa_feedback.F90
+++ b/src/science/casa-cnp/casa_feedback.F90
@@ -53,8 +53,6 @@ SUBROUTINE casa_feedback( ktau, veg, casabiome, casapool, casamet )
 
   integer np,ivt
   real, dimension(mp)  :: ncleafx,npleafx, pleafx, nleafx ! local variables
-  real, dimension(17)                   ::  xnslope
-  data xnslope/0.80,1.00,2.00,1.00,1.00,1.00,0.50,1.00,0.34,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00/
   REAL :: npleafx_coef
 
   ! first initialize
@@ -89,7 +87,7 @@ SUBROUTINE casa_feedback( ktau, veg, casabiome, casapool, casamet )
 
           veg%vcmax(np) = ( casabiome%nintercept(ivt)                          &
                   + casabiome%nslope(ivt) * npleafx_coef * ncleafx(np)         &
-                  / casabiome%sla(ivt) ) * xnslope(ivt) * 1.0e-6
+                  / casabiome%sla(ivt) ) * 1.0e-6
 
        END IF