/[MITgcm]/MITgcm/model/inc/SURFACE.h
ViewVC logotype

Contents of /MITgcm/model/inc/SURFACE.h

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

Revision 1.13 - (show annotations) (download)
Thu Dec 8 15:44:33 2005 UTC (18 years, 5 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint58b_post, checkpoint57y_post, checkpoint58, checkpoint58f_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint58t_post, checkpoint58m_post, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, mitgcm_mapl_00, checkpoint58r_post, checkpoint58n_post, checkpoint58k_post, checkpoint58l_post, checkpoint58g_post, checkpoint58h_post, checkpoint58j_post, checkpoint58i_post, checkpoint58c_post, checkpoint58s_post
Changes since 1.12: +7 -1 lines
File MIME type: text/plain 
First step for a NLFS adjoint
o initially suppress rStar (new flag DISABLE_RSTAR_CODE)
o new init. routines for calc_r_star, calc_surf_dr
o still need to deal with ini_masks_etc
o testreport seemed happy
1 C $Header: /u/gcmpack/MITgcm/model/inc/SURFACE.h,v 1.12 2004/07/06 00:51:35 jmc Exp $
2 C $Name: $
3 C
4 CBOP
5 C !ROUTINE: SURFACE.h
6 C !INTERFACE:
7 C include SURFACE.h
8 C !DESCRIPTION: \bv
9 C *==========================================================*
10 C | SURFACE.h
11 C | o Header file defining surface-related model varaibles
12 C *==========================================================*
13 C | Contains variables relative to the surface position
14 C | that are held fixed in linear free-surface formulation
15 C | but can vary with time with a non-linear free-surface.
16 C *==========================================================*
17 C \ev
18 CEOP
19
20 C-- COMMON /SURF_FIXED/ fixed surface arrays (Real)
21 C Bo_surf :: Boyancy|1/rho [ocean|atmos] at surface level [=g|alpha(p_o)]
22 C recip_Bo :: 1/Bo_surf
23 C topoZ :: topographic height [m] (used mainly for atmosphere)
24 C phi0surf :: starting point for integrating phi_Hyd
25 COMMON /SURF_FIXED/ Bo_surf, recip_Bo, topoZ, phi0surf
26 _RL Bo_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
27 _RL recip_Bo(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
28 _RS topoZ (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
29 _RS phi0surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
30
31 C-- COMMON /SURF_INDEX/ Common block for surface related index
32 C ksurfC :: vertical index of the surface tracer cell
33 C ksurfW :: vertical index of the surface U point
34 C ksurfS :: vertical index of the surface V point
35 C IMPORTANT: ksurfC,W,S = Nr+1 where the fluid column is empty (continent)
36 COMMON /SURF_INDEX/ ksurfC, ksurfW, ksurfS
37 INTEGER ksurfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
38 INTEGER ksurfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
39 INTEGER ksurfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
40
41 #ifdef EXACT_CONSERV
42 C etaHnm1 :: surface r-anomaly, etaH, at previous time level
43 C dEtaHdt :: time derivative of total column height [r_unit/s = w unit]
44 C PmEpR :: keep the fresh water input (=-EmPmR) of the previous time step
45 COMMON /EXACT_ETA_LOCAL/ etaHnm1, dEtaHdt, PmEpR
46 _RL etaHnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
47 _RL dEtaHdt(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
48 _RS PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
49 #endif
50
51 #ifdef NONLIN_FRSURF
52 C-- COMMON /SURF_CHANGE/ transient variables used for Non-Lin Free-Surf
53 C hFac_surfC :: New thickness factor of the surface level
54 C center (Tracer point)
55 C hFac_surfW :: idem, West interface (U point)
56 C hFac_surfS :: idem, South interface (V point)
57 COMMON /SURF_CHANGE/
58 & hFac_surfC, hFac_surfW, hFac_surfS
59 _RS hFac_surfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
60 _RS hFac_surfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
61 _RS hFac_surfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
62
63 C Local variables in common block
64 C Rmin_surf :: minimum r_value of the free surface position
65 C that satisfy the hFacInf criteria
66 COMMON /LOCAL_CALC_SURF_DR/ Rmin_surf
67 _RL Rmin_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
68
69 C-- COMMON /RSTAR_CHANGE/ transient variables used with r* Coordinate
70 C rStarFacC :: = dr/dr* = ratio of r-thickness / r*-thickness = h^n / H
71 C rStarFacW :: same but for West face
72 C rStarFacS :: same but for South face
73 C rStarExpC :: column expansion factor = h^n+1/h^n , Centered
74 C rStarExpW :: column expansion factor = h^n+1/h^n , Western face
75 C rStarExpS :: column expansion factor = h^n+1/h^n , Southern face
76 C rStarDhCDt:: relative time derivative of h_Center = d.eta/dt / H
77 C rStarDhWDt:: relative time derivative of h_West_face (u.point)
78 C rStarDhSDt:: relative time derivative of h_South_face (v.point)
79 COMMON /RSTAR_CHANGE/
80 & rStarFacC, rStarFacW, rStarFacS,
81 & rStarExpC, rStarExpW, rStarExpS,
82 & rStarDhCDt,rStarDhWDt,rStarDhSDt
83 _RL rStarFacC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84 _RL rStarFacW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
85 _RL rStarFacS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86 _RL rStarExpC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87 _RL rStarExpW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88 _RL rStarExpS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
89 _RL rStarDhCDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90 _RL rStarDhWDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91 _RL rStarDhSDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92
93 C-- COMMON /RSTAR_FIXED/ fixed thickness ratio ( r* discretization )
94 C h0FacC :: initial (and fixed in time) hFacC factor
95 C h0FacW :: initial (and fixed in time) hFacW factor
96 C h0FacS :: initial (and fixed in time) hFacS factor
97 COMMON /RSTAR_FIXED/
98 & h0FacC, h0FacW, h0FacS
99 _RS h0FacC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
100 _RS h0FacW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
101 _RS h0FacS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
102
103 #endif /* NONLIN_FRSURF */
  ViewVC Help
Powered by ViewVC 1.1.22