/[MITgcm]/MITgcm/pkg/seaice/seaice_tracer_phys.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_tracer_phys.F

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

revision 1.4 by jmc, Fri Jan 13 14:19:25 2012 UTC revision 1.10 by jmc, Sun Aug 11 02:31:12 2013 UTC
# Line 5  C $Name$ Line 5  C $Name$
5    
6  CStartOfInterface  CStartOfInterface
7        SUBROUTINE SEAICE_TRACER_PHYS( myTime, myIter, myThid )        SUBROUTINE SEAICE_TRACER_PHYS( myTime, myIter, myThid )
8  C     /=======================================================\  C     *=======================================================*
9  C     | SUBROUTINE seaice_tracer_phys                         |  C     | SUBROUTINE seaice_tracer_phys
10  C     | o Time step SItr/SItrEFF as a result of               |  C     | o Time step SItr/SItrEFF as a result of
11  C     |   seaice thermodynamics and specific tracer physics   |  C     |   seaice thermodynamics and specific tracer physics
12  C     \=======================================================/  C     *=======================================================*
13        IMPLICIT NONE        IMPLICIT NONE
14    
15  C     === Global variables ===  C     === Global variables ===
16  #include "SIZE.h"  #include "SIZE.h"
17  #include "EEPARAMS.h"  #include "EEPARAMS.h"
18    #include "PARAMS.h"
19  #include "FFIELDS.h"  #include "FFIELDS.h"
20  #include "DYNVARS.h"  #include "DYNVARS.h"
21  #include "SEAICE_SIZE.h"  #include "SEAICE_SIZE.h"
22  #include "SEAICE.h"  #include "SEAICE.h"
23  #include "SEAICE_PARAMS.h"  #include "SEAICE_PARAMS.h"
24  #include "SEAICE_TRACER.h"  #include "SEAICE_TRACER.h"
25    #ifdef ALLOW_SALT_PLUME
26    # include "SALT_PLUME.h"
27    #endif
28    
29  C     === Routine arguments ===  C     === Routine arguments ===
30  C     INPUT:  C     INPUT:
# Line 55  cgf eventually I will need to handle the Line 59  cgf eventually I will need to handle the
59    
60        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
61        DO bi=myBxLo(myThid),myBxHi(myThid)        DO bi=myBxLo(myThid),myBxHi(myThid)
62        DO iTr=1,SItrMaxNum        DO iTr=1,SItrNumInUse
63    
64  c 0) set ice-ocean and ice-snow exchange values  c 0) set ice-ocean and ice-snow exchange values
65  c =============================================  c =============================================
66        DO J=1,sNy        DO J=1,sNy
67         DO I=1,sNx         DO I=1,sNx
68          SItrFromOcean(i,j)=0. _d 0          SItrFromOcean(i,j)=SItrFromOcean0(iTr)
69          SItrFromFlood(i,j)=0. _d 0          SItrFromFlood(i,j)=SItrFromFlood0(iTr)
70          SItrExpand(i,j)=0. _d 0          SItrExpand(i,j)=SItrExpand0(iTr)
71         ENDDO         ENDDO
72        ENDDO        ENDDO
       if (SItrName(iTr).EQ.'age') then  
 c age tracer: no age in ocean, or effect from ice cover changes  
       elseif (SItrName(iTr).EQ.'salinity') then  
