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

Contents of /MITgcm/pkg/therm_seaice/ice_extrastart.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, 5 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_EXTRASTART
4 C !INTERFACE:
5 SUBROUTINE ICE_EXTRASTART(i,j, bi, bj, myThid, qleft,
6 & esurp,fresh,fsalt,Tf,Ta, compact )
7 C *==========================================================*
8 C | S/R ICE_EXTRASTART
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 Ta ! air temperature
39 _RL compact ! fraction of grid area covered in ice
40 _RL qleft
41
42 C !LOCAL VARIABLES:
43 C == Local variables ==
44 INTEGER i,j
45 _RL qi0av, hi0new,compact2, hice_old
46 _RL qicen (2) ! ice enthalpy (J m-3)
47 _RL hlyr ! individual ice layer thickness (m)
48 _RL hnew(2) ! new ice layer thickness (m)
49 _RL hi ! ice height
50 _RL hs ! snow height
51
52 #ifdef ALLOW_THERM_SEAICE
53
54 hice_old=iceHeight(i,j,bi,bj)
55 hi=iceHeight(i,j,bi,bj)
56 hs=snowHeight(i,j,bi,bj)
57 qicen(1)=Qice1(i,j,bi,bj)
58 qicen(2)=Qice2(i,j,bi,bj)
59 hnew(1)=0.d0
60 hnew(2)=0.d0
61 C start ice
62 qi0av=-rhoi*Lfresh
63 hi0new=-esurp*deltaTtracer/qi0av
64 c redistribute ice over portion of cell if too thin
65 if (hi0new.lt.hicemin) then
66 compact2=(1.d0-compact)*hi0new/hicemin
67 hi0new=hicemin
68 else
69 compact2=(1.d0-compact)
70 endif
71 if ((compact+compact2).gt.1.d0) then
72 print*,'QQQQQQQQQ ice mask > 1, ice_extrastart'
73 stop
74 endif
75 hi=(compact2*hi0new+compact*hi)/
76 & (compact+compact2)
77 hs=(compact2*0.d0+compact*hs)/
78 & (compact+compact2)
79 hlyr=hi/2
80 hnew(2)=hi-hice_old
81 call new_layers_winton(hs,hi,hlyr,hnew,qicen)
82 fresh=(compact2*(-(rhoi*hi0new)/deltaTtracer)+
83 & compact*fresh)/(compact+compact2)
84 fsalt=(compact2*(-(rhoi*hi0new*saltice )
85 & /35.d0/deltaTtracer) + compact*fsalt)/
86 & (compact+compact2)
87 qleft=(-compact2*esurp+compact*qleft)/(compact+compact2)
88 iceMask(i,j,bi,bj)=compact+compact2
89 iceHeight(i,j,bi,bj)=hi
90 snowHeight(i,j,bi,bj)=hs
91 Qice1(i,j,bi,bj)=qicen(1)
92 Qice2(i,j,bi,bj)=qicen(2)
93 cBB
94 cBB print*,'**** ice start **** i,j', i,j
95 cBB print*,'qi1, qi2',Qice1(i,j,bi,bj), Qice2(i,j,bi,bj), Qi0av
96 cBB print*,'hi, fresh',hi0new, fresh
97 cBB print*,'*******************'
98 cBB print*,' '
99 #endif
100
101 RETURN
102 END
  ViewVC Help
Powered by ViewVC 1.1.22