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

Diff of /MITgcm/pkg/exf/exf_getffields.F

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

revision 1.14 by heimbach, Tue Jun 24 16:07:32 2003 UTC revision 1.33 by heimbach, Fri Oct 27 05:16:54 2006 UTC
# Line 1  Line 1 
1  c $Header$  C
2    C $Header$
3    C $Name$
4    
5  #include "EXF_CPPOPTIONS.h"  #include "EXF_OPTIONS.h"
6    
7        subroutine exf_getffields( mycurrenttime, mycurrentiter, mythid )        subroutine exf_getffields( mytime, myiter, mythid )
8    
9  c     ==================================================================  c     ==================================================================
10  c     SUBROUTINE exf_getffields  c     SUBROUTINE exf_getffields
# Line 11  c Line 13  c
13  c     o Read-in atmospheric state and/or surface fluxes from files.  c     o Read-in atmospheric state and/or surface fluxes from files.
14  c  c
15  c       heimbach@mit.edu, 23-May-2003 totally re-structured  c       heimbach@mit.edu, 23-May-2003 totally re-structured
16    c       5-Aug-2003: added USE_EXF_INTERPOLATION for arbitrary input grid
17  c  c
18  c     ==================================================================  c     ==================================================================
19  c     SUBROUTINE exf_getffields  c     SUBROUTINE exf_getffields
# Line 30  c     == global variables == Line 33  c     == global variables ==
33  #include "exf_fields.h"  #include "exf_fields.h"
34  #include "exf_constants.h"  #include "exf_constants.h"
35    
36  #if (defined (ALLOW_ADJOINT_RUN) || \  #ifdef ALLOW_AUTODIFF
      defined (ALLOW_TANGENTLINEAR_RUN) || \  
      defined (ALLOW_ECCO_OPTIMIZATION))  
37  # include "ctrl.h"  # include "ctrl.h"
38  # include "ctrl_dummy.h"  # include "ctrl_dummy.h"
39  #endif  #endif
# Line 40  c     == global variables == Line 41  c     == global variables ==
41  c     == routine arguments ==  c     == routine arguments ==
42    
43        integer mythid        integer mythid
44        integer mycurrentiter        integer myiter
45        _RL     mycurrenttime        _RL     mytime
46    
47  c     == local variables ==  c     == local variables ==
48    
49          integer i, j, bi, bj, interp_method
50          parameter(interp_method=1)
51    
52  c     == end of interface ==  c     == end of interface ==
53    
54  c--   read forcing fields from files and temporal interpolation  c--   read forcing fields from files and temporal interpolation
55    
56  #ifdef ALLOW_ATM_WIND  c     Zonal and meridional wind stress.
57    #ifdef USE_EXF_INTERPOLATION
58  c     Zonal wind.        call exf_set_uv(
59         &     ustressfile, ustressstartdate, ustressperiod,
60         &     ustressstartdate1, ustressstartdate2,
61         &     exf_inscal_ustress, ustress, ustress0, ustress1, ustressmask,
62         &     ustress_lon0, ustress_lon_inc, ustress_lat0, ustress_lat_inc,
63         &     ustress_nlon, ustress_nlat,
64         &     ustress_exfremo_intercept, ustress_exfremo_slope,
65         &     vstressfile, vstressstartdate, vstressperiod,
66         &     vstressstartdate1, vstressstartdate2,
67         &     exf_inscal_vstress, vstress, vstress0, vstress1, vstressmask,
68         &     vstress_lon0, vstress_lon_inc, vstress_lat0, vstress_lat_inc,
69         &     vstress_nlon, vstress_nlat,
70         &     vstress_exfremo_intercept, vstress_exfremo_slope,
71         &     mytime, myiter, mythid )
72    #else /* ifndef USE_EXF_INTERPOLATION */
73        call exf_set_gen(        call exf_set_gen(
74       &     uwindfile, uwindstartdate, uwindperiod,       &     ustressfile, ustressstartdate, ustressperiod,
75       &     exf_inscal_uwind,       &     ustressstartdate1, ustressstartdate2,
76       &     uwind, uwind0, uwind1, uwindmask,       &     exf_inscal_ustress,
77       &     mycurrenttime, mycurrentiter, mythid )       &     ustress_exfremo_intercept, ustress_exfremo_slope,
78         &     ustress, ustress0, ustress1, ustressmask,
79         &     mytime, myiter, mythid )
80          call exf_set_gen(
81         &     vstressfile, vstressstartdate, vstressperiod,
82         &     ustressstartdate1, ustressstartdate2,
83         &     exf_inscal_vstress,
84         &     vstress_exfremo_intercept, vstress_exfremo_slope,
85         &     vstress, vstress0, vstress1, vstressmask,
86         &     mytime, myiter, mythid )
87    #endif /* USE_EXF_INTERPOLATION */
88    
89  c     Meridional wind.  c--   wind speed
90        call exf_set_gen(        call exf_set_gen(
91       &     vwindfile, vwindstartdate, vwindperiod,       &     wspeedfile, wspeedstartdate, wspeedperiod,
92       &     exf_inscal_vwind,       &     wspeedstartdate1, wspeedstartdate2,
93       &     vwind, vwind0, vwind1, vwindmask,       &     exf_inscal_wspeed,
94       &     mycurrenttime, mycurrentiter, mythid )       &     wspeed_exfremo_intercept, wspeed_exfremo_slope,
95         &     wspeed, wspeed0, wspeed1, wspeedmask,
96  #ifdef ALLOW_UWIND_CONTROL  #ifdef USE_EXF_INTERPOLATION
97        call ctrl_get_gen (       &     wspeed_lon0, wspeed_lon_inc,
98       &     xx_uwind_file, xx_uwindstartdate, xx_uwindperiod,       &     wspeed_lat0, wspeed_lat_inc,
99       &     maskw, uwind, xx_uwind0, xx_uwind1, xx_uwind_dummy,       &     wspeed_nlon, wspeed_nlat, xC, yC, interp_method,
      &     mycurrenttime, mycurrentiter, mythid )  
100  #endif  #endif
101         &     mytime, myiter, mythid )
102    
103  #ifdef ALLOW_VWIND_CONTROL  #ifdef ALLOW_ATM_WIND
       call ctrl_get_gen (  
      &     xx_vwind_file, xx_vwindstartdate, xx_vwindperiod,  
      &     masks, vwind, xx_vwind0, xx_vwind1, xx_vwind_dummy,  
      &     mycurrenttime, mycurrentiter, mythid )  
 #endif  
