/[MITgcm]/MITgcm/pkg/therm_seaice/ice_start.F
ViewVC logotype

Contents of /MITgcm/pkg/therm_seaice/ice_start.F

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

Revision 1.1 - (show annotations) (download)
Thu Nov 21 19:11:42 2002 UTC (21 years, 6 months ago) by cheisey
Branch: MAIN
CVS Tags: checkpoint47b_post, checkpoint51f_post, checkpoint48i_post, checkpoint51l_post, checkpoint47a_post, checkpoint51k_post, checkpoint48d_post, checkpoint50b_post, checkpoint47i_post, checkpoint51o_post, checkpoint48g_post, branchpoint-genmake2, checkpoint50c_pre, checkpoint50, checkpoint51j_post, branch-exfmods-tag, checkpoint51q_post, checkpoint47e_post, checkpoint50f_post, checkpoint50a_post, checkpoint48e_post, checkpoint47c_post, checkpoint50f_pre, checkpoint52a_pre, checkpoint48b_post, checkpoint47j_post, checkpoint47d_pre, checkpoint50d_pre, checkpoint47h_post, checkpoint51d_post, checkpoint51, checkpoint51r_post, checkpoint48c_pre, checkpoint52, checkpoint52b_pre, checkpoint48c_post, checkpoint50d_post, checkpoint51o_pre, checkpoint47f_post, checkpoint51t_post, checkpoint51b_pre, checkpoint52a_post, checkpoint51i_post, checkpoint50e_post, checkpoint47g_post, checkpoint50h_post, checkpoint50c_post, checkpoint51a_post, checkpoint51n_pre, checkpoint47d_post, checkpoint50e_pre, checkpoint50b_pre, checkpoint48d_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, checkpoint51e_post, checkpoint51b_post, checkpoint48a_post, checkpoint51h_pre, checkpoint48f_post, checkpoint51i_pre, checkpoint51l_pre, checkpoint50g_post, checkpoint51u_post, checkpoint51c_post, checkpoint51g_post, checkpoint51m_post, ecco_c52_e35, checkpoint48, checkpoint49, checkpoint51f_pre, checkpoint51s_post, checkpoint48h_post
Branch point for: checkpoint51n_branch, branch-nonh, tg2-branch, branch-genmake2, branch-exfmods-curt 
Two packages:  bulk_force (Bulk forcing)
and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
Earlier integration from Stephaine Dutkiewicz
and Patrick Heimbach.

Two ifdef statements for compile time,
ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE

Two switches in data.pkg to turn on at run-time:

cat data.pkg
# Packages
 &PACKAGES
 useBulkForce=.TRUE.,
 useThermSeaIce=.TRUE.,
 &

WARNING:  useSEAICE and useThermSEAICE are mutually exclusive.

The bulk package requires an additional parameter file
with two namelists, data.ice and data.blk.

c ADAPTED FROM:
c LANL CICE.v2.0.2
c-----------------------------------------------------------------------
c.. thermodynamics (vertical physics) based on M. Winton 3-layer model
c.. See Bitz, C. M. and W. H. Lipscomb, 1999:  "An energy-conserving
c..       thermodynamic sea ice model for climate study."  J. Geophys.
c..       Res., 104, 15669 - 15677.
c..     Winton, M., 1999:  "A reformulated three-layer sea ice model."
c..       Submitted to J. Atmos. Ocean. Technol.

c.. authors Elizabeth C. Hunke and William Lipscomb
c..         Fluid Dynamics Group, Los Alamos National Laboratory
c-----------------------------------------------------------------------
1 #include "CPP_OPTIONS.h"
2
3 C !ROUTINE: ICE_START
4 C !INTERFACE:
5 SUBROUTINE ICE_START(i,j, bi, bj, myThid,
6 & esurp,fresh,fsalt,Tf, compact )
7 C *==========================================================*
8 C | S/R ICE_START
9 C | o Check water temperature and start ice layer if
10 C | appropriate.
11 C *==========================================================*
12 C \ev
13
14 C !USES:
15 IMPLICIT NONE
16 C == Global variables ==
17 #include "SIZE.h"
18 #include "DYNVARS.h"
19 #include "EEPARAMS.h"
20 #include "PARAMS.h"
21 #include "GRID.h"
22 cswdice -- add ---
23 #ifdef ALLOW_THERM_SEAICE
24 #include "ICE.h"
25 #include "ICE_DIAGS.h"
26 #include "BULKF_ICE_CONSTANTS.h"
27 #endif
28 cswdice -- end add ---
29
30 C !INPUT/OUTPUT PARAMETERS:
31 C == Routine Arguments ==
32 INTEGER bi,bj
33 INTEGER myThid
34 _RL esurp ! energy available for freezing
35 _RL fresh ! fresh water flux
36 _RL fsalt ! salt water flux
37 _RL Tf ! freezing temperature
38 _RL compact ! fraction of grid area covered in ice
39
40 C !LOCAL VARIABLES:
41 C == Local variables ==
42 INTEGER i,j
43 _RL qi0av, hi0new
44
45 #ifdef ALLOW_THERM_SEAICE
46
47 C start ice
48 compact=0.d0
49 qi0av=-rhoi*Lfresh
50 hi0new=-esurp*deltaTtracer/qi0av
51 c redistribute ice over portion of cell if too thin
52 if (hi0new.lt.hicemin) then
53 compact=hi0new/hicemin
54 iceHeight(i,j,bi,bj)=hicemin
55 else
56 compact=1.d0
57 iceHeight(i,j,bi,bj)=hi0new
58 endif
59 c theta(i,j,1,bi,bj)=Tf
60 Tsrf(i,j,bi,bj)=Tf
61 Tice1(i,j,bi,bj)=Tf
62 Tice2(i,j,bi,bj)=Tf
63 cQQ Qice1(i,j,bi,bj)=qi0av
64 cQQ Qice2(i,j,bi,bj)=qi0av
65 Qice1(i,j,bi,bj) = -cpwater*Tmlt1 +
66 & cpice *(Tmlt1-Tice1(i,j,bi,bj)) +
67 $ Lfresh*(1.-Tmlt1/Tice1(i,j,bi,bj))
68 Qice2(i,j,bi,bj) = -cpice *Tice2(i,j,bi,bj) + Lfresh
69 fresh=-(rhoi*iceheight(i,j,bi,bj))/deltaTtracer
70 cQQ fsalt=-(rhoi*iceheight(i,j,bi,bj)*saltice )
71 cQQ & *0.001/deltaTtracer
72 fsalt=-(rhoi*iceheight(i,j,bi,bj)*saltice )
73 & /35/deltaTtracer
74 iceMask(i,j,bi,bj)=compact
75 cBB
76 cBB print*,'**** ice start **** i,j', i,j
77 cBB print*,'qi1, qi2',Qice1(i,j,bi,bj), Qice2(i,j,bi,bj), Qi0av
78 cBB print*,'hi, fresh',hi0new, fresh
79 cBB print*,'*******************'
80 cBB print*,' '
81 #endif
82
83 RETURN
84 END
  ViewVC Help
Powered by ViewVC 1.1.22