/[MITgcm]/MITgcm/pkg/fizhi/fizhi_tendency_apply.F
ViewVC logotype

Annotation of /MITgcm/pkg/fizhi/fizhi_tendency_apply.F

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

Revision 1.4 - (hide annotations) (download)
Mon Jun 7 18:11:37 2004 UTC (20 years ago) by molod
Branch: MAIN
CVS Tags: checkpoint53d_post, checkpoint54a_pre, checkpoint54a_post, checkpoint54b_post, checkpoint54, checkpoint53g_post, checkpoint53f_post, checkpoint54c_post
Changes since 1.3: +5 -5 lines 
Accomodate new common block names
1 molod 1.4 C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_tendency_apply.F,v 1.3 2004/05/05 00:39:21 edhill Exp $
2 edhill 1.3 C $Name: $
3    
4 molod 1.1 subroutine fizhi_tendency_apply_u(iMin, iMax, jMin, jMax,
5     . bi,bj,kLev,myTime,myThid)
6     C=======================================================================
7     C Routine: fizhi_tendency_apply_u
8     C Interpolate tendencies from physics grid to dynamics grid and
9     C add fizhi tendency terms to U tendency.
10     C
11     C INPUT:
12     C iMin - Working range of tile for applying forcing.
13     C iMax
14     C jMin
15     C jMax
16     C kLev
17     C
18     C Notes: Routine works for one level at a time
19     C Assumes that U and V tendencies are already on C-Grid
20     C=======================================================================
21     implicit none
22    
23     #include "CPP_OPTIONS.h"
24     #include "SIZE.h"
25     #include "GRID.h"
26     #include "EEPARAMS.h"
27     #include "DYNVARS.h"
28     #include "fizhi_SIZE.h"
29 molod 1.4 #include "fizhi_land_SIZE.h"
30 molod 1.1 #include "fizhi_coms.h"
31    
32     integer iMin, iMax, jMin, jMax, kLev, bi, bj, myThid
33     _RL myTime
34    
35     integer i, j
36    
37     do j=jMin,jMax
38     do i=iMin,iMax
39     gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) +
40     . maskC(i,j,kLev,bi,bj) * guphy(i,j,kLev,bi,bj)
41     enddo
42     enddo
43    
44     return
45     end
46     subroutine fizhi_tendency_apply_v(iMin, iMax, jMin, jMax,
47     . bi,bj,kLev,myTime,myThid)
48     C=======================================================================
49     C Routine: fizhi_tendency_apply_v
50     C Interpolate tendencies from physics grid to dynamics grid and
51     C add fizhi tendency terms to V tendency.
52     C
53     C INPUT:
54     C iMin - Working range of tile for applying forcing.
55     C iMax
56     C jMin
57     C jMax
58     C kLev
59     C
60     C Notes: Routine works for one level at a time
61     C Assumes that U and V tendencies are already on C-Grid
62     C=======================================================================
63     implicit none
64    
65     #include "CPP_OPTIONS.h"
66     #include "SIZE.h"
67     #include "GRID.h"
68     #include "EEPARAMS.h"
69     #include "DYNVARS.h"
70     #include "fizhi_SIZE.h"
71 molod 1.4 #include "fizhi_land_SIZE.h"
72 molod 1.1 #include "fizhi_coms.h"
73    
74     integer iMin, iMax, jMin, jMax, kLev, bi, bj, myThid
75     _RL myTime
76    
77     integer i, j
78    
79     do j=jMin,jMax
80     do i=iMin,iMax
81     gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj) +
82     . maskC(i,j,kLev,bi,bj) * gvphy(i,j,kLev,bi,bj)
83     enddo
84     enddo
85    
86     return
87     end
88     subroutine fizhi_tendency_apply_t(iMin, iMax, jMin, jMax,
89     . bi,bj,kLev,myTime,myThid)
90     C=======================================================================
91     C Routine: fizhi_tendency_apply_t
92     C Interpolate tendencies from physics grid to dynamics grid and
93     C add fizhi tendency terms to T (theta) tendency.
94     C
95     C INPUT:
96     C iMin - Working range of tile for applying forcing.
97     C iMax
98     C jMin
99     C jMax
100     C kLev
101     C
102     C Notes: Routine works for one level at a time
103     C=======================================================================
104     implicit none
105    
106     #include "CPP_OPTIONS.h"
107     #include "SIZE.h"
108     #include "GRID.h"
109     #include "EEPARAMS.h"
110     #include "DYNVARS.h"
111     #include "fizhi_SIZE.h"
112 molod 1.4 #include "fizhi_land_SIZE.h"
113 molod 1.1 #include "fizhi_coms.h"
114    
115     integer iMin, iMax, jMin, jMax, kLev, bi, bj, myThid
116     _RL myTime
117    
118     integer i, j
119    
120 molod 1.2 do j=jMin,jMax
121     do i=iMin,iMax
122 molod 1.1 gT(i,j,kLev,bi,bj) = maskC(i,j,kLev,bi,bj)
123     . *( gT(i,j,kLev,bi,bj) + gthphy(i,j,kLev,bi,bj) )
124     enddo
125     enddo
126    
127     return
128     end
129     subroutine fizhi_tendency_apply_s(iMin, iMax, jMin, jMax,
130     . bi,bj,kLev,myTime,myThid)
131     C=======================================================================
132     C Routine: fizhi_tendency_apply_s
133     C Interpolate tendencies from physics grid to dynamics grid and
134     C add fizhi tendency terms to S tendency.
135     C
136     C INPUT:
137     C iMin - Working range of tile for applying forcing.
138     C iMax
139     C jMin
140     C jMax
141     C kLev
142     C
143     C Notes: Routine works for one level at a time
144     C=======================================================================
145     implicit none
146    
147     #include "CPP_OPTIONS.h"
148     #include "SIZE.h"
149     #include "GRID.h"
150     #include "EEPARAMS.h"
151     #include "DYNVARS.h"
152     #include "fizhi_SIZE.h"
153 molod 1.4 #include "fizhi_land_SIZE.h"
154 molod 1.1 #include "fizhi_coms.h"
155    
156     integer iMin, iMax, jMin, jMax, kLev, bi, bj, myThid
157     _RL myTime
158    
159     integer i, j
160    
161 molod 1.2 do j=jMin,jMax
162     do i=iMin,iMax
163 molod 1.1 gS(i,j,kLev,bi,bj) = maskC(i,j,kLev,bi,bj)
164     . *( gS(i,j,kLev,bi,bj) + gsphy(i,j,kLev,bi,bj) )
165     enddo
166     enddo
167    
168     return
169     end
  ViewVC Help
Powered by ViewVC 1.1.22