173 |
c "ice concentration" tracer: no specific process |
c "ice concentration" tracer: no specific process |
174 |
elseif (SItrName(iTr).EQ.'ridge') then |
elseif (SItrName(iTr).EQ.'ridge') then |
175 |
c simple, made up, ice surface roughness index prototype |
c simple, made up, ice surface roughness index prototype |
|
#ifndef SEAICE_GROWTH_LEGACY |
|
176 |
DO J=1,sNy |
DO J=1,sNy |
177 |
DO I=1,sNx |
DO I=1,sNx |
178 |
c ridging increases roughness |
c ridging increases roughness |
186 |
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 |
187 |
ENDDO |
ENDDO |
188 |
ENDDO |
ENDDO |
|
#endif |
|
189 |
endif |
endif |
190 |
c 3) ice-ocean tracer exchange/mapping to external variables |
c 3) ice-ocean tracer exchange/mapping to external variables |
191 |
c ========================================================== |
c ========================================================== |
192 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
193 |
if (SItrMate(iTr).EQ.'HEFF') then |
if (SItrMate(iTr).EQ.'HEFF') then |
194 |
WRITE(diagName,'(A4,I2.2,A2)') 'SItr',iTr,'FX' |
WRITE(diagName,'(A4,I2.2,A2)') 'SItr',iTr,'Fx' |
195 |
tmpscal1=-ONE/SEAICE_deltaTtherm*SEAICE_rhoIce |
tmpscal1=-ONE/SEAICE_deltaTtherm*SEAICE_rhoIce |
196 |
CALL DIAGNOSTICS_SCALE_FILL(SItrBucket(1-OLx,1-OLy,bi,bj,iTr), |
CALL DIAGNOSTICS_SCALE_FILL(SItrBucket(1-OLx,1-OLy,bi,bj,iTr), |
197 |
& tmpscal1, 1, diagName,0,1,2,bi,bj,myThid) |
& tmpscal1, 1, diagName,0,1,2,bi,bj,myThid) |