/[MITgcm]/MITgcm/pkg/thsice/thsice_ini_vars.F
ViewVC logotype

Contents of /MITgcm/pkg/thsice/thsice_ini_vars.F

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

Revision 1.2 - (show annotations) (download)
Wed Apr 7 23:40:34 2004 UTC (20 years, 2 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint52n_post, checkpoint53d_post, checkpoint54a_pre, checkpoint55c_post, checkpoint54e_post, checkpoint54a_post, checkpoint53c_post, checkpoint55d_pre, checkpoint55j_post, checkpoint56b_post, checkpoint55h_post, checkpoint53b_post, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint54d_post, checkpoint56c_post, checkpoint52m_post, checkpoint55, checkpoint53a_post, checkpoint57a_post, checkpoint54, checkpoint54f_post, checkpoint55g_post, checkpoint55f_post, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, checkpoint53g_post, checkpoint55e_post, checkpoint53f_post, checkpoint55a_post, checkpoint53d_pre, checkpoint54c_post, checkpoint56a_post, checkpoint55d_post
Changes since 1.1: +57 -29 lines 
major changes in pkg/thsice: allows atmospheric model (AIM) to use thsice.
- split thsice_therm.F in 2 S/R: thsice_solve4temp.F & thsice_calc_thickn.F
- move most of the ocean & bulk_force interface in thsice_main.F
- add a "slab ocean" component to be used with atmospheric model
1 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_ini_vars.F,v 1.1 2003/11/23 01:20:13 jmc Exp $
2 C $Name: $
3
4 #include "THSICE_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: THSICE_INI_VARS
8 C !INTERFACE:
9 SUBROUTINE THSICE_INI_VARS( myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | S/R THSICE_INI_VARS
14 C | o initialize THermo_SeaICE variables
15 C *==========================================================*
16 C \ev
17
18 C !USES:
19 IMPLICIT NONE
20
21 C === Global variables ===
22 #include "SIZE.h"
23 #include "EEPARAMS.h"
24 #include "PARAMS.h"
25 #include "GRID.h"
26 c #include "DYNVARS.h"
27 #include "THSICE_PARAMS.h"
28 #include "THSICE_VARS.h"
29
30 C !INPUT/OUTPUT PARAMETERS:
31 C == Routine arguments ==
32 C myThid - Number of this instance of INI_RBCS
33 INTEGER myThid
34 CEOP
35
36 #ifdef ALLOW_THSICE
37 C == Local variables ==
38 C bi,bj :: Loop counters
39 C i,j :: Loop counters
40 INTEGER bi, bj
41 INTEGER I, J
42 INTEGER prec
43 CHARACTER*(MAX_LEN_FNAM) fn
44 _RL v2Loc
45
46
47 c set up ice arrays to zero if starting ice
48 DO bj = myByLo(myThid), myByHi(myThid)
49 DO bi = myBxLo(myThid), myBxHi(myThid)
50 C- state variables :
51 DO j=1-Oly,sNy+Oly
52 DO i=1-Olx,sNx+Olx
53 iceMask(i,j,bi,bj) = 0. _d 0
54 iceHeight(i,j,bi,bj)= 0. _d 0
55 snowHeight(i,j,bi,bj)=0. _d 0
56 Tsrf(i,j,bi,bj) = 0. _d 0
57 Tice1(i,j,bi,bj) = 0. _d 0
58 Tice2(i,j,bi,bj) = 0. _d 0
59 Qice1(i,j,bi,bj) = 0. _d 0
60 Qice2(i,j,bi,bj) = 0. _d 0
61 snowAge(i,j,bi,bj) = 0. _d 0
62 ENDDO
63 ENDDO
64 C- fluxes :
65 DO j=1-Oly,sNy+Oly
66 DO i=1-Olx,sNx+Olx
67 sHeating(i,j,bi,bj) = 0. _d 0
68 flxCndBt(i,j,bi,bj) = 0. _d 0
69 snowPrc(i,j,bi,bj) = 0. _d 0
70 #ifdef COUPLE_MODEL
71 dFdT(i,j,bi,bj) = 0. _d 0
72 #endif
73 ENDDO
74 ENDDO
75 C- oceanic mixed layer state :
76 v2Loc = vMxL_default*vMxL_default
77 DO j=1-Oly,sNy+Oly
78 DO i=1-Olx,sNx+Olx
79 hOceMxL(i,j,bi,bj) = hMxL_default
80 tOceMxL(i,j,bi,bj) = 0. _d 0
81 sOceMxL(i,j,bi,bj) = sMxL_default
82 v2ocMxL(i,j,bi,bj) = v2Loc
83 ENDDO
84 ENDDO
85 ENDDO
86 ENDDO
87
88 IF ( startIceModel.LE.0 .AND. nIter0.NE.0 ) THEN
89 C-- Read ice pickup fields
90 _BARRIER
91 prec=precFloat64
92 CALL THSICE_READ_PICKUP( prec, nIter0, myThid )
93
94 _EXCH_XY_R8(iceMask,myThid)
95 _EXCH_XY_R8(iceHeight, myThid)
96 _EXCH_XY_R8(snowHeight,myThid)
97 _EXCH_XY_R8(Tsrf, myThid)
98 _EXCH_XY_R8(Tice1, myThid)
99 _EXCH_XY_R8(Tice2, myThid)
100 _EXCH_XY_R8(Qice1, myThid)
101 _EXCH_XY_R8(Qice2, myThid)
102 _EXCH_XY_R8(snowAge,myThid)
103 #ifdef COUPLE_MODEL
104 _EXCH_XY_R8(dFdT, myThid)
105 #endif
106 ENDIF
107
108 #endif /* ALLOW_THSICE */
109
110 RETURN
111 END
  ViewVC Help
Powered by ViewVC 1.1.22