Fixes for PTM

Bug fixes based on global ALM simulation with PTM.
E3SM-Project · May 23, 2017 · bc031d0 · bc031d0
1 parent 2ca6667
commit bc031d0
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 4 deletions.
diff --git a/components/clm/src/external_models/emi/src/ExternalModelPTMMod.F90 b/components/clm/src/external_models/emi/src/ExternalModelPTMMod.F90
@@ -1567,11 +1567,12 @@ subroutine EM_PTM_TBased_Solve(this, dt, nstep, l2e_list, e2l_list, bounds_clump
                 dz_1d(idx)             = col_z(c,j)*2.d0
                 num_snow_layer_1d(idx) = -col_snl(c)
 
+                ! Save data for internal SoE auxvars
+                snow_water_1d(idx)                         = l2e_h2osno(c)
+
                 if (col_snl(c) /= 0) then
                    ! Save data for boundary SoE auxvars
                    sabg_soil_1d((c-begc)*nlevgrnd + j) = l2e_frac_snow_eff(c)*l2e_sabg_lyr(c,j)
-                   ! Save data for internal SoE auxvars
-                   snow_water_1d(idx)                         = l2e_h2osno(c)
                 else
                    ! Save data for internal SoE auxvars
                    tsurf_tuning_factor_1d(idx) = col_dz(c,j) / &
@@ -1721,7 +1722,7 @@ subroutine EM_PTM_TBased_Solve(this, dt, nstep, l2e_list, e2l_list, bounds_clump
           l = col_landunit(c)
           if ((col_active(c)==1) .and. (lun_lakpoi(l) == 0) .and. (lun_urbpoi(l) == 0)) then
              if (j >= col_snl(c)+1) then
-                e2l_tsnow(c,j-1) = temperature_1d(idx)
+                e2l_tsnow(c,j) = temperature_1d(idx)
              endif
           endif
        enddo

diff --git a/components/clm/src/external_models/mpp b/components/clm/src/external_models/mpp