uncompress_field.F90

Go to the documentation of this file.

SUBROUTINE uncompress_field(KLONG,PSEUIL,PFIELD_IN,PFIELD_OUT)

IMPLICIT NONE
 
INTEGER*4, INTENT(IN) :: KLONG
REAL, INTENT(IN) :: PSEUIL
REAL, DIMENSION(:), INTENT(IN) :: PFIELD_IN
REAL, DIMENSION(:), INTENT(OUT) :: PFIELD_OUT
INTEGER :: ICPT, I, K

icpt = 0

i = 1

pfield_out(:) = 0.
!
! boucle sur les colonnes
DO 
           
  ! si on a dépassé la dernière colonne, on sort de la boucle
  IF (icpt>=klong) EXIT

  ! si la valeur est valide
  IF (pfield_in(i)<pseuil) THEN

    ! on la met dans lwrite à l'indice icpt
    icpt = icpt + 1
    pfield_out(icpt) = pfield_in(i)
!!!!!!!!!!!!!!!!!!!!test temporary: to remove after
    IF (mod(pfield_out(icpt),100.)==0) pfield_out(icpt)=0.
!!!!!!!!!!!!!!!!!!!!!! to remove after
    ! on incrémente i
    i = i+1

  ELSE

    !ideb = icpt + 1
    DO k = 1,nint(pfield_in(i)-(pseuil))
      icpt = icpt + 1
      IF (icpt>klong) EXIT
      pfield_out(icpt) = 0.
    ENDDO

    i = i+1

  ENDIF

ENDDO

END SUBROUTINE uncompress_field