From b08f5482cfd6b257ea73eb452bcad5eaa7e99e59 Mon Sep 17 00:00:00 2001
From: Erik Kluzek <erik@ucar.edu>
Date: Thu, 18 Oct 2018 16:51:49 -0600
Subject: [PATCH] Add changes so that single precision history output is not an
 option

---
 cime_config/namelist_definition_mosart.xml | 6 +++---
 src/riverroute/RtmHistFile.F90             | 2 +-
 src/riverroute/RtmIO.F90                   | 2 +-
 3 files changed, 5 insertions(+), 5 deletions(-)

diff --git a/cime_config/namelist_definition_mosart.xml b/cime_config/namelist_definition_mosart.xml
index d10f7d1..00b3deb 100644
--- a/cime_config/namelist_definition_mosart.xml
+++ b/cime_config/namelist_definition_mosart.xml
@@ -261,13 +261,13 @@
     <type>integer(6)</type> 
     <category>history</category>
     <group>mosart_inparm</group> 
-    <valid_values>1,2</valid_values>
+    <valid_values>1</valid_values>
     <values>
-      <value>2</value>
+      <value>1</value>
     </values>
     <desc>
       Per tape series  history file density (i.e. output precision) 
-      1=double precision, 2=single precision
+      1=double precision, 2=single precision (NOT working)
     </desc>
   </entry>
 
diff --git a/src/riverroute/RtmHistFile.F90 b/src/riverroute/RtmHistFile.F90
index 9f0ad69..7a3ba55 100644
--- a/src/riverroute/RtmHistFile.F90
+++ b/src/riverroute/RtmHistFile.F90
@@ -42,7 +42,7 @@ module RtmHistFile
 !
   integer :: ni
   integer, public :: &
-       rtmhist_ndens(max_tapes) = 2         ! namelist: output density of netcdf history files
+       rtmhist_ndens(max_tapes) = 1         ! namelist: output density of netcdf history files
   integer, public :: &
        rtmhist_mfilt(max_tapes) = 30        ! namelist: number of time samples per tape
   integer, public :: &
diff --git a/src/riverroute/RtmIO.F90 b/src/riverroute/RtmIO.F90
index 6bfae8c..18be51b 100644
--- a/src/riverroute/RtmIO.F90
+++ b/src/riverroute/RtmIO.F90
@@ -1827,7 +1827,7 @@ subroutine ncd_io_real_var1(varname, data, dim1name, &
           call pio_setframe(ncid,vardesc, int(nt,kind=PIO_Offset_kind))
        end if
        if(xtype == ncd_float) then
-          call pio_write_darray(ncid, vardesc, iodesc_plus%iodesc, real(data, kind=r4), status, fillval=real(spval,kind=r4))
+          call shr_sys_abort( subname//' error: Attempt to write out single-precision data which is current NOT implemented (see issue #18)' )
        else
           call pio_write_darray(ncid, vardesc, iodesc_plus%iodesc, data, status, fillval=spval)
        endif