<p><b>laura@ucar.edu</b> 2012-11-16 16:13:15 -0700 (Fri, 16 Nov 2012)</p><p>move the initialization of sea-ice points at the bottom of the subroutine physics_initialize_real so that initialization of some fields at polar latitudes are not overwritten<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-16 23:10:09 UTC (rev 2311)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_initialize_real.F        2012-11-16 23:13:15 UTC (rev 2312)
@@ -80,9 +80,6 @@
call physics_init_sst(mesh,fg)
endif
-!initialize seaice points:
- call physics_init_seaice(mesh,fg)
-
!initialization of the surface background albedo: interpolation of the monthly values to the
!initial date:
initial_date = trim(config_start_time)
@@ -126,6 +123,9 @@
!initialization of soil layers properties:
call init_soil_layers(mesh,fg,dminfo)
+!initialize seaice points:
+call physics_init_seaice(mesh,fg)
+
!define xland over land and ocean:
do iCell = 1, nCellsSolve
if(landmask(iCell) .eq. 1 .or. (landmask(iCell).eq.0 .and. seaice(iCell).eq.1._RKIND)) then
</font>
</pre>