Interpolation of lat, lon and azimuth [0.5 d]

[Ian Culverwell]

https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/common/Makefile.am https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/common/ropp_pp.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/common/ropp_pp_diag2roprof.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/common/ropp_pp_interpolate_latlonaz.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/common/ropp_pp_read_config.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/common/ropp_pp_types.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/icorr/ropp_pp_fit_model_refraction.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/icorr/ropp_pp_fit_model_refraction_new.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/icorr/ropp_pp_ionospheric_correction.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/tools/ropp_pp_invert_tool.f90 https://trac.romsaf.org/ropp/changeset/5427/ropp_src/branches/dev/Share/dmi_trunk_9.0/ropp_pp/tools/ropp_pp_occ_tool.f90

Original changesets:

Interpolation of lat, lon and azimuth: https://trac.romsaf.org/ropp/changeset/4607

[Ian Culverwell]

Algorithm not clearly documented, e.g.

  DO i=1,nx     
...
   IF(use_azm) THEN  ! doesn't see to make sense to me, if azm = bearing wrt due N
    razm = azmtp(i)*deg2rad
    vvec(1,i) = -1.0_wp*cos(rlat)*cos(rlon)*cos(razm) - sin(rlon)*sin(razm)
    vvec(2,i) = -1.0_wp*cos(rlat)*sin(rlon)*cos(razm) + cos(rlon)*sin(razm)
    vvec(3,i) = sin(rlat)*cos(razm)
   ENDIF
  ENDDO

I had to replace

  USE ropp_pp

by

  USE ropp_pp, not_this => ropp_pp_interpolate_latlonaz

No idea why DMI didn't need to do this.

Incorporated in ROPP-11 at r6787.

[Ian Culverwell]

Tests out OK, so closing ticket.