--- MITgcm/pkg/seaice/seaice_tracer_phys.F 2012/01/13 14:19:25 1.4 +++ MITgcm/pkg/seaice/seaice_tracer_phys.F 2012/02/03 13:34:32 1.5 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/seaice/seaice_tracer_phys.F,v 1.4 2012/01/13 14:19:25 jmc Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/seaice/seaice_tracer_phys.F,v 1.5 2012/02/03 13:34:32 gforget Exp $ C $Name: $ #include "SEAICE_OPTIONS.h" @@ -245,9 +245,7 @@ DIAGarray(I,J,1+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr) DIAGarray(I,J,2+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr)*HEFFpost 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,2) - elseif (SItrName(iTr).EQ.'salinity') then + if (SItrName(iTr).EQ.'salinity') then DIAGarray(I,J,3+(iTr-1)*5) = HSALT(i,j,bi,bj)/SEAICE_rhoIce elseif (SItrName(iTr).EQ.'one') then DIAGarray(I,J,3+(iTr-1)*5) = HEFFpost @@ -262,10 +260,6 @@ AREApost=SItrAREA(i,j,bi,bj,3) DIAGarray(I,J,1+(iTr-1)*5) = SItracer(i,j,bi,bj,iTr) 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 ENDDO ENDDO endif