changed precision of coordinates from real to double only for regridded SMC type grids.

Author: ukmo-kitstokes

model/src/w3metamd.F90

@@ -258,8 +258,9 @@ SUBROUTINE META_LIST_APPEND_D(LIST, ATTNAME, DVAL)
 
     TYPE(META_LIST_T), INTENT(INOUT) :: LIST
     CHARACTER(*), INTENT(IN)         :: ATTNAME
-    INTEGER, PARAMETER               :: dp = SELECTED_REAL_KIND(15)  
-    REAL(dp), INTENT(IN)             :: DVAL
+    !INTEGER, PARAMETER               :: dp = SELECTED_REAL_KIND(15)  
+    !REAL(dp), INTENT(IN)             :: DVAL
+    DOUBLE PRECISION, INTENT(IN)     :: DVAL
     !/ ------------------------------------------------------------------- /
     !/ Local parameters
     !/

model/src/w3ounfmetamd.F90

@@ -1271,8 +1271,9 @@ SUBROUTINE GET_ATTVAL_TYPE(BUF, ILINE, ATTV, ATT_TYPE)
     !
     REAL :: R
     INTEGER :: I, IERR
-    INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(15)
-    REAL(dp) :: D    
+    !INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(15)
+    !REAL(dp) :: D    
+    DOUBLE PRECISION :: D
 
     ! Get attribute and type (default to "c" if no type set)
     ATT_TYPE = 'c'
@@ -2337,8 +2338,9 @@ SUBROUTINE WRITE_FREEFORM_META_LIST(NCID, VARID, METALIST, ERR)
     !/
     INTEGER :: I, IVAL
     REAL    :: RVAL
-    INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(15)
-    REAL(dp) :: DVAL
+    !INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(15)
+    !REAL(dp) :: DVAL
+    DOUBLE PRECISION :: DVAL
     TYPE(META_PAIR_T), POINTER :: P
 
     IF(METALIST%N .EQ. 0) RETURN

model/src/ww3_ounf.F90

@@ -2303,15 +2303,19 @@ SUBROUTINE W3EXNC ( NX, NY, IX1, IXN, IY1, IYN, NSEA,             &
                   !LON(1) = X0D
                   !LAT(1) = Y0D
                   DO i=1,NXO
-                    lon(i)=DBLE(X0D+SXD*DBLE(i-1))
+                    lon(i)=REAL(X0D+SXD*DBLE(i-1))
+                    ! KS: Regridded SMC data - use double precision
+                    dbllon(i)=DBLE(X0D+SXD*DBLE(i-1))
 #endif
 #ifdef W3_RTD
                     LON2DEQ(i,:) = lon(i)
 #endif
 #ifdef W3_SMC
                   END DO
                   DO i=1,NYO
-                    lat(i)=DBLE(Y0D+SYD*DBLE(i-1))
+                    lat(i)=REAL(Y0D+SYD*DBLE(i-1))
+                    ! KS: Regridded SMC data - use double precision
+                    dbllat(i)=DBLE(Y0D+SYD*DBLE(i-1))
 #endif
 #ifdef W3_RTD
                     LAT2DEQ(:,i) = lat(i)
@@ -2350,10 +2354,10 @@ SUBROUTINE W3EXNC ( NX, NY, IX1, IXN, IY1, IYN, NSEA,             &
                 !LON(1) = X0D
                 !LAT(1) = Y0D
                 DO I=1,NX
-                  LON(I)=DBLE(X0D+SXD*DBLE(I-1))
+                  LON(I)=REAL(X0D+SXD*DBLE(I-1))
                 END DO
                 DO I=1,NY
-                  LAT(I)=DBLE(Y0D+SYD*DBLE(I-1))
+                  LAT(I)=REAL(Y0D+SYD*DBLE(I-1))
                 END DO
 #ifdef W3_RTD
                 IF ( RTDL ) THEN
@@ -2475,16 +2479,22 @@ SUBROUTINE W3EXNC ( NX, NY, IX1, IXN, IY1, IYN, NSEA,             &
             IF (GTYPE.EQ.RLGTYPE .OR. GTYPE.EQ.SMCTYPE) THEN
               IF(SMCGRD) THEN ! CB: shelter original code from SMC grid
 #ifdef W3_SMC
-                IRET=NF90_PUT_VAR(NCID,VARID(1),LON(:))
-                CALL CHECK_ERR(IRET)
-                IRET=NF90_PUT_VAR(NCID,VARID(2),LAT(:))
-                CALL CHECK_ERR(IRET)
                 IF(SMCOTYPE .EQ. 1) THEN
+                  IRET=NF90_PUT_VAR(NCID,VARID(1),LON(:))
+                  CALL CHECK_ERR(IRET)
+                  IRET=NF90_PUT_VAR(NCID,VARID(2),LAT(:))
+                  CALL CHECK_ERR(IRET)
+
                   ! For type 1 SCM file also put lat/lons and cell sizes:
                   IRET=NF90_PUT_VAR(NCID,VARID(5),SMCCX)
                   CALL CHECK_ERR(IRET)
                   IRET=NF90_PUT_VAR(NCID,VARID(6),SMCCY)
                   CALL CHECK_ERR(IRET)
+                ELSE ! KS: Regridded SMC data - use double precision
+                  IRET=NF90_PUT_VAR(NCID,VARID(1),dbllon(:))
+                  CALL CHECK_ERR(IRET)
+                  IRET=NF90_PUT_VAR(NCID,VARID(2),dbllat(:))
+                  CALL CHECK_ERR(IRET)
                 ENDIF
 #endif
               ELSE ! SMCGRD