C $Header: /u/gcmpack/MITgcm/pkg/gchem/gchem_calc_tendency.F,v 1.6 2017/12/29 19:40:06 jmc Exp $ C $Name: $ #include "GCHEM_OPTIONS.h" CBOP C !ROUTINE: GCHEM_CALC_TENDENCY C !INTERFACE: ========================================================== SUBROUTINE GCHEM_CALC_TENDENCY( I myTime, myIter, myThid ) C !DESCRIPTION: C In the case of GCHEM_SEPARATE_FORCING not defined, C this subroutine computes the tracer tendencies due to a C bio-geogchemistry or ecosystem model and stores them on an array C gchemTendency, that will be incorporated into regular timestepping in C in ptracers_intergrate.F C The current example uses the CFC package, but his is meant to C be replaced by anything that the user provides. C !USES: =============================================================== IMPLICIT NONE #include "SIZE.h" #include "GRID.h" #include "DYNVARS.h" #include "EEPARAMS.h" #include "PARAMS.h" #include "PTRACERS_SIZE.h" #include "PTRACERS_PARAMS.h" #include "PTRACERS_FIELDS.h" #include "GCHEM.h" #ifdef ALLOW_CFC # include "CFC_SIZE.h" #endif #ifdef ALLOW_SPOIL # include "SPOIL_SIZE.h" #endif #include "GCHEM_SIZE.h" #include "GCHEM_FIELDS.h" C !INPUT PARAMETERS: =================================================== C myTime :: current time C myIter :: current timestep C myThid :: thread number _RL myTime INTEGER myIter, myThid CEOP #ifdef GCHEM_ADD2TR_TENDENCY C !LOCAL VARIABLES: ==================================================== C i,j :: loop indices C iTr :: ptracer number C jTr :: corresponding tracer tendency number C bi,bj :: tile indices C k :: vertical level INTEGER i,j INTEGER iTr, jTr INTEGER bi,bj,iMin,iMax,jMin,jMax,k C gchemTendency is re-initialized here DO jTr = 1, GCHEM_tendTr_num DO bj=myByLo(myThid),myByHi(myThid) DO bi=myBxLo(myThid),myBxHi(myThid) DO k = 1, Nr DO j = 1-OLy, sNy+OLy DO i = 1-OLx, sNx+OLx gchemTendency(i,j,k,bi,bj,jTr) = 0. _d 0 ENDDO ENDDO ENDDO ENDDO ENDDO ENDDO C C Here is the place for code to compute bio-geochemical C tendency terms (sometimes referred to as source-minus-sink C terms). The tendencies are stored on gchemTendency, as show C in the CFC-example. C C loop over tiles DO bj=myByLo(myThid),myByHi(myThid) DO bi=myBxLo(myThid),myBxHi(myThid) C define horizontal loop ranges iMin = 1 iMax = sNx jMin = 1 jMax = sNy c DO jTr = 1, GCHEM_tendTr_num c iTr = jTr + gchem_sepFTr_num c DO k = 1, Nr c DO j = 1-OLy, sNy+OLy c DO i = 1-OLx, sNx+OLx c gchemTendency(i,j,k,bi,bj,jTr) = your specific model c ENDDO c ENDDO c ENDDO c ENDDO C------------------------ C chemical forcing | C------------------------ #ifdef ALLOW_CFC IF ( useCFC ) THEN iTr = CFC_pTr_i1 jTr = iTr - gchem_sepFTr_num CALL CFC11_FORCING( I pTracer (1-OLx,1-OLy,1,bi,bj,iTr), U gchemTendency(1-OLx,1-OLy,1,bi,bj,jTr), I bi, bj, iMin, iMax, jMin, jMax, I myTime, myIter, myThid ) iTr = CFC_pTr_i1 + 1 jTr = iTr - gchem_sepFTr_num CALL CFC12_FORCING( I pTracer (1-OLx,1-OLy,1,bi,bj,iTr), U gchemTendency(1-OLx,1-OLy,1,bi,bj,jTr), I bi, bj, iMin, iMax, jMin, jMax, I myTime, myIter, myThid ) ENDIF #endif /* ALLOW_CFC */ C end of tile-loops ENDDO ENDDO #endif /* GCHEM_ADD2TR_TENDENCY */ RETURN END