/[MITgcm]/MITgcm/pkg/seaice/SEAICE_PARAMS.h
ViewVC logotype

Annotation of /MITgcm/pkg/seaice/SEAICE_PARAMS.h

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

Revision 1.5 - (hide annotations) (download)
Tue Feb 18 05:33:55 2003 UTC (21 years, 4 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint48f_post, checkpoint48g_post
Changes since 1.4: +24 -16 lines
File MIME type: text/plain 
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
1 heimbach 1.2 C $Header:
2    
3     #ifdef ALLOW_SEAICE
4    
5     C /==========================================================\
6     C | SEAICE_PARAMS.h |
7     C | o Basic parameter header for sea ice model. |
8     C \==========================================================/
9    
10     C-- COMMON /SEAICE_PARM_L/ Logical parameters of sea ice model.
11     C
12 dimitri 1.4 C SEAICEwriteState - If true, write sea ice state to file;
13     C default is false.
14     C
15     C SEAICEuseDYNAMICS - If false, do not use dynamics;
16     C default is to use dynamics.
17     C
18     C SEAICEuseLSR - If true, use LSR rather than ADI solver;
19     C default is ADI solver.
20     C
21     C SEAICEwindOnCgrid - By default wind files, uwind and vwind,
22     C are defined on North-East B-grid U and V points.
23     C When this flag is set, wind files are defined on
24     C South-West C-grid U and V points.
25     C
26     LOGICAL SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseLSR,
27     & SEAICEwindOnCgrid
28 heimbach 1.2 COMMON /SEAICE_PARM_L/
29 dimitri 1.4 & SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseLSR,
30     & SEAICEwindOnCgrid
31 heimbach 1.2
32     C-- COMMON /SEAICE_PARM_I/ Integer valued parameters of sea ice model.
33 dimitri 1.3 C LAD - time stepping used for sea-ice advection:
34     C 1 = LEAPFROG, 2 = BACKWARD EULER.
35 heimbach 1.2 C IMAX_TICE - number of iterations for ice heat budget 10
36     C NPSEUDO - number of pseudo-timesteps used in dynsolver
37     C DO PSEUDO-TIMESTEPS TO OBTAIN AN ACCURATE VISCOUS-PLASTIC SOLUTION
38     C 10 PSEUDO-TIMESTEPS OR MORE ARE SUGGESTED FOR HIGH-RESOLUTION (~10KM)
39     C 1 PSEUDO-TIMESTEP CAN BE USED FOR LOW-RESOLUTION GLOBAL MODELING
40     C
41     INTEGER LAD, IMAX_TICE, NPSEUDO
42     COMMON /SEAICE_PARM_I/ LAD, IMAX_TICE, NPSEUDO
43    
44     C-- COMMON /SEAICE_PARM_C/ Character valued sea ice model parameters.
45 dimitri 1.5 C uwindFile - File containing uwind
46     C vwindFile - File containing vwind
47     C atempFile - File containing atemp
48     C aqhFile - File containing aqh
49     C lwdownFile - File containing lwdown
50     C swdownFile - File containing swdown
51     C precipFile - File containing precip
52     C evapFile - File containing evap
53     C runoffFile - File containing runoffF
54     C HeffFile - File containing initial sea-ice thickness
55     C !!! NOTE !!! By default, initial sea-ice conditions are
56     C HEFF=1 and AREA=1. This means that there can be large
57     C artificial fluxes of heat and freshwater in the surface layer
58     C during the first model time step. To avoid this it is a good
59     C idea to specify initial HEFF.
60 heimbach 1.2 C
61 dimitri 1.4 CHARACTER*(MAX_LEN_FNAM) uwindFile
62     CHARACTER*(MAX_LEN_FNAM) vwindFile
63     CHARACTER*(MAX_LEN_FNAM) atempFile
64     CHARACTER*(MAX_LEN_FNAM) aqhFile
65 dimitri 1.5 CHARACTER*(MAX_LEN_FNAM) lwdownFile
66     CHARACTER*(MAX_LEN_FNAM) swdownFile
67 dimitri 1.4 CHARACTER*(MAX_LEN_FNAM) precipFile
68 heimbach 1.2 CHARACTER*(MAX_LEN_FNAM) evapFile
69     CHARACTER*(MAX_LEN_FNAM) runoffFile
70 dimitri 1.5 CHARACTER*(MAX_LEN_FNAM) HeffFile
71 dimitri 1.4 COMMON /SEAICE_PARM_C/ uwindFile, vwindFile, atempFile, aqhFile,
72 dimitri 1.5 & lwdownFile, swdownFile, precipFile, evapFile, runoffFile,
73     & HeffFile
74 heimbach 1.2
75     C-- COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model.
76     C SEAICE_dumpFreq - SEAICE dump frequency. (s)
77     C SEAICE_taveFreq - SEAICE time-averaging frequency. (s)
78     C SEAICE_drag - air-ice drag coefficient
79     C SEAICE_waterDrag - water-ice drag coefficient * water density
80     C SEAICE_albedo - summer albedo
81     C SEAICE_drySnowAlb - dry snow albedo
82     C SEAICE_wetSnowAlb - wet snow albedo
83     C SEAICE_waterAlbedo - water albedo
84     C SEAICE_strength - sea-ice strength
85     C SEAICE_sensHeat - sensible heat transfer coefficient
86     C SEAICE_latentWater - latent heat transfer coefficient for water
87     C SEAICE_latentIce - latent heat transfer coefficient for ice
88     C SEAICE_iceConduct - sea-ice conductivity
89     C SEAICE_snowConduct - snow conductivity
90     C SEAICE_emissivity - Stefan-Boltzman constant * emissivity
91     C SEAICE_snowThick - cutoff snow thickness
92     C SEAICE_shortwave - penetration shortwave radiation factor
93     C LSR_ERROR - sets accuracy of LSR solver
94     C DIFF1 - parameter used in advect.F
95     C A22 - parameter used in growth.F
96     C HO - parameter used in growth.F
97     C WindForcingStart - Time of first wind forcing record (s)
98     C WindForcingEnd - Time of last wind forcing record (s)
99     C WindForcingPeriod - Period between wind forcing records (s)
100     C FluxForcingStart - Time of first flux forcing record (s)
101     C FluxForcingEnd - Time of last flux forcing record (s)
102     C FluxForcingPeriod - Period between flux forcing records (s)
103     C SSTForcingStart - Time of first SST forcing record (s)
104     C SSTForcingEnd - Time of last SST forcing record (s)
105     C SSTForcingPeriod - Period between SST forcing records (s)
106     C SSSForcingStart - Time of first SSS forcing record (s)
107     C SSSForcingEnd - Time of last SSS forcing record (s)
108     C SSSForcingPeriod - Period between SSS forcing records (s)
109     C StartingYear - Starting year of integration
110     C EndingYear - Ending year of integration
111     C
112     _RL SEAICE_dumpFreq, SEAICE_taveFreq
113     _RL SEAICE_drag, SEAICE_waterDrag, SEAICE_albedo
114     _RL SEAICE_drySnowAlb, SEAICE_wetSnowAlb
115     _RL SEAICE_waterAlbedo, SEAICE_strength
116     _RL SEAICE_sensHeat, SEAICE_latentWater, SEAICE_latentIce
117     _RL SEAICE_iceConduct, SEAICE_snowConduct, SEAICE_emissivity
118     _RL SEAICE_snowThick, SEAICE_shortwave
119     _RL LSR_ERROR,DIFF1,A22,HO
120     _RL WindForcingStart, WindForcingEnd, WindForcingPeriod
121     _RL FluxForcingStart, FluxForcingEnd, FluxForcingPeriod
122     _RL SSTForcingStart, SSTForcingEnd, SSTForcingPeriod
123     _RL SSSForcingStart, SSSForcingEnd, SSSForcingPeriod
124     _RL StartingYear, EndingYear
125     COMMON /SEAICE_PARM_RL/
126     & SEAICE_dumpFreq, SEAICE_taveFreq,
127     & SEAICE_drag, SEAICE_waterDrag, SEAICE_albedo,
128     & SEAICE_drySnowAlb, SEAICE_wetSnowAlb,
129     & SEAICE_waterAlbedo, SEAICE_strength,
130     & SEAICE_sensHeat, SEAICE_latentWater, SEAICE_latentIce,
131     & SEAICE_iceConduct, SEAICE_snowConduct, SEAICE_emissivity,
132     & SEAICE_snowThick, SEAICE_shortwave,
133     & LSR_ERROR, DIFF1, A22, HO,
134     & WindForcingStart, WindForcingEnd, WindForcingPeriod,
135     & FluxForcingStart, FluxForcingEnd, FluxForcingPeriod,
136     & SSTForcingStart, SSTForcingEnd, SSTForcingPeriod,
137     & SSSForcingStart, SSSForcingEnd, SSSForcingPeriod,
138     & StartingYear, EndingYear
139    
140     C-- COMMON /SEAICE_BOUND_RL/ Various bounding values
141     C MAX_HEFF - maximum ice thickness (m)
142 dimitri 1.4 C MIN_ATEMP - minimum air temperature (deg C)
143 dimitri 1.5 C MIN_LWDOWN - minimum downward longwave (W/m^2)
144 heimbach 1.2 C MAX_TICE - maximum ice temperature (deg C)
145     C MIN_TICE - minimum ice temperature (deg C)
146     C
147 dimitri 1.5 _RL MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE
148 heimbach 1.2 COMMON /SEAICE_BOUND_RL/
149 dimitri 1.5 & MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE
150 dimitri 1.3
151     C-- Seaice time step
152     COMMON /SEAICE_R/ SEAICE_DT
153     _RL SEAICE_DT
154    
155     C-- Constants used by sea-ice model
156     _RL ZERO , ONE , TWO
157     parameter ( ZERO = 0.0 _d 0, ONE = 1.0 _d 0, TWO = 2.0 _d 0 )
158     _RL QUART , HALF
159     parameter ( QUART = 0.25 _d 0, HALF = 0.5 _d 0 )
160 heimbach 1.2
161 dimitri 1.5 #endif /* ALLOW_SEAICE */
  ViewVC Help
Powered by ViewVC 1.1.22