pkg/darwin/darwin_exf_init_varia.F

8fbfd1f382 Oliv*0001 #include "DARWIN_OPTIONS.h"
                0002 #include "EXF_OPTIONS.h"
                0003
                0004 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0005 CBOP
                0006 C     !ROUTINE: DARWIN_EXF_INIT_VARIA
                0007
                0008 C     !INTERFACE:
                0009       SUBROUTINE DARWIN_EXF_INIT_VARIA( myThid )
                0010
                0011 C     !DESCRIPTION:
                0012 C     Initialize DARWIN constants
                0013
                0014 C     !USES:
                0015       IMPLICIT NONE
                0016 #include "EEPARAMS.h"
                0017 #include "SIZE.h"
                0018 #include "GRID.h"
                0019 #include "PARAMS.h"
                0020 #ifdef ALLOW_EXF
d9e107d180 Oliv*0021 #include "EXF_INTERP_SIZE.h"
8fbfd1f382 Oliv*0022 #endif
                0023 #include "DARWIN_SIZE.h"
                0024 #include "DARWIN_EXF_PARAMS.h"
                0025 #include "DARWIN_EXF_FIELDS.h"
                0026
                0027 C     !INPUT PARAMETERS:
                0028       INTEGER myThid
                0029 CEOP
                0030
                0031 #ifdef ALLOW_DARWIN
                0032 #ifdef ALLOW_EXF
                0033
                0034 #ifdef ALLOW_RADTRANS
                0035       INTEGER l
                0036 #endif
                0037
                0038 C     these load constant forcing fields
                0039
fecf0f6104 Oliv*0040       call exf_init_fld (
                0041      &     'PAR', PARfile, PARmask,
                0042      &     PARperiod, darwin_inscal_PAR, PARconst,
                0043      &     surfPAR, PAR0, PAR1,
8fbfd1f382 Oliv*0044 #ifdef USE_EXF_INTERPOLATION
                0045      &     PAR_lon0, PAR_lon_inc,
                0046      &     PAR_lat0, PAR_lat_inc,
6ed0cf4850 Oliv*0047      &     PAR_nlon, PAR_nlat, xC, yC,
                0048      &     PAR_interpMethod,
8fbfd1f382 Oliv*0049 #endif
                0050      &     mythid )
                0051
fecf0f6104 Oliv*0052       call exf_init_fld (
                0053      &     'iron', ironfile, ironmask,
                0054      &     ironperiod, darwin_inscal_iron, ironconst,
                0055      &     inputFe, iron0, iron1,
8fbfd1f382 Oliv*0056 #ifdef USE_EXF_INTERPOLATION
                0057      &     iron_lon0, iron_lon_inc,
                0058      &     iron_lat0, iron_lat_inc,
6ed0cf4850 Oliv*0059      &     iron_nlon, iron_nlat, xC, yC,
                0060      &     iron_interpMethod,
8fbfd1f382 Oliv*0061 #endif
                0062      &     mythid )
                0063
fecf0f6104 Oliv*0064       call exf_init_fld (
                0065      &     'ice', icefile, icemask,
                0066      &     iceperiod, darwin_inscal_ice, iceconst,
                0067      &     iceFrac, ice0, ice1,
8fbfd1f382 Oliv*0068 #ifdef USE_EXF_INTERPOLATION
                0069      &     ice_lon0, ice_lon_inc,
                0070      &     ice_lat0, ice_lat_inc,
6ed0cf4850 Oliv*0071      &     ice_nlon, ice_nlat, xC, yC,
                0072      &     ice_interpMethod,
8fbfd1f382 Oliv*0073 #endif
                0074      &     mythid )
                0075
fecf0f6104 Oliv*0076       call exf_init_fld (
                0077      &     'wind', windfile, windmask,
                0078      &     windperiod, darwin_inscal_wind, windconst,
                0079      &     windSpeed, wind0, wind1,
8fbfd1f382 Oliv*0080 #ifdef USE_EXF_INTERPOLATION
                0081      &     wind_lon0, wind_lon_inc,
                0082      &     wind_lat0, wind_lat_inc,
6ed0cf4850 Oliv*0083      &     wind_nlon, wind_nlat, xC, yC,
                0084      &     wind_interpMethod,
8fbfd1f382 Oliv*0085 #endif
                0086      &     mythid )
                0087
