pkg/exf/exf_getffields.F

c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_getffields.F,v 1.2.6.12 2003/06/19 15:31:51 heimbach Exp $

#include "EXF_CPPOPTIONS.h"

      subroutine exf_getffields( mycurrenttime, mycurrentiter, mythid )

c     ==================================================================
c     SUBROUTINE exf_getffields
c     ==================================================================
c
c     o Read-in atmospheric state and/or surface fluxes from files.
c
c       heimbach@mit.edu, 23-May-2003 totally re-structured
c
c     ==================================================================
c     SUBROUTINE exf_getffields
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "GRID.h"

#include "exf_param.h"
#include "exf_fields.h"
#include "exf_constants.h"

#if (defined (ALLOW_ADJOINT_RUN) || \
     defined (ALLOW_TANGENTLINEAR_RUN) || \
     defined (ALLOW_ECCO_OPTIMIZATION))
# include "ctrl.h"
# include "ctrl_dummy.h"
#endif

c     == routine arguments ==

      integer mythid
      integer mycurrentiter
      _RL     mycurrenttime

c     == local variables ==

c     == end of interface ==

c--   read forcing fields from files and temporal interpolation

#ifdef ALLOW_ATM_WIND

c     Zonal wind.
      call exf_set_gen( 
     &     uwindfile, uwindstartdate, uwindperiod, 
     &     exf_inscal_uwind,
     &     uwind, uwind0, uwind1, uwindmask, 
     &     mycurrenttime, mycurrentiter, mythid )

c     Meridional wind.
      call exf_set_gen( 
     &     vwindfile, vwindstartdate, vwindperiod, 
     &     exf_inscal_vwind,
     &     vwind, vwind0, vwind1, vwindmask, 
     &     mycurrenttime, mycurrentiter, mythid )

#ifdef ALLOW_UWIND_CONTROL
      call ctrl_get_gen ( 
     &     xx_uwind_file, xx_uwindstartdate, xx_uwindperiod,
     &     maskw, uwind, xx_uwind0, xx_uwind1, xx_uwind_dummy,
     &     mycurrenttime, mycurrentiter, mythid )
#endif

#ifdef ALLOW_VWIND_CONTROL
      call ctrl_get_gen ( 
     &     xx_vwind_file, xx_vwindstartdate, xx_vwindperiod,
     &     masks, vwind, xx_vwind0, xx_vwind1, xx_vwind_dummy,
     &     mycurrenttime, mycurrentiter, mythid )
#endif

#else  /* ifndef ALLOW_ATM_WIND */

c     Zonal wind stress.
      call exf_set_gen( 
     &     ustressfile, ustressstartdate, ustressperiod, 
     &     exf_inscal_ustress,
     &     ustress, ustress0, ustress1, ustressmask, 
     &     mycurrenttime, mycurrentiter, mythid )

c     Meridional wind stress.
      call exf_set_gen( 
     &     vstressfile, vstressstartdate, vstressperiod, 
     &     exf_inscal_vstress,
     &     vstress, vstress0, vstress1, vstressmask, 
     &     mycurrenttime, mycurrentiter, mythid )

#endif /* ifndef ALLOW_ATM_WIND */

#ifdef ALLOW_ATM_TEMP

c     Atmospheric temperature.
      call exf_set_gen( 
     &     atempfile, atempstartdate, atempperiod, 
     &     exf_inscal_atemp,
     &     atemp, atemp0, atemp1, atempmask, 
     &     mycurrenttime, mycurrentiter, mythid )

c     Atmospheric humidity.
      call exf_set_gen( 
     &     aqhfile, aqhstartdate, aqhperiod, 
     &     exf_inscal_aqh,
     &     aqh, aqh0, aqh1, aqhmask, 
     &     mycurrenttime, mycurrentiter, mythid )

c     Net long wave radiative flux.
      call exf_set_gen( 
     &     lwfluxfile, lwfluxstartdate, lwfluxperiod, 
     &     exf_inscal_lwflux,
     &     lwflux, lwflux0, lwflux1, lwfluxmask, 
     &     mycurrenttime, mycurrentiter, mythid )

c     Precipitation.
      call exf_set_gen( 
     &     precipfile, precipstartdate, precipperiod, 
     &     exf_inscal_precip,
     &     precip, precip0, precip1, precipmask, 
     &     mycurrenttime, mycurrentiter, mythid )

