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

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

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

Revision 1.1 - (hide 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 cheisey 1.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