fecf0f6104 Oliv*0088       call exf_init_fld (
                0089      &     'pCO2', pCO2file, pCO2mask,
                0090      &     pCO2period, darwin_inscal_pCO2, pCO2const,
                0091      &     atmospCO2, pCO20, pCO21,
8fbfd1f382 Oliv*0092 #ifdef USE_EXF_INTERPOLATION
                0093      &     pCO2_lon0, pCO2_lon_inc,
                0094      &     pCO2_lat0, pCO2_lat_inc,
6ed0cf4850 Oliv*0095      &     pCO2_nlon, pCO2_nlat, xC, yC,
                0096      &     pCO2_interpMethod,
                0097 #endif
                0098      &     mythid )
                0099
                0100       call exf_init_fld (
                0101      &     'ventHe3', ventHe3file, ventHe3mask,
                0102      &     ventHe3period, darwin_inscal_ventHe3, ventHe3const,
                0103      &     ventHe3, ventHe30, ventHe31,
                0104 #ifdef USE_EXF_INTERPOLATION
                0105      &     ventHe3_lon0, ventHe3_lon_inc,
                0106      &     ventHe3_lat0, ventHe3_lat_inc,
                0107      &     ventHe3_nlon, ventHe3_nlat, xC, yC,
                0108      &     ventHe3_interpMethod,
8fbfd1f382 Oliv*0109 #endif
                0110      &     mythid )
                0111
e79020cc6e Oliv*0112       call exf_init_fld (
                0113      &     'DOCrunoff', DOCrunofffile, DOCrunoffmask,
6ed0cf4850 Oliv*0114      &     DOCrunoffperiod, darwin_inscal_DOCrunoff, DOCrunoffconst,
e79020cc6e Oliv*0115      &     DOCrunoff, DOCrunoff0, DOCrunoff1,
                0116 #ifdef USE_EXF_INTERPOLATION
                0117      &     DOCrunoff_lon0, DOCrunoff_lon_inc,
                0118      &     DOCrunoff_lat0, DOCrunoff_lat_inc,
                0119      &     DOCrunoff_nlon, DOCrunoff_nlat, xC, yC,
                0120      &     DOCrunoff_interpMethod,
                0121 #endif
                0122      &     mythid )
                0123
                0124       call exf_init_fld (
                0125      &     'DONrunoff', DONrunofffile, DONrunoffmask,
6ed0cf4850 Oliv*0126      &     DONrunoffperiod, darwin_inscal_DONrunoff, DONrunoffconst,
e79020cc6e Oliv*0127      &     DONrunoff, DONrunoff0, DONrunoff1,
                0128 #ifdef USE_EXF_INTERPOLATION
                0129      &     DONrunoff_lon0, DONrunoff_lon_inc,
                0130      &     DONrunoff_lat0, DONrunoff_lat_inc,
                0131      &     DONrunoff_nlon, DONrunoff_nlat, xC, yC,
                0132      &     DONrunoff_interpMethod,
                0133 #endif
                0134      &     mythid )
                0135
                0136       call exf_init_fld (
6ed0cf4850 Oliv*0137      &     'DOPrunoff', DOPrunofffile, DOPrunoffmask,
                0138      &     DOPrunoffperiod, darwin_inscal_DOPrunoff, DOPrunoffconst,
e79020cc6e Oliv*0139      &     DOPrunoff, DOPrunoff0, DOPrunoff1,
                0140 #ifdef USE_EXF_INTERPOLATION
                0141      &     DOPrunoff_lon0, DOPrunoff_lon_inc,
                0142      &     DOPrunoff_lat0, DOPrunoff_lat_inc,
                0143      &     DOPrunoff_nlon, DOPrunoff_nlat, xC, yC,
                0144      &     DOPrunoff_interpMethod,
                0145 #endif
                0146      &     mythid )
                0147
                0148       call exf_init_fld (
                0149      &     'DINrunoff', DINrunofffile, DINrunoffmask,
6ed0cf4850 Oliv*0150      &     DINrunoffperiod, darwin_inscal_DINrunoff, DINrunoffconst,
e79020cc6e Oliv*0151      &     DINrunoff, DINrunoff0, DINrunoff1,
                0152 #ifdef USE_EXF_INTERPOLATION
                0153      &     DINrunoff_lon0, DINrunoff_lon_inc,
                0154      &     DINrunoff_lat0, DINrunoff_lat_inc,
                0155      &     DINrunoff_nlon, DINrunoff_nlat, xC, yC,
                0156      &     DINrunoff_interpMethod,
                0157 #endif
                0158      &     mythid )
                0159
