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

Contents of /MITgcm/pkg/seaice/diffus.F

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

Revision 1.2 - (show annotations) (download)
Tue Nov 12 20:47:27 2002 UTC (21 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint47a_post, checkpoint47
Changes since 1.1: +83 -0 lines 
Merging from release1_p8 branch:
o New package: pkg/seaice
  Sea ice model by D. Menemenlis (JPL) and Jinlun Zhang (Seattle).
  The sea-ice code is based on Hibler (1979-1980).
  Two sea-ice dynamic solvers, ADI and LSR, are included.
  In addition to computing prognostic sea-ice variables and diagnosing
  the forcing/external data fields that drive the ocean model,
  SEAICE_MODEL also sets theta to the freezing point under sea-ice.
  The implied surface heat flux is then stored in variable
  surfaceTendencyTice, which is needed by KPP package (kpp_calc.F and
  kpp_transport_t.F) to diagnose surface buoyancy fluxes and for the
  non-local transport term.  Because this call precedes model
  thermodynamics, temperature under sea-ice may not be "exactly" at
  the freezing point by the time theta is dumped or time-averaged.
1 C $Header:
2
3 #include "SEAICE_OPTIONS.h"
4
5 CStartOfInterface
6 SUBROUTINE diffus( HEFF, DIFFA, HEFFM, DELTT, myThid )
7 C /==========================================================\
8 C | SUBROUTINE diffus |
9 C | o Add diffusion terms to ice mass conservation equations |
10 C |==========================================================|
11 C \==========================================================/
12 IMPLICIT NONE
13
14 C === Global variables ===
15 #include "SIZE.h"
16 #include "EEPARAMS.h"
17 #include "SEAICE_GRID.h"
18
19 C === Routine arguments ===
20 C myThid - Thread no. that called this routine.
21 _RL HEFF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,3,nSx,nSy)
22 _RL HEFFM (1-OLx:sNx+OLx,1-OLy:sNy+OLy, nSx,nSy)
23 _RL DIFFA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
24 _RL DELTT
25 INTEGER myThid
26 CEndOfInterface
27
28 #ifdef ALLOW_SEAICE
29
30 C === Local variables ===
31 C i,j,k,bi,bj - Loop counters
32
33 INTEGER i, j, k, bi, bj
34 _RL DELTXX1, DELTYY1, DELTXX, DELTYY
35 _RL HEFF1 (1-OLx:sNx+OLx, 1-OLy:sNy+OLy)
36
37 DO bj=myByLo(myThid),myByHi(myThid)
38 DO bi=myBxLo(myThid),myBxHi(myThid)
39
40 DO j=1-OLy,sNy+OLy
41 DO i=1-OLx,sNx+OLx
42 HEFF1(I,J)=0.0
43 ENDDO
44 ENDDO
45
46 C NOW DO DIFFUSION WITH NUIT CONVERSION
47 DO j=1,sNy
48 DO i=1,sNx
49 DELTXX1=DELTT*DIFFA(I,J)
50 DELTYY1=DELTT*DIFFA(I,J)
51 DELTXX=DELTXX1/(DXTICE(I,J,bi,bj)*DXTICE(I,J,bi,bj)
52 & *CSTICE(I,J,bi,bj)*CSTICE(I,J,bi,bj))
53 DELTYY=DELTYY1/(DYTICE(I,J,bi,bj)*DYTICE(I,J,bi,bj)
54 & *CSTICE(I,J,bi,bj))
55 HEFF1(I,J)=DELTXX*((HEFF(I+1,J,3,bi,bj)-HEFF(I,J,3,bi,bj))
56 & *HEFFM(I+1,J,bi,bj)
57 & -(HEFF(I,J,3,bi,bj)-HEFF(I-1,J,3,bi,bj))*HEFFM(I-1,J,bi,bj))
58 & +DELTYY*((HEFF(I,J+1,3,bi,bj)-HEFF(I,J,3,bi,bj))
59 & *CSUICE(I,J,bi,bj)*HEFFM(I,J+1,bi,bj)
60 & -(HEFF(I,J,3,bi,bj)-HEFF(I,J-1,3,bi,bj))
61 & *CSUICE(I,J-1,bi,bj)*HEFFM(I,J-1,bi,bj))
62 ENDDO
63 ENDDO
64
65 ENDDO
66 ENDDO
67
68 _EXCH_XY_R8(HEFF1, myThid)
69
70 DO bj=myByLo(myThid),myByHi(myThid)
71 DO bi=myBxLo(myThid),myBxHi(myThid)
72 DO j=1-OLy,sNy+OLy
73 DO i=1-OLx,sNx+OLx
74 HEFF(I,J,3,bi,bj)=HEFF1(I,J)
75 ENDDO
76 ENDDO
77 ENDDO
78 ENDDO
79
80 #endif ALLOW_SEAICE
81
82 RETURN
83 END
  ViewVC Help
Powered by ViewVC 1.1.22