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

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

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

revision 1.1.6.2 by heimbach, Thu Apr 4 11:08:03 2002 UTC revision 1.92 by jmc, Fri Mar 3 00:57:00 2017 UTC
# Line 1  Line 1 
1  c $Header$  C $Header$
2    C $Name$
3    
4  #include "EXF_CPPOPTIONS.h"  #include "EXF_OPTIONS.h"
5    #ifdef ALLOW_EXCH2
6    # include "W2_OPTIONS.h"
7        subroutine exf_readparms(  #endif /* ALLOW_EXCH2 */
8       I                          mythid  
9       &                        )        SUBROUTINE EXF_READPARMS( myThid )
10    
11  c     ==================================================================  C     ==================================================================
12  c     SUBROUTINE exf_readparms  C     SUBROUTINE exf_readparms
13  c     ==================================================================  C     ==================================================================
14  c  C
15  c     o This routine initialises the package that calculates external  C     o This routine initialises the package that calculates external
16  c       forcing fields for a given timestep of the MITgcmUV. Parameters  C       forcing fields for a given timestep of the MITgcmUV. Parameters
17  c       for this package are set in "data.externalforcing". Some additional  C       for this package are set in "data.externalforcing". Some additional
18  c       precompiler switches have to be specified in "EXF_CPPOPTIONS.h".  C       precompiler switches have to be specified in "EXF_OPTIONS.h".
19  c  C
20  c     started: Christian Eckert eckert@mit.edu  30-Jun-1999  C     started: Christian Eckert eckert@mit.edu  30-Jun-1999
21  c  C
22  c     changed: Christian Eckert eckert@mit.edu  11-Jan-2000  C     changed: Christian Eckert eckert@mit.edu  11-Jan-2000
23  c              - Restructured the code in order to create a package  C              - Restructured the code in order to create a package
24  c                for the MITgcmUV.  C                for the MITgcmUV.
25  c              Christian Eckert eckert@mit.edu  12-Feb-2000  C              Christian Eckert eckert@mit.edu  12-Feb-2000
26  c              - Changed Routine names (package prefix: exf_)  C              - Changed Routine names (package prefix: exf_)
27  c     changed: Patrick Heimbach, heimbach@mit.edu  04-May-2000  C     changed: Patrick Heimbach, heimbach@mit.edu  04-May-2000
28  c              - changed the handling of precip and sflux with respect  C              - changed the handling of precip and sflux with respect
29  c                to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP  C                to CPP options ALLOW_BULKFORMULAE and ALLOW_ATM_TEMP
30  c     changed: Ralf.Giering@FastOpt.de 25-Mai-20000  C     changed: Ralf.Giering@FastOpt.de 25-Mai-20000
31  c              - moved relaxation and climatology to extra routines  C              - moved relaxation and climatology to extra routines
32  c              Patrick Heimbach, heimbach@mit.edu  04-May-2000  C              Patrick Heimbach, heimbach@mit.edu  04-May-2000
33  c              - added obcs parameters  C              - added obcs parameters
34  c     changed: Virginie Thierry, vthierry@ucsd.edu 04-June-2001  C     changed: Virginie Thierry, vthierry@ucsd.edu 04-June-2001
35  c              - added new obcs parameters (for each boundaries)  C              - added new obcs parameters (for each boundaries)
36  c     included runoff D. Stammer, Nov. 25, 2001  C     included runoff D. Stammer, Nov. 25, 2001
37  c  C     included pressure forcing. heimbach@mit.edu 05-Nov-2002
38  c     ==================================================================  C     added "repeatPeriod" for cycling of forcing datasets 19-Dec-2002
39  c     SUBROUTINE exf_readparms  C     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
40  c     ==================================================================  C
41    C     ==================================================================
42    C     SUBROUTINE exf_readparms
43    C     ==================================================================
44    
45        implicit none        implicit none
46    
47  c     == global variables ==  C     == global variables ==
48    
49  #include "EEPARAMS.h"  #include "EEPARAMS.h"
50  #include "SIZE.h"  #include "SIZE.h"
51  #include "cal.h"  #include "PARAMS.h"
52  #include "exf.h"  #ifdef USE_EXF_INTERPOLATION
53  #include "exf_param.h"  # ifdef ALLOW_EXCH2
54  #include "exf_constants.h"  # include "W2_EXCH2_SIZE.h"
55    # include "W2_EXCH2_TOPOLOGY.h"
56    # endif /* ALLOW_EXCH2 */
57    # include "SET_GRID.h"
58    #endif /* USE_EXF_INTERPOLATION */
59    c#include "cal.h"
60    #include "EXF_PARAM.h"
61    #include "EXF_CONSTANTS.h"
62    
63    C     == routine arguments ==
64          INTEGER myThid
65    
66    C     == local variables ==
67    #ifdef USE_EXF_INTERPOLATION
68          INTEGER gridNx, gridNy
69          INTEGER j
70          _RL inp_lon0, inp_lat0, inp_dLon, inp_dLat
71    #endif /* USE_EXF_INTERPOLATION */
72          INTEGER iUnit
73          LOGICAL exf_verbose
74          CHARACTER*(2) exf_yftype
75          CHARACTER*(MAX_LEN_MBUF) msgBuf
76    C     == end of interface ==
77    
78    C     Surface flux data.
79          NAMELIST /EXF_NML_01/
80         &      windstressmax,       repeatPeriod,    exf_albedo,
81         &   ocean_emissivity,     ice_emissivity, snow_emissivity,
82         &          exf_iceCd,          exf_iceCe,     exf_iceCh,
83         &   exf_scal_BulkCdn,     climtempfreeze, useExfCheckRange,
84         &      exf_iprec    ,   exf_iprec_obcs  , exf_yftype,
85         &      exf_verbose  ,   exf_debugLev    , exf_monFreq,
86         & useExfYearlyFields,  twoDigitYear,
87         & useStabilityFct_overIce, readStressOnAgrid, readStressOnCgrid,
88         & rotateStressOnAgrid, useAtmWind, useRelativeWind, noNegativeEvap,
89         & select_ZenAlbedo, useExfZenIncoming,
90         & hu, ht, umin, atmrho, atmcp, cen2kel, gravity_mks,
91         & cdrag_1, cdrag_2, cdrag_3, cstanton_1, cstanton_2, cdalton,
92         & flamb, flami, zolmin, zref,
93         & cvapor_fac, cvapor_exp, cvapor_fac_ice, cvapor_exp_ice,
94         & humid_fac, gamma_blk, saltsat, sstExtrapol, psim_fac
95    
96          NAMELIST /EXF_NML_02/
97         &          hfluxfile,     hfluxstartdate1,     hfluxstartdate2,
98         &          hfluxRepCycle,     hfluxperiod,     hfluxStartTime,
99         &          atempfile,     atempstartdate1,     atempstartdate2,
100         &          atempRepCycle,     atempperiod,     atempStartTime,
101         &            aqhfile,       aqhstartdate1,       aqhstartdate2,
102         &            aqhRepCycle,       aqhperiod,       aqhStartTime,
103         &            hs_file,       hs_startdate1,       hs_startdate2,
104         &            hs_RepCycle,       hs_period,       hs_StartTime,
105         &            hl_file,       hl_startdate1,       hl_startdate2,
106         &            hl_RepCycle,       hl_period,       hl_StartTime,
107         &          sfluxfile,     sfluxstartdate1,     sfluxstartdate2,
108         &          sfluxRepCycle,     sfluxperiod,     sfluxStartTime,
109         &           evapfile,      evapstartdate1,      evapstartdate2,
110         &           evapRepCycle,      evapperiod,      evapStartTime,
111         &         precipfile,    precipstartdate1,    precipstartdate2,
112         &         precipRepCycle,    precipperiod,    precipStartTime,
113         &     snowprecipfile, snowprecipstartdate1, snowprecipstartdate2,
114         &     snowprecipRepCycle, snowprecipperiod, snowprecipStartTime,
115         &         runofffile,    runoffstartdate1,    runoffstartdate2,
116         &         runoffRepCycle,    runoffperiod,    runoffStartTime,
117         &      runoftempfile,
118         &        saltflxfile,   saltflxstartdate1,   saltflxstartdate2,
119         &        saltflxRepCycle,   saltflxperiod,   saltflxStartTime,
120         &        ustressfile,   ustressstartdate1,   ustressstartdate2,
121         &        ustressRepCycle,   ustressperiod,   ustressStartTime,
122         &        vstressfile,   vstressstartdate1,   vstressstartdate2,
123         &        vstressRepCycle,   vstressperiod,   vstressStartTime,
124         &          uwindfile,     uwindstartdate1,     uwindstartdate2,
125         &          uwindRepCycle,     uwindperiod,     uwindStartTime,
126         &          vwindfile,     vwindstartdate1,     vwindstartdate2,
127         &          vwindRepCycle,     vwindperiod,     vwindStartTime,
128         &         wspeedfile,    wspeedstartdate1,    wspeedstartdate2,
129         &         wspeedRepCycle,    wspeedperiod,    wspeedStartTime,
130         &         swfluxfile,    swfluxstartdate1,    swfluxstartdate2,
131         &         swfluxRepCycle,    swfluxperiod,    swfluxStartTime,
132         &         lwfluxfile,    lwfluxstartdate1,    lwfluxstartdate2,
133         &         lwfluxRepCycle,    lwfluxperiod,    lwfluxStartTime,
134         &         swdownfile,    swdownstartdate1,    swdownstartdate2,
135         &         swdownRepCycle,    swdownperiod,    swdownStartTime,
136         &         lwdownfile,    lwdownstartdate1,    lwdownstartdate2,
137         &         lwdownRepCycle,    lwdownperiod,    lwdownStartTime,
138         &      apressurefile, apressurestartdate1, apressurestartdate2,
139         &      apressureRepCycle, apressureperiod, apressureStartTime,
140         &       areamaskfile,  areamaskstartdate1,  areamaskstartdate2,
141         &       areamaskRepCycle,  areamaskperiod,  areamaskStartTime,
142         &        climsstfile,   climsststartdate1,   climsststartdate2,
143         &        climsstRepCycle,   climsstperiod,   climsstStartTime,
144         &        climsssfile,   climsssstartdate1,   climsssstartdate2,
145         &        climsssRepCycle,   climsssperiod,   climsssStartTime,
146         &       climustrfile,  climustrstartdate1,  climustrstartdate2,
147         &       climustrRepCycle,  climustrperiod,  climustrStartTime,
148         &       climvstrfile,  climvstrstartdate1,  climvstrstartdate2,
149         &       climvstrRepCycle,  climvstrperiod,  climvstrStartTime,
150         &     areamaskTauRelax,  climsstTauRelax,    climsssTauRelax,
151         &                       climustrTauRelax,   climvstrTauRelax
152    
153          NAMELIST /EXF_NML_03/
154         &   exf_inscal_hflux,  exf_inscal_sflux,    exf_inscal_evap,
155         & exf_inscal_ustress,  exf_inscal_vstress,
156         &   exf_inscal_uwind,  exf_inscal_vwind,    exf_inscal_wspeed,
157         &   exf_inscal_atemp,  exf_offset_atemp,
158         &   exf_inscal_aqh,    exf_inscal_hs,       exf_inscal_hl,
159         &     exf_inscal_sst,  exf_inscal_sss,
160         &  exf_inscal_swflux,  exf_inscal_lwflux,   exf_inscal_precip,
161         &  exf_inscal_runoff,  exf_inscal_apressure, exf_inscal_snowprecip,
162         &  exf_inscal_runoftemp, exf_inscal_saltflx,
163         &  exf_inscal_swdown,  exf_inscal_lwdown,
164         & exf_inscal_climsst, exf_inscal_climsss,
165         & exf_inscal_climustr, exf_inscal_climvstr,
166         &  exf_outscal_hflux,  exf_outscal_ustress, exf_outscal_vstress,
167         & exf_outscal_swflux,  exf_outscal_sst,     exf_outscal_sss,
168         &  exf_outscal_sflux,  exf_outscal_apressure,
169         & exf_inscal_areamask, exf_outscal_areamask,
170         &  hfluxconst, atempconst, aqhconst, hs_const, hl_const,
171         &  sfluxconst, evapconst, precipconst, snowprecipconst,
172         &  runoffconst, runoftempconst, saltflxconst, ustressconst,
173         &  vstressconst, uwindconst, vwindconst, wspeedconst, swfluxconst,
174         &  lwfluxconst, swdownconst, lwdownconst, apressureconst,
175         &  areamaskconst, climsstconst,   climsssconst,
176         &  climustrconst, climvstrconst,
177         &     hflux_exfremo_intercept, hflux_exfremo_slope,
178         &     atemp_exfremo_intercept, atemp_exfremo_slope,
179         &     aqh_exfremo_intercept, aqh_exfremo_slope,
180         &     hs_exfremo_intercept,  hs_exfremo_slope,
181         &     hl_exfremo_intercept,  hl_exfremo_slope,
182         &     sflux_exfremo_intercept, sflux_exfremo_slope,
183         &     evap_exfremo_intercept, evap_exfremo_slope,
184         &     precip_exfremo_intercept, precip_exfremo_slope,
185         &     snowprecip_exfremo_intercept, snowprecip_exfremo_slope,
186         &     runoff_exfremo_intercept, runoff_exfremo_slope,
187         &     runoftemp_exfremo_intercept, runoftemp_exfremo_slope,
188         &     saltflx_exfremo_intercept, saltflx_exfremo_slope,
189         &     ustress_exfremo_intercept, ustress_exfremo_slope,
190         &     vstress_exfremo_intercept, vstress_exfremo_slope,
191         &     uwind_exfremo_intercept, uwind_exfremo_slope,
192         &     vwind_exfremo_intercept, vwind_exfremo_slope,
193         &     wspeed_exfremo_intercept, wspeed_exfremo_slope,
194         &     swflux_exfremo_intercept, swflux_exfremo_slope,
195         &     lwflux_exfremo_intercept, lwflux_exfremo_slope,
196         &     swdown_exfremo_intercept, swdown_exfremo_slope,
197         &     lwdown_exfremo_intercept, lwdown_exfremo_slope,
198         &     apressure_exfremo_intercept, apressure_exfremo_slope,
199         &     areamask_exfremo_intercept, areamask_exfremo_slope,
200         &     climsst_exfremo_intercept, climsst_exfremo_slope,
201         &     climsss_exfremo_intercept, climsss_exfremo_slope,
202         &     climustr_exfremo_intercept, climustr_exfremo_slope,
203         &     climvstr_exfremo_intercept, climvstr_exfremo_slope
204    
205    #ifdef USE_EXF_INTERPOLATION
206          NAMELIST /EXF_NML_04/
207         & ustress_lon0, ustress_lon_inc, ustress_lat0, ustress_lat_inc,
208         & ustress_nlon, ustress_nlat, ustress_interpMethod,
209         & vstress_lon0, vstress_lon_inc, vstress_lat0, vstress_lat_inc,
210         & vstress_nlon, vstress_nlat, vstress_interpMethod,
211         & hflux_lon0, hflux_lon_inc, hflux_lat0, hflux_lat_inc,
212         & hflux_nlon, hflux_nlat, hflux_interpMethod,
213         & sflux_lon0, sflux_lon_inc, sflux_lat0, sflux_lat_inc,
214         & sflux_nlon, sflux_nlat, sflux_interpMethod,
215         & swflux_lon0, swflux_lon_inc, swflux_lat0, swflux_lat_inc,
216         & swflux_nlon, swflux_nlat, swflux_interpMethod,
217         & lwflux_lon0, lwflux_lon_inc, lwflux_lat0, lwflux_lat_inc,
218         & lwflux_nlon, lwflux_nlat, lwflux_interpMethod,
219         & atemp_lon0, atemp_lon_inc, atemp_lat0, atemp_lat_inc,
220         & atemp_nlon, atemp_nlat, atemp_interpMethod,
221         & aqh_lon0, aqh_lon_inc, aqh_lat0, aqh_lat_inc,
222         & aqh_nlon, aqh_nlat, aqh_interpMethod,
223         & hs_lon0, hs_lon_inc, hs_lat0, hs_lat_inc,
224         & hs_nlon, hs_nlat, hs_interpMethod,
225         & hl_lon0, hl_lon_inc, hl_lat0, hl_lat_inc,
226         & hl_nlon, hl_nlat, hl_interpMethod,
227         & evap_lon0, evap_lon_inc, evap_lat0, evap_lat_inc,
228         & evap_nlon, evap_nlat, evap_interpMethod,
229         & precip_lon0, precip_lon_inc, precip_lat0, precip_lat_inc,
230         & precip_nlon, precip_nlat, precip_interpMethod,
231         & runoff_lon0, runoff_lon_inc, runoff_lat0, runoff_lat_inc,
232         & runoff_nlon, runoff_nlat, runoff_interpMethod,
233         & saltflx_lon0, saltflx_lon_inc,
234         & saltflx_lat0, saltflx_lat_inc,
235         & saltflx_nlon, saltflx_nlat, saltflx_interpMethod,
236         & snowprecip_lon0, snowprecip_lon_inc,
237         & snowprecip_lat0, snowprecip_lat_inc,
238         & snowprecip_nlon, snowprecip_nlat, snowprecip_interpMethod,
239         & uwind_lon0, uwind_lon_inc, uwind_lat0, uwind_lat_inc,
240         & uwind_nlon, uwind_nlat, uwind_interpMethod,
241         & vwind_lon0, vwind_lon_inc, vwind_lat0, vwind_lat_inc,
242         & vwind_nlon, vwind_nlat, vwind_interpMethod,
243         & wspeed_lon0, wspeed_lon_inc, wspeed_lat0, wspeed_lat_inc,
244         & wspeed_nlon, wspeed_nlat, wspeed_interpMethod,
245         & swdown_lon0, swdown_lon_inc, swdown_lat0, swdown_lat_inc,
246         & swdown_nlon, swdown_nlat, swdown_interpMethod,
247         & lwdown_lon0, lwdown_lon_inc, lwdown_lat0, lwdown_lat_inc,
248         & lwdown_nlon, lwdown_nlat, lwdown_interpMethod,
249         & apressure_lon0, apressure_lon_inc,
250         & apressure_lat0, apressure_lat_inc,
251         & apressure_nlon, apressure_nlat, apressure_interpMethod,
252         & areamask_lon0, areamask_lon_inc, areamask_lat0, areamask_lat_inc,
253         & areamask_nlon, areamask_nlat, areamask_interpMethod,
254         & climsst_lon0, climsst_lon_inc, climsst_lat0, climsst_lat_inc,
255         & climsst_nlon, climsst_nlat, climsst_interpMethod,
256         & climsss_lon0, climsss_lon_inc,climsss_lat0, climsss_lat_inc,
257         & climsss_nlon, climsss_nlat, climsss_interpMethod,
258         & climustr_lon0, climustr_lon_inc, climustr_lat0, climustr_lat_inc,
259         & climustr_nlon, climustr_nlat, climustr_interpMethod,
260         & climvstr_lon0, climvstr_lon_inc, climvstr_lat0, climvstr_lat_inc,
261         & climvstr_nlon, climvstr_nlat, climvstr_interpMethod,
262         & exf_output_interp
263    #endif /* USE_EXF_INTERPOLATION */
264    
265  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
266  # include "OBCS_OPTIONS.h"        NAMELIST /EXF_NML_OBCS/
267         &    useOBCSYearlyFields,
268         &    obcsNstartdate1,   obcsNstartdate2,   obcsNstartTime,
269         &        obcsNperiod,   obcsNrepCycle,
270         &    obcsSstartdate1,   obcsSstartdate2,   obcsSstartTime,
271         &        obcsSperiod,   obcsSrepCycle,
272         &    obcsEstartdate1,   obcsEstartdate2,   obcsEstartTime,
273         &        obcsEperiod,   obcsErepCycle,
274         &    obcsWstartdate1,   obcsWstartdate2,   obcsWstartTime,
275         &        obcsWperiod,   obcsWrepCycle,
276         &    siobNstartdate1,   siobNstartdate2,   siobNstartTime,
277         &        siobNperiod,   siobNrepCycle,
278         &    siobSstartdate1,   siobSstartdate2,   siobSstartTime,
279         &        siobSperiod,   siobSrepCycle,
280         &    siobEstartdate1,   siobEstartdate2,   siobEstartTime,
281         &        siobEperiod,   siobErepCycle,
282         &    siobWstartdate1,   siobWstartdate2,   siobWstartTime,
283         &        siobWperiod,   siobWrepCycle
284    #endif /* ALLOW_OBCS */
285    
286    #ifdef USE_EXF_INTERPOLATION
287    # ifdef ALLOW_EXCH2
288          gridNx = exch2_mydNx(1)
289          gridNy = exch2_mydNy(1)
290    # else /* ALLOW_EXCH2 */
291          gridNx = Nx
292          gridNy = Ny
293    # endif /* ALLOW_EXCH2 */
294    #endif /* USE_EXF_INTERPOLATION */
295    
296          IF ( .NOT.useEXF ) THEN
297    C-    pkg EXF is not used
298            _BEGIN_MASTER(myThid)
299    C-    Track pkg activation status:
300    C     print a (weak) warning if data.exf is found
301             CALL PACKAGES_UNUSED_MSG( 'useEXF', ' ', ' ' )
302            _END_MASTER(myThid)
303            RETURN
304          ENDIF
305    
306          _BEGIN_MASTER(myThid)
307    
308    C     Set default values.
309    
310    c     exf_verbose        = debugMode
311          exf_verbose        = .FALSE.
312          exf_debugLev       = debugLevel
313          exf_monFreq        = monitorFreq
314          useExfCheckRange   = .TRUE.
315          select_ZenAlbedo   = 0
316          useExfZenIncoming  = .FALSE.
317          readStressOnAgrid  = .FALSE.
318          rotateStressOnAgrid = .FALSE.
319          readStressOnCgrid  = .FALSE.
320    #ifdef ALLOW_ATM_WIND
321          useAtmWind         = .TRUE.
322    #else
323          useAtmWind         = .FALSE.
324  #endif  #endif
325          useRelativeWind    = .FALSE.
326          noNegativeEvap     = .FALSE.
327    
328  c     == routine arguments ==  C-  default value should be set to main model parameter:
329    c     cen2kel     =  celsius2K
330        integer mythid  c     gravity_mks = gravity
331    c     atmcp       =  atm_Cp
332  c     == local variables ==  c     humid_fac   =  atm_Rq     <- default is zero !!!
333    
334        integer i        cen2kel        =      273.150  _d 0
335          gravity_mks    =        9.81   _d 0
336  c     == external ==        atmrho         =        1.200  _d 0
337          atmcp          =     1005.000  _d 0
338        integer  ilnblnk        flamb          =  2500000.000  _d 0
339        external ilnblnk        flami          =   334000.000  _d 0
340          cvapor_fac     =   640380.000  _d 0
341  c     == end of interface ==        cvapor_exp     =     5107.400  _d 0
342          cvapor_fac_ice = 11637800.000  _d 0
343  c     Surface flux data.        cvapor_exp_ice =     5897.800  _d 0
344        namelist /exf_nml/        humid_fac      =        0.606  _d 0
345       &    hfluxstartdate1,    hfluxstartdate2,   hfluxperiod,        gamma_blk      =        0.010  _d 0
346       &    atempstartdate1,    atempstartdate2,   atempperiod,        saltsat        =        0.980  _d 0
347       &      aqhstartdate1,      aqhstartdate2,     aqhperiod,        sstExtrapol    =        0.0    _d 0
348       &    sfluxstartdate1,    sfluxstartdate2,   sfluxperiod,        cdrag_1        =        0.0027000 _d 0
349       &   precipstartdate1,   precipstartdate2,  precipperiod,        cdrag_2        =        0.0001420 _d 0
350       &   runoffstartdate1,   runoffstartdate2,  runoffperiod,        cdrag_3        =        0.0000764 _d 0
351       &  ustressstartdate1,  ustressstartdate2, ustressperiod,        cstanton_1     =        0.0327 _d 0
352       &  vstressstartdate1,  vstressstartdate2, vstressperiod,        cstanton_2     =        0.0180 _d 0
353       &    uwindstartdate1,    uwindstartdate2,   uwindperiod,        cdalton        =        0.0346 _d 0
354       &    vwindstartdate1,    vwindstartdate2,   vwindperiod,        zolmin         =     -100.000  _d 0
355       &   swfluxstartdate1,   swfluxstartdate2,  swfluxperiod,        psim_fac       =        5.000  _d 0
356       &   lwfluxstartdate1,   lwfluxstartdate2,  lwfluxperiod,        zref           =       10.000  _d 0
357       &    obcsNstartdate1,    obcsNstartdate2,   obcsNperiod,        hu             =       10.000  _d 0
358       &    obcsSstartdate1,    obcsSstartdate2,   obcsSperiod,        ht             =        2.000  _d 0
359       &    obcsEstartdate1,    obcsEstartdate2,   obcsEperiod,        umin           =        0.5    _d 0
360       &    obcsWstartdate1,    obcsWstartdate2,   obcsWperiod,        useStabilityFct_overIce = .FALSE.
361       &          hfluxfile,          atempfile,       aqhfile,        exf_iceCd        = 1.63 _d -3
362       &          sfluxfile,         precipfile,    runofffile,        exf_iceCe        = 1.63 _d -3
363       &        ustressfile,        vstressfile,        exf_iceCh        = 1.63 _d -3
364       &          uwindfile,          vwindfile,        exf_albedo       = 0.1 _d 0
365       &         swfluxfile,         lwfluxfile,  C--   this default is chosen to be backward compatible with
366       &          exf_iprec,         exf_yftype  C--   an earlier setting of 5.5 = ocean_emissivity*stefanBoltzmann
367          ocean_emissivity = 5.50 _d -8 / 5.670 _d -8
368        _BEGIN_MASTER(mythid)        ice_emissivity   = 0.95 _d 0
369          snow_emissivity  = 0.95 _d 0
 c     Set default values.  
