<p><b>laura@ucar.edu</b> 2013-02-11 16:42:52 -0700 (Mon, 11 Feb 2013)</p><p>updated the skin temperature over ocean points, including over sea-ice points when an SST update file is used to update the skin temperatures over oceans. Previously, the skin temperatures were not updated to SSTs when sea-ice was present.<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        2013-02-11 22:35:53 UTC (rev 2459)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_initialize_real.F        2013-02-11 23:42:52 UTC (rev 2460)
@@ -124,7 +124,7 @@
call init_soil_layers(mesh,fg,dminfo)
!initialize seaice points:
-call physics_init_seaice(mesh,fg)
+ call physics_init_seaice(mesh,fg)
!define xland over land and ocean:
do iCell = 1, nCellsSolve
@@ -478,8 +478,15 @@
!update the skin temperature and the soil temperature of the first soil layer with the updated
!sea-surface temperatures:
+!change made so that the SSTs read for the surface update file are the same as the skin temperature
+!over the oceans.
+!do iCell = 1, nCells
+! if(landmask(iCell) == 0 .and. xice(iCell) == 0) then
+! tsk(iCell) = sst(iCell)
+! endif
+!enddo
do iCell = 1, nCells
- if(landmask(iCell) == 0 .and. xice(iCell) == 0) then
+ if(landmask(iCell) == 0) then
tsk(iCell) = sst(iCell)
endif
enddo
</font>
</pre>