11d365acff Oliv*0160       call exf_init_fld (
                0161      &     'NO3runoff', NO3runofffile, NO3runoffmask,
                0162      &     NO3runoffperiod, darwin_inscal_NO3runoff, NO3runoffconst,
                0163      &     NO3runoff, NO3runoff0, NO3runoff1,
                0164 #ifdef USE_EXF_INTERPOLATION
                0165      &     NO3runoff_lon0, NO3runoff_lon_inc,
                0166      &     NO3runoff_lat0, NO3runoff_lat_inc,
                0167      &     NO3runoff_nlon, NO3runoff_nlat, xC, yC,
                0168      &     NO3runoff_interpMethod,
                0169 #endif
                0170      &     mythid )
                0171
                0172       call exf_init_fld (
                0173      &     'NO2runoff', NO2runofffile, NO2runoffmask,
                0174      &     NO2runoffperiod, darwin_inscal_NO2runoff, NO2runoffconst,
                0175      &     NO2runoff, NO2runoff0, NO2runoff1,
                0176 #ifdef USE_EXF_INTERPOLATION
                0177      &     NO2runoff_lon0, NO2runoff_lon_inc,
                0178      &     NO2runoff_lat0, NO2runoff_lat_inc,
                0179      &     NO2runoff_nlon, NO2runoff_nlat, xC, yC,
                0180      &     NO2runoff_interpMethod,
                0181 #endif
                0182      &     mythid )
                0183
                0184       call exf_init_fld (
                0185      &     'NH4runoff', NH4runofffile, NH4runoffmask,
                0186      &     NH4runoffperiod, darwin_inscal_NH4runoff, NH4runoffconst,
                0187      &     NH4runoff, NH4runoff0, NH4runoff1,
                0188 #ifdef USE_EXF_INTERPOLATION
                0189      &     NH4runoff_lon0, NH4runoff_lon_inc,
                0190      &     NH4runoff_lat0, NH4runoff_lat_inc,
                0191      &     NH4runoff_nlon, NH4runoff_nlat, xC, yC,
                0192      &     NH4runoff_interpMethod,
                0193 #endif
                0194      &     mythid )
                0195
