/[MITgcm]/MITgcm_contrib/torge/itd/code/SEAICE.h
ViewVC logotype

Diff of /MITgcm_contrib/torge/itd/code/SEAICE.h

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

revision 1.1 by dimitri, Fri Apr 27 22:22:17 2012 UTC revision 1.6 by torge, Fri May 3 18:59:39 2013 UTC
# Line 63  C                    for metric terms in Line 63  C                    for metric terms in
63  #endif /* SEAICE_CGRID */  #endif /* SEAICE_CGRID */
64    
65  C--   Dynamical variables  C--   Dynamical variables
66    #ifdef SEAICE_ITD
67          COMMON/SEAICE_DYNVARS_1/AREA,HEFF,HSNOW,UICE,VICE,
68         &                        AREAITD,HEFFITD,HSNOWITD
69    #else
70        COMMON/SEAICE_DYNVARS_1/AREA,HEFF,HSNOW,UICE,VICE        COMMON/SEAICE_DYNVARS_1/AREA,HEFF,HSNOW,UICE,VICE
71    #endif
72        _RL AREA       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL AREA       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
73        _RL HEFF       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL HEFF       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
74        _RL HSNOW      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL HSNOW      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
75        _RL UICE       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL UICE       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
76        _RL VICE       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL VICE       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
77    #ifdef SEAICE_ITD
78  #ifdef SEAICE_GROWTH_LEGACY        _RL AREAITD    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nITD,nSx,nSy)
79        COMMON/SEAICE_DYNVARS_2/ areaNm1, hEffNm1        _RL HEFFITD    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nITD,nSx,nSy)
80        _RL areaNm1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL HSNOWITD   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nITD,nSx,nSy)
81        _RL hEffNm1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  #endif
 #endif /* SEAICE_GROWTH_LEGACY */  
82    
83  C     uIceC :: average of UICE between last two time steps  C     uIceC :: average of UICE between last two time steps
84  C     vIceC :: average of VICE between last two time steps  C     vIceC :: average of VICE between last two time steps
85        COMMON/SEAICE_DYNVARS_3/        COMMON/SEAICE_DYNVARS_3/
86       &     ETA,ZETA,PRESS, e11, e22, e12,       &     ETA,etaZ,ZETA,PRESS, e11, e22, e12,
87       &     DRAGS,DRAGA,FORCEX,FORCEY,       &     FORCEX,FORCEY,
88       &     uIceC, vIceC, uIceNm1, vIceNm1       &     uIceC, vIceC, uIceNm1, vIceNm1
89        _RL ETA        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL ETA        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90          _RL etaZ       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91        _RL ZETA       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL ZETA       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92  C     ice strength/pressure term  C     ice strength/pressure term
93        _RL PRESS      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL PRESS      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 91  C     strain rate tensor Line 96  C     strain rate tensor
96        _RL e22        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL e22        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
97        _RL e12        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL e12        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
98  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)  
99        _RL FORCEX     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL FORCEX     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
100        _RL FORCEY     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL FORCEY     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
101        _RL uIceC      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL uIceC      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 100  C Line 103  C
103        _RL uIceNm1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL uIceNm1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
104        _RL vIceNm1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL vIceNm1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
105    
106    #if (defined (ALLOW_MEAN_SFLUX_COST_CONTRIBUTION) || defined (ALLOW_SSH_GLOBMEAN_COST_CONTRIBUTION))
107    C--   Dynamical variables
108          COMMON/SEAICE_DYNVARS_COST/ AREAforAtmFW, frWtrAtm
109          _RL AREAforAtmFW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
110          _RL frWtrAtm     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
111    #endif
112    
113  #ifndef SEAICE_CGRID  #ifndef SEAICE_CGRID
114        COMMON/SEAICE_DYNVARS_BGRID/ AMASS, DAIRN        COMMON/SEAICE_DYNVARS_BGRID/ AMASS, DAIRN
115        _RL AMASS      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL AMASS      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 131  C     frWtrIce contains m of freshwater Line 141  C     frWtrIce contains m of freshwater
141  C              that is, ice due to precipitation or snow  C              that is, ice due to precipitation or snow
142  C     frWtrAtm contains freshwater flux from the atmosphere  C     frWtrAtm contains freshwater flux from the atmosphere
143        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  
144        _RL saltWtrIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL saltWtrIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
145        _RL frWtrIce   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL frWtrIce   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
146    
# Line 199  C--   KGEO    Level used as a proxy for Line 205  C--   KGEO    Level used as a proxy for
205  #ifdef ALLOW_SEAICE_COST_EXPORT  #ifdef ALLOW_SEAICE_COST_EXPORT
206        _RL uHeffExportCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL uHeffExportCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
207        _RL vHeffExportCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL vHeffExportCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
208          _RL icevolMeanCell(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
209        COMMON /SEAICE_COST_EXPORT_R/        COMMON /SEAICE_COST_EXPORT_R/
210       &       uHeffExportCell, vHeffExportCell       &       uHeffExportCell, vHeffExportCell,
211         &       icevolMeanCell
212  #endif  #endif
213    
214  C     SWFracB :: fraction of surface Short-Wave radiation reaching  C     SWFracB :: fraction of surface Short-Wave radiation reaching
# Line 209  C                the bottom of ocean sur Line 217  C                the bottom of ocean sur
217        COMMON /SEAICE_SW_R/        COMMON /SEAICE_SW_R/
218       &       SWFracB       &       SWFracB
219    
220    #ifdef SEAICE_ALLOW_JFNK
221    C     diagnostics for the JFNK solver
222          INTEGER totalNewtonIters
223          INTEGER totalNewtonFails
224          INTEGER totalKrylovIters
225          INTEGER totalKrylovFails
226          INTEGER totalJFNKtimeSteps
227          COMMON /SEAICE_JFNK_I/
228         &     totalNewtonIters, totalNewtonFails,
229         &     totalKrylovIters, totalKrylovFails,
230         &     totalJFNKtimeSteps
231          INTEGER nVec
232          PARAMETER ( nVec=2*sNx*sNy )
233          _RL scalarProductMetric( nVec, 1, nSx, nSy )
234          COMMON /SEAICE_JFNK_RL/ scalarProductMetric
235    #endif /* SEAICE_ALLOW_JFNK */
236    
237  CEH3 ;;; Local Variables: ***  CEH3 ;;; Local Variables: ***
238  CEH3 ;;; mode:fortran ***  CEH3 ;;; mode:fortran ***
239  CEH3 ;;; End: ***  CEH3 ;;; End: ***
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22