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

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

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


Revision 1.5 - (hide 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 mlosch 1.5 C $Header: /u/gcmpack/MITgcm/pkg/shelfice/SHELFICE.h,v 1.4 2006/08/14 16:52:46 mlosch Exp $
2 mlosch 1.1 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 mlosch 1.2 C SHELFICEloadAnomalyFile - name of shelfice load anomaly file
20 mlosch 1.1 C SHELFICEDragLinear - linear drag at bottom shelfice (1/s)
21     C SHELFICEDragQuadratic - quadratic drag at bottom shelfice (1/m)
22 mlosch 1.3 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 mlosch 1.1 C SHELFICElatentHeat - latent heat of fusion (J/kg)
27     C useISOMIPTD - use simple ISOMIP thermodynamics
28 mlosch 1.5 C SHELFICEconserve - use conservative form of H&O-thermodynamics
29     C following Jenkins et al. (2001, JPO)
30 mlosch 1.4 C SHELFICEboundaryLayer - turn on vertical merging of cells to for a
31     C boundary layer of drF thickness
32 mlosch 1.1 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 mlosch 1.2 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 mlosch 1.1 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 mlosch 1.3 & SHELFICEheatTransCoeff, SHELFICEsaltTransCoeff,
58     & rhoShelfice, SHELFICEkappa,
59 mlosch 1.1 & SHELFICElatentHeat, recip_SHELFICElatentHeat,
60 mlosch 1.3 & SHELFICEheatCapacity_Cp,
61     & SHELFICEthetaSurface,
62 mlosch 1.1 & SHELFICEDragLinear, SHELFICEDragQuadratic
63     _RL SHELFICE_dumpFreq, SHELFICE_taveFreq
64 mlosch 1.3 _RL SHELFICEheatTransCoeff
65     _RL SHELFICEsaltTransCoeff
66 mlosch 1.1 _RL SHELFICElatentHeat
67 mlosch 1.3 _RL SHELFICEheatCapacity_Cp
68     _RL rhoShelfice
69     _RL SHELFICEkappa
70 mlosch 1.1 _RL recip_SHELFICElatentHeat
71     _RL SHELFICEDragLinear
72     _RL SHELFICEDragQuadratic
73 mlosch 1.3 _RL SHELFICEthetaSurface
74 mlosch 1.1
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 mlosch 1.2 & shelficeLoadAnomaly,
84 mlosch 1.1 & shelficeHeatFlux,
85     & shelfIceFreshWaterFlux
86     _RS R_shelfIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87 mlosch 1.2 _RS shelficeLoadAnomaly (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88 mlosch 1.1 _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 mlosch 1.5 LOGICAL SHELFICEconserve
94 mlosch 1.4 LOGICAL SHELFICEboundaryLayer
95 mlosch 1.1 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 mlosch 1.5 & SHELFICEconserve,
105 mlosch 1.4 & SHELFICEboundaryLayer,
106 mlosch 1.1 & no_slip_shelfice,
107     & SHELFICEwriteState,
108     & SHELFICE_dump_mdsio,
109     & SHELFICE_tave_mdsio,
110     & SHELFICE_dump_mnc,
111     & SHELFICE_tave_mnc
112    
113 mlosch 1.2 CHARACTER*(MAX_LEN_FNAM) SHELFICEloadAnomalyFile
114     COMMON /SHELFICE_PARM_C/
115     & SHELFICEloadAnomalyFile
116    
117 mlosch 1.1 #endif /* ALLOW_SHELFICE */

  ViewVC Help
Powered by ViewVC 1.1.22