104    
105  #else  /* ifndef ALLOW_ATM_WIND */  c     Zonal and meridional wind.
106    #ifdef USE_EXF_INTERPOLATION
107          call exf_set_uv(
108         &     uwindfile, uwindstartdate, uwindperiod,
109         &     uwindstartdate1, uwindstartdate2,
110         &     exf_inscal_uwind, uwind, uwind0, uwind1, uwindmask,
111         &     uwind_lon0, uwind_lon_inc, uwind_lat0, uwind_lat_inc,
112         &     uwind_nlon, uwind_nlat,
113         &     uwind_exfremo_intercept, uwind_exfremo_slope,
114         &     vwindfile, vwindstartdate, vwindperiod,
115         &     vwindstartdate1, vwindstartdate2,
116         &     exf_inscal_vwind, vwind, vwind0, vwind1, vwindmask,
117         &     vwind_lon0, vwind_lon_inc, vwind_lat0, vwind_lat_inc,
118         &     vwind_nlon, vwind_nlat,
119         &     vwind_exfremo_intercept, vwind_exfremo_slope,
120         &     mytime, myiter, mythid )
121    #else /* ifndef USE_EXF_INTERPOLATION */
122          call exf_set_gen(
123         &     uwindfile, uwindstartdate, uwindperiod,
124         &     uwindstartdate1, uwindstartdate2,
125         &     exf_inscal_uwind,
126         &     uwind_exfremo_intercept, uwind_exfremo_slope,
127         &     uwind, uwind0, uwind1, uwindmask,
128         &     mytime, myiter, mythid )
129          call exf_set_gen(
130         &     vwindfile, vwindstartdate, vwindperiod,
131         &     vwindstartdate1, vwindstartdate2,
132         &     exf_inscal_vwind,
133         &     vwind_exfremo_intercept, vwind_exfremo_slope,
134         &     vwind, vwind0, vwind1, vwindmask,
135         &     mytime, myiter, mythid )
136    #endif /* USE_EXF_INTERPOLATION */
137    
138  c     Zonal wind stress.  #endif /* ALLOW_ATM_WIND */
       call exf_set_gen(  
      &     ustressfile, ustressstartdate, ustressperiod,  
      &     exf_inscal_ustress,  
      &     ustress, ustress0, ustress1, ustressmask,  
      &     mycurrenttime, mycurrentiter, mythid )  
