Implemented option to exclude output of the level 1b bending angle [2 d]

[Ian Culverwell]

https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_io/ropp/ropp_io.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_io/ropp/ropp_io_read.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_io/ropp/ropp_io_read_ncdf_get.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_io/tools/eum2bufr_ec.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_io/tools/eum2ropp.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_io/tools/ropp2bufr_ec.f90 (similar/consistent updates to ropp2bufr_mo.f90 has not been made since https://trac.romsaf.org/ropp/changeset/3430)

Original changesets:

Implemented option to exclude output of the level 1b bending angle: https://trac.romsaf.org/ropp/changeset/5362

[Ian Culverwell]

This has been done at r5970. (This includes changes to the eum2ropp man page.) There is no change to the default behaviour of extracting level 1b data with eum2ropp (and, reasonably, there is no option but to keep extracting it with eum2bufr and ropp2bufr). But if we run eum2ropp with the -x option on the default EUM file, the output contains no level 1b data. OK.

Note that other elements of this changeset alter

 pcd = 23 ;

 overall_qual = 0 ;

(i.e. PCD summary, offline, rising and bangle are all set, and the profile is of poor quality) to

 pcd = 6 ;

 overall_qual = 100 ;

(i.e. PCD offline and rising are set, and the profile is of good quality).

This make sense in view of the change in coding from

  data%PCD = 0

  ! Overall nominal / non-nominal
  CALL ncdf_getvar('/quality/overall_quality_ok', readbyte1)
  IF (readbyte1 == 0) data%PCD = IBSET(data%PCD, PCD_summary)

  ! BA nominal (set to same value as overall)
  IF (readbyte1 == 0) data%PCD = IBSET(data%PCD, PCD_bangle) ! set to same value as overall

  ! quality indicator, currently just set to 100% or 0%
  data%overall_qual = 100.0_wp * readbyte1
  data%units%overall_qual = '%'

to

  data%PCD = 0
  data%overall_qual = 100.0_wp
  data%units%overall_qual = '%'

  IF (getlevel1b) THEN

    ! Overall nominal / non-nominal
    tdir = '/quality/'
    CALL ncdf_getvar(TRIM(tdir)//'overall_quality_ok', readbyte1)
    IF (readbyte1 == 0) data%PCD = IBSET(data%PCD, PCD_summary)

    ! BA nominal (set to same value as overall)
    IF (readbyte1 == 0) data%PCD = IBSET(data%PCD, PCD_bangle) ! set to same value as overall

    ! quality indicator, currently just set to 100% or 0%
    data%overall_qual = 100.0_wp * readbyte1

  ENDIF

given that /quality/overall_quality_ok = 0 for this dataset.

Closing ticket.