<p><b>laura@ucar.edu</b> 2013-04-18 10:43:19 -0600 (Thu, 18 Apr 2013)</p><p>Revised the initialization of some variables prior to calling the surface layer scheme. Results are not changed.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_sfclayer.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_sfclayer.F        2013-04-18 16:41:33 UTC (rev 2768)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_sfclayer.F        2013-04-18 16:43:19 UTC (rev 2769)
@@ -137,31 +137,31 @@
     tsk_p(i,j)    = sfc_input    % skintemp % array(i)
     xland_p(i,j)  = sfc_input    % xland    % array(i)       
     !inout variables:
+    br_p(i,j)     = diag_physics % br     % array(i)
+    cpm_p(i,j)    = diag_physics % cpm    % array(i)
+    chs_p(i,j)    = diag_physics % chs    % array(i)
+    chs2_p(i,j)   = diag_physics % chs2   % array(i)
+    cqs2_p(i,j)   = diag_physics % cqs2   % array(i)
+    fh_p(i,j)     = diag_physics % fh     % array(i)
+    fm_p(i,j)     = diag_physics % fm     % array(i)
+    flhc_p(i,j)   = diag_physics % flhc   % array(i)
+    flqc_p(i,j)   = diag_physics % flqc   % array(i)
+    gz1oz0_p(i,j) = diag_physics % gz1oz0 % array(i)
     hfx_p(i,j)    = diag_physics % hfx    % array(i)
     qfx_p(i,j)    = diag_physics % qfx    % array(i)
+    qgh_p(i,j)    = diag_physics % qgh    % array(i)
     qsfc_p(i,j)   = diag_physics % qsfc   % array(i) 
+    lh_p(i,j)     = diag_physics % lh     % array(i)
     mol_p(i,j)    = diag_physics % mol    % array(i) 
+    psim_p(i,j)   = diag_physics % psim   % array(i)
+    psih_p(i,j)   = diag_physics % psih   % array(i)
     regime_p(i,j) = diag_physics % regime % array(i)
+    rmol_p(i,j)   = diag_physics % rmol   % array(i)
     ust_p(i,j)    = diag_physics % ust    % array(i)
+    wspd_p(i,j)   = diag_physics % wspd   % array(i)
     znt_p(i,j)    = diag_physics % znt    % array(i) 
     zol_p(i,j)    = diag_physics % zol    % array(i) 
     !output variables:
-    br_p(i,j)     = 0._RKIND
-    cpm_p(i,j)    = cp
-    chs_p(i,j)    = 0._RKIND
-    chs2_p(i,j)   = 0._RKIND
-    cqs2_p(i,j)   = 0._RKIND
-    flhc_p(i,j)   = 0._RKIND
-    flqc_p(i,j)   = 0._RKIND
-    fh_p(i,j)     = 0._RKIND
-    fm_p(i,j)     = 0._RKIND
-    gz1oz0_p(i,j) = 0._RKIND
-    lh_p(i,j)     = 0._RKIND
-    psim_p(i,j)   = 0._RKIND
-    psih_p(i,j)   = 0._RKIND
-    qgh_p(i,j)    = 0._RKIND
-    rmol_p(i,j)   = 0._RKIND
-    wspd_p(i,j)   = 0._RKIND
     q2_p(i,j)     = 0._RKIND
     t2m_p(i,j)    = 0._RKIND
     th2m_p(i,j)   = 0._RKIND
@@ -188,6 +188,8 @@
     diag_physics % chs    % array(i) = chs_p(i,j)
     diag_physics % chs2   % array(i) = chs2_p(i,j)
     diag_physics % cqs2   % array(i) = cqs2_p(i,j)
+    diag_physics % fh     % array(i) = fh_p(i,j)
+    diag_physics % fm     % array(i) = fm_p(i,j)
     diag_physics % flhc   % array(i) = flhc_p(i,j)
     diag_physics % flqc   % array(i) = flqc_p(i,j)
     diag_physics % gz1oz0 % array(i) = gz1oz0_p(i,j)

</font>
</pre>