LIBTOOLS/conf/config.gfortranR64

+#
+#   Configuration file for PGF (64-bit reals).
+#
+AR      = ar
+ARFLAGS = rv
+#
+CC      = gcc
+CFLAGS  = -O2 -D__hpux -DREAL_8 -DREAL_BIGGER_THAN_INTEGER
+FASTCFLAGS = 
+#
+FC      = gfortran
+VECTFFLAGS =
+CPPFLAGS = -D__hpux -DREAL_8 -DREAL_BIGGER_THAN_INTEGER -Dextend2o -Dg95
+FFLAGS = -w -g -O2 -fdefault-real-8
+#
+LDFLAGS = -L . -l emos$(R64)
+RANLIB  = /bin/true
+CT      = /bin/true
+NPROC   = 1

LIBTOOLS/conf/listing

+>lst.conf
+for i in config.* ; do
+echo $i >>lst.conf
+echo '---------------' >> lst.conf
+cat $i >>lst.conf
+echo '======================================='>> lst.conf
+done
+

LIBTOOLS/lib/COMPRESS/Makefile

+LIBCOMP = liblficomp.a
+#######################################
+DIR_OBJ = ./$(ARCH)
+
+VPATH = src:$(DIR_OBJ)
+INC = -I$(DIR_OBJ)
+
+DIR_CONF:=$(shell pwd|sed -e 's/lib\/.*/conf/')
+
+include $(DIR_CONF)/config.$(ARCH)
+include Rules.$(ARCH)
+
+
+%.o:%.f90
+	$(CPP) $(INC) $(CPPFLAGS) $< >  $(DIR_OBJ)/cpp_$(*F).f90
+	$(F90) $(INC) -c $(F90FLAGS) $(DIR_OBJ)/cpp_$(*F).f90 -o $(DIR_OBJ)/$(*F).o
+	-@mv  *.mod $(DIR_OBJ)/. 2> /dev/null || echo pas de module dans $*.f90
+
+%.o:%.c
+	$(CC) $(INC) $(CFLAGS) $(CPPFLAGS) -c $< -o $(DIR_OBJ)/$(*F).o
+
+$(LIBCOMP) : $(DIR_OBJ)/.dummy $(OBJS)
+	cd $(DIR_OBJ);$(AR) crv $@ $(OBJS)
+
+$(DIR_OBJ)/.dummy :
+	mkdir -p $(DIR_OBJ)
+	@touch $(DIR_OBJ)/.dummy
+
+compress.o : searchgrp.o comppar.o
+decompress.o : searchgrp.o comppar.o
+
+clean:
+	(if [ -d $(DIR_OBJ) ] ; then cd $(DIR_OBJ); rm -f cpp_*.f90 *.o ; fi)
+
+distclean:
+	rm -rf $(DIR_OBJ)
+
+

LIBTOOLS/lib/COMPRESS/Rules.AIX32

+CPPFLAGS += -DBIG_endian
+
+OBJS=ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.AIX64

+CPPFLAGS += -DBIG_endian
+
+OBJS=ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.HPNAGf95

+F90FLAGS += -O2
+CPPFLAGS += -DNAGf95 -DBIG_endian -DNO_UNDERSCORE
+
+OBJS=comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.HPf90

+CFLAGS   += -Ae
+F90FLAGS += -O3
+CPPFLAGS += -DNO_UNDERSCORE -DBIG_endian
+
+OBJS=ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.LXNAGf95

+F77FLAGS +=
+F90FLAGS += -O2
+CPPFLAGS = -DNAGf95 -DLITTLE_endian
+
+OBJS=comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.LXg95

+F77FLAGS +=
+F90FLAGS += -O2
+CPPFLAGS = -DLITTLE_endian
+
+OBJS=comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o ieee_is_nan.o

LIBTOOLS/lib/COMPRESS/Rules.LXgfortran

+F77FLAGS +=
+F90FLAGS += -O2
+CPPFLAGS = -DLITTLE_endian
+
+OBJS=comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o ieee_is_nan.o

LIBTOOLS/lib/COMPRESS/Rules.LXifort

+F77FLAGS +=
+F90FLAGS += -O2
+CPPFLAGS = -DLITTLE_endian
+
+OBJS=comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o ieee_is_nan.o

LIBTOOLS/lib/COMPRESS/Rules.LXpgf90

+CPPFLAGS = -DLITTLE_endian
+
+OBJS=ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.SGI32

