--- ncarg2d/src/libncarg_gks/bwi/argb2ci.f.orig	2019-02-28 08:44:39.000000000 +0900
+++ ncarg2d/src/libncarg_gks/bwi/argb2ci.f	2020-07-25 14:27:05.000000000 +0900
@@ -16,10 +16,10 @@
       integer index, nearest
       integer i
       integer ARGBMASK, RMASK, GMASK, BMASK
-      parameter (ARGBMASK = Z'40000000')
-      parameter (RMASK     = Z'00FF0000')
-      parameter (GMASK     = Z'0000FF00')
-      parameter (BMASK     = Z'000000FF')
+      data ARGBMASK/Z'40000000'/
+      data RMASK   /Z'00FF0000'/
+      data GMASK   /Z'0000FF00'/
+      data BMASK   /Z'000000FF'/
       real r, g, b, dist, mindist
 
       if (iand(index, ARGBMASK).eq.0) then
@@ -31,8 +31,8 @@
       mindist = 2e31
       nearest = 0
       do i=1,mol
-          r = (iand(index, RMASK) / Z'0000FFFF') / 255.
-          g = (iand(index, GMASK) / Z'000000FF') / 255.
+          r = (iand(index, RMASK) / int(Z'0000FFFF')) / 255.
+          g = (iand(index, GMASK) / int(Z'000000FF')) / 255.
           b = (iand(index, BMASK))               / 255.
 
 C         we don't need absolute distance, so forego the sqrt operation...