e79020cc6e Oliv*0196       call exf_init_fld (
                0197      &     'IPrunoff', IPrunofffile, IPrunoffmask,
6ed0cf4850 Oliv*0198      &     IPrunoffperiod, darwin_inscal_IPrunoff, IPrunoffconst,
e79020cc6e Oliv*0199      &     IPrunoff, IPrunoff0, IPrunoff1,
                0200 #ifdef USE_EXF_INTERPOLATION
                0201      &     IPrunoff_lon0, IPrunoff_lon_inc,
                0202      &     IPrunoff_lat0, IPrunoff_lat_inc,
                0203      &     IPrunoff_nlon, IPrunoff_nlat, xC, yC,
                0204      &     IPrunoff_interpMethod,
                0205 #endif
                0206      &     mythid )
                0207
                0208       call exf_init_fld (
                0209      &     'DSirunoff', DSirunofffile, DSirunoffmask,
6ed0cf4850 Oliv*0210      &     DSirunoffperiod, darwin_inscal_DSirunoff, DSirunoffconst,
e79020cc6e Oliv*0211      &     DSirunoff, DSirunoff0, DSirunoff1,
                0212 #ifdef USE_EXF_INTERPOLATION
                0213      &     DSirunoff_lon0, DSirunoff_lon_inc,
                0214      &     DSirunoff_lat0, DSirunoff_lat_inc,
                0215      &     DSirunoff_nlon, DSirunoff_nlat, xC, yC,
                0216      &     DSirunoff_interpMethod,
                0217 #endif
                0218      &     mythid )
                0219
                0220       call exf_init_fld (
                0221      &     'POCrunoff', POCrunofffile, POCrunoffmask,
6ed0cf4850 Oliv*0222      &     POCrunoffperiod, darwin_inscal_POCrunoff, POCrunoffconst,
e79020cc6e Oliv*0223      &     POCrunoff, POCrunoff0, POCrunoff1,
                0224 #ifdef USE_EXF_INTERPOLATION
                0225      &     POCrunoff_lon0, POCrunoff_lon_inc,
                0226      &     POCrunoff_lat0, POCrunoff_lat_inc,
                0227      &     POCrunoff_nlon, POCrunoff_nlat, xC, yC,
                0228      &     POCrunoff_interpMethod,
                0229 #endif
                0230      &     mythid )
                0231
                0232       call exf_init_fld (
                0233      &     'POPrunoff', POPrunofffile, POPrunoffmask,
6ed0cf4850 Oliv*0234      &     POPrunoffperiod, darwin_inscal_POPrunoff, POPrunoffconst,
e79020cc6e Oliv*0235      &     POPrunoff, POPrunoff0, POPrunoff1,
                0236 #ifdef USE_EXF_INTERPOLATION
                0237      &     POPrunoff_lon0, POPrunoff_lon_inc,
                0238      &     POPrunoff_lat0, POPrunoff_lat_inc,
                0239      &     POPrunoff_nlon, POPrunoff_nlat, xC, yC,
                0240      &     POPrunoff_interpMethod,
                0241 #endif
                0242      &     mythid )
                0243
                0244       call exf_init_fld (
                0245      &     'PONrunoff', PONrunofffile, PONrunoffmask,
6ed0cf4850 Oliv*0246      &     PONrunoffperiod, darwin_inscal_PONrunoff, PONrunoffconst,
e79020cc6e Oliv*0247      &     PONrunoff, PONrunoff0, PONrunoff1,
                0248 #ifdef USE_EXF_INTERPOLATION
                0249      &     PONrunoff_lon0, PONrunoff_lon_inc,
                0250      &     PONrunoff_lat0, PONrunoff_lat_inc,
                0251      &     PONrunoff_nlon, PONrunoff_nlat, xC, yC,
                0252      &     PONrunoff_interpMethod,
                0253 #endif
                0254      &     mythid )
                0255
                0256       call exf_init_fld (
                0257      &     'DICrunoff', DICrunofffile, DICrunoffmask,
6ed0cf4850 Oliv*0258      &     DICrunoffperiod, darwin_inscal_DICrunoff, DICrunoffconst,
e79020cc6e Oliv*0259      &     DICrunoff, DICrunoff0, DICrunoff1,
                0260 #ifdef USE_EXF_INTERPOLATION
                0261      &     DICrunoff_lon0, DICrunoff_lon_inc,
                0262      &     DICrunoff_lat0, DICrunoff_lat_inc,
                0263      &     DICrunoff_nlon, DICrunoff_nlat, xC, yC,
                0264      &     DICrunoff_interpMethod,
                0265 #endif
                0266      &     mythid )
                0267
6ed0cf4850 Oliv*0268
8fbfd1f382 Oliv*0269       IF (darwin_loadFieldsEarly) THEN
                0270 C       these load time-varying forcing fields
                0271 C       (usually done in darwin_fields_load, but may be needed for chl)
fecf0f6104 Oliv*0272         call exf_set_fld(
                0273      &     'PAR', PARfile, PARmask,
                0274      &     PARStartTime, PARperiod, PARRepCycle,
8fbfd1f382 Oliv*0275      &     darwin_inscal_PAR,
                0276      &     PAR_exfremo_intercept, PAR_exfremo_slope,
fecf0f6104 Oliv*0277      &     surfPAR, PAR0, PAR1,
8fbfd1f382 Oliv*0278 #ifdef USE_EXF_INTERPOLATION
                0279      &     PAR_lon0, PAR_lon_inc,
                0280      &     PAR_lat0, PAR_lat_inc,
6ed0cf4850 Oliv*0281      &     PAR_nlon, PAR_nlat, xC, yC,
                0282      &     PAR_interpMethod,
8fbfd1f382 Oliv*0283 #endif
                0284      &     starttime, nIter0, mythid )
                0285