+F77LAGS  +=
+F90FLAGS += -O2
+CFLAGS   +=
+CPPFLAGS += -DBIG_endian
+
+OBJS = ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.SGI64

+F77FLAGS +=
+F90FLAGS += -O2
+CPPFLAGS += -DBIG_endian
+
+OBJS = ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.SX5

+F90FLAGS +=
+CPPFLAGS += -DBIG_endian -DSX5
+
+OBJS=ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.SX8

+F90FLAGS +=
+CPPFLAGS += -DBIG_endian -DSX5
+
+OBJS=ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/Rules.VPP

+F77FLAGS +=
+F90FLAGS += 
+CPPFLAGS += -DVPP -DBIG_endian
+
+OBJS=ieee_is_nan.o comppar.o compress.o decompress.o nearestpow2.o searchgrp.o bitbuff.o

LIBTOOLS/lib/COMPRESS/src/bitbuff.c

+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+
+#ifdef VPP
+# include <sys/types.h>
+typedef __uint64_t WORD;
+#else
+#ifdef SX5
+typedef unsigned long uint64_t;
+#else
+# include <inttypes.h>
+#endif
+typedef uint64_t WORD;
+#endif
+
+#define WORDSIZE 64
+
+#ifdef NO_UNDERSCORE
+# define SET_FILLIDX set_fillidx
+# define GET_FILLIDX get_fillidx
+# define FILL_BBUFF  fill_bbuff
+# define SET_EXTRACTIDX set_extractidx
+# define GET_EXTRACTIDX get_extractidx
+# define EXTRACT_BBUFF  extract_bbuff
+#else
+# define SET_FILLIDX set_fillidx_
+# define GET_FILLIDX get_fillidx_
+# define FILL_BBUFF  fill_bbuff_
+# define SET_EXTRACTIDX set_extractidx_
+# define GET_EXTRACTIDX get_extractidx_
+# define EXTRACT_BBUFF  extract_bbuff_
+#endif
+
+int outidx = 0;
+int outbrem = WORDSIZE ;
+
+int inidx = 0;
+int inbrem = WORDSIZE;
+
+void SET_FILLIDX(unsigned *idx, unsigned *bitoffset){
+    inidx  = *idx;
+    inidx += (*bitoffset/WORDSIZE);
+    inbrem = WORDSIZE - (*bitoffset%WORDSIZE);
+}
+
+void GET_FILLIDX(unsigned *idx, unsigned *bitoffset){
+    *idx    = inidx;
+    *bitoffset = WORDSIZE - inbrem;
+}
+
+void FILL_BBUFF(WORD *out, int *n, unsigned *val){
+    /* inidx = index of the current buffer elt to fill */
+    /* inbrem = number of bits remaining on buffer elt out[idx] */
+
+    /* fill buffer out with n low bits of val */
+
+    if (inbrem >= *n){
+	inbrem = inbrem - *n;
+	/* turn to 0 the n bits of out */
+	out[inidx] &= ~(~(~(WORD)0 << *n) << inbrem);
+	/* now set the n bits of out to val */
+	out[inidx] |= (*val & ~(~(WORD)0 << *n)) << inbrem;
+	return;
+    } else {
+	int nex = *n - inbrem; /* number of bits that will be filled later */
+	if (inbrem != 0){
+	    /* turn to 0 the inbrem lower bits of out */
+	    out[inidx] &= (~(WORD)0 << inbrem) ;
+	    /* now set the inbrem lower bits of out with val */
+	    out[inidx] |= ((*val >> nex) & ~(~(WORD)0 << inbrem));
+	}
+	inidx++;
+	inbrem = WORDSIZE;
+	FILL_BBUFF(out, &nex, val);
+    }
+
+}
+
+void SET_EXTRACTIDX(unsigned *idx, unsigned *bitoffset) {
+    outidx = *idx;
+    outidx += (*bitoffset/WORDSIZE);
+    outbrem = WORDSIZE-(*bitoffset%WORDSIZE);
+}
+
+void GET_EXTRACTIDX(unsigned *idx, unsigned *bitoffset){
+    *idx = outidx;
+    *bitoffset = WORDSIZE - outbrem;
+}
+
+    
+void extract_bbuff_rec(WORD *buff, int *n, unsigned *val) {
+    
+    if (outbrem >= *n){
+	outbrem = outbrem - *n;
+	*val = (*val << *n) | (unsigned)((buff[outidx]>>outbrem) & ~(~(WORD)0 << *n));
+	return;
+    } else {
+	int nex = *n - outbrem;
+	if (outbrem != 0){
+	    *val = (*val << outbrem)| (unsigned)(buff[outidx] & ~(~(WORD)0 << outbrem));
+
+	}
+	outidx++;
+	outbrem=WORDSIZE;
+	extract_bbuff_rec(buff,&nex,val);
+    }
+}
+
+void EXTRACT_BBUFF(WORD *buff, int *n, unsigned *val) {
+    
+    unsigned tmpval;
+
+    tmpval=0;
+    extract_bbuff_rec(buff,n,&tmpval);
+    *val = tmpval;
+}
+