#ifdef ALLOW_ATEMP_CONTROL
      call ctrl_get_gen ( 
     &     xx_atemp_file, xx_atempstartdate, xx_atempperiod,
     &     maskc, atemp, xx_atemp0, xx_atemp1, xx_atemp_dummy,
     &     mycurrenttime, mycurrentiter, mythid )
#endif

#ifdef ALLOW_AQH_CONTROL
      call ctrl_get_gen ( 
     &     xx_aqh_file, xx_aqhstartdate, xx_aqhperiod,
     &     maskc, aqh, xx_aqh0, xx_aqh1, xx_aqh_dummy,
     &     mycurrenttime, mycurrentiter, mythid )
#endif

#else  /* ifndef ALLOW_ATM_TEMP */

c     Atmospheric heat flux.
      call exf_set_gen  ( 
     &     hfluxfile, hfluxstartdate, hfluxperiod, exf_inscal_hflux,
     &     hflux, hflux0, hflux1, hfluxmask, 
     &     mycurrenttime, mycurrentiter, mythid )

c     Salt flux.
      call exf_set_gen  ( 
     &     sfluxfile, sfluxstartdate, sfluxperiod, exf_inscal_sflux,
     &     sflux, sflux0, sflux1, sfluxmask, 
     &     mycurrenttime, mycurrentiter, mythid )

#endif /* ifndef ALLOW_ATM_TEMP */

#if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
c     Net short wave radiative flux.
      call exf_set_gen  ( 
     &     swfluxfile, swfluxstartdate, swfluxperiod, exf_inscal_swflux,
     &     swflux, swflux0, swflux1, swfluxmask, 
     &     mycurrenttime, mycurrentiter, mythid )
#endif

#ifdef EXF_READ_EVAP
c     Evaporation
      call exf_set_gen  ( 
     &     evapfile, evapstartdate, evapperiod, exf_inscal_evap,
     &     evap, evap0, evap1, evapmask, 
     &     mycurrenttime, mycurrentiter, mythid )
#endif

#ifdef ALLOW_DOWNWARD_RADIATION

c     Downward shortwave radiation.
      call exf_set_gen  ( 
     &     swdownfile, swdownstartdate, swdownperiod, exf_inscal_swdown,
     &     swdown, swdown0, swdown1, swdownmask, 
     &     mycurrenttime, mycurrentiter, mythid )

c     Downward longwave radiation.
      call exf_set_gen  ( 
     &     lwdownfile, lwdownstartdate, lwdownperiod, exf_inscal_lwdown,
     &     lwdown, lwdown0, lwdown1, lwdownmask, 
     &     mycurrenttime, mycurrentiter, mythid )

#endif

#ifdef ATMOSPHERIC_LOADING
c     Atmos. pressure forcing
      call exf_set_gen  ( 
     &     apressurefile, apressurestartdate, apressureperiod, 
     &     exf_inscal_apressure,
     &     apressure, apressure0, apressure1, apressuremask, 
     &     mycurrenttime, mycurrentiter, mythid )
#endif


      end