fecf0f6104 Oliv*0286         call exf_set_fld(
                0287      &     'iron', ironfile, ironmask,
                0288      &     ironStartTime, ironperiod, ironRepCycle,
8fbfd1f382 Oliv*0289      &     darwin_inscal_iron,
                0290      &     iron_exfremo_intercept, iron_exfremo_slope,
fecf0f6104 Oliv*0291      &     inputFe, iron0, iron1,
8fbfd1f382 Oliv*0292 #ifdef USE_EXF_INTERPOLATION
                0293      &     iron_lon0, iron_lon_inc,
                0294      &     iron_lat0, iron_lat_inc,
6ed0cf4850 Oliv*0295      &     iron_nlon, iron_nlat, xC, yC,
                0296      &     iron_interpMethod,
8fbfd1f382 Oliv*0297 #endif
                0298      &     starttime, nIter0, mythid )
                0299
fecf0f6104 Oliv*0300         call exf_set_fld(
                0301      &     'ice', icefile, icemask,
                0302      &     iceStartTime, iceperiod, iceRepCycle,
8fbfd1f382 Oliv*0303      &     darwin_inscal_ice,
                0304      &     ice_exfremo_intercept, ice_exfremo_slope,
fecf0f6104 Oliv*0305      &     iceFrac, ice0, ice1,
8fbfd1f382 Oliv*0306 #ifdef USE_EXF_INTERPOLATION
                0307      &     ice_lon0, ice_lon_inc,
                0308      &     ice_lat0, ice_lat_inc,
6ed0cf4850 Oliv*0309      &     ice_nlon, ice_nlat, xC, yC,
                0310      &     ice_interpMethod,
8fbfd1f382 Oliv*0311 #endif
                0312      &     starttime, nIter0, mythid )
                0313
fecf0f6104 Oliv*0314         call exf_set_fld(
                0315      &     'wind', windfile, windmask,
                0316      &     windStartTime, windperiod, windRepCycle,
8fbfd1f382 Oliv*0317      &     darwin_inscal_wind,
                0318      &     wind_exfremo_intercept, wind_exfremo_slope,
fecf0f6104 Oliv*0319      &     windSpeed, wind0, wind1,
8fbfd1f382 Oliv*0320 #ifdef USE_EXF_INTERPOLATION
                0321      &     wind_lon0, wind_lon_inc,
                0322      &     wind_lat0, wind_lat_inc,
6ed0cf4850 Oliv*0323      &     wind_nlon, wind_nlat, xC, yC,
                0324      &     wind_interpMethod,
8fbfd1f382 Oliv*0325 #endif
                0326      &     starttime, nIter0, mythid )
                0327
fecf0f6104 Oliv*0328         call exf_set_fld(
                0329      &     'pCO2', pCO2file, pCO2mask,
                0330      &     pCO2StartTime, pCO2period, pCO2RepCycle,
8fbfd1f382 Oliv*0331      &     darwin_inscal_pCO2,
                0332      &     pCO2_exfremo_intercept, pCO2_exfremo_slope,
fecf0f6104 Oliv*0333      &     atmospCO2, pCO20, pCO21,
8fbfd1f382 Oliv*0334 #ifdef USE_EXF_INTERPOLATION
                0335      &     pCO2_lon0, pCO2_lon_inc,
                0336      &     pCO2_lat0, pCO2_lat_inc,
6ed0cf4850 Oliv*0337      &     pCO2_nlon, pCO2_nlat, xC, yC,
                0338      &     pCO2_interpMethod,
                0339 #endif
                0340      &     starttime, nIter0, mythid )
                0341
                0342         call exf_set_fld(
                0343      &     'ventHe3', ventHe3file, ventHe3mask,
                0344      &     ventHe3StartTime, ventHe3period, ventHe3RepCycle,
                0345      &     darwin_inscal_ventHe3,
                0346      &     ventHe3_exfremo_intercept, ventHe3_exfremo_slope,
                0347      &     ventHe3, ventHe30, ventHe31,
                0348 #ifdef USE_EXF_INTERPOLATION
                0349      &     ventHe3_lon0, ventHe3_lon_inc,
                0350      &     ventHe3_lat0, ventHe3_lat_inc,
                0351      &     ventHe3_nlon, ventHe3_nlat, xC, yC,
                0352      &     ventHe3_interpMethod,
8fbfd1f382 Oliv*0353 #endif
                0354      &     starttime, nIter0, mythid )
                0355
