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

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

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

revision 1.44 by heimbach, Sat Mar 5 18:06:06 2011 UTC revision 1.45 by jmc, Wed Apr 20 01:34:27 2011 UTC
# Line 62  C     afx       :: horizontal advective Line 62  C     afx       :: horizontal advective
62  C     afy       :: horizontal advective flux, y direction  C     afy       :: horizontal advective flux, y direction
63  C     gFld      :: tendency of seaice field  C     gFld      :: tendency of seaice field
64  C     xA,yA     :: "areas" of X and Y face of tracer cells  C     xA,yA     :: "areas" of X and Y face of tracer cells
65        INTEGER i, j, bi, bj, iTracer        INTEGER i, j, bi, bj
66        INTEGER ks        INTEGER ks
67        LOGICAL SEAICEmultiDimAdvection        LOGICAL SEAICEmultiDimAdvection
68    #ifdef SEAICE_AGE
69          INTEGER iTracer
70    #endif
71    
72  C-    MPI+MTH: apply exch (sure with exch1) only to array in common block  C-    MPI+MTH: apply exch (sure with exch1) only to array in common block
73        COMMON / SEAICE_ADVDIFF_LOCAL / uc, vc        COMMON / SEAICE_ADVDIFF_LOCAL / uc, vc
# Line 299  C-- Line 302  C--
302            CALL SEAICE_ADVECTION(            CALL SEAICE_ADVECTION(
303       I       GAD_ICEAGETR+iTracer-1, SEAICEadvSchAge,       I       GAD_ICEAGETR+iTracer-1, SEAICEadvSchAge,
304       I       uc(1-OLx,1-OLy,bi,bj), vc(1-OLx,1-OLy,bi,bj),       I       uc(1-OLx,1-OLy,bi,bj), vc(1-OLx,1-OLy,bi,bj),
305       I       uTrans, vTrans, IceAgeTr(1-OLx,1-OLy,bi,bj,iTracer),       I       uTrans, vTrans, IceAgeTr(1-OLx,1-OLy,bi,bj,iTracer),
306       I       recip_heff,       I       recip_heff,
307       O       gFld, afx, afy,       O       gFld, afx, afy,
308       I       bi, bj, myTime, myIter, myThid )       I       bi, bj, myTime, myIter, myThid )
# Line 315  C     now do the "explicit" time step Line 318  C     now do the "explicit" time step
318            DO j=1,sNy            DO j=1,sNy
319             DO i=1,sNx             DO i=1,sNx
320              IceAgeTr(i,j,bi,bj,iTracer) = HEFFM(i,j,bi,bj) * (              IceAgeTr(i,j,bi,bj,iTracer) = HEFFM(i,j,bi,bj) * (
321       &          IceAgeTr(i,j,bi,bj,iTracer)       &          IceAgeTr(i,j,bi,bj,iTracer)
322       &          + SEAICE_deltaTtherm * gFld(i,j)       &          + SEAICE_deltaTtherm * gFld(i,j)
323       &          )       &          )
324             ENDDO             ENDDO
# Line 356  CADJ STORE vc   = comlev1, key = ikey_dy Line 359  CADJ STORE vc   = comlev1, key = ikey_dy
359  #ifdef SEAICE_AGE  #ifdef SEAICE_AGE
360         IF ( SEAICEadvAge ) THEN         IF ( SEAICEadvAge ) THEN
361          DO iTracer = 1, SEAICE_num          DO iTracer = 1, SEAICE_num
362           CALL ADVECT( uc, vc, iceAgeTr(1-Olx,1-Oly,1,1,iTracer),           CALL ADVECT( uc, vc, iceAgeTr(1-Olx,1-Oly,1,1,iTracer),
363       &    fldNm1, HEFFM, myThid )       &    fldNm1, HEFFM, myThid )
364          ENDDO          ENDDO
365         ENDIF         ENDIF
Legend:
Removed from v.1.44  
changed lines
  Added in v.1.45
  ViewVC Help
Powered by ViewVC 1.1.22