/[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.6 - (show annotations) (download)
Tue Feb 18 05:33:54 2003 UTC (21 years, 4 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint50c_post, checkpoint50c_pre, checkpoint48i_post, checkpoint50, checkpoint50d_post, checkpoint50b_pre, checkpoint48f_post, checkpoint48h_post, checkpoint50a_post, checkpoint50f_pre, checkpoint50e_pre, checkpoint50e_post, checkpoint50d_pre, checkpoint49, checkpoint48g_post, checkpoint50b_post
Changes since 1.5: +2 -2 lines 
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
1 c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_getforcing.F,v 1.2.4.4 2003/01/12 08:20:10 dimitri 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_getffields <- this one does almost everything
27 c | | 1. reads in fields, either flux or atmos. state,
28 c | | depending on CPP options
29 c | | 2. If forcing and control is flux, we're already done here
30 c | | 3. If forcing is atmos. state, then
31 c | | (a) if control is atmos. state, then the control variable
32 c | | anomalies are read here
33 c | | * ctrl_getatemp
34 c | | * ctrl_getaqh
35 c | | * ctrl_getuwind
36 c | | * ctrl_getvwind
37 c | | 4. flux fields are interpolated on current time
38 c | | (b) next the flux fields are computed via bulk formulae
39 c | |
40 c | |-- exf_obcs
41 c | If open boundaries are prescribed externally,
42 c | OB values for current time step are read here
43 c | 1. for each boundary (north, south, west, east)
44 c | sliced fields are read for T,S,U,V
45 c | 2. if OB's are part of control vector,
46 c | control variable anomalies are added to OB values
47 c | * ctrl_getobcsn
48 c | * ctrl_getobcss
49 c | * ctrl_getobcsw
50 c | * ctrl_getobcse
51 c |
52 c |-- exf_getsurfacefluxes
53 c | If forcing and control is flux, then the
54 c | control vector anomalies are read here
55 c | * ctrl_getheatflux
56 c | * ctrl_getsaltflux
57 c | * ctrl_getzonstress
58 c | * call ctrl_getmerstress
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 mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
88 c *=================================================================
89 c | SUBROUTINE exf_getforcing
90 c *=================================================================
91 C \ev
92
93 C !USES:
94 implicit none
95
96 c == global variables ==
97
98 C !INPUT/OUTPUT PARAMETERS:
99 c == routine arguments ==
100 integer mythid
101 integer myiter
102 _RL mytime
103
104 C !LOCAL VARIABLES:
105 c == local variables ==
106
107 c == end of interface ==
108 CEOP
109
110 c Get values of climatological fields.
111 call exf_getclim( mytime, myiter, mythid )
112
113 c Get the surface forcing fields.
114 call exf_getffields( mytime, myiter, mythid )
115
116 c Get values of the surface flux anomalies.
117 call exf_getsurfacefluxes( mytime, myiter, mythid )
118
119 c Map the forcing fields onto the corresponding model fields.
120 call exf_mapfields( mythid )
121
122 end
123
  ViewVC Help
Powered by ViewVC 1.1.22