implemented bug fix 2023

decyjphr-org · Apr 12, 2023 · f313e61 · f313e61
1 parent 0636208
commit f313e61
Show file tree

Hide file tree

Showing 3 changed files with 13 additions and 2 deletions.
diff --git 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
@@ -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
@@ -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