/[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.20 - (hide annotations) (download)
Wed Sep 2 19:28:44 2009 UTC (14 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint63p, checkpoint63q, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint62, checkpoint63, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61v, checkpoint61w, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.19: +2 -2 lines 
use right type of S/R to read "Tsrf" (RL type)
1 jmc 1.20 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_ini_vars.F,v 1.19 2009/08/25 21:33:28 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "THSICE_OPTIONS.h"
5 jmc 1.8 #ifdef ALLOW_AIM
6     # include "AIM_OPTIONS.h"
7     #endif
8 jmc 1.1
9 jmc 1.2 CBOP
10     C !ROUTINE: THSICE_INI_VARS
11     C !INTERFACE:
12 jmc 1.1 SUBROUTINE THSICE_INI_VARS( myThid )
13 jmc 1.2
14     C !DESCRIPTION: \bv
15 jmc 1.1 C *==========================================================*
16 jmc 1.2 C | S/R THSICE_INI_VARS
17 jmc 1.1 C | o initialize THermo_SeaICE variables
18     C *==========================================================*
19 jmc 1.2 C \ev
20    
21     C !USES:
22 jmc 1.1 IMPLICIT NONE
23    
24     C === Global variables ===
25     #include "SIZE.h"
26     #include "EEPARAMS.h"
27     #include "PARAMS.h"
28 jmc 1.7 #include "FFIELDS.h"
29 jmc 1.1 #include "THSICE_PARAMS.h"
30 jmc 1.2 #include "THSICE_VARS.h"
31 jmc 1.8 #ifdef ALLOW_AIM
32     # include "AIM_FFIELDS.h"
33     #endif
34 jmc 1.1
35 jmc 1.2 C !INPUT/OUTPUT PARAMETERS:
36 jmc 1.1 C == Routine arguments ==
37 jmc 1.8 C myThid :: My Thread Id. number
38 jmc 1.1 INTEGER myThid
39 jmc 1.2 CEOP
40 jmc 1.1
41     #ifdef ALLOW_THSICE
42     C == Local variables ==
43 jmc 1.2 C bi,bj :: Loop counters
44     C i,j :: Loop counters
45 jmc 1.1 INTEGER bi, bj
46 jmc 1.8 INTEGER i, j
47 jmc 1.1 INTEGER prec
48 jmc 1.5 c CHARACTER*(MAX_LEN_FNAM) fn
49 jmc 1.2 _RL v2Loc
50 dimitri 1.15 _RL Tf
51 jmc 1.1
52     c set up ice arrays to zero if starting ice
53 jmc 1.8 DO bj = myByLo(myThid), myByHi(myThid)
54     DO bi = myBxLo(myThid), myBxHi(myThid)
55 jmc 1.2 C- state variables :
56 jmc 1.6 DO j=1-OLy,sNy+OLy
57     DO i=1-OLx,sNx+OLx
58 jmc 1.2 iceMask(i,j,bi,bj) = 0. _d 0
59     iceHeight(i,j,bi,bj)= 0. _d 0
60 jmc 1.1 snowHeight(i,j,bi,bj)=0. _d 0
61 jmc 1.2 Tsrf(i,j,bi,bj) = 0. _d 0
62     Tice1(i,j,bi,bj) = 0. _d 0
63     Tice2(i,j,bi,bj) = 0. _d 0
64     Qice1(i,j,bi,bj) = 0. _d 0
65     Qice2(i,j,bi,bj) = 0. _d 0
66     snowAge(i,j,bi,bj) = 0. _d 0
67     ENDDO
68     ENDDO
69     C- fluxes :
70 jmc 1.6 DO j=1-OLy,sNy+OLy
71     DO i=1-OLx,sNx+OLx
72 jmc 1.2 sHeating(i,j,bi,bj) = 0. _d 0
73     flxCndBt(i,j,bi,bj) = 0. _d 0
74     snowPrc(i,j,bi,bj) = 0. _d 0
75 jmc 1.4 siceAlb(i,j,bi,bj) = 0. _d 0
76 jmc 1.6 icFlxSW (i,j,bi,bj) = 0. _d 0
77     icFlxAtm(i,j,bi,bj) = 0. _d 0
78     icFrwAtm(i,j,bi,bj) = 0. _d 0
79 jmc 1.11 C- needed when using advection/diffusion:
80 heimbach 1.9 oceFWfx(i,j,bi,bj) = 0. _d 0
81     oceSflx(i,j,bi,bj) = 0. _d 0
82     oceQnet(i,j,bi,bj) = 0. _d 0
83 jmc 1.2 ENDDO
84     ENDDO
85     C- oceanic mixed layer state :
86     v2Loc = vMxL_default*vMxL_default
87 jmc 1.6 DO j=1-OLy,sNy+OLy
88     DO i=1-OLx,sNx+OLx
89 jmc 1.2 hOceMxL(i,j,bi,bj) = hMxL_default
90     tOceMxL(i,j,bi,bj) = 0. _d 0
91     sOceMxL(i,j,bi,bj) = sMxL_default
92     v2ocMxL(i,j,bi,bj) = v2Loc
93 jmc 1.1 ENDDO
94     ENDDO
95 jmc 1.8 #ifdef ALLOW_AIM
96     IF ( useAIM ) THEN
97     C- Mask mixed layer depth : depth is used in thsice slab_ocean
98     C and this mask is used in thsice_advdiff and if coupled
99     DO j=1-OLy,sNy+OLy
100     DO i=1-OLx,sNx+OLx
101     IF ( aim_landFr(i,j,bi,bj).EQ.1. _d 0 )
102     & hOceMxL(i,j,bi,bj) = 0.
103     ENDDO
104     ENDDO
105     ENDIF
106     #endif /* ALLOW_AIM */
107 jmc 1.1 ENDDO
108 jmc 1.8 ENDDO
109 jmc 1.1
110 jmc 1.3 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
111    
112 jmc 1.19 IF ( startIceModel.LE.0 .AND.
113     & ( nIter0.NE.0 .OR. pickupSuff.NE.' ' )
114     & ) THEN
115 jmc 1.1 C-- Read ice pickup fields
116     prec=precFloat64
117 jmc 1.2 CALL THSICE_READ_PICKUP( prec, nIter0, myThid )
118 jmc 1.1
119 jmc 1.3 ELSE
120     C-- Read initial conditions:
121 dimitri 1.15 IF ( thSIceFract_InitFile .NE. ' ' ) THEN
122     CALL READ_REC_XY_RL(thSIceFract_InitFile,iceMask,1,0,myThid)
123     ENDIF
124 dimitri 1.12 IF ( thSIceThick_InitFile .NE. ' ' ) THEN
125     CALL READ_REC_XY_RL(thSIceThick_InitFile,iceHeight,1,0,myThid)
126     ENDIF
127 jmc 1.3 IF ( thSIceSnowH_InitFile .NE. ' ' ) THEN
128     CALL READ_REC_XY_RL(thSIceSnowH_InitFile,snowHeight,1,0,myThid)
129     ENDIF
130     IF ( thSIceSnowA_InitFile .NE. ' ' ) THEN
131     CALL READ_REC_XY_RL(thSIceSnowA_InitFile,snowAge,1,0,myThid)
132     ENDIF
133     IF ( thSIceEnthp_InitFile .NE. ' ' ) THEN
134     CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice1,1,0,myThid)
135     CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice2,2,0,myThid)
136 dimitri 1.15 ENDIF
137     IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN
138 jmc 1.20 CALL READ_REC_XY_RL(thSIceTsurf_InitFile,Tsrf,1,0,myThid)
139 dimitri 1.15 ENDIF
140     IF ( thSIceEnthp_InitFile .EQ. ' ' ) THEN
141     C- enthalpy of new ice in J/kg, taken from thsice_extend.F with Tf beeing
142     C the freezing Temp of seawater computed from a fixed salinity (31.5 psu)
143 dimitri 1.16 C Tf = -mu_Tf*salinity = -1.70 deg C; Qice1 ~ 3.2e5; Qice2 ~ 3.4e5
144 dimitri 1.15 Tf = -1.70 _d 0
145 dimitri 1.12 DO bj = myByLo(myThid), myByHi(myThid)
146     DO bi = myBxLo(myThid), myBxHi(myThid)
147     DO j=1-OLy,sNy+OLy
148     DO i=1-OLx,sNx+OLx
149     IF (iceMask(i,j,bi,bj) .NE. 0. _d 0) THEN
150 dimitri 1.15 Qice1(i,j,bi,bj) = -cpWater*Tmlt1
151     & + cpIce *(Tmlt1-Tf) + Lfresh*(1. _d 0-Tmlt1/Tf)
152     Qice2(i,j,bi,bj) = -cpIce *Tf + Lfresh
153 dimitri 1.12 ENDIF
154     ENDDO
155     ENDDO
156     ENDDO
157     ENDDO
158 jmc 1.3 ENDIF
159     ENDIF
160    
161 jmc 1.18 _EXCH_XY_RL(iceMask,myThid)
162     _EXCH_XY_RL(iceHeight, myThid)
163     _EXCH_XY_RL(snowHeight,myThid)
164     _EXCH_XY_RL(Tsrf, myThid)
165     _EXCH_XY_RL(Tice1, myThid)
166     _EXCH_XY_RL(Tice2, myThid)
167     _EXCH_XY_RL(Qice1, myThid)
168     _EXCH_XY_RL(Qice2, myThid)
169     _EXCH_XY_RL(snowAge,myThid)
170 jmc 1.1
171 jmc 1.7 C-- Initialise Sea-Ice Loading for SeaIce-Dynamics :
172 dimitri 1.15 IF ( useSEAICE ) THEN
173     DO bj = myByLo(myThid), myByHi(myThid)
174     DO bi = myBxLo(myThid), myBxHi(myThid)
175     DO j=1-OLy,sNy+OLy
176     DO i=1-OLx,sNx+OLx
177     sIceLoad(i,j,bi,bj) = ( snowHeight(i,j,bi,bj)*rhos
178     & + iceHeight(i,j,bi,bj)*rhoi
179     & )*iceMask(i,j,bi,bj)
180 jmc 1.7 ENDDO
181 dimitri 1.15 ENDDO
182     ENDDO
183     ENDDO
184     ENDIF
185 jmc 1.7
186 jmc 1.1 #endif /* ALLOW_THSICE */
187    
188     RETURN
189     END
  ViewVC Help
Powered by ViewVC 1.1.22