73  c salinity tracer:  c salinity tracer:
74          if ( (SItrName(iTr).EQ.'salinity').AND.
75         &      (SItrFromOceanFrac(iTr).GT.ZERO) ) then
76         DO J=1,sNy         DO J=1,sNy
77          DO I=1,sNx          DO I=1,sNx
78           SItrFromOcean(i,j)=SIsal0           SItrFromOcean(i,j)=SItrFromOceanFrac(iTr)*salt(I,j,ks,bi,bj)
79  #ifdef SEAICE_VARIABLE_SALINITY           SItrFromFlood(i,j)=SItrFromFloodFrac(iTr)*salt(I,j,ks,bi,bj)
          if (SIsalFRAC.GT.0.)  
      &   SItrFromOcean(i,j)=SIsalFRAC*salt(I,j,ks,bi,bj)  
 #endif  
 c as of now, flooding implies no salt extraction from ocean  
         ENDDO  
        ENDDO  
       elseif (SItrName(iTr).EQ.'one') then  
 c "ice concentration" tracer that should remain .EQ.1.  
        DO J=1,sNy  
         DO I=1,sNx  
          SItrFromOcean(i,j)=1. _d 0  
          SItrFromFlood(i,j)=1. _d 0  
80          ENDDO          ENDDO
81         ENDDO         ENDDO
82        endif        endif
# Line 183  c salinity tracer: no specific process Line 174  c salinity tracer: no specific process
174  c "ice concentration" tracer: no specific process  c "ice concentration" tracer: no specific process
175        elseif (SItrName(iTr).EQ.'ridge') then        elseif (SItrName(iTr).EQ.'ridge') then
176  c simple, made up, ice surface roughness index prototype  c simple, made up, ice surface roughness index prototype
 #ifndef SEAICE_GROWTH_LEGACY  
177         DO J=1,sNy         DO J=1,sNy
178          DO I=1,sNx          DO I=1,sNx
179  c ridging increases roughness  c ridging increases roughness
# Line 197  c ice melt reduces ridges/roughness Line 187  c ice melt reduces ridges/roughness
187            SItracer(i,j,bi,bj,iTr)=SItracer(i,j,bi,bj,iTr)*tmpscal1            SItracer(i,j,bi,bj,iTr)=SItracer(i,j,bi,bj,iTr)*tmpscal1
188          ENDDO          ENDDO
189         ENDDO         ENDDO
 #endif  
