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

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

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

Revision 1.2 - (hide annotations) (download)
Wed Apr 7 23:40:34 2004 UTC (20 years, 1 month 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 jmc 1.2 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_ini_vars.F,v 1.1 2003/11/23 01:20:13 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "THSICE_OPTIONS.h"
5    
6 jmc 1.2 CBOP
7     C !ROUTINE: THSICE_INI_VARS
8     C !INTERFACE:
9 jmc 1.1 SUBROUTINE THSICE_INI_VARS( myThid )
10 jmc 1.2
11     C !DESCRIPTION: \bv
12 jmc 1.1 C *==========================================================*
13 jmc 1.2 C | S/R THSICE_INI_VARS
14 jmc 1.1 C | o initialize THermo_SeaICE variables
15     C *==========================================================*
16 jmc 1.2 C \ev
17    
18     C !USES:
19 jmc 1.1 IMPLICIT NONE
20    
21     C === Global variables ===
22     #include "SIZE.h"
23     #include "EEPARAMS.h"
24     #include "PARAMS.h"
25     #include "GRID.h"
26 jmc 1.2 c #include "DYNVARS.h"
27 jmc 1.1 #include "THSICE_PARAMS.h"
28 jmc 1.2 #include "THSICE_VARS.h"
29 jmc 1.1
30 jmc 1.2 C !INPUT/OUTPUT PARAMETERS:
31 jmc 1.1 C == Routine arguments ==
32     C myThid - Number of this instance of INI_RBCS
33     INTEGER myThid
34 jmc 1.2 CEOP
35 jmc 1.1
36     #ifdef ALLOW_THSICE
37     C == Local variables ==
38 jmc 1.2 C bi,bj :: Loop counters
39     C i,j :: Loop counters
40 jmc 1.1 INTEGER bi, bj
41     INTEGER I, J
42     INTEGER prec
43     CHARACTER*(MAX_LEN_FNAM) fn
44 jmc 1.2 _RL v2Loc
45 jmc 1.1
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 jmc 1.2 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 jmc 1.1 snowHeight(i,j,bi,bj)=0. _d 0
56 jmc 1.2 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 jmc 1.1 #ifdef COUPLE_MODEL
71 jmc 1.2 dFdT(i,j,bi,bj) = 0. _d 0
72 jmc 1.1 #endif
73 jmc 1.2 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 jmc 1.1 ENDDO
84     ENDDO
85     ENDDO
86     ENDDO
87    
88 jmc 1.2 IF ( startIceModel.LE.0 .AND. nIter0.NE.0 ) THEN
89 jmc 1.1 C-- Read ice pickup fields
90     _BARRIER
91     prec=precFloat64
92 jmc 1.2 CALL THSICE_READ_PICKUP( prec, nIter0, myThid )
93 jmc 1.1
94 jmc 1.2 _EXCH_XY_R8(iceMask,myThid)
95 jmc 1.1 _EXCH_XY_R8(iceHeight, myThid)
96 jmc 1.2 _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 jmc 1.1 #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