Ready to test CO2 changes

Author: Steve Goldhaber

src/physics/cam/ghg_data.F90

@@ -136,7 +136,7 @@ subroutine trcmix(name, lchnk, ncol, clat, pmid, q)
    real(r8) dlat            ! latitude in degrees
    real(r8) ptrop           ! pressure level of tropopause
    real(r8) pratio          ! pressure divided by ptrop
-   real(r8) trop_mmr        ! tropospheric mass mixing ratio
+   real(r8) trop_mmr(pcols) ! tropospheric mass mixing ratio
    real(r8) scale           ! pressure scale height
 !-----------------------------------------------------------------------
 
@@ -150,12 +150,12 @@ subroutine trcmix(name, lchnk, ncol, clat, pmid, q)
 
    else if (name == 'CO2') then
 
-      q = chem_surfvals_co2_rad(lchnk, ncol)
+      q = chem_surfvals_co2_rad()
 
    else if (name == 'CH4') then
 
       ! set tropospheric mass mixing ratios
-      trop_mmr = rmwch4 * chem_surfvals_get('CH4VMR', lchnk, ncol)
+      trop_mmr(:ncol) = rmwch4 * chem_surfvals_get('CH4VMR', lchnk, ncol)
 
       do k = 1,pver
          do i = 1,ncol
@@ -172,18 +172,18 @@ subroutine trcmix(name, lchnk, ncol, clat, pmid, q)
 
             ! determine output mass mixing ratios
             if (pmid(i,k) >= ptrop) then
-               q(i,k) = trop_mmr
+               q(i,k) = trop_mmr(i)
             else
                pratio = pmid(i,k)/ptrop
-               q(i,k) = trop_mmr * (pratio)**scale
+               q(i,k) = trop_mmr(i) * (pratio)**scale
             end if
          end do
       end do
 
    else if (name == 'N2O') then
 
       ! set tropospheric mass mixing ratios
-      trop_mmr = rmwn2o * chem_surfvals_get('N2OVMR', lchnk, ncol)
+      trop_mmr(:ncol) = rmwn2o * chem_surfvals_get('N2OVMR', lchnk, ncol)
 
       do k = 1,pver
          do i = 1,ncol
@@ -200,18 +200,18 @@ subroutine trcmix(name, lchnk, ncol, clat, pmid, q)
 
             ! determine output mass mixing ratios
             if (pmid(i,k) >= ptrop) then
-               q(i,k) = trop_mmr
+               q(i,k) = trop_mmr(i)
             else
                pratio = pmid(i,k)/ptrop
-               q(i,k) = trop_mmr * (pratio)**scale
+               q(i,k) = trop_mmr(i) * (pratio)**scale
             end if
          end do
       end do
 
    else if (name == 'CFC11') then
 
       ! set tropospheric mass mixing ratios
-      trop_mmr = rmwf11 * chem_surfvals_get('F11VMR', lchnk, ncol)
+      trop_mmr(:ncol) = rmwf11 * chem_surfvals_get('F11VMR', lchnk, ncol)
 
       do k = 1,pver
          do i = 1,ncol
@@ -228,18 +228,18 @@ subroutine trcmix(name, lchnk, ncol, clat, pmid, q)
 
             ! determine output mass mixing ratios
             if (pmid(i,k) >= ptrop) then
-               q(i,k) = trop_mmr
+               q(i,k) = trop_mmr(i)
             else
                pratio = pmid(i,k)/ptrop
-               q(i,k) = trop_mmr * (pratio)**scale
+               q(i,k) = trop_mmr(i) * (pratio)**scale
             end if
          end do
       end do
 
    else if (name == 'CFC12') then
 
       ! set tropospheric mass mixing ratios
-      trop_mmr = rmwf12 * chem_surfvals_get('F12VMR', lchnk, ncol)
+      trop_mmr(:ncol) = rmwf12 * chem_surfvals_get('F12VMR', lchnk, ncol)
 
       do k = 1,pver
          do i = 1,ncol
@@ -256,10 +256,10 @@ subroutine trcmix(name, lchnk, ncol, clat, pmid, q)
 
             ! determine output mass mixing ratios
             if (pmid(i,k) >= ptrop) then
-               q(i,k) = trop_mmr
+               q(i,k) = trop_mmr(i)
             else
                pratio = pmid(i,k)/ptrop
-               q(i,k) = trop_mmr * (pratio)**scale
+               q(i,k) = trop_mmr(i) * (pratio)**scale
             end if
          end do
       end do