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

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

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

Revision 1.4 - (show annotations) (download)
Sat Dec 28 10:11:11 2002 UTC (21 years, 5 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint47j_post, checkpoint48d_pre, checkpoint47f_post, checkpoint48d_post, checkpoint48a_post, checkpoint48e_post, checkpoint47i_post, checkpoint47h_post, checkpoint48c_post, checkpoint48, checkpoint47g_post, checkpoint48b_post, checkpoint48c_pre
Branch point for: c24_e25_ice, ecco-branch
Changes since 1.3: +8 -9 lines 
checkpoint47f_post
Merging from release1_p10:
o modifications for using pkg/exf with pkg/seaice
  - pkg/seaice CPP options SEAICE_EXTERNAL_FORCING
    and SEAICE_EXTERNAL_FLUXES
  - pkg/exf CPP options EXF_READ_EVAP and
    EXF_NO_BULK_COMPUTATIONS
  - usage examples are Experiments 8 and 9 in
    verification/lab_sea/README
  - verification/lab_sea default experiment now uses
    pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf
1 C $Header:
2
3 #include "SEAICE_OPTIONS.h"
4
5 CStartOfInterface
6 SUBROUTINE groatb( A22, myThid )
7 C /==========================================================\
8 C | SUBROUTINE groatb |
9 C | o Calculate ice growth |
10 C |==========================================================|
11 C \==========================================================/
12 IMPLICIT NONE
13
14 C === Global variables ===
15 #include "SIZE.h"
16 #include "EEPARAMS.h"
17 #include "PARAMS.h"
18 #include "DYNVARS.h"
19 #include "FFIELDS.h"
20 #include "SEAICE.h"
21 #include "SEAICE_FFIELDS.h"
22
23 C === Routine arguments ===
24 C myThid - Thread no. that called this routine.
25 _RL A22
26 INTEGER myThid
27 CEndOfInterface
28
29 #ifdef ALLOW_SEAICE
30
31 C === Local variables ===
32 C i,j,bi,bj - Loop counters
33
34 INTEGER i, j, bi, bj
35 INTEGER KOPEN
36 _RL U1, V1
37
38 _RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy)
39 _RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy)
40
41 C if this is not done before then do it here
42 _EXCH_XY_R8(UWIND, myThid)
43 _EXCH_XY_R8(VWIND, myThid)
44
45 c$taf loop = parallel
46 DO bj=myByLo(myThid),myByHi(myThid)
47 c$taf loop = parallel
48 DO bi=myBxLo(myThid),myBxHi(myThid)
49
50 C DETERMINE AMOUNT OF OPEN WATER AND ICE THICKNESS
51 DO J=1,sNy
52 DO I=1,sNx
53 AREA(I,J,2,bi,bj)=MAX(A22,AREA(I,J,2,bi,bj))
54 FHEFF(I,J,bi,bj)=0.0 _d 0
55 HICE(I,J)=HEFF(I,J,2,bi,bj)/AREA(I,J,2,bi,bj)
56 ENDDO
57 ENDDO
58
59 C NOW DETERMINE MIXED LAYER TEMPERATURE
60 DO J=1,sNy
61 DO I=1,sNx
62 TMIX(I,J,bi,bj)=theta(I,J,1,bi,bj)+273.16 _d +00
63 #ifdef SEAICE_DEBUG
64 TMIX(I,J,bi,bj)=MAX(TMIX(I,J,bi,bj),271.2 _d +00)
65 #endif SEAICE_DEBUG
66 ENDDO
67 ENDDO
68
69 DO J=1,sNy
70 DO I=1,sNx
71 U1=0.25 _d 0*(UWIND(I-1,J-1,bi,bj)+UWIND(I-1,J,bi,bj)
72 1 +UWIND(I,J-1,bi,bj)+UWIND(I,J,bi,bj))
73 V1=0.25 _d 0*(VWIND(I-1,J-1,bi,bj)+VWIND(I-1,J,bi,bj)
74 1 +VWIND(I,J-1,bi,bj)+VWIND(I,J,bi,bj))
75 UG(I,J)=SQRT(U1**2+V1**2)
76 ENDDO
77 ENDDO
78
79 C NOW DETERMINE GROWTH RATES
80 C FIRST DO OPEN WATER
81 KOPEN=-1
82 CALL BUDGET(UG, TMIX, HICE, FO, KOPEN, bi, bj)
83 C NOW DO ICE
84 KOPEN=1
85 CALL BUDGET(UG, TICE, HICE, FICE, KOPEN, bi, bj)
86
87 ENDDO
88 ENDDO
89
90 #endif ALLOW_SEAICE
91
92 RETURN
93 END
  ViewVC Help
Powered by ViewVC 1.1.22