/[MITgcm]/MITgcm/pkg/shelfice/SHELFICE.h
ViewVC logotype

Contents of /MITgcm/pkg/shelfice/SHELFICE.h

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


Revision 1.5 - (show annotations) (download)
Wed Oct 10 09:26:39 2007 UTC (16 years, 11 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint60, checkpoint61, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59k, checkpoint59j, checkpoint61b, checkpoint61c, checkpoint61a
Changes since 1.4: +5 -1 lines
File MIME type: text/plain
add option for convservative form of Hellmer&Olbers (1989) thermodynamics
- requires a little reorganization, which affects the testreport results
  of isomip.htd (only 7 digits of cg2d agree) => update experiment
  The results are extremely sensitive to details of the algorithm
  (e.g., order in which sums and differences are taken, etc) so that
  replacing a term with a (theoretically) exactly equal other term
  causes the big difference; possibly, because the solution of a
  quadratic equation is involed. I am not happy about this, but the
  present formulation appears to be a little more robust.
- add a few comments

1 C $Header: /u/gcmpack/MITgcm/pkg/shelfice/SHELFICE.h,v 1.4 2006/08/14 16:52:46 mlosch Exp $
2 C $Name: $
3
4 #ifdef ALLOW_SHELFICE
5
6 CBOP
7 C !ROUTINE: SHELFICE.h
8
9 C !DESCRIPTION: \bv
10 C /==========================================================\
11 C | SHELFICE.h |
12 C | o Basic header thermodnynamic shelf ice package. |
13 C | Contains all SHELFICE field declarations. |
14 C \==========================================================/
15
16 C-----------------------------------------------------------------------
17 C
18 C-- Constants that can be set in data.shelfice
19 C SHELFICEloadAnomalyFile - name of shelfice load anomaly file
20 C SHELFICEDragLinear - linear drag at bottom shelfice (1/s)
21 C SHELFICEDragQuadratic - quadratic drag at bottom shelfice (1/m)
22 C SHELFICEheatTransCoeff - heat transfer coefficient that determines
23 C heat flux into shelfice (m/s)
24 C SHELFICEsaltTransCoeff - salinity transfer coefficient that determines
25 C salt flux into shelfice (m/s)
26 C SHELFICElatentHeat - latent heat of fusion (J/kg)
27 C useISOMIPTD - use simple ISOMIP thermodynamics
28 C SHELFICEconserve - use conservative form of H&O-thermodynamics
29 C following Jenkins et al. (2001, JPO)
30 C SHELFICEboundaryLayer - turn on vertical merging of cells to for a
31 C boundary layer of drF thickness
32 C no_slip_shelfice - set slip conditions for shelfice separately,
33 C (by default the same as no_slip_bottom)
34 C SHELFICEwriteState - enable output
35 C SHELFICE_dump_mnc - use netcdf for snapshot output
36 C SHELFICE_tave_mnc - use netcdf for time-averaged output
37 C SHELFICE_dumpFreq - analoguous to dumpFreq (= default)
38 C SHELFICE_taveFreq - analoguous to taveFreq (= default)
39 C
40 C-- Fields
41 C ktopC - index of the top "wet cell" (2D)
42 C R_shelfIce - shelfice topography [m]
43 C shelficeLoadAnomaly - pressure load anomaly of shelfice [Pa]
44 C shelficeHeatFlux - upward heat flux [W/m^2]
45 C shelficeFreshWaterFlux - upward fresh water flux (virt. salt flux) [m/s]
46 C shelficeForcingT - analogue of surfaceForcingT
47 C shelficeForcingS - analogue of surfaceForcingS
48 C-----------------------------------------------------------------------
49 C \ev
50 CEOP
51
52 COMMON /SHELFICE_PARMS_I/ kTopC
53 INTEGER kTopC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
54
55 COMMON /SHELFICE_PARMS_R/
56 & SHELFICE_dumpFreq, SHELFICE_taveFreq,
57 & SHELFICEheatTransCoeff, SHELFICEsaltTransCoeff,
58 & rhoShelfice, SHELFICEkappa,
59 & SHELFICElatentHeat, recip_SHELFICElatentHeat,
60 & SHELFICEheatCapacity_Cp,
61 & SHELFICEthetaSurface,
62 & SHELFICEDragLinear, SHELFICEDragQuadratic
63 _RL SHELFICE_dumpFreq, SHELFICE_taveFreq
64 _RL SHELFICEheatTransCoeff
65 _RL SHELFICEsaltTransCoeff
66 _RL SHELFICElatentHeat
67 _RL SHELFICEheatCapacity_Cp
68 _RL rhoShelfice
69 _RL SHELFICEkappa
70 _RL recip_SHELFICElatentHeat
71 _RL SHELFICEDragLinear
72 _RL SHELFICEDragQuadratic
73 _RL SHELFICEthetaSurface
74
75 COMMON /SHELFICE_FIELDS_RL/
76 & shelficeForcingT,
77 & shelficeForcingS
78 _RL shelficeForcingT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
79 _RL shelficeForcingS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
80
81 COMMON /SHELFICE_FIELDS_RS/
82 & R_shelfIce,
83 & shelficeLoadAnomaly,
84 & shelficeHeatFlux,
85 & shelfIceFreshWaterFlux
86 _RS R_shelfIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87 _RS shelficeLoadAnomaly (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88 _RS shelficeHeatFlux (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
89 _RS shelficeFreshWaterFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90
91 LOGICAL SHELFICEisOn
92 LOGICAL useISOMIPTD
93 LOGICAL SHELFICEconserve
94 LOGICAL SHELFICEboundaryLayer
95 LOGICAL no_slip_shelfice
96 LOGICAL SHELFICEwriteState
97 LOGICAL SHELFICE_dump_mdsio
98 LOGICAL SHELFICE_tave_mdsio
99 LOGICAL SHELFICE_dump_mnc
100 LOGICAL SHELFICE_tave_mnc
101 COMMON /SHELFICE_PARMS_L/
102 & SHELFICEisOn,
103 & useISOMIPTD,
104 & SHELFICEconserve,
105 & SHELFICEboundaryLayer,
106 & no_slip_shelfice,
107 & SHELFICEwriteState,
108 & SHELFICE_dump_mdsio,
109 & SHELFICE_tave_mdsio,
110 & SHELFICE_dump_mnc,
111 & SHELFICE_tave_mnc
112
113 CHARACTER*(MAX_LEN_FNAM) SHELFICEloadAnomalyFile
114 COMMON /SHELFICE_PARM_C/
115 & SHELFICEloadAnomalyFile
116
117 #endif /* ALLOW_SHELFICE */

  ViewVC Help
Powered by ViewVC 1.1.22