<p><b>mhoffman@lanl.gov</b> 2013-03-27 21:30:50 -0600 (Wed, 27 Mar 2013)</p><p>BRANCH COMMIT - land ice<br>
<br>
Adding a calculation of layerThickness during initialization. Otherwise the call made to vertical_remap that happens in land_ice_diagnostic_solve during init messes up the values of tracers.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_mpas_core.F
===================================================================
--- branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_mpas_core.F 2013-03-27 20:36:54 UTC (rev 2675)
+++ branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_mpas_core.F 2013-03-28 03:30:50 UTC (rev 2676)
@@ -196,7 +196,7 @@
err = 0
state => block % state % time_levs(1) % state ! initial state
-
+
! Call init routines ====
call land_ice_setup_vertical_coords(mesh, err_tmp)
err = ior(err, err_tmp)
@@ -227,9 +227,16 @@
! Initialize state ====
-
call mpas_timer_start("initial state calculation")
+ ! We need to assign the layerThickness before calling land_ice_diagnostic_solve because
+ ! it will perform vertical remapping which needs layer thicknesses to map from.
+ ! Vertical remapping is not necessary (and perhaps not desired) at this point
+ ! in the code, but it is currently part of land_ice_diagnostic_solve.
+ do iCell = 1, mesh % nCells
+ state % layerThickness % array(:,iCell) = mesh % layerThicknessFractions % array * state % thickness % array(iCell)
+ enddo
+
call land_ice_diagnostic_solve(mesh, state, err)
! vertexMask and edgeMask needs halo updates before they can be used. Halo updates should occur outside of block loops.
</font>
</pre>