Merge branch 'hotfix_v3.1'

This merge corrects a bug in the extrapolation of RH
below the surface for MPAS-Atmosphere.

Author: mgduda

README.md

@@ -1,4 +1,4 @@
-MPAS-v3.0
+MPAS-v3.1
 ====
 
 The Model for Prediction Across Scales (MPAS) is a collaborative project for

src/core_atmosphere/Registry.xml

@@ -1,5 +1,5 @@
 <?xml version="1.0"?>
-<registry model="mpas" core="atmosphere" version="3.0">
+<registry model="mpas" core="atmosphere" version="3.1">
 
 <!-- **************************************************************************************** -->
 <!-- ************************************** Dimensions ************************************** -->

src/core_init_atmosphere/Registry.xml

@@ -1,5 +1,5 @@
 <?xml version="1.0"?>
-<registry model="mpas" core="init_atmosphere" version="3.0">
+<registry model="mpas" core="init_atmosphere" version="3.1">
 
 <!-- **************************************************************************************** -->
 <!-- ************************************** Dimensions ************************************** -->

src/core_init_atmosphere/mpas_init_atm_cases.F

@@ -4004,7 +4004,7 @@ subroutine init_atm_case_gfs(block, mesh, nCells, nEdges, nVertLevels, fg, state
             sorted_arr(2,k) = rh_fg(k,iCell)
          end do
          call mpas_quicksort(nfglevels_actual, sorted_arr)
-         do k=1,nVertLevels
+         do k=nVertLevels,1,-1
             target_z = 0.5 * (zgrid(k,iCell) + zgrid(k+1,iCell))
 !           scalars(index_qv,k,iCell) = vertical_interp(target_z, nfglevels_actual, sorted_arr, order=1, extrap=0)
             scalars(index_qv,k,iCell) = vertical_interp(target_z, nfglevels_actual-1, &

src/core_landice/Registry.xml

@@ -1,5 +1,5 @@
 <?xml version="1.0"?>
-<registry model="mpas" core="landice" version="3.0">
+<registry model="mpas" core="landice" version="3.1">
 
 
 <!-- ======================================================================= -->

src/core_ocean/Registry.xml

@@ -1,5 +1,5 @@
 <?xml version="1.0"?>
-<registry model="mpas" core="ocean" version="3.0">
+<registry model="mpas" core="ocean" version="3.1">
 
 	<dims>
 		<dim name="nCells" units="unitless"

src/core_sw/Registry.xml

@@ -1,5 +1,5 @@
 <?xml version="1.0"?>
-<registry model="mpas" core="sw" version="3.0">
+<registry model="mpas" core="sw" version="3.1">
 	<dims>
 		<dim name="nCells"/>
 		<dim name="nEdges"/>