updating version of clubb

Author: bstephens82

Skx_module.F90

@@ -75,7 +75,7 @@ subroutine Skx_func( nz, ngrdcol, xp2, xp3, &
     !$acc parallel loop gang vector collapse(2) default(present)
     do k = 1, nz
       do i = 1, ngrdcol
-        Skx(i,k) = xp3(i,k) / ( ( xp2(i,k) + Skx_denom_tol ) * sqrt( xp2(i,k) + Skx_denom_tol ) )
+        Skx(i,k) = xp3(i,k) * sqrt(( xp2(i,k) + Skx_denom_tol )**(-3))
       end do
     end do
     !$acc end parallel loop

adg1_adg2_3d_luhar_pdf.F90

@@ -147,7 +147,7 @@ subroutine ADG1_pdf_driver( nz, ngrdcol,                              & ! In
 
     integer :: j  ! Loop index
 
-    !$acc declare create( w_1_n, w_2_n )
+    !$acc enter data create( w_1_n, w_2_n )
 
     ! Calculate the mixture fraction and the PDF component means and variances
     ! of w.
@@ -204,6 +204,8 @@ subroutine ADG1_pdf_driver( nz, ngrdcol,                              & ! In
        enddo ! i=1, sclr_dim
     endif ! l_scalar_calc
 
+    !$acc exit data delete( w_1_n, w_2_n )
+
     return
 
   end subroutine ADG1_pdf_driver
@@ -1358,7 +1360,7 @@ subroutine backsolve_Luhar_params( Sk_max, Skx,          & ! In
   end subroutine backsolve_Luhar_params
 
   !=============================================================================
-  pure function max_cubic_root( a_coef, b_coef, c_coef, d_coef ) &
+  function max_cubic_root( a_coef, b_coef, c_coef, d_coef ) &
   result( max_root )
 
     ! Description: