63 |
#endif /* SEAICE_CGRID */ |
#endif /* SEAICE_CGRID */ |
64 |
|
|
65 |
C-- Dynamical variables |
C-- Dynamical variables |
66 |
|
CToM<<< |
67 |
|
#ifdef SEAICE_ITD |
68 |
|
COMMON/SEAICE_DYNVARS_1/AREA,HEFF,HSNOW,UICE,VICE, |
69 |
|
& AREAITD,HEFFITD,HSNOWITD |
70 |
|
#else |
71 |
COMMON/SEAICE_DYNVARS_1/AREA,HEFF,HSNOW,UICE,VICE |
COMMON/SEAICE_DYNVARS_1/AREA,HEFF,HSNOW,UICE,VICE |
72 |
|
#endif |
73 |
|
C>>>ToM |
74 |
_RL AREA (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL AREA (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
75 |
_RL HEFF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL HEFF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
76 |
_RL HSNOW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL HSNOW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
77 |
_RL UICE (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL UICE (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
78 |
_RL VICE (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL VICE (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
79 |
|
CToM<<< |
80 |
|
#ifdef SEAICE_ITD |
81 |
|
_RL AREAITD (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nITD,nSx,nSy) |
82 |
|
_RL HEFFITD (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nITD,nSx,nSy) |
83 |
|
_RL HSNOWITD (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nITD,nSx,nSy) |
84 |
|
#endif |
85 |
|
C>>>ToM |
86 |
|
|
87 |
#ifdef SEAICE_GROWTH_LEGACY |
#ifdef SEAICE_GROWTH_LEGACY |
88 |
COMMON/SEAICE_DYNVARS_2/ areaNm1, hEffNm1 |
COMMON/SEAICE_DYNVARS_2/ areaNm1, hEffNm1 |
93 |
C uIceC :: average of UICE between last two time steps |
C uIceC :: average of UICE between last two time steps |
94 |
C vIceC :: average of VICE between last two time steps |
C vIceC :: average of VICE between last two time steps |
95 |
COMMON/SEAICE_DYNVARS_3/ |
COMMON/SEAICE_DYNVARS_3/ |
96 |
& ETA,ZETA,PRESS, e11, e22, e12, |
& ETA,etaZ,ZETA,PRESS, e11, e22, e12, |
97 |
& DRAGS,DRAGA,FORCEX,FORCEY, |
& FORCEX,FORCEY, |
98 |
& uIceC, vIceC, uIceNm1, vIceNm1 |
& uIceC, vIceC, uIceNm1, vIceNm1 |
99 |
_RL ETA (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ETA (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
100 |
|
_RL etaZ (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
101 |
_RL ZETA (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ZETA (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
102 |
C ice strength/pressure term |
C ice strength/pressure term |
103 |
_RL PRESS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL PRESS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
106 |
_RL e22 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL e22 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
107 |
_RL e12 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL e12 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
108 |
C |
C |
|
_RL DRAGS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RL DRAGA (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
109 |
_RL FORCEX (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL FORCEX (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
110 |
_RL FORCEY (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL FORCEY (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
111 |
_RL uIceC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL uIceC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
113 |
_RL uIceNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL uIceNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
114 |
_RL vIceNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL vIceNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
115 |
|
|
116 |
|
#if (defined (ALLOW_MEAN_SFLUX_COST_CONTRIBUTION) || defined (ALLOW_SSH_GLOBMEAN_COST_CONTRIBUTION)) |
117 |
|
C-- Dynamical variables |
118 |
|
COMMON/SEAICE_DYNVARS_COST/ AREAforAtmFW, frWtrAtm |
119 |
|
_RL AREAforAtmFW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
120 |
|
_RL frWtrAtm (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
121 |
|
#endif |
122 |
|
|
123 |
#ifndef SEAICE_CGRID |
#ifndef SEAICE_CGRID |
124 |
COMMON/SEAICE_DYNVARS_BGRID/ AMASS, DAIRN |
COMMON/SEAICE_DYNVARS_BGRID/ AMASS, DAIRN |
125 |
_RL AMASS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL AMASS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
151 |
C that is, ice due to precipitation or snow |
C that is, ice due to precipitation or snow |
152 |
C frWtrAtm contains freshwater flux from the atmosphere |
C frWtrAtm contains freshwater flux from the atmosphere |
153 |
COMMON/ICEFLUX/ saltWtrIce, frWtrIce |
COMMON/ICEFLUX/ saltWtrIce, frWtrIce |
|
#ifdef ALLOW_MEAN_SFLUX_COST_CONTRIBUTION |
|
|
& , frWtrAtm |
|
|
_RL frWtrAtm (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
#endif |
|
154 |
_RL saltWtrIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL saltWtrIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
155 |
_RL frWtrIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL frWtrIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
156 |
|
|
215 |
#ifdef ALLOW_SEAICE_COST_EXPORT |
#ifdef ALLOW_SEAICE_COST_EXPORT |
216 |
_RL uHeffExportCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL uHeffExportCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
217 |
_RL vHeffExportCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL vHeffExportCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
218 |
|
_RL icevolMeanCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
219 |
COMMON /SEAICE_COST_EXPORT_R/ |
COMMON /SEAICE_COST_EXPORT_R/ |
220 |
& uHeffExportCell, vHeffExportCell |
& uHeffExportCell, vHeffExportCell, |
221 |
|
& icevolMeanCell |
222 |
#endif |
#endif |
223 |
|
|
224 |
C SWFracB :: fraction of surface Short-Wave radiation reaching |
C SWFracB :: fraction of surface Short-Wave radiation reaching |
227 |
COMMON /SEAICE_SW_R/ |
COMMON /SEAICE_SW_R/ |
228 |
& SWFracB |
& SWFracB |
229 |
|
|
230 |
|
#ifdef SEAICE_ALLOW_JFNK |
231 |
|
C diagnostics for the JFNK solver |
232 |
|
INTEGER totalNewtonIters |
233 |
|
INTEGER totalNewtonFails |
234 |
|
INTEGER totalKrylovIters |
235 |
|
INTEGER totalKrylovFails |
236 |
|
INTEGER totalJFNKtimeSteps |
237 |
|
COMMON /SEAICE_JFNK_I/ |
238 |
|
& totalNewtonIters, totalNewtonFails, |
239 |
|
& totalKrylovIters, totalKrylovFails, |
240 |
|
& totalJFNKtimeSteps |
241 |
|
#endif /* SEAICE_ALLOW_JFNK */ |
242 |
|
|
243 |
CEH3 ;;; Local Variables: *** |
CEH3 ;;; Local Variables: *** |
244 |
CEH3 ;;; mode:fortran *** |
CEH3 ;;; mode:fortran *** |
245 |
CEH3 ;;; End: *** |
CEH3 ;;; End: *** |