e79020cc6e Oliv*0356         call exf_set_fld(
                0357      &     'DOCrunoff', DOCrunofffile, DOCrunoffmask,
                0358      &     DOCrunoffStartTime, DOCrunoffperiod, DOCrunoffRepCycle,
                0359      &     darwin_inscal_DOCrunoff,
                0360      &     DOCrunoff_exfremo_intercept, DOCrunoff_exfremo_slope,
                0361      &     DOCrunoff, DOCrunoff0, DOCrunoff1,
                0362 #ifdef USE_EXF_INTERPOLATION
                0363      &     DOCrunoff_lon0, DOCrunoff_lon_inc,
                0364      &     DOCrunoff_lat0, DOCrunoff_lat_inc,
                0365      &     DOCrunoff_nlon, DOCrunoff_nlat, xC, yC,
                0366      &     DOCrunoff_interpMethod,
                0367 #endif
                0368      &     starttime, nIter0, mythid )
                0369
                0370         call exf_set_fld(
                0371      &     'DONrunoff', DONrunofffile, DONrunoffmask,
                0372      &     DONrunoffStartTime, DONrunoffperiod, DONrunoffRepCycle,
                0373      &     darwin_inscal_DONrunoff,
                0374      &     DONrunoff_exfremo_intercept, DONrunoff_exfremo_slope,
                0375      &     DONrunoff, DONrunoff0, DONrunoff1,
                0376 #ifdef USE_EXF_INTERPOLATION
                0377      &     DONrunoff_lon0, DONrunoff_lon_inc,
                0378      &     DONrunoff_lat0, DONrunoff_lat_inc,
                0379      &     DONrunoff_nlon, DONrunoff_nlat, xC, yC,
                0380      &     DONrunoff_interpMethod,
                0381 #endif
                0382      &     starttime, nIter0, mythid )
                0383
                0384         call exf_set_fld(
                0385      &     'DOPrunoff', DOPrunofffile, DOPrunoffmask,
                0386      &     DOPrunoffStartTime, DOPrunoffperiod, DOPrunoffRepCycle,
                0387      &     darwin_inscal_DOPrunoff,
                0388      &     DOPrunoff_exfremo_intercept, DOPrunoff_exfremo_slope,
                0389      &     DOPrunoff, DOPrunoff0, DOPrunoff1,
                0390 #ifdef USE_EXF_INTERPOLATION
                0391      &     DOPrunoff_lon0, DOPrunoff_lon_inc,
                0392      &     DOPrunoff_lat0, DOPrunoff_lat_inc,
                0393      &     DOPrunoff_nlon, DOPrunoff_nlat, xC, yC,
                0394      &     DOPrunoff_interpMethod,
                0395 #endif
                0396      &     starttime, nIter0, mythid )
                0397
                0398         call exf_set_fld(
                0399      &     'DINrunoff', DINrunofffile, DINrunoffmask,
                0400      &     DINrunoffStartTime, DINrunoffperiod, DINrunoffRepCycle,
                0401      &     darwin_inscal_DINrunoff,
                0402      &     DINrunoff_exfremo_intercept, DINrunoff_exfremo_slope,
                0403      &     DINrunoff, DINrunoff0, DINrunoff1,
                0404 #ifdef USE_EXF_INTERPOLATION
                0405      &     DINrunoff_lon0, DINrunoff_lon_inc,
                0406      &     DINrunoff_lat0, DINrunoff_lat_inc,
                0407      &     DINrunoff_nlon, DINrunoff_nlat, xC, yC,
                0408      &     DINrunoff_interpMethod,
                0409 #endif
                0410      &     starttime, nIter0, mythid )
                0411