370    
371  c     Calendar data.  C     Calendar data.
372        hfluxstartdate1    = 0        hfluxstartdate1    = 0
373        hfluxstartdate2    = 0        hfluxstartdate2    = 0
374        hfluxperiod        = 0.0 _d 0        hfluxperiod        = 0.0 _d 0
375          hfluxconst         = 0.0 _d 0
376          hflux_exfremo_intercept = 0.0 _d 0
377          hflux_exfremo_slope = 0.0 _d 0
378    
379        atempstartdate1    = 0        atempstartdate1    = 0
380        atempstartdate2    = 0        atempstartdate2    = 0
381        atempperiod        = 0.0 _d 0        atempperiod        = 0.0 _d 0
382          atempconst         = celsius2K
383          atemp_exfremo_intercept = 0.0 _d 0
384          atemp_exfremo_slope = 0.0 _d 0
385    
386        aqhstartdate1      = 0        aqhstartdate1      = 0
387        aqhstartdate2      = 0        aqhstartdate2      = 0
388        aqhperiod          = 0.0 _d 0        aqhperiod          = 0.0 _d 0
389          aqhconst           = 0.0 _d 0
390          aqh_exfremo_intercept = 0.0 _d 0
391          aqh_exfremo_slope  = 0.0 _d 0
392    
393          hs_startdate1      = 0
394          hs_startdate2      = 0
395          hs_period          = 0.0 _d 0
396          hs_const           = 0.0 _d 0
397          hs_exfremo_intercept = 0.0 _d 0
398          hs_exfremo_slope   = 0.0 _d 0
399    
400          hl_startdate1      = 0
401          hl_startdate2      = 0
402          hl_period          = 0.0 _d 0
403          hl_const           = 0.0 _d 0
404          hl_exfremo_intercept = 0.0 _d 0
405          hl_exfremo_slope   = 0.0 _d 0
406    
407        sfluxstartdate1    = 0        sfluxstartdate1    = 0
408        sfluxstartdate2    = 0        sfluxstartdate2    = 0
409        sfluxperiod        = 0.0 _d 0        sfluxperiod        = 0.0 _d 0
410          sfluxconst         = 0.0 _d 0
411          sflux_exfremo_intercept = 0.0 _d 0
412          sflux_exfremo_slope = 0.0 _d 0
413    
414          evapstartdate1   = 0
415          evapstartdate2   = 0
416          evapperiod       = 0.0 _d 0
417          evapconst        = 0.0 _d 0
418          evap_exfremo_intercept = 0.0 _d 0
419          evap_exfremo_slope = 0.0 _d 0
420    
421        precipstartdate1   = 0        precipstartdate1   = 0
422        precipstartdate2   = 0        precipstartdate2   = 0
423        precipperiod       = 0.0 _d 0        precipperiod       = 0.0 _d 0
424          precipconst        = 0.0 _d 0
425          precip_exfremo_intercept = 0.0 _d 0
426          precip_exfremo_slope = 0.0 _d 0
427    
428          snowprecipstartdate1   = 0
429          snowprecipstartdate2   = 0
430          snowprecipperiod       = 0.0 _d 0
431          snowprecipconst        = 0.0 _d 0
432          snowprecip_exfremo_intercept = 0.0 _d 0
433          snowprecip_exfremo_slope = 0.0 _d 0
434    
435        runoffstartdate1   = 0        runoffstartdate1   = 0
436        runoffstartdate2   = 0        runoffstartdate2   = 0
437        runoffperiod       = 0.0 _d 0        runoffperiod       = 0.0 _d 0
438          runoffconst        = 0.0 _d 0
439          runoff_exfremo_intercept = 0.0 _d 0
440          runoff_exfremo_slope = 0.0 _d 0
441    
442          runoftempconst              = 0.0 _d 0
443          runoftemp_exfremo_intercept = 0.0 _d 0
444          runoftemp_exfremo_slope     = 0.0 _d 0
445    
446          saltflxstartdate1  = 0
447          saltflxstartdate2  = 0
448          saltflxperiod      = 0.0 _d 0
449          saltflxconst       = 0.0 _d 0
450          saltflx_exfremo_intercept = 0.0 _d 0
451          saltflx_exfremo_slope = 0.0 _d 0
452    
453        ustressstartdate1  = 0        ustressstartdate1  = 0
454        ustressstartdate2  = 0        ustressstartdate2  = 0
455        ustressperiod      = 0.0 _d 0        ustressperiod      = 0.0 _d 0
456          ustressconst       = 0.0 _d 0
457          ustress_exfremo_intercept = 0.0 _d 0
458          ustress_exfremo_slope = 0.0 _d 0
459    
460        vstressstartdate1  = 0        vstressstartdate1  = 0
461        vstressstartdate2  = 0        vstressstartdate2  = 0
462        vstressperiod      = 0.0 _d 0        vstressperiod      = 0.0 _d 0
463          vstressconst       = 0.0 _d 0
464          vstress_exfremo_intercept = 0.0 _d 0
465          vstress_exfremo_slope = 0.0 _d 0
466    
467        uwindstartdate1    = 0        uwindstartdate1    = 0
468        uwindstartdate2    = 0        uwindstartdate2    = 0
469        uwindperiod        = 0.0 _d 0        uwindperiod        = 0.0 _d 0
470          uwindconst         = 0.0 _d 0
471          uwind_exfremo_intercept = 0.0 _d 0
472          uwind_exfremo_slope = 0.0 _d 0
473    
474        vwindstartdate1    = 0        vwindstartdate1    = 0
475        vwindstartdate2    = 0        vwindstartdate2    = 0
476        vwindperiod        = 0.0 _d 0        vwindperiod        = 0.0 _d 0
477          vwindconst         = 0.0 _d 0
478          vwind_exfremo_intercept = 0.0 _d 0
479          vwind_exfremo_slope = 0.0 _d 0
480    
481          wspeedstartdate1    = 0
482          wspeedstartdate2    = 0
483          wspeedperiod        = 0.0 _d 0
484          wspeedconst         = 0.0 _d 0
485          wspeed_exfremo_intercept = 0.0 _d 0
486          wspeed_exfremo_slope = 0.0 _d 0
487    
488        swfluxstartdate1   = 0        swfluxstartdate1   = 0
489        swfluxstartdate2   = 0        swfluxstartdate2   = 0
490        swfluxperiod       = 0.0 _d 0        swfluxperiod       = 0.0 _d 0
491          swfluxconst        = 0.0 _d 0
492          swflux_exfremo_intercept = 0.0 _d 0
493          swflux_exfremo_slope = 0.0 _d 0
494    
495        lwfluxstartdate1   = 0        lwfluxstartdate1   = 0
496        lwfluxstartdate2   = 0        lwfluxstartdate2   = 0
497        lwfluxperiod       = 0.0 _d 0        lwfluxperiod       = 0.0 _d 0
498          lwfluxconst        = 0.0 _d 0
499          lwflux_exfremo_intercept = 0.0 _d 0
500          lwflux_exfremo_slope = 0.0 _d 0
501    
502          swdownstartdate1   = 0
503          swdownstartdate2   = 0
504          swdownperiod       = 0.0 _d 0
505          swdownconst        = 0.0 _d 0
506          swdown_exfremo_intercept = 0.0 _d 0
507          swdown_exfremo_slope = 0.0 _d 0
508    
509          lwdownstartdate1   = 0
510          lwdownstartdate2   = 0
511          lwdownperiod       = 0.0 _d 0
512          lwdownconst        = 0.0 _d 0
513          lwdown_exfremo_intercept = 0.0 _d 0
514          lwdown_exfremo_slope = 0.0 _d 0
515    
516          apressurestartdate1    = 0
517          apressurestartdate2    = 0
518          apressureperiod        = 0.0 _d 0
519          apressureconst         = 0.0 _d 0
520          apressure_exfremo_intercept = 0.0 _d 0
521          apressure_exfremo_slope = 0.0 _d 0
522    
523          areamaskstartdate1    = 0
524          areamaskstartdate2    = 0
525          areamaskperiod        = 0.0 _d 0
526          areamaskTauRelax      = 0.0 _d 0
527          areamaskconst         = 0.0 _d 0
528          areamask_exfremo_intercept = 0. _d 0
529          areamask_exfremo_slope = 0. _d 0
530    
531          climsststartdate1  = 0
532          climsststartdate2  = 0
533          climsstperiod      = 0
534          climsstTauRelax    = 0.0 _d 0
535          climsstconst         = 0.0 _d 0
536          climsst_exfremo_intercept = 0.0 _d 0
537          climsst_exfremo_slope = 0.0 _d 0
538    
539          climsssstartdate1  = 0
540          climsssstartdate2  = 0
541          climsssperiod      = 0
542          climsssTauRelax    = 0.0 _d 0
543          climsssconst         = 0.0 _d 0
544          climsss_exfremo_intercept = 0.0 _d 0
545          climsss_exfremo_slope = 0.0 _d 0
546    
547          climustrstartdate1  = 0
548          climustrstartdate2  = 0
549          climustrperiod      = 0
550          climustrTauRelax    = 0.0 _d 0
551          climustrconst         = 0.0 _d 0
552          climustr_exfremo_intercept = 0.0 _d 0
553          climustr_exfremo_slope = 0.0 _d 0
554    
555          climvstrstartdate1  = 0
556          climvstrstartdate2  = 0
557          climvstrperiod      = 0
558          climvstrTauRelax    = 0.0 _d 0
559          climvstrconst         = 0.0 _d 0
560          climvstr_exfremo_intercept = 0.0 _d 0
561          climvstr_exfremo_slope = 0.0 _d 0
562    
563          useOBCSYearlyFields = .FALSE.
564        obcsNstartdate1    = 0        obcsNstartdate1    = 0
565        obcsNstartdate2    = 0        obcsNstartdate2    = 0
566        obcsNperiod        = 0.0 _d 0        obcsNperiod        = 0.0 _d 0
   
