From b5709af577fbbb5aeb3ca3b8f072bb812513f60f Mon Sep 17 00:00:00 2001 From: Enrico Fucile Date: Thu, 28 May 2015 09:56:33 +0100 Subject: [PATCH] ECC-129 fixed --- examples/F90/bufr_read_scatterometer.f90 | 9 +- examples/F90/bufr_read_scatterometer.sh | 4046 +++++++++++++++++++++- fortran/eccodes_f90_tail.f90 | 10 +- fortran/grib_f90_tail.f90 | 6 +- 4 files changed, 4060 insertions(+), 11 deletions(-) diff --git a/examples/F90/bufr_read_scatterometer.f90 b/examples/F90/bufr_read_scatterometer.f90 index 1a896ec39..56b41a9e3 100644 --- a/examples/F90/bufr_read_scatterometer.f90 +++ b/examples/F90/bufr_read_scatterometer.f90 @@ -24,8 +24,9 @@ integer :: ifile integer :: iret integer :: ibufr integer :: i, count=0 -integer(kind=4) :: numObs +integer(kind=4) :: numObs,ii real(kind=8), dimension(:), allocatable :: latVal,lonVal,bscatterVal +real(kind=8), dimension(:), allocatable :: year call codes_open_file(ifile,'../../data/bufr/asca_139.bufr','r') @@ -59,6 +60,12 @@ real(kind=8), dimension(:), allocatable :: latVal,lonVal,bscatterVal !Get longitude (for all the subsets) call codes_get(ibufr,'latitude',lonVal); + allocate(year(numObs)) + call codes_get(ibufr,'year',year); + do ii= 1, size(year) + write(*,*) 'year(',ii,')=',year(ii) + enddo + !Get backScatter for beam two. We use an access by condition for this key. !(for all the subsets) call codes_get(ibufr,'/beamIdentifier=2/backscatter',bscatterVal); diff --git a/examples/F90/bufr_read_scatterometer.sh b/examples/F90/bufr_read_scatterometer.sh index 310838bef..24a994538 100755 --- a/examples/F90/bufr_read_scatterometer.sh +++ b/examples/F90/bufr_read_scatterometer.sh @@ -16,6 +16,7 @@ label="bufr_read_scatterometer_f" #Define tmp file fTmp=${label}.tmp.txt +ref=${label}.tmp.ref rm -f $fTmp | true #We check "asca_139.bufr". The path is @@ -26,10 +27,4049 @@ REDIRECT=/dev/null #Write the key values into a file ${examples_dir}/f_bufr_read_scatterometer 2> $REDIRECT > $fTmp -#TODO: check the results +# check the results -#cat $fTmp +cat > $ref <