11d365acff Oliv*0412         call exf_set_fld(
                0413      &     'NO3runoff', NO3runofffile, NO3runoffmask,
                0414      &     NO3runoffStartTime, NO3runoffperiod, NO3runoffRepCycle,
                0415      &     darwin_inscal_NO3runoff,
                0416      &     NO3runoff_exfremo_intercept, NO3runoff_exfremo_slope,
                0417      &     NO3runoff, NO3runoff0, NO3runoff1,
                0418 #ifdef USE_EXF_INTERPOLATION
                0419      &     NO3runoff_lon0, NO3runoff_lon_inc,
                0420      &     NO3runoff_lat0, NO3runoff_lat_inc,
                0421      &     NO3runoff_nlon, NO3runoff_nlat, xC, yC,
                0422      &     NO3runoff_interpMethod,
                0423 #endif
                0424      &     starttime, nIter0, mythid )
                0425
                0426         call exf_set_fld(
                0427      &     'NO2runoff', NO2runofffile, NO2runoffmask,
                0428      &     NO2runoffStartTime, NO2runoffperiod, NO2runoffRepCycle,
                0429      &     darwin_inscal_NO2runoff,
                0430      &     NO2runoff_exfremo_intercept, NO2runoff_exfremo_slope,
                0431      &     NO2runoff, NO2runoff0, NO2runoff1,
                0432 #ifdef USE_EXF_INTERPOLATION
                0433      &     NO2runoff_lon0, NO2runoff_lon_inc,
                0434      &     NO2runoff_lat0, NO2runoff_lat_inc,
                0435      &     NO2runoff_nlon, NO2runoff_nlat, xC, yC,
                0436      &     NO2runoff_interpMethod,
                0437 #endif
                0438      &     starttime, nIter0, mythid )
                0439
                0440         call exf_set_fld(
                0441      &     'NH4runoff', NH4runofffile, NH4runoffmask,
                0442      &     NH4runoffStartTime, NH4runoffperiod, NH4runoffRepCycle,
                0443      &     darwin_inscal_NH4runoff,
                0444      &     NH4runoff_exfremo_intercept, NH4runoff_exfremo_slope,
                0445      &     NH4runoff, NH4runoff0, NH4runoff1,
                0446 #ifdef USE_EXF_INTERPOLATION
                0447      &     NH4runoff_lon0, NH4runoff_lon_inc,
                0448      &     NH4runoff_lat0, NH4runoff_lat_inc,
                0449      &     NH4runoff_nlon, NH4runoff_nlat, xC, yC,
                0450      &     NH4runoff_interpMethod,
                0451 #endif
                0452      &     starttime, nIter0, mythid )
                0453
e79020cc6e Oliv*0454         call exf_set_fld(
                0455      &     'IPrunoff', IPrunofffile, IPrunoffmask,
                0456      &     IPrunoffStartTime, IPrunoffperiod, IPrunoffRepCycle,
                0457      &     darwin_inscal_IPrunoff,
                0458      &     IPrunoff_exfremo_intercept, IPrunoff_exfremo_slope,
                0459      &     IPrunoff, IPrunoff0, IPrunoff1,
                0460 #ifdef USE_EXF_INTERPOLATION
                0461      &     IPrunoff_lon0, IPrunoff_lon_inc,
                0462      &     IPrunoff_lat0, IPrunoff_lat_inc,
                0463      &     IPrunoff_nlon, IPrunoff_nlat, xC, yC,
                0464      &     IPrunoff_interpMethod,
                0465 #endif
                0466      &     starttime, nIter0, mythid )
                0467
                0468         call exf_set_fld(
                0469      &     'DSirunoff', DSirunofffile, DSirunoffmask,
                0470      &     DSirunoffStartTime, DSirunoffperiod, DSirunoffRepCycle,
                0471      &     darwin_inscal_DSirunoff,
                0472      &     DSirunoff_exfremo_intercept, DSirunoff_exfremo_slope,
                0473      &     DSirunoff, DSirunoff0, DSirunoff1,
                0474 #ifdef USE_EXF_INTERPOLATION
                0475      &     DSirunoff_lon0, DSirunoff_lon_inc,
                0476      &     DSirunoff_lat0, DSirunoff_lat_inc,
                0477      &     DSirunoff_nlon, DSirunoff_nlat, xC, yC,
                0478      &     DSirunoff_interpMethod,
                0479 #endif
                0480      &     starttime, nIter0, mythid )
                0481
                0482         call exf_set_fld(
                0483      &     'POCrunoff', POCrunofffile, POCrunoffmask,
                0484      &     POCrunoffStartTime, POCrunoffperiod, POCrunoffRepCycle,
                0485      &     darwin_inscal_POCrunoff,
                0486      &     POCrunoff_exfremo_intercept, POCrunoff_exfremo_slope,
                0487      &     POCrunoff, POCrunoff0, POCrunoff1,
                0488 #ifdef USE_EXF_INTERPOLATION
                0489      &     POCrunoff_lon0, POCrunoff_lon_inc,
                0490      &     POCrunoff_lat0, POCrunoff_lat_inc,
                0491      &     POCrunoff_nlon, POCrunoff_nlat, xC, yC,
                0492      &     POCrunoff_interpMethod,
                0493 #endif
                0494      &     starttime, nIter0, mythid )
                0495
                0496         call exf_set_fld(
                0497      &     'POPrunoff', POPrunofffile, POPrunoffmask,
                0498      &     POPrunoffStartTime, POPrunoffperiod, POPrunoffRepCycle,
                0499      &     darwin_inscal_POPrunoff,
                0500      &     POPrunoff_exfremo_intercept, POPrunoff_exfremo_slope,
                0501      &     POPrunoff, POPrunoff0, POPrunoff1,
                0502 #ifdef USE_EXF_INTERPOLATION
                0503      &     POPrunoff_lon0, POPrunoff_lon_inc,
                0504      &     POPrunoff_lat0, POPrunoff_lat_inc,
                0505      &     POPrunoff_nlon, POPrunoff_nlat, xC, yC,
                0506      &     POPrunoff_interpMethod,
                0507 #endif
                0508      &     starttime, nIter0, mythid )
                0509
                0510         call exf_set_fld(
                0511      &     'PONrunoff', PONrunofffile, PONrunoffmask,
                0512      &     PONrunoffStartTime, PONrunoffperiod, PONrunoffRepCycle,
                0513      &     darwin_inscal_PONrunoff,
                0514      &     PONrunoff_exfremo_intercept, PONrunoff_exfremo_slope,
                0515      &     PONrunoff, PONrunoff0, PONrunoff1,
                0516 #ifdef USE_EXF_INTERPOLATION
                0517      &     PONrunoff_lon0, PONrunoff_lon_inc,
                0518      &     PONrunoff_lat0, PONrunoff_lat_inc,
                0519      &     PONrunoff_nlon, PONrunoff_nlat, xC, yC,
                0520      &     PONrunoff_interpMethod,
                0521 #endif
                0522      &     starttime, nIter0, mythid )
                0523
                0524         call exf_set_fld(
                0525      &     'DICrunoff', DICrunofffile, DICrunoffmask,
                0526      &     DICrunoffStartTime, DICrunoffperiod, DICrunoffRepCycle,
                0527      &     darwin_inscal_DICrunoff,
                0528      &     DICrunoff_exfremo_intercept, DICrunoff_exfremo_slope,
                0529      &     DICrunoff, DICrunoff0, DICrunoff1,
                0530 #ifdef USE_EXF_INTERPOLATION
                0531      &     DICrunoff_lon0, DICrunoff_lon_inc,
                0532      &     DICrunoff_lat0, DICrunoff_lat_inc,
                0533      &     DICrunoff_nlon, DICrunoff_nlat, xC, yC,
                0534      &     DICrunoff_interpMethod,
                0535 #endif
                0536      &     starttime, nIter0, mythid )
                0537