567        obcsSstartdate1    = 0        obcsSstartdate1    = 0
568        obcsSstartdate2    = 0        obcsSstartdate2    = 0
569        obcsSperiod        = 0.0 _d 0        obcsSperiod        = 0.0 _d 0
   
570        obcsEstartdate1    = 0        obcsEstartdate1    = 0
571        obcsEstartdate2    = 0        obcsEstartdate2    = 0
572        obcsEperiod        = 0.0 _d 0        obcsEperiod        = 0.0 _d 0
   
573        obcsWstartdate1    = 0        obcsWstartdate1    = 0
574        obcsWstartdate2    = 0        obcsWstartdate2    = 0
575        obcsWperiod        = 0.0 _d 0        obcsWperiod        = 0.0 _d 0
576    
577  c     Data files.        siobNstartdate1    = UNSET_I
578          siobNstartdate2    = UNSET_I
579          siobNperiod        = UNSET_RL
580          siobSstartdate1    = UNSET_I
581          siobSstartdate2    = UNSET_I
582          siobSperiod        = UNSET_RL
583          siobEstartdate1    = UNSET_I
584          siobEstartdate2    = UNSET_I
585          siobEperiod        = UNSET_RL
586          siobWstartdate1    = UNSET_I
587          siobWstartdate2    = UNSET_I
588          siobWperiod        = UNSET_RL
589    
590          repeatPeriod       = 0.0 _d 0
591          windstressmax      = 2.0 _d 0
592    
593          exf_scal_BulkCdn   = 1.0  _d 0
594    
595    C     Initialise freezing temperature of sea water
596          climtempfreeze     = -1.9 _d 0
597    
598    C     Data files.
599        hfluxfile          = ' '        hfluxfile          = ' '
600        atempfile          = ' '        atempfile          = ' '
601        aqhfile            = ' '        aqhfile            = ' '
602          hs_file            = ' '
603          hl_file            = ' '
604          evapfile           = ' '
605        precipfile         = ' '        precipfile         = ' '
606          snowprecipfile     = ' '
607        sfluxfile          = ' '        sfluxfile          = ' '
608        runofffile         = ' '        runofffile         = ' '
609          runoftempfile      = ' '
610          saltflxfile        = ' '
611        ustressfile        = ' '        ustressfile        = ' '
612        vstressfile        = ' '        vstressfile        = ' '
613        uwindfile          = ' '        uwindfile          = ' '
614        vwindfile          = ' '        vwindfile          = ' '
615          wspeedfile         = ' '
616        swfluxfile         = ' '        swfluxfile         = ' '
617        lwfluxfile         = ' '        lwfluxfile         = ' '
618          swdownfile         = ' '
619  c     Initialise the date arrays.        lwdownfile         = ' '
620        do i = 1,4        apressurefile      = ' '
621           hfluxstartdate(i)    = 0        areamaskfile       = ' '
622           atempstartdate(i)    = 0        climsstfile        = ' '
623           aqhstartdate(i)      = 0        climsssfile        = ' '
624           precipstartdate(i)   = 0        climustrfile       = ' '
625           sfluxstartdate(i)    = 0        climvstrfile       = ' '
626           runoffstartdate(i)   = 0  
627           ustressstartdate(i)  = 0  C     Start Time.
628           vstressstartdate(i)  = 0        hfluxStartTime     = UNSET_RL
629           uwindstartdate(i)    = 0        atempStartTime     = UNSET_RL
630           vwindstartdate(i)    = 0        aqhStartTime       = UNSET_RL
631           swfluxstartdate(i)   = 0        hs_StartTime       = UNSET_RL
632           lwfluxstartdate(i)   = 0        hl_StartTime       = UNSET_RL
633           obcsNstartdate(i)    = 0        evapStartTime      = UNSET_RL
634           obcsSstartdate(i)    = 0        precipStartTime    = UNSET_RL
635           obcsEstartdate(i)    = 0        snowprecipStartTime= UNSET_RL
636           obcsWstartdate(i)    = 0        sfluxStartTime     = UNSET_RL
637        enddo        runoffStartTime    = UNSET_RL
638          saltflxStartTime   = UNSET_RL
639  c     Initialise file type and field precision        ustressStartTime   = UNSET_RL
640        exf_iprec       = 32        vstressStartTime   = UNSET_RL
641        exf_yftype      = 'RL'        uwindStartTime     = UNSET_RL
642          vwindStartTime     = UNSET_RL
643  c     Check for the availability of the right calendar version.        wspeedStartTime    = UNSET_RL
644        if ( calendarversion .ne. usescalendarversion ) then        swfluxStartTime    = UNSET_RL
645           print*,' exf_readparms: You are not using the appropriate'        lwfluxStartTime    = UNSET_RL
646           print*,'           version of the calendar package.'        swdownStartTime    = UNSET_RL
647           print*        lwdownStartTime    = UNSET_RL
648           print*,' You are using Calendar version: ', calendarversion        apressureStartTime = UNSET_RL
649           print*,' Please use    Calendar version: ', usescalendarversion        areamaskStartTime  = UNSET_RL
650           stop ' stopped in exf_readparms.'        climsstStartTime   = UNSET_RL
651        endif        climsssStartTime   = UNSET_RL
652          climustrStartTime  = UNSET_RL
653  c     Next, read the forcing data file.        climvstrStartTime  = UNSET_RL
654        call nml_filter( 'data.exf', scrunit1, myThid )        obcsNstartTime     = UNSET_RL
655        if (scrunit1 .eq. 0) then        obcsSstartTime     = UNSET_RL
656           stop 'exf_readparms: reading namelist failed'        obcsEstartTime     = UNSET_RL
657        end if        obcsWstartTime     = UNSET_RL
658        read(  scrunit1, nml = exf_nml )        siobNstartTime     = UNSET_RL
659        close( scrunit1 )        siobSstartTime     = UNSET_RL
660          siobEstartTime     = UNSET_RL
661  c     Complete the start date specifications for the forcing        siobWstartTime     = UNSET_RL
662  c     fields to get a complete calendar date array.  
663    C     Initialise file type and field precision
664  c     check for consistency        exf_iprec            = 32
665          exf_iprec_obcs       = UNSET_I
666        if (.NOT.        exf_yftype           = 'RL'
667       &     (exf_iprec .EQ. 32 .OR. exf_iprec .EQ. 64)        useExfYearlyFields   = .FALSE.
668       &     ) then        twoDigitYear         = .FALSE.
669           stop 'stop in exf_readparms: value of exf_iprec not allowed'  
670        else if (.NOT.  C     Input scaling factors.
671       &        (exf_yftype .EQ. 'RS' .OR.        exf_inscal_hflux     =  1. _d 0
672       &        exf_yftype .EQ. 'RL')        exf_inscal_sflux     =  1. _d 0
673       &        ) then        exf_inscal_ustress   =  1. _d 0
674           stop 'stop in exf_readparms: value of exf_yftype not allowed'        exf_inscal_vstress   =  1. _d 0
675        end if        exf_inscal_uwind     =  1. _d 0
676          exf_inscal_vwind     =  1. _d 0
677  #ifdef ALLOW_RUNOFF        exf_inscal_wspeed    =  1. _d 0
678        call cal_FullDate(  runoffstartdate1,  runoffstartdate2,        exf_inscal_swflux    =  1. _d 0
679       &                    runoffstartdate,            mythid )        exf_inscal_lwflux    =  1. _d 0
680  #endif        exf_inscal_precip    =  1. _d 0
681          exf_inscal_snowprecip=  1. _d 0
682  #ifdef ALLOW_BULKFORMULAE        exf_inscal_sst       =  1. _d 0
683          exf_inscal_sss       =  1. _d 0
684  #ifdef ALLOW_ATM_TEMP        exf_inscal_atemp     =  1. _d 0
685        call cal_FullDate(   atempstartdate1,   atempstartdate2,        exf_offset_atemp     =  0. _d 0
686       &                     atempstartdate,            mythid )        exf_inscal_aqh       =  1. _d 0
687        call cal_FullDate(     aqhstartdate1,     aqhstartdate2,        exf_inscal_hs        =  1. _d 0
688       &                       aqhstartdate,            mythid )        exf_inscal_hl        =  1. _d 0
689        call cal_FullDate(  swfluxstartdate1,  swfluxstartdate2,        exf_inscal_evap      =  1. _d 0
690       &                    swfluxstartdate,            mythid )        exf_inscal_apressure =  1. _d 0
691        call cal_FullDate(  lwfluxstartdate1,  lwfluxstartdate2,        exf_inscal_runoff    =  1. _d 0
692       &                    lwfluxstartdate,            mythid )        exf_inscal_runoftemp =  1. _d 0
693        call cal_FullDate(  precipstartdate1,  precipstartdate2,        exf_inscal_saltflx   =  1. _d 0
694       &                    precipstartdate,            mythid )        exf_inscal_swdown    =  1. _d 0
695  #else        exf_inscal_lwdown    =  1. _d 0
696        call cal_FullDate(   hfluxstartdate1,   hfluxstartdate2,        exf_inscal_climsst   =  1. _d 0
697       &                     hfluxstartdate,            mythid )        exf_inscal_climsss   =  1. _d 0
698        call cal_FullDate(   sfluxstartdate1,  sfluxstartdate2,        exf_inscal_climustr  =  1. _d 0
699       &                     sfluxstartdate,           mythid )        exf_inscal_climvstr  =  1. _d 0
700  #ifdef ALLOW_KPP        exf_inscal_areamask  =  1. _d 0
701        call cal_FullDate(  swfluxstartdate1,  swfluxstartdate2,  
702       &                    swfluxstartdate,            mythid )  C     Output scaling factors.
703  #endif        exf_outscal_hflux    =  1. _d 0
704          exf_outscal_sflux    =  1. _d 0
705  #endif        exf_outscal_ustress  =  1. _d 0
706          exf_outscal_vstress  =  1. _d 0
707  #ifdef ALLOW_ATM_WIND        exf_outscal_swflux   =  1. _d 0
708        call cal_FullDate(   uwindstartdate1,   uwindstartdate2,        exf_outscal_sst      =  1. _d 0
709       &                     uwindstartdate,            mythid )        exf_outscal_sss      =  1. _d 0
710        call cal_FullDate(   vwindstartdate1,   vwindstartdate2,        exf_outscal_apressure=  1. _d 0
711       &                     vwindstartdate,            mythid )        exf_outscal_areamask =  1. _d 0
712  #else  
713        call cal_FullDate( ustressstartdate1, ustressstartdate2,  #ifdef USE_EXF_INTERPOLATION
714       &                   ustressstartdate,            mythid )  C--   set default input location to match (in case of simple Lat-Lonp grid)
715        call cal_FullDate( vstressstartdate1, vstressstartdate2,  C     model grid cell-center position (leading to trivial interpolation)
716       &                   vstressstartdate,            mythid )        inp_lon0 = xgOrigin + delX(1)*exf_half
717  #endif        inp_lat0 = ygOrigin + delY(1)*exf_half
718          inp_dLon = delX(1)
719  #else        inp_dLat = delY(1)
720        call cal_FullDate(   hfluxstartdate1,  hfluxstartdate2,  
721       &                     hfluxstartdate,           mythid )        ustress_lon0   = inp_lon0
722        call cal_FullDate(   sfluxstartdate1,  sfluxstartdate2,        uwind_lon0     = inp_lon0
723       &                     sfluxstartdate,           mythid )        vstress_lon0   = inp_lon0
724        call cal_FullDate( ustressstartdate1, ustressstartdate2,        hflux_lon0     = inp_lon0
725       &                   ustressstartdate,            mythid )        sflux_lon0     = inp_lon0
726        call cal_FullDate( vstressstartdate1, vstressstartdate2,        swflux_lon0    = inp_lon0
727       &                   vstressstartdate,            mythid )        runoff_lon0    = inp_lon0
728  #ifdef ALLOW_KPP        saltflx_lon0   = inp_lon0
729        call cal_FullDate(  swfluxstartdate1,  swfluxstartdate2,        atemp_lon0     = inp_lon0
730       &                    swfluxstartdate,            mythid )        aqh_lon0       = inp_lon0
731  #endif        hs_lon0        = inp_lon0
732          hl_lon0        = inp_lon0
733  #endif        evap_lon0      = inp_lon0
734          precip_lon0    = inp_lon0
735          snowprecip_lon0= inp_lon0
736          vwind_lon0     = inp_lon0
737          wspeed_lon0    = inp_lon0
738          lwflux_lon0    = inp_lon0
739          swdown_lon0    = inp_lon0
740          lwdown_lon0    = inp_lon0
741          apressure_lon0 = inp_lon0
742          areamask_lon0  = inp_lon0
743          vstress_lat0   = inp_lat0
744          vwind_lat0     = inp_lat0
745          wspeed_lat0    = inp_lat0
746          ustress_lat0   = inp_lat0
747          hflux_lat0     = inp_lat0
748          sflux_lat0     = inp_lat0
749          runoff_lat0    = inp_lat0
750          saltflx_lat0   = inp_lat0
751          swflux_lat0    = inp_lat0
752          atemp_lat0     = inp_lat0
753          aqh_lat0       = inp_lat0
754          hs_lat0        = inp_lat0
755          hl_lat0        = inp_lat0
756          evap_lat0      = inp_lat0
757          precip_lat0    = inp_lat0
758          snowprecip_lat0= inp_lat0
759          uwind_lat0     = inp_lat0
760          lwflux_lat0    = inp_lat0
761          swdown_lat0    = inp_lat0
762          lwdown_lat0    = inp_lat0
763          apressure_lat0 = inp_lat0
764          areamask_lat0  = inp_lat0
765          ustress_nlon   = gridNx
766          ustress_nlat   = gridNy
767          vstress_nlon   = gridNx
768          vstress_nlat   = gridNy
769          hflux_nlon     = gridNx
770          hflux_nlat     = gridNy
771          sflux_nlon     = gridNx
772          sflux_nlat     = gridNy
773          swflux_nlon    = gridNx
774          swflux_nlat    = gridNy
775          runoff_nlon    = gridNx
776          runoff_nlat    = gridNy
777          saltflx_nlon   = gridNx
778          saltflx_nlat   = gridNy
779          atemp_nlon     = gridNx
780          atemp_nlat     = gridNy
781          aqh_nlon       = gridNx
782          aqh_nlat       = gridNy
783          hs_nlon        = gridNx
784          hs_nlat        = gridNy
785          hl_nlon        = gridNx
786          hl_nlat        = gridNy
787          evap_nlon      = gridNx
788          evap_nlat      = gridNy
789          precip_nlon    = gridNx
790          precip_nlat    = gridNy
791          snowprecip_nlon= gridNx
792          snowprecip_nlat= gridNy
793          uwind_nlon     = gridNx
794          uwind_nlat     = gridNy
795          vwind_nlon     = gridNx
796          vwind_nlat     = gridNy
797          wspeed_nlon    = gridNx
798          wspeed_nlat    = gridNy
799          lwflux_nlon    = gridNx
800          lwflux_nlat    = gridNy
801          swdown_nlon    = gridNx
802          swdown_nlat    = gridNy
803          lwdown_nlon    = gridNx
804          lwdown_nlat    = gridNy
805          apressure_nlon = gridNx
806          apressure_nlat = gridNy
807          areamask_nlon  = gridNx
808          areamask_nlat  = gridNy
809          ustress_lon_inc   = inp_dLon
810          vstress_lon_inc   = inp_dLon
811          hflux_lon_inc     = inp_dLon
812          sflux_lon_inc     = inp_dLon
813          swflux_lon_inc    = inp_dLon
814          runoff_lon_inc    = inp_dLon
815          saltflx_lon_inc   = inp_dLon
816          atemp_lon_inc     = inp_dLon
817          aqh_lon_inc       = inp_dLon
818          hs_lon_inc        = inp_dLon
819          hl_lon_inc        = inp_dLon
820          evap_lon_inc      = inp_dLon
821          precip_lon_inc    = inp_dLon
822          snowprecip_lon_inc= inp_dLon
823          uwind_lon_inc     = inp_dLon
824          vwind_lon_inc     = inp_dLon
825          wspeed_lon_inc    = inp_dLon
826          lwflux_lon_inc    = inp_dLon
827          swdown_lon_inc    = inp_dLon
828          lwdown_lon_inc    = inp_dLon
829          apressure_lon_inc = inp_dLon
830          areamask_lon_inc  = inp_dLon
831    
832          climsst_lon0    = inp_lon0
833          climsss_lon0    = inp_lon0
834          climustr_lon0   = inp_lon0
835          climvstr_lon0   = inp_lon0
836          climsst_lat0    = inp_lat0
837          climsss_lat0    = inp_lat0
838          climustr_lat0   = inp_lat0
839          climvstr_lat0   = inp_lat0
840          climsst_nlon    = gridNx
841          climsst_nlat    = gridNy
842          climsss_nlon    = gridNx
843          climsss_nlat    = gridNy
844          climustr_nlon   = gridNx
845          climustr_nlat   = gridNy
846          climvstr_nlon   = gridNx
847          climvstr_nlat   = gridNy
848          climsst_lon_inc = inp_dLon
849          climsss_lon_inc = inp_dLon
850          climustr_lon_inc= inp_dLon
851          climvstr_lon_inc= inp_dLon
852    
853          DO j=1,MAX_LAT_INC
854            IF (j.LT.gridNy) THEN
855              inp_dLat = (delY(j) + delY(j+1))*exf_half
856            ELSE
857              inp_dLat = 0.
858            ENDIF
859            ustress_lat_inc(j)   = inp_dLat
860            vstress_lat_inc(j)   = inp_dLat
861            hflux_lat_inc(j)     = inp_dLat
862            sflux_lat_inc(j)     = inp_dLat
863            swflux_lat_inc(j)    = inp_dLat
864            runoff_lat_inc(j)    = inp_dLat
865            saltflx_lat_inc(j)   = inp_dLat
866            atemp_lat_inc(j)     = inp_dLat
867            aqh_lat_inc(j)       = inp_dLat
868            hs_lat_inc(j)        = inp_dLat
869            hl_lat_inc(j)        = inp_dLat
870            evap_lat_inc(j)      = inp_dLat
871            precip_lat_inc(j)    = inp_dLat
872            snowprecip_lat_inc(j)= inp_dLat
873            uwind_lat_inc(j)     = inp_dLat
874            vwind_lat_inc(j)     = inp_dLat
875            wspeed_lat_inc(j)    = inp_dLat
876            lwflux_lat_inc(j)    = inp_dLat
877            swdown_lat_inc(j)    = inp_dLat
878            lwdown_lat_inc(j)    = inp_dLat
879            apressure_lat_inc(j) = inp_dLat
880            areamask_lat_inc(j)  = inp_dLat
881            climsst_lat_inc(j)   = inp_dLat
882            climsss_lat_inc(j)   = inp_dLat
883            climustr_lat_inc(j)  = inp_dLat
884            climvstr_lat_inc(j)  = inp_dLat
885          ENDDO
886    
887          ustress_interpMethod   = 12
888          vstress_interpMethod   = 22
889          hflux_interpMethod     =  1
890          sflux_interpMethod     =  1
891          swflux_interpMethod    =  1
892          runoff_interpMethod    =  1
893          saltflx_interpMethod   =  1
894          atemp_interpMethod     =  1
895          aqh_interpMethod       =  1
896          hs_interpMethod        =  1
897          hl_interpMethod        =  1
898          evap_interpMethod      =  1
899          precip_interpMethod    =  1
900          snowprecip_interpMethod=  1
901          uwind_interpMethod     = 12
902          vwind_interpMethod     = 22
903          wspeed_interpMethod    =  1
904          lwflux_interpMethod    =  1
905          swdown_interpMethod    =  1
906          lwdown_interpMethod    =  1
907          apressure_interpMethod =  1
908          areamask_interpMethod  =  1
909          climsst_interpMethod   =  2
910          climsss_interpMethod   =  2
911          climustr_interpMethod  = 12
912          climvstr_interpMethod  = 22
913    
914          exf_output_interp = .FALSE.
915    #endif /* USE_EXF_INTERPOLATION */
916    
917    C--   Next, read pkg/exf parameter file.
918          WRITE(msgBuf,'(A)') 'EXF_READPARMS: opening data.exf'
919          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
920         &                    SQUEEZE_RIGHT, myThid )
921    
922          CALL OPEN_COPY_DATA_FILE(
923         I                          'data.exf', 'EXF_READPARMS',
924         O                          iUnit,
925         I                          myThid )
926    
927          WRITE(msgBuf,'(A)')
928         &     'EXF_READPARMS: reading EXF_NML_01'
929          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
930         &                    SQUEEZE_RIGHT, myThid )
931          READ(  iUnit, nml = EXF_NML_01 )
932    C-    Set default fldRepeatCycle to repeatPeriod
933          hfluxRepCycle      = repeatPeriod
934          atempRepCycle      = repeatPeriod
935          aqhRepCycle        = repeatPeriod
936          hs_RepCycle        = repeatPeriod
937          hl_RepCycle        = repeatPeriod
938          evapRepCycle       = repeatPeriod
939          precipRepCycle     = repeatPeriod
940          snowprecipRepCycle = repeatPeriod
941          sfluxRepCycle      = repeatPeriod
942          runoffRepCycle     = repeatPeriod
943          saltflxRepCycle    = repeatPeriod
944          ustressRepCycle    = repeatPeriod
945          vstressRepCycle    = repeatPeriod
946          uwindRepCycle      = repeatPeriod
947          vwindRepCycle      = repeatPeriod
948          wspeedRepCycle     = repeatPeriod
949          swfluxRepCycle     = repeatPeriod
950          lwfluxRepCycle     = repeatPeriod
951          swdownRepCycle     = repeatPeriod
952          lwdownRepCycle     = repeatPeriod
953          apressureRepCycle  = repeatPeriod
954          areamaskRepCycle   = repeatPeriod
955          climsstRepCycle    = repeatPeriod
956          climsssRepCycle    = repeatPeriod
957          climustrRepCycle   = repeatPeriod
958          climvstrRepCycle   = repeatPeriod
959    C-
960          obcsNrepCycle      = repeatPeriod
961          obcsSrepCycle      = repeatPeriod
962          obcsErepCycle      = repeatPeriod
963          obcsWrepCycle      = repeatPeriod
964          siobNrepCycle      = UNSET_RL
965          siobSrepCycle      = UNSET_RL
966          siobErepCycle      = UNSET_RL
967          siobWrepCycle      = UNSET_RL
968    
969          WRITE(msgBuf,'(A)')
970         &     'EXF_READPARMS: reading EXF_NML_02'
971          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
972         &                    SQUEEZE_RIGHT, myThid )
973          READ(  iUnit, nml = EXF_NML_02 )
974          WRITE(msgBuf,'(A)')
975         &     'EXF_READPARMS: reading EXF_NML_03'
976          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
977         &                    SQUEEZE_RIGHT, myThid )
978          READ(  iUnit, nml = EXF_NML_03 )
979    #ifdef USE_EXF_INTERPOLATION
980          WRITE(msgBuf,'(A)')
981         &     'EXF_READPARMS: reading EXF_NML_04'
982          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
983         &                    SQUEEZE_RIGHT, myThid )
984          READ(  iUnit, nml = EXF_NML_04 )
985    #endif /* USE_EXF_INTERPOLATION */
986    
987  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
988  #ifdef ALLOW_OBCS_NORTH        IF ( useOBCS ) THEN
989        call cal_FullDate(  obcsNstartdate1,  obcsNstartdate2,         WRITE(msgBuf,'(A)')
990       &                    obcsNstartdate,           mythid )       &      'EXF_READPARMS: reading EXF_NML_OBCS'
991  #endif         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
992  #ifdef ALLOW_OBCS_SOUTH       &                     SQUEEZE_RIGHT, myThid )
993        call cal_FullDate(  obcsSstartdate1,  obcsSstartdate2,         READ(  iUnit, nml = EXF_NML_OBCS )
994       &                    obcsSstartdate,           mythid )        ENDIF
995  #endif         IF(siobNstartdate1.EQ.UNSET_I ) siobNstartdate1=obcsNstartdate1
996  #ifdef ALLOW_OBCS_EAST         IF(siobNstartdate2.EQ.UNSET_I ) siobNstartdate2=obcsNstartdate2
997        call cal_FullDate(  obcsEstartdate1,  obcsEstartdate2,         IF(siobNperiod    .EQ.UNSET_RL) siobNperiod    =obcsNperiod
998       &                    obcsEstartdate,           mythid )         IF(siobNrepCycle  .EQ.UNSET_RL) siobNrepCycle  =obcsNrepCycle
999           IF(siobSstartdate1.EQ.UNSET_I ) siobSstartdate1=obcsSstartdate1
1000           IF(siobSstartdate2.EQ.UNSET_I ) siobSstartdate2=obcsSstartdate2
1001           IF(siobSperiod    .EQ.UNSET_RL) siobSperiod    =obcsSperiod
1002           IF(siobSrepCycle  .EQ.UNSET_RL) siobSrepCycle  =obcsSrepCycle
1003           IF(siobEstartdate1.EQ.UNSET_I ) siobEstartdate1=obcsEstartdate1
1004           IF(siobEstartdate2.EQ.UNSET_I ) siobEstartdate2=obcsEstartdate2
1005           IF(siobEperiod    .EQ.UNSET_RL) siobEperiod    =obcsEperiod
1006           IF(siobErepCycle  .EQ.UNSET_RL) siobErepCycle  =obcsErepCycle
1007           IF(siobWstartdate1.EQ.UNSET_I ) siobWstartdate1=obcsWstartdate1
1008           IF(siobWstartdate2.EQ.UNSET_I ) siobWstartdate2=obcsWstartdate2
1009           IF(siobWperiod    .EQ.UNSET_RL) siobWperiod    =obcsWperiod
1010           IF(siobWrepCycle  .EQ.UNSET_RL) siobWrepCycle  =obcsWrepCycle
1011    
1012           IF(exf_iprec_obcs .EQ. UNSET_I) exf_iprec_obcs =exf_iprec
1013    #endif /* ALLOW_OBCS */
1014    
1015          WRITE(msgBuf,'(A)')
1016         &     'EXF_READPARMS: finished reading data.exf'
1017          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1018         &                    SQUEEZE_RIGHT, myThid )
1019    
1020          CLOSE( iUnit )
1021    
1022    C--   Retired parameters
1023          IF ( exf_yftype.NE.'RL' ) THEN
1024           STOP 'S/R EXF_READPARAMS: value of exf_yftype not allowed'
1025          ENDIF
1026    
1027    C--   Derive other parameters:
1028          IF ( exf_verbose ) exf_debugLev = MAX( exf_debugLev, debLevD )
1029          hq = ht
1030          stressIsOnCgrid = readStressOnCgrid
1031    #if ( defined (ALLOW_BULKFORMULAE) )
1032          IF ( useAtmWind ) stressIsOnCgrid = .FALSE.
1033  #endif  #endif
1034  #ifdef ALLOW_OBCS_WEST  #ifdef USE_EXF_INTERPOLATION
1035        call cal_FullDate(  obcsWstartdate1,  obcsWstartdate2,        IF ( (ustress_interpMethod.GE.1 .AND. ustressfile.NE.' ') .OR.
1036       &                    obcsWstartdate,           mythid )       &     (vstress_interpMethod.GE.1 .AND. vstressfile.NE.' ') )
1037         &   stressIsOnCgrid = .FALSE.
1038    #endif /* USE_EXF_INTERPOLATION */
1039    
1040          useExfZenAlbedo = select_ZenAlbedo.GE.1
1041         &            .AND. select_ZenAlbedo.LE.3
1042    
1043    C--   Overwrite tauThetaClimRelax but stop if already set.
1044    C- Note: need this, even if EXF option ALLOW_CLIMSST_RELAXATION is undef;
1045    C        this prevents to apply relaxation towards potentially wrong SST since,
1046    C        with EXF, we skip the update of loaded SST in EXTERNAL_FIELDS_LOAD.
1047    C- Note2: let s see whether we can put this back under ifdef
1048    C        ALLOW_CLIMSST_RELAXATION, but always call EXTERNAL_FIELDS_LOAD.
1049    C        If ALLOW_CLIMSST_RELAXATION is undef, clim.relaxation could still
1050    C        be done outside of exf.
1051    #ifdef ALLOW_CLIMSST_RELAXATION
1052          IF ( tauThetaClimRelax.NE.0. _d 0 ) THEN
1053            WRITE(msgBuf,'(2A)') 'EXF_READPARMS: ',
1054         &   'with EXF, cannot use "tauThetaClimRelax" in "data"'
1055            CALL PRINT_ERROR( msgBuf, myThid )
1056            WRITE(msgBuf,'(2A)') 'since SST relax. is handled by EXF',
1057         &   ' (data.exf, "climsstTauRelax")'
1058            CALL PRINT_ERROR( msgBuf, myThid )
1059            STOP 'ABNORMAL END: S/R EXF_READPARMS'
1060          ENDIF
1061          tauThetaClimRelax = climsstTauRelax
1062  #endif  #endif
1063    
1064    #ifdef ALLOW_CLIMSSS_RELAXATION
1065    C--   Overwrite tauSaltClimRelax but stop if already set.
1066          IF ( tauSaltClimRelax.NE.0. _d 0 ) THEN
1067            WRITE(msgBuf,'(2A)') 'EXF_READPARMS: ',
1068         &   'with EXF, cannot use "tauSaltClimRelax" in "data"'
1069            CALL PRINT_ERROR( msgBuf, myThid )
1070            WRITE(msgBuf,'(2A)') 'since SSS relax. is handled by EXF',
1071         &   ' (data.exf, "climsssTauRelax")'
1072            CALL PRINT_ERROR( msgBuf, myThid )
1073            STOP 'ABNORMAL END: S/R EXF_READPARMS'
1074          ENDIF
1075          tauSaltClimRelax = climsssTauRelax
1076  #endif  #endif
1077    
1078        _END_MASTER( mythid )  C     Complete the start date specifications for the forcing
1079    C     fields to get a complete calendar date array.
1080    C     => moved to EXF_INIT_FIXED
1081    
1082          _END_MASTER( myThid )
1083        _BARRIER        _BARRIER
1084    
1085  c--   Summarize the External forcing's setup.        RETURN
1086        call exf_summary( mythid )        END
   
   
 c--   set climatology parameters  
       call exf_clim_readparms( mythid )  
   
 c--   summarize climatologic forcing configuration  
       call exf_clim_summary( mythid )  
   
       end  
Legend:
Removed from v.1.1.6.2  
changed lines
  Added in v.1.92
  ViewVC Help
Powered by ViewVC 1.1.22