<p><b>laura@ucar.edu</b> 2012-08-06 11:04:50 -0600 (Mon, 06 Aug 2012)</p><p>moved the initialization of the arrays sst and xice so that they are not overwritten when both data are present.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-08-06 15:42:25 UTC (rev 2087)
+++ branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-08-06 17:04:50 UTC (rev 2088)
@@ -113,7 +113,7 @@
 
       else if (config_test_case == 8 ) then
 
-         write(0,*) ' real-data surface (SST) update test case '
+         write(0,*) 'real-data surface (SST) update test case '
          block_ptr =&gt; domain % blocklist
          do while (associated(block_ptr))
             call init_atm_test_case_sfc(domain, domain % dminfo, block_ptr % mesh, block_ptr % fg, block_ptr % state % time_levs(1) % state, &amp;
@@ -4412,10 +4412,9 @@
 
             !initialization of sea-surface temperature (SST) and sea-ice fraction (XICE) arrays,
             !prior to reading the input data:
-            fg % sst  % array (1:grid%nCells) = 0.0
-            fg % xice % array (1:grid%nCells) = 0.0
 
             if (index(field % field, 'SKINTEMP') /= 0 .or. index(field % field, 'SST') /= 0) then
+               fg % sst  % array (1:grid%nCells) = 0.0_RKIND
 
                ! Interpolation routines use real(kind=RKIND), so copy from default real array
                allocate(slab_r8(field % nx, field % ny))
@@ -4480,6 +4479,7 @@
                deallocate(field % slab)
 
             else if (index(field % field, 'SEAICE') /= 0) then
+               fg % xice % array (1:grid%nCells) = 0.0_RKIND
 
                ! Interpolation routines use real(kind=RKIND), so copy from default real array
                allocate(slab_r8(field % nx, field % ny))

</font>
</pre>