/[MITgcm]/MITgcm/model/src/ini_cori.F
ViewVC logotype

Contents of /MITgcm/model/src/ini_cori.F

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


Revision 1.17 - (show annotations) (download)
Sun Nov 23 04:19:11 2003 UTC (20 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint52l_pre, hrcube4, hrcube5, checkpoint52d_pre, checkpoint52j_pre, checkpoint54d_post, checkpoint54e_post, checkpoint52l_post, checkpoint52k_post, checkpoint54, checkpoint53, checkpoint52f_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint52e_pre, checkpoint52e_post, checkpoint53d_post, checkpoint54b_post, checkpoint52m_post, checkpoint52b_post, checkpoint52c_post, checkpoint52f_pre, checkpoint54a_pre, checkpoint53c_post, checkpoint54a_post, checkpoint53a_post, checkpoint52d_post, checkpoint53g_post, checkpoint52i_post, checkpoint52h_pre, checkpoint53f_post, checkpoint52j_post, branch-netcdf, checkpoint52n_post, checkpoint53b_pre, checkpoint53b_post, checkpoint53d_pre, checkpoint54c_post
Branch point for: netcdf-sm0
Changes since 1.16: +2 -1 lines
#include "PACKAGES_CONFIG.h" was missing.

1 C $Header: /u/gcmpack/MITgcm/model/src/ini_cori.F,v 1.16 2003/11/03 22:40:22 edhill Exp $
2 C $Name: $
3
4 #include "PACKAGES_CONFIG.h"
5 #include "CPP_OPTIONS.h"
6
7 CBOP
8 C !ROUTINE: INI_CORI
9 C !INTERFACE:
10 SUBROUTINE INI_CORI( myThid )
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | SUBROUTINE INI_CORI
14 C | o Initialise coriolis term.
15 C *==========================================================*
16 C \ev
17
18 C !USES:
19 IMPLICIT NONE
20 C == Global variables ==
21 #include "SIZE.h"
22 #include "EEPARAMS.h"
23 #include "PARAMS.h"
24 #include "GRID.h"
25 #include "DYNVARS.h"
26
27 C !INPUT/OUTPUT PARAMETERS:
28 C == Routine arguments ==
29 C myThid - Number of this instance of INI_CORI
30 INTEGER myThid
31
32 C !LOCAL VARIABLES:
33 C == Local variables ==
34 C bi,bj - Loop counters
35 C I,J,K
36 C facGrid - Factor for grid to meter conversion
37 INTEGER bi, bj
38 INTEGER I, J, K
39 _RL facGrid
40 CEOP
41
42 C-- Initialise coriolis parameter
43 IF ( useConstantF ) THEN
44 C o Constant F case
45 DO bj = myByLo(myThid), myByHi(myThid)
46 DO bi = myBxLo(myThid), myBxHi(myThid)
47 DO K=1,Nr
48 DO J=1-Oly,sNy+Oly
49 DO I=1-Olx,sNx+Olx
50 fCori(i,j,bi,bj)=f0
51 fCoriG(i,j,bi,bj)=f0
52 fCoriCos(i,j,bi,bj)=0.
53 ENDDO
54 ENDDO
55 ENDDO
56 ENDDO
57 ENDDO
58 ELSEIF ( useBetaPlaneF ) THEN
59 C o Beta plane case
60 facGrid = 1. _d 0
61 IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere
62 DO bj = myByLo(myThid), myByHi(myThid)
63 DO bi = myBxLo(myThid), myBxHi(myThid)
64 DO K=1,Nr
65 DO J=1-Oly,sNy+Oly
66 DO I=1-Olx,sNx+Olx
67 fCori(i,j,bi,bj)=f0+beta*_yC(i,j,bi,bj)*facGrid
68 fCoriG(i,j,bi,bj)=f0+beta*yG(i,j,bi,bj)*facGrid
69 fCoriCos(i,j,bi,bj)=0.
70 ENDDO
71 ENDDO
72 ENDDO
73 ENDDO
74 ENDDO
75 ELSEIF ( useSphereF ) THEN
76 C o Spherical case
77 C Note in this case we assume yC is in degrees.
78 DO bj = myByLo(myThid), myByHi(myThid)
79 DO bi = myBxLo(myThid), myBxHi(myThid)
80 DO K=1,Nr
81 DO J=1-Oly,sNy+Oly
82 DO I=1-Olx,sNx+Olx
83 fCori(i,j,bi,bj)=
84 & 2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad)
85 fCoriG(i,j,bi,bj)=
86 & 2. _d 0*omega*sin(yG(i,j,bi,bj)*deg2rad)
87 fCoriCos(i,j,bi,bj)=
88 & 2. _d 0*omega*cos(_yC(i,j,bi,bj)*deg2rad)
89 ENDDO
90 ENDDO
91 ENDDO
92 ENDDO
93 ENDDO
94 ELSE
95 C o Special custom form
96 DO bj = myByLo(myThid), myByHi(myThid)
97 DO bi = myBxLo(myThid), myBxHi(myThid)
98 DO K=1,Nr
99 DO J=1-Oly,sNy+Oly
100 DO I=1-Olx,sNx+Olx
101 fCori(i,j,bi,bj)=0.
102 fCoriG(i,j,bi,bj)=0.
103 fCoriCos(i,j,bi,bj)=0.
104 ENDDO
105 ENDDO
106 ENDDO
107 ENDDO
108 ENDDO
109 ENDIF
110 C
111 c _EXCH_XY_R4(fCori,myThid)
112 c CALL EXCH_Z_XY_RS(fCoriG,myThid)
113 c _EXCH_XY_R4(fCoriCos,myThid)
114
115 #ifdef ALLOW_MONITOR
116 CALL MON_PRINTSTATS_RS(1,fCori,'fCori',myThid)
117 CALL MON_PRINTSTATS_RS(1,fCoriG,'fCoriG',myThid)
118 CALL MON_PRINTSTATS_RS(1,fCoriCos,'fCoriCos',myThid)
119 #endif
120
121 RETURN
122 END

  ViewVC Help
Powered by ViewVC 1.1.22