139    
140  c     Meridional wind stress.  c     Atmospheric heat flux.
141        call exf_set_gen(        call exf_set_gen  (
142       &     vstressfile, vstressstartdate, vstressperiod,       &     hfluxfile, hfluxstartdate, hfluxperiod,
143       &     exf_inscal_vstress,       &     hfluxstartdate1, hfluxstartdate2,
144       &     vstress, vstress0, vstress1, vstressmask,       &     exf_inscal_hflux,
145       &     mycurrenttime, mycurrentiter, mythid )       &     hflux_exfremo_intercept, hflux_exfremo_slope,
146         &     hflux, hflux0, hflux1, hfluxmask,
147    #ifdef USE_EXF_INTERPOLATION
148         &     hflux_lon0, hflux_lon_inc, hflux_lat0, hflux_lat_inc,
149         &     hflux_nlon, hflux_nlat, xC, yC, interp_method,
150    #endif
151         &     mytime, myiter, mythid )
152    
153  #endif /* ifndef ALLOW_ATM_WIND */  c     Salt flux.
154          call exf_set_gen  (
155         &     sfluxfile, sfluxstartdate, sfluxperiod,
156         &     sfluxstartdate1, sfluxstartdate2,
157         &     exf_inscal_sflux,
158         &     sflux_exfremo_intercept, sflux_exfremo_slope,
159         &     sflux, sflux0, sflux1, sfluxmask,
160    #ifdef USE_EXF_INTERPOLATION
161         &     sflux_lon0, sflux_lon_inc, sflux_lat0, sflux_lat_inc,
162         &     sflux_nlon, sflux_nlat, xC, yC, interp_method,
163    #endif
164         &     mytime, myiter, mythid )
165    
166  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
167    
168  c     Atmospheric temperature.  c     Atmospheric temperature.
169        call exf_set_gen(        call exf_set_gen(
170       &     atempfile, atempstartdate, atempperiod,       &     atempfile, atempstartdate, atempperiod,
171         &     atempstartdate1, atempstartdate2,
172       &     exf_inscal_atemp,       &     exf_inscal_atemp,
173         &     atemp_exfremo_intercept, atemp_exfremo_slope,
174       &     atemp, atemp0, atemp1, atempmask,       &     atemp, atemp0, atemp1, atempmask,
175       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
176         &     atemp_lon0, atemp_lon_inc, atemp_lat0, atemp_lat_inc,
177         &     atemp_nlon, atemp_nlat, xC, yC, interp_method,
178    #endif
179         &     mytime, myiter, mythid )
180          do bj = mybylo(mythid),mybyhi(mythid)
181           do bi = mybxlo(mythid),mybxhi(mythid)
182            do j = 1,sny
183             do i = 1,snx
184              atemp(i,j,bi,bj) = atemp(i,j,bi,bj) + exf_offset_atemp
185             enddo
186            enddo
187           enddo
188          enddo
189    
190  c     Atmospheric humidity.  c     Atmospheric humidity.
191        call exf_set_gen(        call exf_set_gen(
192       &     aqhfile, aqhstartdate, aqhperiod,       &     aqhfile, aqhstartdate, aqhperiod,
193         &     aqhstartdate1, aqhstartdate2,
194       &     exf_inscal_aqh,       &     exf_inscal_aqh,
195         &     aqh_exfremo_intercept, aqh_exfremo_slope,
196       &     aqh, aqh0, aqh1, aqhmask,       &     aqh, aqh0, aqh1, aqhmask,
197       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
198         &     aqh_lon0, aqh_lon_inc, aqh_lat0, aqh_lat_inc,
199         &     aqh_nlon, aqh_nlat, xC, yC, interp_method,
200    #endif
201         &     mytime, myiter, mythid )
202    
203  c     Net long wave radiative flux.  c     Net long wave radiative flux.
204        call exf_set_gen(        call exf_set_gen(
205       &     lwfluxfile, lwfluxstartdate, lwfluxperiod,       &     lwfluxfile, lwfluxstartdate, lwfluxperiod,
206         &     lwfluxstartdate1, lwfluxstartdate2,
207       &     exf_inscal_lwflux,       &     exf_inscal_lwflux,
208         &     lwflux_exfremo_intercept, lwflux_exfremo_slope,
209       &     lwflux, lwflux0, lwflux1, lwfluxmask,       &     lwflux, lwflux0, lwflux1, lwfluxmask,
210       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
211         &     lwflux_lon0, lwflux_lon_inc, lwflux_lat0, lwflux_lat_inc,
212         &     lwflux_nlon, lwflux_nlat, xC, yC, interp_method,
213    #endif
214         &     mytime, myiter, mythid )
215    
216  c     Precipitation.  c     Precipitation.
217        call exf_set_gen(        call exf_set_gen(
218       &     precipfile, precipstartdate, precipperiod,       &     precipfile, precipstartdate, precipperiod,
219         &     precipstartdate1, precipstartdate2,
220       &     exf_inscal_precip,       &     exf_inscal_precip,
221         &     precip_exfremo_intercept, precip_exfremo_slope,
222       &     precip, precip0, precip1, precipmask,       &     precip, precip0, precip1, precipmask,
223       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
224         &     precip_lon0, precip_lon_inc, precip_lat0, precip_lat_inc,
225  #ifdef ALLOW_ATEMP_CONTROL       &     precip_nlon, precip_nlat, xC, yC, interp_method,
       call ctrl_get_gen (  
      &     xx_atemp_file, xx_atempstartdate, xx_atempperiod,  
      &     maskc, atemp, xx_atemp0, xx_atemp1, xx_atemp_dummy,  
      &     mycurrenttime, mycurrentiter, mythid )  
226  #endif  #endif
227         &     mytime, myiter, mythid )
228    
229  #ifdef ALLOW_AQH_CONTROL  c     Snow.
230        call ctrl_get_gen (        call exf_set_gen(
231       &     xx_aqh_file, xx_aqhstartdate, xx_aqhperiod,       &     snowprecipfile, snowprecipstartdate, snowprecipperiod,
232       &     maskc, aqh, xx_aqh0, xx_aqh1, xx_aqh_dummy,       &     snowprecipstartdate1, snowprecipstartdate2,
233       &     mycurrenttime, mycurrentiter, mythid )       &     exf_inscal_snowprecip,
234         &     snowprecip_exfremo_intercept, snowprecip_exfremo_slope,
235         &     snowprecip, snowprecip0, snowprecip1, snowprecipmask,
236    #ifdef USE_EXF_INTERPOLATION
237         &     snowprecip_lon0, snowprecip_lon_inc,
238         &     snowprecip_lat0, snowprecip_lat_inc,
239         &     snowprecip_nlon, snowprecip_nlat, xC, yC, interp_method,
240  #endif  #endif
241         &     mytime, myiter, mythid )
242    
243  #else  /* ifndef ALLOW_ATM_TEMP */  #endif /* 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 */  
244    
245  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)  #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
246  c     Net short wave radiative flux.  c     Net short wave radiative flux.
247        call exf_set_gen  (        call exf_set_gen  (
248       &     swfluxfile, swfluxstartdate, swfluxperiod, exf_inscal_swflux,       &     swfluxfile, swfluxstartdate, swfluxperiod,
249         &     swfluxstartdate1, swfluxstartdate2,
250         &     exf_inscal_swflux,
251         &     swflux_exfremo_intercept, swflux_exfremo_slope,
252       &     swflux, swflux0, swflux1, swfluxmask,       &     swflux, swflux0, swflux1, swfluxmask,
253       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
254         &     swflux_lon0, swflux_lon_inc, swflux_lat0, swflux_lat_inc,
255         &     swflux_nlon, swflux_nlat, xC, yC, interp_method,
256    #endif
257         &     mytime, myiter, mythid )
258  #endif  #endif
259    
260  #ifdef EXF_READ_EVAP  #ifdef EXF_READ_EVAP
261  c     Evaporation  c     Evaporation
262        call exf_set_gen  (        call exf_set_gen  (
263       &     evapfile, evapstartdate, evapperiod, exf_inscal_evap,       &     evapfile, evapstartdate, evapperiod,
264         &     evapstartdate1, evapstartdate2,
265         &     exf_inscal_evap,
266         &     evap_exfremo_intercept, evap_exfremo_slope,
267       &     evap, evap0, evap1, evapmask,       &     evap, evap0, evap1, evapmask,
268       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
269         &     evap_lon0, evap_lon_inc, evap_lat0, evap_lat_inc,
270         &     evap_nlon, evap_nlat, xC, yC, interp_method,
271    #endif
272         &     mytime, myiter, mythid )
273  #endif  #endif
274    
275  #ifdef ALLOW_DOWNWARD_RADIATION  #ifdef ALLOW_DOWNWARD_RADIATION
276    
277  c     Downward shortwave radiation.  c     Downward shortwave radiation.
278        call exf_set_gen  (        call exf_set_gen  (
279       &     swdownfile, swdownstartdate, swdownperiod, exf_inscal_swdown,       &     swdownfile, swdownstartdate, swdownperiod,
280         &     swdownstartdate1, swdownstartdate2,
281         &     exf_inscal_swdown,
282         &     swdown_exfremo_intercept, swdown_exfremo_slope,
283       &     swdown, swdown0, swdown1, swdownmask,       &     swdown, swdown0, swdown1, swdownmask,
284       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
285         &     swdown_lon0, swdown_lon_inc, swdown_lat0, swdown_lat_inc,
286         &     swdown_nlon, swdown_nlat, xC, yC, interp_method,
287    #endif
288         &     mytime, myiter, mythid )
289    
290  c     Downward longwave radiation.  c     Downward longwave radiation.
291        call exf_set_gen  (        call exf_set_gen  (
292       &     lwdownfile, lwdownstartdate, lwdownperiod, exf_inscal_lwdown,       &     lwdownfile, lwdownstartdate, lwdownperiod,
293         &     lwdownstartdate1, lwdownstartdate2,
294         &     exf_inscal_lwdown,
295         &     lwdown_exfremo_intercept, lwdown_exfremo_slope,
296       &     lwdown, lwdown0, lwdown1, lwdownmask,       &     lwdown, lwdown0, lwdown1, lwdownmask,
297       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
298         &     lwdown_lon0, lwdown_lon_inc, lwdown_lat0, lwdown_lat_inc,
299         &     lwdown_nlon, lwdown_nlat, xC, yC, interp_method,
300    #endif
301         &     mytime, myiter, mythid )
302    
303  #endif  #endif
304    
# Line 193  c     Downward longwave radiation. Line 306  c     Downward longwave radiation.
306  c     Atmos. pressure forcing  c     Atmos. pressure forcing
307        call exf_set_gen  (        call exf_set_gen  (
308       &     apressurefile, apressurestartdate, apressureperiod,       &     apressurefile, apressurestartdate, apressureperiod,
309         &     apressurestartdate1, apressurestartdate2,
310       &     exf_inscal_apressure,       &     exf_inscal_apressure,
311         &     apressure_exfremo_intercept, apressure_exfremo_slope,
312       &     apressure, apressure0, apressure1, apressuremask,       &     apressure, apressure0, apressure1, apressuremask,
313       &     mycurrenttime, mycurrentiter, mythid )  #ifdef USE_EXF_INTERPOLATION
314         &     apressure_lon0, apressure_lon_inc, apressure_lat0,
315         &     apressure_lat_inc, apressure_nlon, apressure_nlat, xC, yC,
316         &     interp_method,
317    #endif
318         &     mytime, myiter, mythid )
319  #endif  #endif
320    
321    c-- Control variables for atmos. state
322    
323    #ifdef ALLOW_ATEMP_CONTROL
324          call ctrl_get_gen (
325         &     xx_atemp_file, xx_atempstartdate, xx_atempperiod,
326         &     maskc, atemp, xx_atemp0, xx_atemp1, xx_atemp_dummy,
327         &     xx_atemp_remo_intercept, xx_atemp_remo_slope,
328         &     mytime, myiter, mythid )
329    #endif
330    
331    #ifdef ALLOW_AQH_CONTROL
332          call ctrl_get_gen (
333         &     xx_aqh_file, xx_aqhstartdate, xx_aqhperiod,
334         &     maskc, aqh, xx_aqh0, xx_aqh1, xx_aqh_dummy,
335         &     xx_aqh_remo_intercept, xx_aqh_remo_slope,
336         &     mytime, myiter, mythid )
337    #endif
338    
339    #ifdef ALLOW_PRECIP_CONTROL
340          call ctrl_get_gen (
341         &     xx_precip_file, xx_precipstartdate, xx_precipperiod,
342         &     maskc, precip, xx_precip0, xx_precip1, xx_precip_dummy,
343         &     xx_precip_remo_intercept, xx_precip_remo_slope,
344         &     mytime, myiter, mythid )
345    #endif
346    
347    #ifdef ALLOW_SWFLUX_CONTROL
348          call ctrl_get_gen (
349         &     xx_swflux_file, xx_swfluxstartdate, xx_swfluxperiod,
350         &     maskc, swflux, xx_swflux0, xx_swflux1, xx_swflux_dummy,
351         &     xx_swflux_remo_intercept, xx_swflux_remo_slope,
352         &     mytime, myiter, mythid )
353    #endif
354    
355    #ifdef ALLOW_SWDOWN_CONTROL
356          call ctrl_get_gen (
357         &     xx_swdown_file, xx_swdownstartdate, xx_swdownperiod,
358         &     maskc, swdown, xx_swdown0, xx_swdown1, xx_swdown_dummy,
359         &     xx_swdown_remo_intercept, xx_swdown_remo_slope,
360         &     mytime, myiter, mythid )
361    #endif
362    
363    #ifdef ALLOW_LWFLUX_CONTROL
364          call ctrl_get_gen (
365         &     xx_lwflux_file, xx_lwfluxstartdate, xx_lwfluxperiod,
366         &     maskc, lwflux, xx_lwflux0, xx_lwflux1, xx_lwflux_dummy,
367         &     xx_lwflux_remo_intercept, xx_lwflux_remo_slope,
368         &     mytime, myiter, mythid )
369    #endif
370    
371    #ifdef ALLOW_LWDOWN_CONTROL
372          call ctrl_get_gen (
373         &     xx_lwdown_file, xx_lwdownstartdate, xx_lwdownperiod,
374         &     maskc, lwdown, xx_lwdown0, xx_lwdown1, xx_lwdown_dummy,
375         &     xx_lwdown_remo_intercept, xx_lwdown_remo_slope,
376         &     mytime, myiter, mythid )
377    #endif
378    
379    #ifdef ALLOW_EVAP_CONTROL
380          call ctrl_get_gen (
381         &     xx_evap_file, xx_evapstartdate, xx_evapperiod,
382         &     maskc, evap, xx_evap0, xx_evap1, xx_evap_dummy,
383         &     xx_evap_remo_intercept, xx_evap_remo_slope,
384         &     mytime, myiter, mythid )
385    #endif
386    
387    #ifdef ALLOW_SNOWPRECIP_CONTROL
388          call ctrl_get_gen (
389         &     xx_snowprecip_file, xx_snowprecipstartdate,
390         &     xx_snowprecipperiod,
391         &     maskc, snowprecip, xx_snowprecip0, xx_snowprecip1,
392         &     xx_snowprecip_dummy,
393         &     xx_snowprecip_remo_intercept, xx_snowprecip_remo_slope,
394         &     mytime, myiter, mythid )
395    #endif
396    
397    #ifdef ALLOW_APRESSURE_CONTROL
398          call ctrl_get_gen (
399         &     xx_apressure_file, xx_apressurestartdate,
400         &     xx_apressureperiod,
401         &     maskc, apressure, xx_apressure0, xx_apressure1,
402         &     xx_apressure_dummy,
403         &     xx_apressure_remo_intercept, xx_apressure_remo_slope,
404         &     mytime, myiter, mythid )
405    #endif
406    
407    #ifdef ALLOW_UWIND_CONTROL
408          call ctrl_get_gen (
409         &     xx_uwind_file, xx_uwindstartdate, xx_uwindperiod,
410         &     maskc, uwind, xx_uwind0, xx_uwind1, xx_uwind_dummy,
411         &     xx_uwind_remo_intercept, xx_uwind_remo_slope,
412         &     mytime, myiter, mythid )
413    #endif /* ALLOW_UWIND_CONTROL */
414    
415    #ifdef ALLOW_VWIND_CONTROL
416          call ctrl_get_gen (
417         &     xx_vwind_file, xx_vwindstartdate, xx_vwindperiod,
418         &     maskc, vwind, xx_vwind0, xx_vwind1, xx_vwind_dummy,
419         &     xx_vwind_remo_intercept, xx_vwind_remo_slope,
420         &     mytime, myiter, mythid )
421    #endif /* ALLOW_VWIND_CONTROL */
422    
423        end        end
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22