/[MITgcm]/MITgcm/pkg/ctrl/ctrl_getobcse.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_getobcse.F

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

revision 1.3 by heimbach, Fri Jul 18 21:10:16 2003 UTC revision 1.4 by heimbach, Sun Oct 26 00:58:03 2003 UTC
# Line 148  cgg    surface level. This velocity is n Line 148  cgg    surface level. This velocity is n
148  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
149  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
150                        utop = tmpfldyz(j,k,bi,bj)*                        utop = tmpfldyz(j,k,bi,bj)*
151       &                maskW(i+ip1,j,k,bi,bj) * delZ(k) + utop       &                maskW(i+ip1,j,k,bi,bj) * delR(k) + utop
152  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
153                        if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then                        if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then
154                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro
# Line 156  cgg    Add the barotropic velocity compo Line 156  cgg    Add the barotropic velocity compo
156                      enddo                      enddo
157  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
158                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)
159       &                                      - utop / delZ(1)       &                                      - utop / delR(1)
160                    enddo                    enddo
161                  enddo                  enddo
162                enddo                enddo
# Line 183  cgg    surface level. This velocity is n Line 183  cgg    surface level. This velocity is n
183  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
184  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
185                        utop = tmpfldyz(j,k,bi,bj)*                        utop = tmpfldyz(j,k,bi,bj)*
186       &                maskS(i,j,k,bi,bj) * delZ(k) + utop       &                maskS(i,j,k,bi,bj) * delR(k) + utop
187  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
188                        if (maskS(i,j,k,bi,bj) .ne. 0.) then                        if (maskS(i,j,k,bi,bj) .ne. 0.) then
189                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro
# Line 191  cgg    Add the barotropic velocity compo Line 191  cgg    Add the barotropic velocity compo
191                      enddo                      enddo
192  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
193                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)
194       &                                      - utop / delZ(1)       &                                      - utop / delR(1)
195                    enddo                    enddo
196                  enddo                  enddo
197                enddo                enddo
# Line 262  cgg    surface level. This velocity is n Line 262  cgg    surface level. This velocity is n
262  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
263  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
264                        utop = tmpfldyz(j,k,bi,bj)*                        utop = tmpfldyz(j,k,bi,bj)*
265       &                maskW(i+ip1,j,k,bi,bj) * delZ(k) + utop       &                maskW(i+ip1,j,k,bi,bj) * delR(k) + utop
266  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
267                        if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then                        if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then
268                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro
# Line 270  cgg    Add the barotropic velocity compo Line 270  cgg    Add the barotropic velocity compo
270                      enddo                      enddo
271  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
272                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)
273       &                                      - utop / delZ(1)       &                                      - utop / delR(1)
274                    enddo                    enddo
275                  enddo                  enddo
276                enddo                enddo
# Line 297  cgg    surface level. This velocity is n Line 297  cgg    surface level. This velocity is n
297  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
298  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
299                        utop = tmpfldyz(j,k,bi,bj)*                        utop = tmpfldyz(j,k,bi,bj)*
300       &                maskS(i,j,k,bi,bj) * delZ(k) + utop       &                maskS(i,j,k,bi,bj) * delR(k) + utop
301  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
302                        if (maskS(i,j,k,bi,bj) .ne. 0.) then                        if (maskS(i,j,k,bi,bj) .ne. 0.) then
303                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro
# Line 305  cgg    Add the barotropic velocity compo Line 305  cgg    Add the barotropic velocity compo
305                      enddo                      enddo
306  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
307                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)
308       &                                      - utop / delZ(1)       &                                      - utop / delR(1)
309                    enddo                    enddo
310                  enddo                  enddo
311                enddo                enddo
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22