pkg/exf/exf_getclim.F

c $Header: /u/gcmpack/development/heimbach/ecco_env/pkg/exf/exf_getclim.F,v 1.1 2001/02/02 19:43:46 heimbach Exp $

#include "EXF_CPPOPTIONS.h"


      subroutine exf_getclim(
     I                        mytime,
     I                        myiter,
     I                        mythid
     &                      )

c     ==================================================================
c     SUBROUTINE exf_getclim
c     ==================================================================
c
c     o Get the climatogy fields for the current time step. The switches
c       for the inclusion of the individual forcing components have to
c       be set in EXF_CPPOPTIONS.h .
c
c       A note on surface fluxes:
c
c       The MITgcmUV's vertical coordinate z is positive upward.
c       This implies that a positive flux is out of the ocean
c       model. However, the wind stress forcing is not treated
c       this way. A positive zonal wind stress accelerates the
c       model ocean towards the east.
c
c     started: Ralf.Giering@FastOpt.de 25-Mai-2000
c
c     ==================================================================
c     SUBROUTINE exf_getclim
c     ==================================================================

      implicit none

c     == global variables ==
#include "EEPARAMS.h"
#include "SIZE.h"
#include "exf_fields.h"
#include "exf_clim_fields.h"

c     == routine arguments ==

c     mythid - thread number for this instance of the routine.

      integer mythid
      integer myiter
      _RL     mytime

c     == local variables ==

c     == end of interface ==

#ifdef ALLOW_CLIMTEMP_RELAXATION
c     Get values of climatological temperature fields.
      call exf_set_climtemp ( climtemp
     I                      , mytime, myiter, mythid )
c--   Update the tile edges.
      _EXCH_XY_R8(climtemp, mythid)
#else
c     Do not include relaxation to climatological temperature fields.
#endif

#ifdef ALLOW_CLIMSALT_RELAXATION
c     Get values of climatological salinity fields.
      call exf_set_climsalt ( climsalt
     I                      , mytime, myiter, mythid )
c--   Update the tile edges.
      _EXCH_XY_R8(climsalt, mythid)
#else
c     Do not include relaxation to climatological salinity fields.
#endif

#ifdef ALLOW_CLIMSST_RELAXATION
c     Get values of climatological sst fields.
      call exf_set_climsst ( climsst
     I                     , mytime, myiter, mythid )
c--   Update the tile edges.
      _EXCH_XY_R8(climsst, mythid)
#else
c     Do not include relaxation to climatological sst fields.
#endif

#ifdef ALLOW_CLIMSSS_RELAXATION
c     Get values of climatological sst fields.
      call exf_set_climsss ( climsss
     I                     , mytime, myiter, mythid )
c--   Update the tile edges.
      _EXCH_XY_R8(climsss, mythid)
#else
c     Do not include relaxation to climatological sst fields.
#endif

      end
1	c $Header: /u/gcmpack/development/heimbach/ecco_env/pkg/exf/exf_getclim.F,v 1.1 2001/02/02 19:43:46 heimbach Exp $
2
3	#include "EXF_CPPOPTIONS.h"
4
5
6	subroutine exf_getclim(
7	I mytime,
8	I myiter,
9	I mythid
10	& )
11
12	c ==================================================================
13	c SUBROUTINE exf_getclim
14	c ==================================================================
15	c
16	c o Get the climatogy fields for the current time step. The switches
17	c for the inclusion of the individual forcing components have to
18	c be set in EXF_CPPOPTIONS.h .
19	c
20	c A note on surface fluxes:
21	c
22	c The MITgcmUV's vertical coordinate z is positive upward.
23	c This implies that a positive flux is out of the ocean
24	c model. However, the wind stress forcing is not treated
25	c this way. A positive zonal wind stress accelerates the
26	c model ocean towards the east.
27	c
28	c started: Ralf.Giering@FastOpt.de 25-Mai-2000
29	c
30	c ==================================================================
31	c SUBROUTINE exf_getclim
32	c ==================================================================
33
34	implicit none
35
36	c == global variables ==
37	#include "EEPARAMS.h"
38	#include "SIZE.h"
39	#include "exf_fields.h"
40	#include "exf_clim_fields.h"
41
42	c == routine arguments ==
43
44	c mythid - thread number for this instance of the routine.
45
46	integer mythid
47	integer myiter
48	_RL mytime
49
50	c == local variables ==
51
52	c == end of interface ==
53
54	#ifdef ALLOW_CLIMTEMP_RELAXATION
55	c Get values of climatological temperature fields.
56	call exf_set_climtemp ( climtemp
57	I , mytime, myiter, mythid )
58	c-- Update the tile edges.
59	_EXCH_XY_R8(climtemp, mythid)
60	#else
61	c Do not include relaxation to climatological temperature fields.
62	#endif
63
64	#ifdef ALLOW_CLIMSALT_RELAXATION
65	c Get values of climatological salinity fields.
66	call exf_set_climsalt ( climsalt
67	I , mytime, myiter, mythid )
68	c-- Update the tile edges.
69	_EXCH_XY_R8(climsalt, mythid)
70	#else
71	c Do not include relaxation to climatological salinity fields.
72	#endif
73
74	#ifdef ALLOW_CLIMSST_RELAXATION
75	c Get values of climatological sst fields.
76	call exf_set_climsst ( climsst
77	I , mytime, myiter, mythid )
78	c-- Update the tile edges.
79	_EXCH_XY_R8(climsst, mythid)
80	#else
81	c Do not include relaxation to climatological sst fields.
82	#endif
83
84	#ifdef ALLOW_CLIMSSS_RELAXATION
85	c Get values of climatological sst fields.
86	call exf_set_climsss ( climsss
87	I , mytime, myiter, mythid )
88	c-- Update the tile edges.
89	_EXCH_XY_R8(climsss, mythid)
90	#else
91	c Do not include relaxation to climatological sst fields.
92	#endif
93
94	end