/[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.8 by gforget, Mon Mar 5 14:52:16 2012 UTC revision 1.10 by jmc, Sun Aug 11 02:31:12 2013 UTC
# Line 15  C     *================================= Line 15  C     *=================================
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"
# Line 173  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 187  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    
202        if ( (SItrName(iTr).EQ.'salinity').AND.        if ( (SItrName(iTr).EQ.'salinity').AND.
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22