LIBTOOLS/lib/COMPRESS/src/comppar.f90

+!-----------------------------------------------------------------
+!--------------- special set of characters for RCS information
+!-----------------------------------------------------------------
+! $Source$ $Revision$ $Date$
+!-----------------------------------------------------------------
+!-----------------------------------------------------------------
+!-----------------------------------------------------------------
+MODULE MODD_COMPPAR
+IMPLICIT NONE 
+! Debug mode : set LPDEBUG to .TRUE.
+LOGICAL,PARAMETER :: LPDEBUG = .FALSE.
+
+
+! contains coding parameters for (de)compress routines
+
+INTEGER,PARAMETER :: JPCSTENCOD = 1 ! constant array 
+INTEGER,PARAMETER :: JPSOPENCOD = 2 ! second order packing 
+INTEGER,PARAMETER :: JPEXTENCOD = 3 ! second order packing with min/max values excluded
+
+! Extended code when JPEXTENCOD enabled
+!
+! BE CAREFUL : 3 bits are reserved for coding this code => max value is 7
+INTEGER,PARAMETER :: JPCONST   = 0 ! constant value array
+INTEGER,PARAMETER :: JPNORM    = 1 ! same as JPSOPENCOD
+INTEGER,PARAMETER :: JPMINEXCL = 2 ! Min value is isolated
+INTEGER,PARAMETER :: JPMAXEXCL = 3 ! Max value is isolated
+INTEGER,PARAMETER :: JPMINMAXEXCL = 4 ! Min&Max values are isolated
+INTEGER,PARAMETER :: JP2VAL       = 5 ! 2 different values in array
+INTEGER,PARAMETER :: JP3VAL       = 6 ! 3 different values in array
+INTEGER,PARAMETER :: JPOTHER      = 7 ! for future use
+INTEGER,PARAMETER :: JPLOG        = 8
+END MODULE MODD_COMPPAR

LIBTOOLS/lib/COMPRESS/src/compress.f90

