pkg/cfc/cfc11_surfforcing.F

8ac664a04c Step*0001 #include "GCHEM_OPTIONS.h"
                0002
6cebf9fc9a Jean*0003 CBOP
                0004 C     !ROUTINE: CFC11_SURFFORCING
                0005 C     !INTERFACE:
                0006       SUBROUTINE CFC11_SURFFORCING(
                0007      I                          pTr_CFC11, AtmosCFC11,
                0008      O                          fluxCFC11,
                0009      I                          bi, bj, iMin, iMax, jMin, jMax,
                0010      I                          myTime, myIter, myThid )
                0011
                0012 C     !DESCRIPTION:
                0013 C     *==========================================================*
                0014 C     | SUBROUTINE CFC11_SURFFORCING
                0015 C     | o Calculate the CFC11 air-sea fluxes
                0016 C     *==========================================================*
                0017
                0018 C     !USES:
8ac664a04c Step*0019       IMPLICIT NONE
                0020
                0021 C     == GLobal variables ==
                0022 #include "SIZE.h"
                0023 #include "EEPARAMS.h"
                0024 #include "PARAMS.h"
                0025 #include "GRID.h"
3247fa4f0b Jean*0026 #include "DYNVARS.h"
8ac664a04c Step*0027 #include "CFC.h"
                0028
6cebf9fc9a Jean*0029 C     !INPUT/OUTPUT PARAMETERS:
                0030 C     pTr_CFC11  :: ocean CFC11 concentration
                0031 C     AtmosCFC11 :: atmospheric CFC11 field
                0032 C     fluxCFC11  :: air-sea CFC11 fluxes
                0033 C     bi, bj     :: current tile indices
                0034 C     iMin,iMax  :: computation domain, 1rst index bounds
                0035 C     jMin,jMax  :: computation domain, 2nd  index bounds
                0036 C     myTime     :: current time in simulation
                0037 C     myIter     :: current iteration number
                0038 C     myThid     :: my Thread Id number
                0039       _RL  pTr_CFC11 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
                0040       _RL  AtmosCFC11(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0041       _RL  fluxCFC11 (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0042       INTEGER bi, bj
                0043       INTEGER iMin, iMax, jMin, jMax
8ac664a04c Step*0044       _RL myTime
6cebf9fc9a Jean*0045       INTEGER myIter, myThid
                0046 CEOP
8ac664a04c Step*0047
                0048 #ifdef ALLOW_PTRACERS
                0049 #ifdef ALLOW_CFC
                0050
6cebf9fc9a Jean*0051 C     !LOCAL VARIABLES:
                0052 C i, j, k - Loop counters
                0053       INTEGER i,j,k
8ac664a04c Step*0054 C Solubility relation coefficients
                0055       _RL SchmidtNocfc11(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0056       _RL SolCFC11(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
3247fa4f0b Jean*0057 c     _RL cfc11sat(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
8ac664a04c Step*0058       _RL Kwexch(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0059       _RL Csat(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0060       _RL ttemp, ttemp2
                0061
6cebf9fc9a Jean*0062       k=1
8ac664a04c Step*0063
                0064 C calculate SCHMIDT NO. for O2
                0065         DO j=jMin,jMax
6cebf9fc9a Jean*0066          DO i=iMin,iMax
                0067            IF ( maskC(i,j,k,bi,bj).NE.zeroRS ) THEN
8ac664a04c Step*0068 C calculate SCHMIDT NO. for CFC11
6cebf9fc9a Jean*0069              SchmidtNocfc11(i,j) =
                0070      &            sca_11_1
8ac664a04c Step*0071      &          + sca_11_2 * theta(i,j,k,bi,bj)
6cebf9fc9a Jean*0072      &          + sca_11_3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
                0073      &          + sca_11_4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
8ac664a04c Step*0074      &                *theta(i,j,k,bi,bj)
                0075 c calculate solubility for CFC11
6cebf9fc9a Jean*0076              ttemp=( theta(i,j,k,bi,bj) + 273.16 _d 0)* 0.01 _d 0
                0077              ttemp2=( B3_11 * ttemp + B2_11 )*ttemp + B1_11
                0078              SolCFC11(i,j) =
                0079      &           EXP ( A1_11
8ac664a04c Step*0080      &           +  A2_11 / ttemp
                0081      &           +  A3_11 * log( ttemp )
                0082      &           +  A4_11 * ttemp * ttemp
                0083      &           +  Salt(i,j,k,bi,bj)* ttemp2 )
                0084 c     conversion from mol/(l * atm) to mol/(m^3 * atm)
6cebf9fc9a Jean*0085              SolCFC11(i,j) = 1000. _d 0 * SolCFC11(i,j)
8ac664a04c Step*0086 c     conversion from mol/(m^3 * atm) to mol/(m3 * pptv)
6cebf9fc9a Jean*0087              SolCFC11(i,j) = 1. _d -12 * SolCFC11(i,j)
8ac664a04c Step*0088 C Determine surface flux (Fcfc11)
6cebf9fc9a Jean*0089              Csat(i,j) = SolCFC11(i,j)*AtmosP(i,j,bi,bj)
                0090      &                     *AtmosCFC11(i,j)
                0091              Kwexch(i,j) = (1. _d 0 - fIce(i,j,bi,bj))
                0092      &                   * pisVel(i,j,bi,bj)
                0093      &                   / SQRT(SchmidtNoCFC11(i,j)/660. _d 0)
                0094              fluxCFC11(i,j) =
                0095      &           Kwexch(i,j)*(Csat(i,j) - pTr_CFC11(i,j,1))
8ac664a04c Step*0096
6cebf9fc9a Jean*0097            ELSE
                0098              fluxCFC11(i,j) = 0. _d 0
                0099            ENDIF
8ac664a04c Step*0100
                0101          ENDDO
6cebf9fc9a Jean*0102         ENDDO
8ac664a04c Step*0103
6cebf9fc9a Jean*0104 #endif /* ALLOW_CFC */
                0105 #endif /* ALLOW_PTRACERS */
8ac664a04c Step*0106
6cebf9fc9a Jean*0107       RETURN
                0108       END