/[MITgcm]/MITgcm/pkg/thsice/THSICE_VARS.h
ViewVC logotype

Contents of /MITgcm/pkg/thsice/THSICE_VARS.h

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


Revision 1.7 - (show annotations) (download)
Thu May 25 18:03:24 2006 UTC (18 years ago) by jmc
Branch: MAIN
CVS Tags: checkpoint58l_post, mitgcm_mapl_00, checkpoint58u_post, checkpoint58w_post, checkpoint58r_post, checkpoint58n_post, checkpoint58x_post, checkpoint58t_post, checkpoint58h_post, checkpoint58q_post, checkpoint58j_post, checkpoint58f_post, checkpoint58i_post, checkpoint58g_post, checkpoint58o_post, checkpoint58k_post, checkpoint58v_post, checkpoint58s_post, checkpoint58p_post, checkpoint58m_post
Changes since 1.6: +14 -3 lines
File MIME type: text/plain
- put i,j loops inside S/R: THSICE_ALBEDO, THSICE_SOLVE4TEMP, THSICE_EXTEND
   and THSICE_CALC_THICKN
- split thsice_step_fwd.F in 2 S/R: thsice_step_temp.F & thsice_step_fwd.F

1 C $Header: /u/gcmpack/MITgcm/pkg/thsice/THSICE_VARS.h,v 1.6 2006/05/14 14:06:13 jmc Exp $
2 C $Name: $
3
4 #ifdef ALLOW_THSICE
5
6 C !ROUTINE: THSICE_VARS.h
7 C -------------------------------
8 C THSICE_VARS.h
9 C variable for thermodynamics - Sea-Ice model
10 C -------------------------------
11
12 C-- COMMON /THSICE_VARS/ state variables of sea-ice model :
13 C iceMask :: sea-ice fraction: no ice=0, grid all ice 1 []
14 C iceHeight :: depth of ice layer [m]
15 C snowHeight:: depth of snow layer [m]
16 C Tsrf :: temperature at surface [oC]
17 C Tice1 :: temperature of ice layer 1 [oC]
18 C Tice2 :: temperature of ice layer 2 [oC]
19 C Qice1 :: enthalpy of ice layer 1 [J/kg]
20 C Qice2 :: enthalpy of ice layer 2 [J/kg]
21 C snowAge :: snow age [s]
22 COMMON /THSICE_VARS/
23 & iceMask, iceHeight, snowHeight,
24 & Tsrf, Tice1, Tice2,
25 & Qice1, Qice2, snowAge
26
27 _RL iceMask(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
28 _RL iceHeight (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
29 _RL snowHeight(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
30 _RL Tsrf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
31 _RL Tice1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
32 _RL Tice2 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
33 _RL Qice1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
34 _RL Qice2 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
35 _RL snowAge(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
36
37 C-- COMMON /THSICE_FLUX/ hold fluxes needed for thsice model
38 C sHeating :: surf heating left to melt snow or ice (= Atmos-Conduction)
39 C flxCndBt :: heat flux conducted through the ice to bottom surface
40 C snowPrc :: snow precipitation [kg/m2/s]
41 C siceAlb :: area weighted sea-ice albedo [0-1]
42 C dFdT :: heat deriveative for coupled model
43 C atmospheric fluxes (change along the time-stepping):
44 C icFlxSW :: short-wave heat flux (+=down) over sea-ice
45 C (downward SW / net SW @ surface / net SW below sea-ice)
46 C icFlxAtm :: Atmospheric surf. heat flux over sea-ice [W/m2] (+=down)
47 C (over sea-ice only / weighted by ice-fraction)
48 C icFrwAtm :: fresh-water flux (E-P) from the atmosphere [kg/m2/s] (+=up)
49 C ( ice Evap only / ice E-P / ice - ocean weighted E-P )
50 C oceQnet :: net heat flux to the ocean (+=down) [W/m2]
51 C oceQsw :: net short-wave that enter the ocean (+=down) [W/m2]
52 C oceFWfx :: net fresh water flux to the ocean (+=down) [kg/m2]
53 C oceSflx :: net salt flux to the ocean (+=down) [psu.kg/m2]
54 COMMON / THSICE_FLUX /
55 c & oceQsw, oceQnet, oceFWfx, oceSflx,
56 & sHeating, flxCndBt,
57 & snowPrc, siceAlb,
58 #ifdef COUPLE_MODEL
59 & dFdT,
60 #endif
61 & icFlxSW, icFlxAtm, icFrwAtm
62
63 c _RL oceQnet(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
64 c _RL oceQsw (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
65 c _RL oceFWfx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
66 c _RL oceSflx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
67 _RL sHeating(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
68 _RL flxCndBt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
69 _RL snowPrc (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
70 _RL siceAlb (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
71 #ifdef COUPLE_MODEL
72 _RL dFdT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
73 #endif
74 _RL icFlxSW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
75 _RL icFlxAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
76 _RL icFrwAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
77
78 C-- COMMON / THSICE_OCEMXLAYER / oceanic mixed layer state
79 C hOceMxL :: thickness of the ocean mixed layer [m]
80 C tOceMxL :: temperature in the ocean mixed layer [oC]
81 C sOceMxL :: salinity in the ocean mixed layer [psu]
82 C v2ocMxL :: velocity (square) in the mixed layer [m2/s2]
83 COMMON / THSICE_OCEMXLAYER /
84 & hOceMxL, tOceMxL, sOceMxL, v2ocMxL
85
86 _RL hOceMxL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87 _RL tOceMxL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88 _RL sOceMxL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
89 _RL v2ocMxL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90
91 #ifdef ALLOW_AUTODIFF_TAMC
92 INTEGER iicekey
93 #endif
94
95 #endif /* ALLOW_THSICE */
96
97 CEH3 ;;; Local Variables: ***
98 CEH3 ;;; mode:fortran ***
99 CEH3 ;;; End: ***

  ViewVC Help
Powered by ViewVC 1.1.22