8fbfd1f382 Oliv*0538         _EXCH_XY_RL( surfPAR, myThid )
                0539         _EXCH_XY_RL( inputFe, myThid )
                0540         _EXCH_XY_RL( iceFrac, myThid )
                0541         _EXCH_XY_RL( windSpeed, myThid )
                0542         _EXCH_XY_RL( atmospCO2, myThid )
6ed0cf4850 Oliv*0543         _EXCH_XY_RL( ventHe3, myThid )
e79020cc6e Oliv*0544         _EXCH_XY_RL( DOCrunoff, myThid )
                0545         _EXCH_XY_RL( DONrunoff, myThid )
                0546         _EXCH_XY_RL( DOPrunoff, myThid )
                0547         _EXCH_XY_RL( DINrunoff, myThid )
11d365acff Oliv*0548         _EXCH_XY_RL( NO3runoff, myThid )
                0549         _EXCH_XY_RL( NO2runoff, myThid )
                0550         _EXCH_XY_RL( NH4runoff, myThid )
e79020cc6e Oliv*0551         _EXCH_XY_RL( IPrunoff, myThid )
                0552         _EXCH_XY_RL( DSirunoff, myThid )
                0553         _EXCH_XY_RL( POCrunoff, myThid )
                0554         _EXCH_XY_RL( POPrunoff, myThid )
                0555         _EXCH_XY_RL( PONrunoff, myThid )
                0556         _EXCH_XY_RL( DICrunoff, myThid )
                0557
8fbfd1f382 Oliv*0558       ENDIF
                0559
                0560 #endif /* ALLOW_EXF */
                0561 #endif /* ALLOW_DARWIN */
                0562
                0563       RETURN
                0564       END