<p><b>laura@ucar.edu</b> 2012-08-06 11:04:50 -0600 (Mon, 06 Aug 2012)</p><p>moved the initialization of the arrays sst and xice so that they are not overwritten when both data are present.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-08-06 15:42:25 UTC (rev 2087)
+++ branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-08-06 17:04:50 UTC (rev 2088)
@@ -113,7 +113,7 @@
else if (config_test_case == 8 ) then
- write(0,*) ' real-data surface (SST) update test case '
+ write(0,*) 'real-data surface (SST) update test case '
block_ptr => domain % blocklist
do while (associated(block_ptr))
call init_atm_test_case_sfc(domain, domain % dminfo, block_ptr % mesh, block_ptr % fg, block_ptr % state % time_levs(1) % state, &
@@ -4412,10 +4412,9 @@
!initialization of sea-surface temperature (SST) and sea-ice fraction (XICE) arrays,
!prior to reading the input data:
- fg % sst % array (1:grid%nCells) = 0.0
- fg % xice % array (1:grid%nCells) = 0.0
if (index(field % field, 'SKINTEMP') /= 0 .or. index(field % field, 'SST') /= 0) then
+ fg % sst % array (1:grid%nCells) = 0.0_RKIND
! Interpolation routines use real(kind=RKIND), so copy from default real array
allocate(slab_r8(field % nx, field % ny))
@@ -4480,6 +4479,7 @@
deallocate(field % slab)
else if (index(field % field, 'SEAICE') /= 0) then
+ fg % xice % array (1:grid%nCells) = 0.0_RKIND
! Interpolation routines use real(kind=RKIND), so copy from default real array
allocate(slab_r8(field % nx, field % ny))
</font>
</pre>