Declare variables passed to w3nco library as double precission

physics/aerinterp.F90

@@ -98,7 +98,7 @@ END SUBROUTINE read_aerdata
 !
 !**********************************************************************
       SUBROUTINE read_aerdataf ( me, master, iflip, idate, FHOUR, errmsg, errflg)
-      use machine, only: kind_phys, kind_io4, kind_io8
+      use machine, only: kind_phys, kind_dbl_prec
       use aerclm_def
 
 !--- in/out
@@ -111,7 +111,7 @@ SUBROUTINE read_aerdataf ( me, master, iflip, idate, FHOUR, errmsg, errflg)
       logical      :: file_exist
       integer  IDAT(8),JDAT(8)
       real(kind=kind_phys) rjday
-      real(8) RINC(5)
+      real(kind=kind_dbl_prec) rinc(5)
       integer jdow, jdoy, jday
 
       integer, allocatable  :: invardims(:)
@@ -218,7 +218,7 @@ END SUBROUTINE setindxaer
       SUBROUTINE aerinterpol( me,master,nthrds,npts,IDATE,FHOUR,iflip, jindx1,jindx2, &
                              ddy,iindx1,iindx2,ddx,lev,prsl,aerout)
 !
-      use machine, only: kind_phys, kind_io4, kind_io8
+      use machine, only: kind_phys, kind_dbl_prec
       use aerclm_def
 
       implicit none
