FV3 Bundle
ObsOpBaseQG.h
Go to the documentation of this file.
1 /*
2  * (C) Copyright 2017-2018 UCAR
3  *
4  * This software is licensed under the terms of the Apache Licence Version 2.0
5  * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
6  */
7 
8 #ifndef QG_MODEL_OBSOPBASEQG_H_
9 #define QG_MODEL_OBSOPBASEQG_H_
10 
11 #include <map>
12 #include <string>
13 
14 #include <boost/noncopyable.hpp>
15 
16 #include "eckit/config/Configuration.h"
17 #include "model/ObsSpaceQG.h"
18 #include "oops/base/Variables.h"
19 #include "oops/util/abor1_cpp.h"
20 #include "oops/util/Printable.h"
21 
22 namespace qg {
23 class GomQG;
24 class ObsBias;
25 class ObsVecQG;
26 
27 // -----------------------------------------------------------------------------
28 /// Base class for observation operators
29 
30 class ObsOpBaseQG : public util::Printable,
31  private boost::noncopyable {
32  public:
34  virtual ~ObsOpBaseQG() {}
35 
36 /// Obs Operator
37  virtual void simulateObs(const GomQG &, ObsVecQG &, const ObsBias &) const = 0;
38 
39 /// Other
40  virtual const oops::Variables & variables() const = 0; // Required from Model
41 
42  private:
43  virtual void print(std::ostream &) const = 0;
44 };
45 
46 // -----------------------------------------------------------------------------
47 
48 /// Obs Operator Factory
49 class ObsOpFactory {
50  public:
51  static ObsOpBaseQG * create(const ObsSpaceQG &, const eckit::Configuration &);
52  virtual ~ObsOpFactory() { getMakers().clear(); }
53  protected:
54  explicit ObsOpFactory(const std::string &);
55  private:
56  virtual ObsOpBaseQG * make(const ObsSpaceQG &, const eckit::Configuration &) = 0;
57  static std::map < std::string, ObsOpFactory * > & getMakers() {
58  static std::map < std::string, ObsOpFactory * > makers_;
59  return makers_;
60  }
61 };
62 
63 // -----------------------------------------------------------------------------
64 
65 template<class T>
66 class ObsOpMaker : public ObsOpFactory {
67  virtual ObsOpBaseQG * make(const ObsSpaceQG & odb, const eckit::Configuration & conf)
68  { return new T(odb, conf); }
69  public:
70  explicit ObsOpMaker(const std::string & name) : ObsOpFactory(name) {}
71 };
72 
73 // -----------------------------------------------------------------------------
74 
75 } // namespace qg
76 
77 #endif // QG_MODEL_OBSOPBASEQG_H_
GomQG class to handle local model values for QG model.
Definition: GomQG.h:31
virtual void simulateObs(const GomQG &, ObsVecQG &, const ObsBias &) const =0
Obs Operator.
ObsOpMaker(const std::string &name)
Definition: ObsOpBaseQG.h:70
Base class for observation operators.
Definition: ObsOpBaseQG.h:30
virtual ObsOpBaseQG * make(const ObsSpaceQG &, const eckit::Configuration &)=0
Definition: conf.py:1
virtual ~ObsOpFactory()
Definition: ObsOpBaseQG.h:52
Obs Operator Factory.
Definition: ObsOpBaseQG.h:49
static std::map< std::string, ObsOpFactory *> & getMakers()
Definition: ObsOpBaseQG.h:57
virtual ObsOpBaseQG * make(const ObsSpaceQG &odb, const eckit::Configuration &conf)
Definition: ObsOpBaseQG.h:67
virtual const oops::Variables & variables() const =0
Other.
character(len=32) name
static ObsOpBaseQG * create(const ObsSpaceQG &, const eckit::Configuration &)
Definition: ObsOpBaseQG.cc:28
virtual void print(std::ostream &) const =0
Class to handle observation bias parameters.
************************************************************************GNU Lesser General Public License **This file is part of the GFDL Flexible Modeling System(FMS). ! *! *FMS is free software without even the implied warranty of MERCHANTABILITY or *FITNESS FOR A PARTICULAR PURPOSE See the GNU General Public License *for more details **You should have received a copy of the GNU Lesser General Public *License along with FMS If see< http:! ***********************************************************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! MPP_TRANSMIT !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! subroutine MPP_TRANSMIT_(put_data, put_len, to_pe, get_data, get_len, from_pe, block, tag, recv_request, send_request)!a message-passing routine intended to be reminiscent equally of both MPI and SHMEM!put_data and get_data are contiguous MPP_TYPE_ arrays!at each call, your put_data array is put to to_pe 's get_data! your get_data array is got from from_pe 's put_data!i.e we assume that typically(e.g updating halo regions) each PE performs a put _and_ a get!special PE designations:! NULL_PE:to disable a put or a get(e.g at boundaries)! ANY_PE:if remote PE for the put or get is to be unspecific! ALL_PES:broadcast and collect operations(collect not yet implemented)!ideally we would not pass length, but this f77-style call performs better(arrays passed by address, not descriptor)!further, this permits< length > contiguous words from an array of any rank to be passed(avoiding f90 rank conformance check)!caller is responsible for completion checks(mpp_sync_self) before and after integer, intent(in) ::put_len, to_pe, get_len, from_pe MPP_TYPE_, intent(in) ::put_data(*) MPP_TYPE_, intent(out) ::get_data(*) logical, intent(in), optional ::block integer, intent(in), optional ::tag integer, intent(out), optional ::recv_request, send_request logical ::block_comm integer ::i MPP_TYPE_, allocatable, save ::local_data(:) !local copy used by non-parallel code(no SHMEM or MPI) integer ::comm_tag integer ::rsize if(.NOT.module_is_initialized) call mpp_error(FATAL, 'MPP_TRANSMIT:You must first call mpp_init.') if(to_pe.EQ.NULL_PE .AND. from_pe.EQ.NULL_PE) return block_comm=.true. if(PRESENT(block)) block_comm=block if(debug) then call SYSTEM_CLOCK(tick) write(stdout_unit,'(a, i18, a, i6, a, 2i6, 2i8)')&'T=', tick, ' PE=', pe, ' MPP_TRANSMIT begin:to_pe, from_pe, put_len, get_len=', to_pe, from_pe, put_len, get_len end if comm_tag=DEFAULT_TAG if(present(tag)) comm_tag=tag!do put first and then get if(to_pe.GE.0 .AND. to_pe.LT.npes) then!use non-blocking sends if(debug .and.(current_clock.NE.0)) call SYSTEM_CLOCK(start_tick)!z1l:truly non-blocking send.! if(request(to_pe).NE.MPI_REQUEST_NULL) then !only one message from pe-> to_pe in queue *PE waiting for to_pe ! call error else get_len so only do gets but you cannot have a pure get with MPI call a get means do a wait to ensure put on remote PE is complete error call increase mpp_nml request_multiply call MPP_TRANSMIT get_len end if return end subroutine MPP_TRANSMIT_ ! MPP_BROADCAST ! subroutine but that doesn t allow !broadcast to a subset of PEs This version and mpp_transmit will remain !backward compatible intent(inout) & T
ObsVecQG class to handle vectors in observation space for QG model.
Definition: ObsVecQG.h:27
virtual ~ObsOpBaseQG()
Definition: ObsOpBaseQG.h:34
Wrapper around ObsHelpQG, mostly to hide the factory.
Definition: ObsSpaceQG.h:35
ObsOpFactory(const std::string &)
Definition: ObsOpBaseQG.cc:18
The namespace for the qg model.