/[MITgcm]/MITgcm/pkg/exf/exf_getforcing.F
ViewVC logotype

Contents of /MITgcm/pkg/exf/exf_getforcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.3 - (show annotations) (download)
Thu Feb 7 20:00:09 2002 UTC (22 years, 4 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint44e_post, checkpoint46l_post, checkpoint46g_pre, checkpoint46f_post, checkpoint44f_post, checkpoint46b_post, checkpoint46l_pre, chkpt44d_post, checkpoint44e_pre, checkpoint46d_pre, checkpoint45d_post, checkpoint46j_pre, checkpoint44h_pre, checkpoint46a_post, checkpoint46j_post, checkpoint46k_post, chkpt44c_pre, checkpoint45a_post, checkpoint44g_post, checkpoint46e_pre, checkpoint45b_post, checkpoint46b_pre, release1_final_v1, checkpoint46c_pre, checkpoint46, checkpoint44b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint45c_post, checkpoint44h_post, checkpoint46g_post, checkpoint46i_post, checkpoint46c_post, checkpoint46e_post, checkpoint45, checkpoint46h_post, chkpt44c_post, checkpoint44f_pre, checkpoint46d_post
Branch point for: release1_final
Changes since 1.2: +9 -18 lines 
o merge of relevant stuff from the ecco-branch:
  - genmake: removed $S64 overwrite for case SunOS
  - pkg/exf: update and corrections for field swapping and obcs
  - pkg/ecco: parameter lists for the_model_main, the_main_loop
              harmonized between ECCO and MITgcm
  - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice
                  updated checkpointing_lev... lists for obcs
  - model/src: minor changes in forward_step, plot_field
               added directive for divided adjoint in the_main_loop
  - pkg/mdsio: added mdsio_gl_slice
1 c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_getforcing.F,v 1.2.6.1 2002/02/06 23:18:53 heimbach Exp $
2
3 #include "CPP_OPTIONS.h"
4
5 CBOI
6 C
7 C !TITLE: EXTERNAL FORCING
8 C !AUTHORS: mitgcm developers ( support@mitgcm.org )
9 C !AFFILIATION: Massachussetts Institute of Technology
10 C !DATE:
11 C !INTRODUCTION: External forcing package
12 c \bv
13 c * The external forcing package, in conjunction with the
14 c calendar package, enables the handling of realistic forcing
15 c fields of differing temporal forcing patterns.
16 c * It comprises climatological restoring and relaxation
17 c * Bulk formulae are implemented to convert atmospheric fields
18 c to surface fluxes.
19 c
20 C !CALLING SEQUENCE:
21 c ...
22 c exf_getforcing
23 c |
24 c |-- exf_getclim (get climatological fields used e.g. for relax.)
25 c |
26 c |-- exf_getsurfacefluxes
27 c | |
28 c | |-- exf_getffields <- this one does almost everything
29 c | | 1. reads in fields, either flux or atmos. state,
30 c | | depending on CPP options
31 c | | 2. If forcing and control is flux, we're already done here
32 c | | 3. If forcing is atmos. state, then
33 c | | (a) if control is atmos. state, then the control variable
34 c | | anomalies are read here
35 c | | * ctrl_getatemp
36 c | | * ctrl_getaqh
37 c | | * ctrl_getuwind
38 c | | * ctrl_getvwind
39 c | | 4. flux fields are interpolated on current time
40 c | | (b) next the flux fields are computed via bulk formulae
41 c | | 5. if forcing and control is flux, then the control vector
42 c | | anomalies are read here
43 c | | * ctrl_getheatflux
44 c | | * ctrl_getsaltflux
45 c | | * ctrl_getzonstress
46 c | | * call ctrl_getmerstress
47 c | |
48 c | |-- exf_obcs
49 c | | If open boundaries are prescribed externally,
50 c | | OB values for current time step are read here
51 c | | 1. for each boundary (north, south, west, east)
52 c | | sliced fields are read for T,S,U,V
53 c | | 2. if OB's are part of control vector,
54 c | | control variable anomalies are added to OB values
55 c | | * ctrl_getobcsn
56 c | | * ctrl_getobcss
57 c | | * ctrl_getobcsw
58 c | | * ctrl_getobcse
59 c | |
60 c | |-- exf_mapfields
61 c | | Forcing fields from exf package are mapped onto
62 c | | mitgcm forcing arrays.
63 c | | Mapping enables a runtime rescaling of fields
64 c
65 c \ev
66 CEOI
67
68 CBOP
69 C !ROUTINE: exf_getforcing
70 C !INTERFACE:
71 subroutine exf_getforcing( mytime, myiter, mythid )
72
73 C !DESCRIPTION: \bv
74 c *=================================================================
75 c | SUBROUTINE exf_getforcing
76 c *=================================================================
77 c o Get the forcing fields for the current time step. The switches
78 c for the inclusion of the individual forcing components have to
79 c be set in EXF_CPPOPTIONS.h (or ECCO_CPPOPTIONS.h).
80 c A note on surface fluxes:
81 c The MITgcmUV's vertical coordinate z is positive upward.
82 c This implies that a positive flux is out of the ocean
83 c model. However, the wind stress forcing is not treated
84 c this way. A positive zonal wind stress accelerates the
85 c model ocean towards the east.
86 c started: eckert@mit.edu, heimbach@mit.edu, ralf@ocean.mit.edu
87 c *=================================================================
88 c | SUBROUTINE exf_getforcing
89 c *=================================================================
90 C \ev
91
92 C !USES:
93 implicit none
94
95 c == global variables ==
96
97 C !INPUT/OUTPUT PARAMETERS:
98 c == routine arguments ==
99 integer mythid
100 integer myiter
101 _RL mytime
102
103 C !LOCAL VARIABLES:
104 c == local variables ==
105
106 c == end of interface ==
107 CEOP
108
109 c Get values of climatological fields.
110 call exf_getclim( mytime, myiter, mythid )
111
112 c Get the surface forcing fields.
113 call exf_getffields( mytime, myiter, mythid )
114
115 c Get values of the surface flux anomalies.
116 call exf_getsurfacefluxes( mytime, myiter, mythid )
117
118 c Map the forcing fields onto the corresponding model fields.
119 call exf_mapfields( mythid )
120
121 end
122
  ViewVC Help
Powered by ViewVC 1.1.22