[Dart-dev] [10478] DART/trunk/models/cam/model_mod.f90: Fixed error in convert_vert, where new_array(3) was not set when vert_coord='pressure'.

[nancy at ucar.edu]

Revision: 10478
Author:   raeder
Date:     2016-07-08 13:57:33 -0600 (Fri, 08 Jul 2016)
Log Message:
-----------
Fixed error in convert_vert, where new_array(3) was not set when vert_coord='pressure'.

Modified Paths:
--------------
    DART/trunk/models/cam/model_mod.f90

-------------- next part --------------
Modified: DART/trunk/models/cam/model_mod.f90
===================================================================
--- DART/trunk/models/cam/model_mod.f90	2016-07-07 22:47:33 UTC (rev 10477)
+++ DART/trunk/models/cam/model_mod.f90	2016-07-08 19:57:33 UTC (rev 10478)
@@ -5212,7 +5212,7 @@
       endif
    enddo findzero
    write(string1, '(A)') 'invalid value for namelist "highest_obs_pressure_Pa"'
-   write(string2, '(A)') 'value is too large (and so located too low in atmosphere)'
+   write(string2, '(A)') 'value is too large (located out of the pure pressure levels of the atmosphere)'
    write(string3, '(A,F9.3,A)') 'must specify a value located above pressure ', p_col(lowest_ok), ' Pascals'
    call error_handler(E_ERR, 'set_highest_obs_limit', string1, source, revision, revdate, &
                       text2=string2, text3=string3)
@@ -5798,6 +5798,7 @@
    old_pressure = (1.0_r8 - frac) * p_col(bot_lev) + frac * p_col(top_lev)
 
    if (vert_coord == 'pressure') then
+      new_array(3) = old_pressure
       new_which = VERTISPRESSURE
    else if (vert_coord == 'log_invP') then
       new_array(3) = scale_height(p_surface=p_surf, p_above=old_pressure)