◆ count_distinct_var6d()

integer function volgrid6d_var_class::count_distinct_var6d	(	type(volgrid6d_var), dimension(:), intent(in)	vect,
		logical, dimension(:), intent(in), optional	mask,
		logical, intent(in), optional	back
	)

private

conta gli elementi distinti in vect

Definizione alla linea 597 del file volgrid6d_var_class.F90.

 100 CONTINUE
 
 rewind(un)
 ALLOCATE(conv_fwd(n),stat=stallo)
 IF (stallo /= 0) THEN
   CALL l4f_log(l4f_fatal,"allocating memory")
   CALL raise_fatal_error()
 ENDIF
 
 conv_fwd(:) = vg6d_v7d_var_conv_miss
 CALL import_var_conv(un, conv_fwd)
 CLOSE(un)
 
 ! backward, bufr to grib
 un = open_package_file('vargrib2bufr.csv', filetype_data)
 ! use the same file for now
 !un = open_package_file('varbufr2grib.csv', filetype_data)
 n=0
 DO WHILE(.true.)
   READ(un,*,end=300)
   n = n + 1
 ENDDO
 
 300 CONTINUE
 
 rewind(un)
 ALLOCATE(conv_bwd(n),stat=stallo)
 IF (stallo /= 0) THEN
   CALL l4f_log(l4f_fatal,"allocating memory")
   CALL raise_fatal_error()
 end if
 
 conv_bwd(:) = vg6d_v7d_var_conv_miss
 CALL import_var_conv(un, conv_bwd)
 DO i = 1, n
   conv_bwd(i)%c_func%a = 1./conv_bwd(i)%c_func%a
   conv_bwd(i)%c_func%b = - conv_bwd(i)%c_func%b
 ENDDO
 CLOSE(un)
 
 CONTAINS
 
 SUBROUTINE import_var_conv(un, conv_type)
 INTEGER, INTENT(in) :: un
 TYPE(vg6d_v7d_var_conv), INTENT(out) :: conv_type(:)
 
 INTEGER :: i
 TYPE(csv_record) :: csv
 CHARACTER(len=1024) :: line
 CHARACTER(len=10) :: btable
 INTEGER :: centre, category, number, discipline
 
 DO i = 1, SIZE(conv_type)
   READ(un,'(A)',end=200)line
   CALL init(csv, line)
   CALL csv_record_getfield(csv, btable)
   CALL csv_record_getfield(csv) ! skip fields for description and unit,
   CALL csv_record_getfield(csv) ! they correspond to grib information, not bufr Btable
   CALL init(conv_type(i)%v7d_var, btable=btable)
 
   CALL csv_record_getfield(csv, centre)
   CALL csv_record_getfield(csv, category)
   CALL csv_record_getfield(csv, number)
   CALL csv_record_getfield(csv, discipline)
   CALL init(conv_type(i)%vg6d_var, centre=centre, category=category, &
    number=number, discipline=discipline) ! controllare l'ordine
 
   CALL csv_record_getfield(csv, conv_type(i)%c_func%a)
   CALL csv_record_getfield(csv, conv_type(i)%c_func%b)