FV3 Bundle
StateWriter.h
Go to the documentation of this file.
1 /*
2  * (C) Copyright 2009-2016 ECMWF.
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  * In applying this licence, ECMWF does not waive the privileges and immunities
7  * granted to it by virtue of its status as an intergovernmental organisation nor
8  * does it submit to any jurisdiction.
9  */
10 
11 #ifndef OOPS_BASE_STATEWRITER_H_
12 #define OOPS_BASE_STATEWRITER_H_
13 
14 #include "oops/base/PostBase.h"
15 
16 #include "eckit/config/LocalConfiguration.h"
17 #include "oops/util/DateTime.h"
18 
19 namespace oops {
20 
21 /// Handles writing-out of forecast fields.
22 /*!
23  * Write out forecast fields.
24  */
25 
26 template <typename FLDS> class StateWriter : public PostBase<FLDS> {
27  public:
28  StateWriter(const util::DateTime & start, const eckit::Configuration & conf):
29  PostBase<FLDS>(start, conf), ppConfig_(conf) {}
30  ~StateWriter() {}
31 
32  private:
33  const eckit::LocalConfiguration ppConfig_;
34  void doProcessing(const FLDS & xx) override {xx.write(ppConfig_);}
35 };
36 
37 } // namespace oops
38 
39 #endif // OOPS_BASE_STATEWRITER_H_
PostBase.h
oops::StateWriter
Handles writing-out of forecast fields.
Definition: StateWriter.h:26
conf
Definition: conf.py:1
oops::StateWriter::StateWriter
StateWriter(const util::DateTime &start, const eckit::Configuration &conf)
Definition: StateWriter.h:28
oops::PostBase
Handles post-processing of model fields.
Definition: PostBase.h:33
oops
The namespace for the main oops code.
Definition: ErrorCovarianceFV3JEDI.h:24
oops::StateWriter::ppConfig_
const eckit::LocalConfiguration ppConfig_
Definition: StateWriter.h:33
start
************************************************************************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:! ***********************************************************************subroutine READ_RECORD_CORE_(unit, field, nwords, data, start, axsiz) integer, intent(in) ::unit type(fieldtype), intent(in) ::field integer, intent(in) ::nwords MPP_TYPE_, intent(inout) ::data(nwords) integer, intent(in) ::start(:), axsiz(:) integer(SHORT_KIND) ::i2vals(nwords)!rab used in conjunction with transfer intrinsic to determine size of a variable integer(KIND=1) ::one_byte(8) integer ::word_sz!#ifdef __sgi integer(INT_KIND) ::ivals(nwords) real(FLOAT_KIND) ::rvals(nwords)!#else! integer ::ivals(nwords)! real ::rvals(nwords)!#endif real(DOUBLE_KIND) ::r8vals(nwords) pointer(ptr1, i2vals) pointer(ptr2, ivals) pointer(ptr3, rvals) pointer(ptr4, r8vals) if(mpp_io_stack_size< nwords) call mpp_io_set_stack_size(nwords) call mpp_error(FATAL, 'MPP_READ currently requires use_netCDF option') end subroutine READ_RECORD_CORE_ subroutine READ_RECORD_(unit, field, nwords, data, time_level, domain, position, tile_count, start_in, axsiz_in)!routine that is finally called by all mpp_read routines to perform the read!a non-netCDF record contains:! field ID! a set of 4 coordinates(is:ie, js:je) giving the data subdomain! a timelevel and a timestamp(=NULLTIME if field is static)! 3D real data(stored as 1D)!if you are using direct access I/O, the RECL argument to OPEN must be large enough for the above!in a global direct access file, record position on PE is given by %record.!Treatment of timestamp:! We assume that static fields have been passed without a timestamp.! Here that is converted into a timestamp of NULLTIME.! For non-netCDF fields, field is treated no differently, but is written! with a timestamp of NULLTIME. There is no check in the code to prevent! the user from repeatedly writing a static field. integer, intent(in) ::unit, nwords type(fieldtype), intent(in) ::field MPP_TYPE_, intent(inout) ::data(nwords) integer, intent(in), optional ::time_level type(domain2D), intent(in), optional ::domain integer, intent(in), optional ::position, tile_count integer, intent(in), optional ::start_in(:), axsiz_in(:) integer, dimension(size(field%axes(:))) ::start, axsiz integer ::tlevel !, subdomain(4) integer ::i, error, is, ie, js, je, isg, ieg, jsg, jeg type(domain2d), pointer ::io_domain=> tlevel if(PRESENT(start_in) .AND. PRESENT(axsiz_in)) then if(size(start(! the data domain and compute domain must refer to the subdomain being passed ! In this ! since that attempts to gather all data on PE start
Definition: mpp_read_2Ddecomp.h:197
oops::StateWriter::doProcessing
void doProcessing(const FLDS &xx) override
Actual processing.
Definition: StateWriter.h:34