190        endif        endif
191  c 3) ice-ocean tracer exchange/mapping to external variables  c 3) ice-ocean tracer exchange/mapping to external variables
192  c ==========================================================  c ==========================================================
193  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
194        if (SItrMate(iTr).EQ.'HEFF') then        IF ( useDiagnostics .AND. SItrMate(iTr).EQ.'HEFF') THEN
195          WRITE(diagName,'(A4,I2.2,A2)') 'SItr',iTr,'FX'          WRITE(diagName,'(A4,I2.2,A2)') 'SItr',iTr,'Fx'
196          tmpscal1=-ONE/SEAICE_deltaTtherm*SEAICE_rhoIce          tmpscal1=-ONE/SEAICE_deltaTtherm*SEAICE_rhoIce
197          CALL DIAGNOSTICS_SCALE_FILL(SItrBucket(1-oLx,1-oLy,bi,bj,iTr),          CALL DIAGNOSTICS_SCALE_FILL(SItrBucket(1-OLx,1-OLy,bi,bj,iTr),
198       &   tmpscal1, 1, diagName,0,1,2,bi,bj,myThid)       &   tmpscal1, 1, diagName,0,1,2,bi,bj,myThid)
199        endif        ENDIF
200  #endif  #endif
201        if (SItrName(iTr).EQ.'age') then  
202  c age tracer: not passed to ocean        if ( (SItrName(iTr).EQ.'salinity').AND.
203        elseif (SItrName(iTr).EQ.'salinity') then       &     (SEAICE_salinityTracer) ) then
204  c salinity tracer: salt flux  c salinity tracer: salt flux
205  c      DO J=1,sNy          DO J=1,sNy
206  c       DO I=1,sNx           DO I=1,sNx
207  c        saltFlux(I,J,bi,bj) = - SItrBucket(i,j,bi,bj,iTr)            saltFlux(I,J,bi,bj) = - SItrBucket(i,j,bi,bj,iTr)
208  c    &     *HEFFM(I,J,bi,bj)/SEAICE_deltaTtherm*SEAICE_rhoIce       &      *HEFFM(I,J,bi,bj)/SEAICE_deltaTtherm*SEAICE_rhoIce
209  c note: at this point of the time step, that is the correct sign  c note: at this point of the time step, that is the correct sign
210  c        saltPlumeFlux(I,J,bi,bj) = ...  #ifdef ALLOW_SALT_PLUME
211  c       ENDDO  c should work for both constant and variable ice salinity -- to be tested
212  c      ENDDO            saltPlumeFlux(I,J,bi,bj) = MAX(ZERO,saltFlux(I,J,bi,bj))
213        elseif (SItrName(iTr).EQ.'one') then       &      *SPsalFRAC*(salt(I,j,ks,bi,bj)-SItrFromOcean(i,j))
214  c "ice concentration" tracer: not passed to ocean  #endif
215             ENDDO
216            ENDDO
217        endif        endif
218    
219        DO J=1,sNy        DO J=1,sNy
220         DO I=1,sNx         DO I=1,sNx
221  #ifdef ALLOW_SITRACER_DEBUG_DIAG  #ifdef ALLOW_SITRACER_DEBUG_DIAG
# Line 234  c empty bucket Line 226  c empty bucket
226          SItrBucket(i,j,bi,bj,iTr)=0. _d 0          SItrBucket(i,j,bi,bj,iTr)=0. _d 0
227         ENDDO         ENDDO
228        ENDDO        ENDDO
229    
230  c TAF? elseif (SItrMate(iTr).EQ.'AREA') then  c TAF? elseif (SItrMate(iTr).EQ.'AREA') then
231    
232  c 4) diagnostics  c 4) diagnostics
233  c ==============  c ==============
234  #ifdef ALLOW_SITRACER_DEBUG_DIAG  #ifdef ALLOW_SITRACER_DEBUG_DIAG
# Line 245  c ============== Line 239  c ==============
239          DIAGarray(I,J,1+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)          DIAGarray(I,J,1+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)
240          DIAGarray(I,J,2+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)*HEFFpost          DIAGarray(I,J,2+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)*HEFFpost
241  c DIAGarray(:,:,3) is the term of comparison for DIAGarray(:,:,2)  c DIAGarray(:,:,3) is the term of comparison for DIAGarray(:,:,2)
242          if (SItrName(iTr).EQ.'age') then          if (SItrName(iTr).EQ.'salinity') then
           DIAGarray(I,J,3+(iTr-1)*5) = IceAgeTr(i,j,bi,bj,2)  
         elseif (SItrName(iTr).EQ.'salinity') then  
243            DIAGarray(I,J,3+(iTr-1)*5) = HSALT(i,j,bi,bj)/SEAICE_rhoIce            DIAGarray(I,J,3+(iTr-1)*5) = HSALT(i,j,bi,bj)/SEAICE_rhoIce
244          elseif (SItrName(iTr).EQ.'one') then          elseif (SItrName(iTr).EQ.'one') then
245            DIAGarray(I,J,3+(iTr-1)*5) = HEFFpost            DIAGarray(I,J,3+(iTr-1)*5) = HEFFpost
# Line 262  c DIAGarray(:,:,5) is the tracer flux fr Line 254  c DIAGarray(:,:,5) is the tracer flux fr
254          AREApost=SItrAREA(i,j,bi,bj,3)          AREApost=SItrAREA(i,j,bi,bj,3)
255          DIAGarray(I,J,1+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)          DIAGarray(I,J,1+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)
256          DIAGarray(I,J,2+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)*AREApost          DIAGarray(I,J,2+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)*AREApost
 c DIAGarray(:,:,3) is the term of comparison for DIAGarray(:,:,2)  
         if (SItrName(iTr).EQ.'age') then  
           DIAGarray(I,J,3+(iTr-1)*5) = IceAgeTr(i,j,bi,bj,1)  
         endif  
257         ENDDO         ENDDO
258        ENDDO        ENDDO
259        endif        endif
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22