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

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

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

revision 1.1 by mlosch, Thu Feb 16 10:41:48 2006 UTC revision 1.2 by mlosch, Wed Mar 15 19:49:04 2006 UTC
# Line 48  CADJ STORE heff  = comlev1, key = ikey_d Line 48  CADJ STORE heff  = comlev1, key = ikey_d
48  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
49    
50    
51    C--   This would be the natural way to do diffusion (explicitly)
52    C     For now we stick to the modified Eulerian time step
53    CML      DO bj=myByLo(myThid),myByHi(myThid)
54    CML       DO bi=myBxLo(myThid),myBxHi(myThid)
55    CML        CALL GAD_DIFF_X(bi,bj,k,xA,diff1,localT,df,myThid)
56    CML        DO j=1-Oly,sNy+Oly
57    CML         DO i=1-Olx,sNx+Olx
58    CML          fZon(i,j) = fZon(i,j) + df(i,j)
59    CML         ENDDO
60    CML        ENDDO
61    CML        CALL GAD_DIFF_Y(bi,bj,k,yA,diff1,localT,df,myThid)
62    CML        DO j=1-Oly,sNy+Oly
63    CML         DO i=1-Olx,sNx+Olx
64    CML          fMer(i,j) = fMer(i,j) + df(i,j)
65    CML         ENDDO
66    CML        ENDDO
67    CMLC--   Divergence of fluxes: update scalar field
68    CML        DO j=1-Oly,sNy+Oly-1
69    CML         DO i=1-Olx,sNx+Olx-1
70    CML          HEFF(i,j,1,bi,bj)=HEFF(i,j,1,bi,bj) + DELTT *
71    CML     &         maskC(i,j,kSurface,bi,bj)*recip_rA(i,j,bi,bj)
72    CML     &       *( (fZon(i+1,j)-fZon(i,j))
73    CML     &         +(fMer(i,j+1)-fMer(i,j))
74    CML     &                       )
75    CML     &         )
76    CML         ENDDO
77    CML        ENDDO
78    CML       ENDDO
79    CML      ENDDO
80    
81  C NOW DO DIFFUSION ON H(I,J,3)  C NOW DO DIFFUSION ON H(I,J,3)
82  C NOW CALCULATE DIFFUSION COEF ROUGHLY  C NOW CALCULATE DIFFUSION COEF ROUGHLY
83        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22