<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>