+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
+!MNH_LIC for details. version 1.
+!-----------------------------------------------------------------
+SUBROUTINE COMPRESS_FIELD(XTAB,KX,KY,KNBTOT,KNBUSE)
+USE MODD_COMPPAR
+USE MODE_SEARCHGRP
+
+USE,INTRINSIC :: IEEE_ARITHMETIC
+
+IMPLICIT NONE 
+
+REAL,PARAMETER :: PPFLOATMIN = 2.0**(-126)
+
+INTEGER, INTENT(IN) :: KX,KY
+!INTEGER, INTENT(IN) :: KNBLEV
+INTEGER, INTENT(IN) :: KNBTOT
+REAL(KIND=8),DIMENSION(KNBTOT),INTENT(INOUT) :: XTAB
+
+INTEGER, INTENT(OUT) :: KNBUSE
+
+INTEGER :: INBLEV
+INTEGER,DIMENSION(:), ALLOCATABLE :: ITAB
+REAL :: XMIN,XMAX
+TYPE(SOP_t) :: SOPRES
+INTEGER :: IND1, IND2
+INTEGER :: GELT,IBE
+INTEGER :: ILEVNBELT
+INTEGER :: NBITCOD
+INTEGER :: II, JI, JJ
+INTEGER :: BITOFFSET
+INTEGER :: GRPIDX,GRPOFF,IDXSAVE,OFFSAVE
+INTEGER :: nbgroupmod
+INTEGER :: IEXTCOD
+CHARACTER(LEN=8),PARAMETER :: KEYWORD='COMPRESS'
+REAL,DIMENSION(KNBTOT) :: XWORKTAB
+LOGICAL :: LUPREAL,LNAN
+logical :: gnansupport
+
+ILEVNBELT = KX*KY
+LUPREAL = .FALSE.
+LNAN    = .FALSE.
+
+if ( IEEE_SUPPORT_NAN( xtab(1)) ) then
+  gnansupport=.true.
+else
+  gnansupport=.false.
+end if
+
+! Check for NAN and change Upper and Lower bound according to 32bits real limits.
+DO JI=1,KNBTOT
+  IF ( gnansupport .and. IEEE_IS_NAN(XTAB(JI)) ) THEN
+    XTAB(JI)=0.
+    LNAN = .TRUE.
+  ELSE IF (ABS(XTAB(JI)) > HUGE(1.0_4)) THEN
+    XTAB(JI) = SIGN(REAL(HUGE(1.0_4)/1.1,8),XTAB(JI))
+    LUPREAL = .TRUE.
+  ELSEIF (ABS(XTAB(JI)) < TINY(1.0_4)) THEN
+    XTAB(JI) = 0.
+  END IF
+END DO
+
+XMIN=MINVAL(XTAB)
+XMAX=MAXVAL(XTAB)
+PRINT *,'MINVAL,MAXVAL= ',XMIN,XMAX
+IF (LNAN) PRINT *,"==================> NAN values DETECTED : set to 0.0"
+IF (LUPREAL) PRINT *,"==================> OVERFLOW values DETECTED : set to ",HUGE(1.0_4)/1.1
+
+! Convert 64 bits real to 32 bits real
+XWORKTAB(:) = XTAB(:)
+!
+! BEWARE : Now XTAB is overwritten. 
+!          XWORKTAB contains the 32 bits floating point data.
+!
+CALL SET_FILLIDX(0,0)
+! store 8 characters header string in buffer
+DO II=1,LEN(KEYWORD)
+  CALL FILL_BBUFF(XTAB,8,ICHAR(KEYWORD(II:II)))
+END DO
+
+! is whole array XTAB64 a constant field ?
+
+IF (xmin == xmax) THEN
+  PRINT *,"--------> CONSTANT ARRAY !"
+  CALL FILL_BBUFF(XTAB,32,JPCSTENCOD)
+  CALL FILL_BBUFF(XTAB,32,KNBTOT)
+  CALL FILL_BBUFF(XTAB,32,xmin)
+  CALL GET_FILLIDX(KNBUSE,BITOFFSET)
+  KNBUSE=KNBUSE+1
+  RETURN
+END IF
+
+
+INBLEV = KNBTOT/(ILEVNBELT)
+IF (KNBTOT /= (INBLEV*ILEVNBELT)) THEN
+  PRINT *,'Pb in COMPRESS_FIELD : KNBTOT must be a multiple of KX*KY'
+  STOP
+END IF
+
+
+
+ALLOCATE(ITAB(ILEVNBELT))
+CALL INI_SOPDATA(SOPRES)
+
+CALL FILL_BBUFF(XTAB,32,JPEXTENCOD)
+CALL FILL_BBUFF(XTAB,32,KNBTOT)
+CALL FILL_BBUFF(XTAB,32,KX)
+CALL FILL_BBUFF(XTAB,32,KY)
+
+DO JI=1,INBLEV
+  IND1=(JI-1)*ILEVNBELT+1
+  IND2=JI*ILEVNBELT
+  IF (LPDEBUG) PRINT *,"---- Compressing Level ",JI," ----"
+  CALL COMP_FOPEXT(XWORKTAB(IND1:IND2),ITAB,IEXTCOD)
+  IF (IEXTCOD /= JPCONST) THEN
+    CALL INVERTCOL(ITAB,KX,KY)
+    CALL RECSEARCH(ITAB,SOPRES)
+    GELT = MAXVAL(SOPRES%IEND(1:SOPRES%NBGRP)-SOPRES%IBEG(1:SOPRES%NBGRP)+1)
+    IBE = FMINBITS_IN_WORD(GELT)
+    CALL GET_FILLIDX(GRPIDX,GRPOFF) ! save the idx/offset for future NBGRP modification
+    CALL FILL_BBUFF(XTAB,32,SOPRES%NBGRP)
+    CALL FILL_BBUFF(XTAB,5,IBE)
+  
+    NBGROUPMOD = SOPRES%NBGRP
+    DO II=1,SOPRES%NBGRP
+      GELT = SOPRES%IEND(II)-SOPRES%IBEG(II)+1
+      nbitcod  = FMINBITS_IN_WORD(SOPRES%VALMAX(II)-SOPRES%VALMIN(II))
+      !    PRINT *, 'Groupe',II,'(',GELT,')',':',SOPRES%IBEG(II),SOPRES%IEND(II),&
+      !         &'MIN,MAX=',SOPRES%VALMIN(II),SOPRES%VALMAX(II),&
+      !         &'(',SOPRES%VALMAX(II)-SOPRES%VALMIN(II),'/',&
+      !         &nbitcod,')'
+      IF (nbitcod >= 16) THEN
+        PRINT *,'-----> ERREUR FATALE : Groupe',II,'codage sur ',nbitcod,'bits'
+      END IF
+      IF (GELT > 1) THEN
+        ! Plus d'un element dans le groupe
+        IF ((17*GELT) < (17+4+IBE+nbitcod*GELT)) THEN
+          ! on prefere GELT groupes de 1 elt
+          DO JJ=SOPRES%IBEG(II),SOPRES%IEND(II)
+            ! 1 seul elt par groupe
+            CALL FILL_BBUFF(XTAB,1,1)
+            CALL FILL_BBUFF(XTAB,16,ITAB(JJ))
+          END DO
+          NBGROUPMOD = NBGROUPMOD+GELT-1
+        ELSE
+          CALL FILL_BBUFF(XTAB,1,0)
+          CALL FILL_BBUFF(XTAB,16,SOPRES%VALMIN(II))
+          CALL FILL_BBUFF(XTAB,4,nbitcod)
+          CALL FILL_BBUFF(XTAB,IBE,GELT)
+          IF (nbitcod > 0) THEN
+            DO JJ=SOPRES%IBEG(II),SOPRES%IEND(II)
+              ! stockage des GELT carts/VALMIN
+              CALL FILL_BBUFF(XTAB,nbitcod,ITAB(JJ)-SOPRES%VALMIN(II))
+            END DO
+          END IF
+        END IF
+      ELSE
+        ! 1 seul elt dans groupe
+        CALL FILL_BBUFF(XTAB,1,1)
+        CALL FILL_BBUFF(XTAB,16,SOPRES%VALMIN(II))
+      END IF
+    END DO
+    IF (NBGROUPMOD > SOPRES%NBGRP) THEN
+      ! we must change the number of elements 
+      CALL GET_FILLIDX(IDXSAVE,OFFSAVE) ! save the current idx/offset
+      CALL SET_FILLIDX(GRPIDX,GRPOFF)   
+      CALL FILL_BBUFF(XTAB,32,NBGROUPMOD)
+      CALL SET_FILLIDX(IDXSAVE,OFFSAVE) ! restore the current idx/offset
+    END IF
+  END IF
+END DO
+
+CALL GET_FILLIDX(IDXSAVE,OFFSAVE)
+KNBUSE=IDXSAVE+1
+
+DEALLOCATE(ITAB)
+
+CONTAINS 
+
+SUBROUTINE COMP_FOPEXT(PTAB,KTAB,KEXTCOD)
+REAL,    DIMENSION(:), INTENT(IN) :: PTAB
+INTEGER, DIMENSION(:), INTENT(OUT):: KTAB 
+INTEGER,               INTENT(OUT):: KEXTCOD
+
+LOGICAL,DIMENSION(SIZE(PTAB)) :: GMASK
+REAL,DIMENSION(SIZE(PTAB)) :: PTABWORK
+REAL :: XMIN1,XMAX1,XRANGE1
+REAL :: XMIN2,XMAX2,XRANGE2
+REAL :: XREF,XMAX,XCOEFF
+INTEGER :: INTRANGE
+INTEGER :: INDCOR   ! correction d'index pour la supression du min
+LOGICAL :: GMINEXCL,GMAXEXCL,GLOG
+INTEGER :: IEXTCOD2
+REAL, PARAMETER :: XUNDEF = 999.
+REAL, PARAMETER :: XUNDEFSURF =  1.E+20
+
+
+!! G. TANGUY avril 2010 : on change la valeur indfinie 999. a une valeur
+!indfinie plus grande que sera de faon certaine le max du champ s'il est
+!present. POur a on travaille dans le tableau de travail PTABWORK
+PTABWORK=PTAB
+WHERE(PTABWORK == XUNDEF)
+        PTABWORK=XUNDEFSURF
+END WHERE
+
+XMIN1=MINVAL(PTABWORK(:))
+XMAX1=MAXVAL(PTABWORK(:))
+XRANGE1=XMAX1-XMIN1
+IF (LPDEBUG) PRINT *,"XMIN1,XMAX1,XRANGE1 = ",XMIN1,XMAX1,XRANGE1
+
+IF (XRANGE1 > 0.) THEN
+  XMIN2=MINVAL(PTABWORK,MASK=PTABWORK>XMIN1)
+  XMAX2=MAXVAL(PTABWORK,MASK=PTABWORK<XMAX1)
+  XRANGE2 = XMAX2-XMIN2
+  IF (LPDEBUG) PRINT *,"XMIN2,XMAX2,XRANGE2 = ",XMIN2,XMAX2,XRANGE2
+  IF (XRANGE2 > 0.) THEN
+    GLOG     = .FALSE.
+    GMINEXCL = .FALSE.
+    GMAXEXCL = .FALSE.
+    GMASK(:) = .TRUE.
+    INDCOR = 0
+    KEXTCOD = JPNORM
+    INTRANGE=65535
+    XREF = XMIN1
+    XMAX = XMAX1
+
+    ! Check for range between 0 and 1 to convert to LOG values
+    IF (XMIN1 >= 0. .AND. XMAX1 < 1.) THEN
+      IF ((XMAX2/XMIN2)>10.) THEN 
+        GLOG = .TRUE.
+        KEXTCOD = JPOTHER
+        IEXTCOD2 = JPLOG
+        INTRANGE=INTRANGE-1
+        INDCOR = 1           ! On reserve la valeur 0 dans tous les cas
+        IF (XMIN1 == 0.0) THEN
+          XREF = LOG(XMIN2)
+          WHERE (PTABWORK < XMIN2)
+            KTAB  = 0
+            GMASK = .FALSE.
+          END WHERE
+        ELSE
+          XREF = LOG(XMIN1)
+        END IF
+        XMAX1 = LOG(XMAX1)
+        XMAX  = XMAX1
+        XMAX2 = LOG(XMAX2)
+        XRANGE2 = XMAX2 - XREF
+        IF (LPDEBUG) PRINT *,"EXTENCOD,  LOG conversion enabled : XMIN1, XREF, XMAX1, XMAX2 =",&
+             &XMIN1,XREF,XMAX1,XMAX2
+      END IF
+    ELSE
+      ! Check for MIN value exclusion
+      IF (XMIN1 == XUNDEFSURF .OR. (XMIN2-XMIN1) > XRANGE2) THEN
+        ! Min value excluded 
+        GMINEXCL = .TRUE.
+        XREF=XMIN2
+        INTRANGE=INTRANGE-1
+        INDCOR = 1
+        WHERE (PTABWORK < XMIN2)
+          KTAB = 0
+          GMASK = .FALSE.
+        END WHERE
+        IF (LPDEBUG) PRINT *,"EXTENCOD,     Min value isolated :",XMIN1
+        KEXTCOD = JPMINEXCL
+        IF (XMIN1 == XUNDEFSURF) THEN 
+                XMIN1=XUNDEF
+        END IF
+      END IF
+      ! Check for MAX value exclusion
+      IF (XMAX1 == XUNDEFSURF .OR. (XMAX1-XMAX2) > XRANGE2) THEN
+        ! Max value excluded
+        GMAXEXCL = .TRUE.
+        XMAX=XMAX2
+        INTRANGE=INTRANGE-1
+        WHERE (PTABWORK > XMAX2)
+          KTAB = 65535
+          GMASK = .FALSE.
+        END WHERE
+        
+        IF (GMINEXCL) THEN
+          KEXTCOD = JPMINMAXEXCL ! Min et Max exclus
+          IF (LPDEBUG) PRINT *,"EXTENCOD, and Max value isolated :",XMAX1
+        ELSE
+          KEXTCOD = JPMAXEXCL ! Max exclus
+          IF (LPDEBUG) PRINT *,"EXTENCOD, Max value isolated :",XMAX1
+        END IF
+        ! avril 2010 : on remet la valeur indefine de mesonh 999.
+        IF (XMAX1 == XUNDEFSURF) THEN 
+                XMAX1=XUNDEF
+        END IF
+      END IF
+    END IF
+    !
+    XCOEFF=(XMAX-XREF)/INTRANGE
+    IF (XCOEFF < PPFLOATMIN) THEN
+      XCOEFF = PPFLOATMIN
+      PRINT *, "very low range DATA : XCOEFF set to",XCOEFF
+    END IF
+    IF (LPDEBUG) PRINT *,"XCOEFF = ",XCOEFF
+    IF (GLOG) THEN
+      WHERE(GMASK)
+        KTAB = INDCOR + NINT((LOG(PTABWORK)-XREF)/XCOEFF)
+      END WHERE
+    ELSE
+      WHERE(GMASK)
+        KTAB = INDCOR + NINT((PTABWORK(:)-XREF)/XCOEFF)
+      END WHERE
+    END IF
+    IF (LPDEBUG) PRINT *,"KEXTCOD = ",KEXTCOD
+    CALL FILL_BBUFF(XTAB,3,KEXTCOD)
+    IF (GLOG)     CALL FILL_BBUFF(XTAB,3,IEXTCOD2)
+    IF (GMINEXCL) CALL FILL_BBUFF(XTAB,32,XMIN1)
+    IF (GMAXEXCL) CALL FILL_BBUFF(XTAB,32,XMAX1)
+    CALL FILL_BBUFF(XTAB,32,XREF)
+    CALL FILL_BBUFF(XTAB,32,XCOEFF)
+  ELSE
+    IF (XRANGE2 < 0.) THEN
+      ! only 2 values in PTAB array
+      !
+      ! KTAB(i)= 0 if PTAB(i)==XMIN1
+      !          1 if PTAB(i)==XMAX1
+      !
+      IF (LPDEBUG) PRINT *,"EXTENCOD, 2 values in array :",XMIN1,XMAX1
+        IF (XMAX1 == XUNDEFSURF) THEN 
+                XMAX1=XUNDEF
+        END IF
+        IF (XMIN1 == XUNDEFSURF) THEN 
+                XMIN1=XUNDEF
+        END IF
+      KEXTCOD = JP2VAL
+      CALL FILL_BBUFF(XTAB,3,KEXTCOD)
+      CALL FILL_BBUFF(XTAB,32,XMIN1)
+      CALL FILL_BBUFF(XTAB,32,XMAX1)
+      WHERE (PTABWORK < XMAX1)
+        KTAB = 0
+      ELSEWHERE
+        KTAB = 1
+      END WHERE
+    ELSE
+      ! XRANGE2 == 0. <==> XMIN2=XMAX2 
+      ! 3 values in PTAB array :
+      !
+      !          0 if PTAB(i)==XMIN1      ! KTAB(i)= 1 if PTAB(i)==XMIN2(=XMAX2)
+      !          2 if PTAB(i)==XMAX1
+      !
+      IF (LPDEBUG) PRINT *,"EXTENCOD, 3 values in array :",XMIN1,XMIN2,XMAX1
+        IF (XMAX1 == XUNDEFSURF) THEN 
+                XMAX1=XUNDEF
+        END IF
+        IF (XMIN1 == XUNDEFSURF) THEN 
+                XMIN1=XUNDEF
+        END IF
+
+      KEXTCOD = JP3VAL
+      CALL FILL_BBUFF(XTAB,3,KEXTCOD)
+      CALL FILL_BBUFF(XTAB,32,XMIN1)
+      CALL FILL_BBUFF(XTAB,32,XMIN2)
+      CALL FILL_BBUFF(XTAB,32,XMAX1)
+      WHERE (PTABWORK < XMIN2)
+        KTAB = 0
+      ELSEWHERE
+        KTAB = 1
+      END WHERE
+      WHERE (PTABWORK > XMIN2) KTAB = 2
+    END IF
+    
+  END IF
+ELSE
+        IF (XMIN1 == XUNDEFSURF) THEN 
+                XMIN1=XUNDEF
+        END IF
+
+  ! Constant array found : save its 32 bits real value.
+  KEXTCOD=JPCONST
+  CALL FILL_BBUFF(XTAB,3,KEXTCOD)
+  CALL FILL_BBUFF(XTAB,32,XMIN1)
+  IF (LPDEBUG) PRINT *,"EXTENCOD, constant array : ",XMIN1
+END IF
+END SUBROUTINE COMP_FOPEXT
+
+END SUBROUTINE COMPRESS_FIELD