1	heimbach	1.14	c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_getffields.F,v 1.2.6.12 2003/06/19 15:31:51 heimbach Exp $
2	heimbach	1.1
3			#include "EXF_CPPOPTIONS.h"
4
5	heimbach	1.13	subroutine exf_getffields( mycurrenttime, mycurrentiter, mythid )
6	heimbach	1.1
7			c ==================================================================
8	heimbach	1.13	c SUBROUTINE exf_getffields
9	heimbach	1.1	c ==================================================================
10			c
11	dimitri	1.8	c o Read-in atmospheric state and/or surface fluxes from files.
12			c
13	heimbach	1.13	c heimbach@mit.edu, 23-May-2003 totally re-structured
14	dimitri	1.7	c
15	heimbach	1.1	c ==================================================================
16	heimbach	1.13	c SUBROUTINE exf_getffields
17	heimbach	1.1	c ==================================================================
18
19			implicit none
20
21			c == global variables ==
22
23			#include "EEPARAMS.h"
24			#include "SIZE.h"
25			#include "PARAMS.h"
26			#include "DYNVARS.h"
27			#include "GRID.h"
28
29	dimitri	1.9	#include "exf_param.h"
30	heimbach	1.1	#include "exf_fields.h"
31			#include "exf_constants.h"
32
33	heimbach	1.14	#if (defined (ALLOW_ADJOINT_RUN) \|\| \
34			defined (ALLOW_TANGENTLINEAR_RUN) \|\| \
35			defined (ALLOW_ECCO_OPTIMIZATION))
36			# include "ctrl.h"
37			# include "ctrl_dummy.h"
38			#endif
39
40	heimbach	1.1	c == routine arguments ==
41
42			integer mythid
43			integer mycurrentiter
44			_RL mycurrenttime
45
46			c == local variables ==
47
48			c == end of interface ==
49
50	dimitri	1.8	c-- read forcing fields from files and temporal interpolation
51
52			#ifdef ALLOW_ATM_WIND
53
54			c Zonal wind.
55	heimbach	1.13	call exf_set_gen(
56			& uwindfile, uwindstartdate, uwindperiod,
57			& exf_inscal_uwind,
58			& uwind, uwind0, uwind1, uwindmask,
59			& mycurrenttime, mycurrentiter, mythid )
60	dimitri	1.8
61			c Meridional wind.
62	heimbach	1.13	call exf_set_gen(
63			& vwindfile, vwindstartdate, vwindperiod,
64			& exf_inscal_vwind,
65			& vwind, vwind0, vwind1, vwindmask,
66			& mycurrenttime, mycurrentiter, mythid )
67	dimitri	1.8
68			#ifdef ALLOW_UWIND_CONTROL
69	heimbach	1.14	call ctrl_get_gen (
70			& xx_uwind_file, xx_uwindstartdate, xx_uwindperiod,
71			& maskw, uwind, xx_uwind0, xx_uwind1, xx_uwind_dummy,
72			& mycurrenttime, mycurrentiter, mythid )
73	dimitri	1.8	#endif
74
75			#ifdef ALLOW_VWIND_CONTROL
76	heimbach	1.14	call ctrl_get_gen (
77			& xx_vwind_file, xx_vwindstartdate, xx_vwindperiod,
78			& masks, vwind, xx_vwind0, xx_vwind1, xx_vwind_dummy,
79			& mycurrenttime, mycurrentiter, mythid )
80	dimitri	1.8	#endif
81
82			#else /* ifndef ALLOW_ATM_WIND */
83
84			c Zonal wind stress.
85	heimbach	1.13	call exf_set_gen(
86			& ustressfile, ustressstartdate, ustressperiod,
87			& exf_inscal_ustress,
88			& ustress, ustress0, ustress1, ustressmask,
89			& mycurrenttime, mycurrentiter, mythid )
90	dimitri	1.8
91			c Meridional wind stress.
92	heimbach	1.13	call exf_set_gen(
93			& vstressfile, vstressstartdate, vstressperiod,
94			& exf_inscal_vstress,
95			& vstress, vstress0, vstress1, vstressmask,
96			& mycurrenttime, mycurrentiter, mythid )
97	dimitri	1.8
98			#endif /* ifndef ALLOW_ATM_WIND */
99
100	heimbach	1.1	#ifdef ALLOW_ATM_TEMP
101	dimitri	1.8
102	heimbach	1.1	c Atmospheric temperature.
103	heimbach	1.13	call exf_set_gen(
104			& atempfile, atempstartdate, atempperiod,
105			& exf_inscal_atemp,
106			& atemp, atemp0, atemp1, atempmask,
107			& mycurrenttime, mycurrentiter, mythid )
108	heimbach	1.1
109			c Atmospheric humidity.
110	heimbach	1.13	call exf_set_gen(
111			& aqhfile, aqhstartdate, aqhperiod,
112			& exf_inscal_aqh,
113			& aqh, aqh0, aqh1, aqhmask,
114			& mycurrenttime, mycurrentiter, mythid )
115	heimbach	1.1
116			c Net long wave radiative flux.
117	heimbach	1.13	call exf_set_gen(
118			& lwfluxfile, lwfluxstartdate, lwfluxperiod,
119			& exf_inscal_lwflux,
120			& lwflux, lwflux0, lwflux1, lwfluxmask,
121			& mycurrenttime, mycurrentiter, mythid )
122	heimbach	1.1
123			c Precipitation.
124	heimbach	1.13	call exf_set_gen(
125			& precipfile, precipstartdate, precipperiod,
126			& exf_inscal_precip,
127			& precip, precip0, precip1, precipmask,
128			& mycurrenttime, mycurrentiter, mythid )
129	heimbach	1.1
130	heimbach	1.3	#ifdef ALLOW_ATEMP_CONTROL
131	heimbach	1.14	call ctrl_get_gen (
132			& xx_atemp_file, xx_atempstartdate, xx_atempperiod,
133			& maskc, atemp, xx_atemp0, xx_atemp1, xx_atemp_dummy,
134			& mycurrenttime, mycurrentiter, mythid )
135	heimbach	1.3	#endif
136
137			#ifdef ALLOW_AQH_CONTROL
138	heimbach	1.14	call ctrl_get_gen (
139			& xx_aqh_file, xx_aqhstartdate, xx_aqhperiod,
140			& maskc, aqh, xx_aqh0, xx_aqh1, xx_aqh_dummy,
141			& mycurrenttime, mycurrentiter, mythid )
142	heimbach	1.3	#endif
143	heimbach	1.1
144	dimitri	1.8	#else /* ifndef ALLOW_ATM_TEMP */
145	heimbach	1.1
146			c Atmospheric heat flux.
147	heimbach	1.13	call exf_set_gen (
148			& hfluxfile, hfluxstartdate, hfluxperiod, exf_inscal_hflux,
149			& hflux, hflux0, hflux1, hfluxmask,
150			& mycurrenttime, mycurrentiter, mythid )
151	heimbach	1.1
152			c Salt flux.
153	heimbach	1.13	call exf_set_gen (
154			& sfluxfile, sfluxstartdate, sfluxperiod, exf_inscal_sflux,
155			& sflux, sflux0, sflux1, sfluxmask,
156			& mycurrenttime, mycurrentiter, mythid )
157	dimitri	1.8
158			#endif /* ifndef ALLOW_ATM_TEMP */
159	heimbach	1.1
160	dimitri	1.8	#if defined(ALLOW_ATM_TEMP) \|\| defined(SHORTWAVE_HEATING)
161	heimbach	1.1	c Net short wave radiative flux.
162	heimbach	1.13	call exf_set_gen (
163			& swfluxfile, swfluxstartdate, swfluxperiod, exf_inscal_swflux,
164			& swflux, swflux0, swflux1, swfluxmask,
165			& mycurrenttime, mycurrentiter, mythid )
166	heimbach	1.3	#endif
167
168	dimitri	1.8	#ifdef EXF_READ_EVAP
169			c Evaporation
170	heimbach	1.13	call exf_set_gen (
171			& evapfile, evapstartdate, evapperiod, exf_inscal_evap,
172			& evap, evap0, evap1, evapmask,
173			& mycurrenttime, mycurrentiter, mythid )
174	heimbach	1.3	#endif
175	heimbach	1.1
176	dimitri	1.8	#ifdef ALLOW_DOWNWARD_RADIATION
177	heimbach	1.1
178	dimitri	1.8	c Downward shortwave radiation.
179	heimbach	1.13	call exf_set_gen (
180			& swdownfile, swdownstartdate, swdownperiod, exf_inscal_swdown,
181			& swdown, swdown0, swdown1, swdownmask,
182			& mycurrenttime, mycurrentiter, mythid )
183	heimbach	1.1
184	dimitri	1.8	c Downward longwave radiation.
185	heimbach	1.13	call exf_set_gen (
186			& lwdownfile, lwdownstartdate, lwdownperiod, exf_inscal_lwdown,
187			& lwdown, lwdown0, lwdown1, lwdownmask,
188			& mycurrenttime, mycurrentiter, mythid )
189	heimbach	1.1
190	heimbach	1.3	#endif
191	heimbach	1.1
192	heimbach	1.12	#ifdef ATMOSPHERIC_LOADING
193			c Atmos. pressure forcing
194	heimbach	1.13	call exf_set_gen (
195			& apressurefile, apressurestartdate, apressureperiod,
196			& exf_inscal_apressure,
197			& apressure, apressure0, apressure1, apressuremask,
198			& mycurrenttime, mycurrentiter, mythid )
199	heimbach	1.3	#endif
200
201	heimbach	1.1
202			end