diff --git a/code/initcld.f b/code/initcld.f
index 5d0934a..0e48256 100644
--- a/code/initcld.f
+++ b/code/initcld.f
@@ -81,8 +81,8 @@ subroutine initcld
       d_nuc  = a0 / sqrt(2.)
 
 c  Contact parameter of water ice on dust ( m=cos(theta) )
-!     mteta  = 0.95        !Franck's number
-      mteta  = 0.975       !Franck recommends - newer value
+      mteta  = 0.95        !Franck's number
+!      mteta  = 0.975       !Franck recommends - newer value
 
 c  nu_eff = exp( dev^2) - 1
 c  Standard deviation of the dust distribution
@@ -189,6 +189,11 @@ subroutine initcld
         enddo
       enddo
 
+      print*,'qexti 1 = ',qexti_dst(1,:)
+      print*,'qscati 1 = ',qscati_dst(1,:)
+      print*,'gi 1 = ',gi_dst(1,:)
+
+
       close(60)
       close(61)
       close(62)
diff --git a/code/optci.f90 b/code/optci.f90
index d49eefc..70c4499 100644
--- a/code/optci.f90
+++ b/code/optci.f90
@@ -59,6 +59,7 @@ SUBROUTINE OPTCI(DTAUI,TAUCUMI,CO2I,PLEV,PFGASREF,TGASREF,       &
       real*8  :: QREFV(L_LEVELS+1)
       real*8  :: TAUREF(L_LEVELS+1)
       real*8  :: TAUREF_save(L_LEVELS+1)
+      real*8  :: TAUREFCLD_save(L_LEVELS+1)
 
 !     For clouds
       real*8  :: Qxicld(L_LEVELS+1,L_NSPECTI)
@@ -97,6 +98,7 @@ SUBROUTINE OPTCI(DTAUI,TAUCUMI,CO2I,PLEV,PFGASREF,TGASREF,       &
 
       do k=1,L_LEVELS+1
         tauref_save(k) = tauref(k)
+        taurefcld_save(k) = taurefcld(k)
       end do
 
 !  Determine the total gas opacity throughout the column, for each
@@ -270,6 +272,7 @@ SUBROUTINE OPTCI(DTAUI,TAUCUMI,CO2I,PLEV,PFGASREF,TGASREF,       &
 
       do k=1,L_LEVELS+1
         tauref(k) = tauref_save(k)
+        taurefcld(k) = taurefcld_save(k)
       end do
 
       RETURN
diff --git a/code/optcv.f90 b/code/optcv.f90
index d468315..a80451c 100644
--- a/code/optcv.f90
+++ b/code/optcv.f90
@@ -47,6 +47,7 @@ SUBROUTINE OPTCV(DTAUV,TAUV,TAUCUMV,CO2V,PLEV,PFGASREF,          &
       real*8  :: Qextref(L_LEVELS+1)
       real*8  :: TAUREF(L_LEVELS+1)
       real*8  :: TAUREF_save(L_LEVELS+1)
+      real*8  :: TAUREFCLD_save(L_LEVELS+1)
 
 !     For clouds
       real*8  :: Qxvcld(L_LEVELS+1,L_NSPECTV)
@@ -85,6 +86,7 @@ SUBROUTINE OPTCV(DTAUV,TAUV,TAUCUMV,CO2V,PLEV,PFGASREF,          &
 
       do k=1,L_LEVELS+1
         tauref_save(k) = tauref(k)
+        taurefcld_save(k) = taurefcld(k)
       end do
 
 !  Determine the total gas opacity throughout the column, for each
@@ -275,6 +277,7 @@ SUBROUTINE OPTCV(DTAUV,TAUV,TAUCUMV,CO2V,PLEV,PFGASREF,          &
 
       do k=1,L_LEVELS+1
         tauref(k) = tauref_save(k)
+        taurefcld(k) = taurefcld_save(k)
       end do
 
       RETURN