@@ -238,6 +238,7 @@ SUBROUTINE aerinterpol( me,master,nthrds,npts,IDATE,FHOUR,iflip, jindx1,jindx2,
       real(kind=kind_phys) aerpm(npts,levsaer,ntrcaer)
       real(kind=kind_phys) prsl(npts,lev), aerpres(npts,levsaer)
       real(kind=kind_phys) rjday
+      real(kind=kind_dbl_prec) rinc(5)
       integer jdow, jdoy, jday
 !
       IDAT = 0

physics/h2ointerp.f90

@@ -131,7 +131,7 @@ subroutine h2ointerpol(me,npts,idate,fhour,jindx1,jindx2,h2oplout,ddy)
 !
 ! May 2015 Shrinivas Moorthi - Prepare for H2O interpolation
 !
-      use machine , only : kind_phys
+      use machine , only : kind_phys, kind_dbl_prec
       use h2o_def
       implicit none
       integer             j,j1,j2,l,npts,nc,n1,n2
@@ -145,6 +145,7 @@ subroutine h2ointerpol(me,npts,idate,fhour,jindx1,jindx2,h2oplout,ddy)
       real(kind=kind_phys) ddy(npts)
       real(kind=kind_phys) h2oplout(npts,levh2o,h2o_coeff)
       real(kind=kind_phys) rjday
+      real(kind=kind_dbl_prec) rinc(5)
       integer              jdow, jdoy, jday
 !
       idat    = 0

physics/iccninterp.F90

@@ -129,7 +129,7 @@ END SUBROUTINE setindxci
       SUBROUTINE ciinterpol(me,npts,IDATE,FHOUR,jindx1,jindx2,ddy, &
                  iindx1,iindx2,ddx,lev, prsl, ciplout,ccnout)
 !
-      USE MACHINE,  ONLY : kind_phys
+      USE MACHINE,  ONLY : kind_phys, kind_dbl_prec
       use iccn_def
       implicit none
       integer   i1,i2, iday,j,j1,j2,l,npts,nc,n1,n2,lev,k,i
@@ -145,6 +145,7 @@ SUBROUTINE ciinterpol(me,npts,IDATE,FHOUR,jindx1,jindx2,ddy, &
       real(kind=kind_phys) ccnout(npts,lev),ccnpm(npts,kcipl)
       real(kind=kind_phys) cipres(npts,kcipl), prsl(npts,lev)
       real(kind=kind_phys) rjday
+      real(kind=kind_dbl_prec) rinc(5)
       integer jdow, jdoy, jday
 !
       IDAT=0

physics/ozinterp.f90

@@ -135,7 +135,7 @@ END SUBROUTINE setindxoz
 !
       SUBROUTINE ozinterpol(me,npts,IDATE,FHOUR,jindx1,jindx2,ozplout,ddy)
 !
-      USE MACHINE,  ONLY : kind_phys
+      USE MACHINE,  ONLY : kind_phys, kind_dbl_prec
       USE OZNE_DEF
       implicit none
       integer             iday,j,j1,j2,l,npts,nc,n1,n2
@@ -148,6 +148,7 @@ SUBROUTINE ozinterpol(me,npts,IDATE,FHOUR,jindx1,jindx2,ozplout,ddy)
       real(kind=kind_phys) DDY(npts)
       real(kind=kind_phys) ozplout(npts,levozp,oz_coeff)
       real(kind=kind_phys) rjday
+      real(kind=kind_dbl_prec) rinc(5)
       integer jdow, jdoy, jday
 !
       IDAT=0

physics/sfcsub.F

@@ -2747,7 +2747,7 @@ subroutine fixrdg(lugb,idim,jdim,fngrib,                          &
       real (kind=kind_io8) gdata(idim*jdim)
       logical gaus
       real (kind=kind_io8) blno,blto
-      real (kind=kind_io8), allocatable :: data8(:)
+      real (kind=kind_dbl_prec), allocatable :: data8(:)
 !
       logical*1, allocatable :: lbms(:)
 !
@@ -6969,7 +6969,7 @@ subroutine clima(lugb,iy,im,id,ih,fh,len,lsoil,slmskl,slmskw,     &
      &,                gaus, blno, blto, me,lprnt,iprnt, fnalbc2, ialb  &
      &,                tile_num_ch, i_index, j_index)
 !
-      use machine , only : kind_io8,kind_io4
+      use machine , only : kind_io8,kind_io4, kind_dbl_prec
       implicit none
       character(len=*), intent(in) :: tile_num_ch
       integer, intent(in) :: i_index(len), j_index(len)
@@ -7029,7 +7029,7 @@ subroutine clima(lugb,iy,im,id,ih,fh,len,lsoil,slmskl,slmskw,     &
       data dayhf/ 15.5, 45.0, 74.5,105.0,135.5,166.0,
      &           196.5,227.5,258.0,288.5,319.0,349.5,380.5/
 !
-      real (kind=kind_io8) fha(5)
+      real (kind=kind_dbl_prec) fha(5)
       integer ida(8),jda(8),ivtyp, kpd7
 !
       real (kind=kind_io8), allocatable :: tsf(:,:),sno(:,:),
@@ -8296,7 +8296,7 @@ subroutine fixrdc(lugb,fngrib,kpds5,kpds7,mon,slmask,             &
 !
       real (kind=kind_io8) gdata(len), slmask(len)
       real (kind=kind_io8), allocatable :: data(:,:), rslmsk(:,:)
-      real (kind=kind_io8), allocatable :: data8(:)
+      real (kind=kind_dbl_prec), allocatable :: data8(:)
       real (kind=kind_io8), allocatable :: rlngrb(:), rltgrb(:)
 !
       logical lmask, yr2kc, gaus, ijordr
@@ -8471,7 +8471,7 @@ subroutine fixrda(lugb,fngrib,kpds5,slmask,                       &
 !
       real (kind=kind_io8) gdata(len), slmask(len)
       real (kind=kind_io8), allocatable :: data(:,:),rslmsk(:,:)
-      real (kind=kind_io8), allocatable :: data8(:)
+      real (kind=kind_dbl_prec), allocatable :: data8(:)
       real (kind=kind_io8), allocatable :: rlngrb(:), rltgrb(:)
 !
       logical lmask, yr2kc, gaus, ijordr
@@ -8492,7 +8492,7 @@ subroutine fixrda(lugb,fngrib,kpds5,slmask,                       &
       integer mjday(12)
       data mjday/31,28,31,30,31,30,31,31,30,31,30,31/
 !
-      real (kind=kind_io8) fha(5)
+      real (kind=kind_dbl_prec) fha(5)
       integer ida(8),jda(8)
 !
       allocate(data8(1:mdata))