<p><b>laura@ucar.edu</b> 2012-11-19 10:38:23 -0700 (Mon, 19 Nov 2012)</p><p>an other update to the initialization of surface variables over sea-ice points<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_initialize_real.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_initialize_real.F        2012-11-19 16:47:52 UTC (rev 2313)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_initialize_real.F        2012-11-19 17:38:23 UTC (rev 2314)
@@ -460,7 +460,7 @@
  integer:: iCell,nCells
  integer,dimension(:),pointer:: landmask
 
- real(kind=RKIND),dimension(:),pointer  :: sst,tsk
+ real(kind=RKIND),dimension(:),pointer  :: sst,tsk,xice
  real(kind=RKIND),dimension(:,:),pointer:: tslb
 
 !--------------------------------------------------------------------------------------------------
@@ -474,14 +474,16 @@
  sst  =&gt; input % sst      % array
  tsk  =&gt; input % skintemp % array
  tslb =&gt; input % tslb     % array
+ xice =&gt; input % xice     % array
 
 !update the skin temperature and the soil temperature of the first soil layer with the updated
 !sea-surface temperatures:
  do iCell = 1, nCells
-    if(landmask(iCell) == 0) then
+    if(landmask(iCell) == 0 .and. xice(iCell) == 0) then
        tsk(iCell) = sst(iCell)
     endif
  enddo
+
  write(0,*) '--- end subroutine physics_update_sst:'
 
  end subroutine physics_init_sst

</font>
</pre>