<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 => input % sst % array
tsk => input % skintemp % array
tslb => input % tslb % array
+ xice => 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>