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

Diff of /MITgcm/pkg/ctrl/ctrl_getobcss.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 149  cgg    surface level. This velocity is n Line 149  cgg    surface level. This velocity is n
149  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
150  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
151                        vtop = tmpfldxz(i,k,bi,bj)*                        vtop = tmpfldxz(i,k,bi,bj)*
152       &                maskS(i,j+jp1,k,bi,bj) * delZ(k) + vtop       &                maskS(i,j+jp1,k,bi,bj) * delR(k) + vtop
153  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
154                        if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then                        if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
155                          tmpfldxz(i,k,bi,bj) = tmpfldxz(i,k,bi,bj)+ vbaro                          tmpfldxz(i,k,bi,bj) = tmpfldxz(i,k,bi,bj)+ vbaro
# Line 157  cgg    Add the barotropic velocity compo Line 157  cgg    Add the barotropic velocity compo
157                      enddo                      enddo
158  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
159                      tmpfldxz(i,1,bi,bj) = tmpfldxz(i,1,bi,bj)                      tmpfldxz(i,1,bi,bj) = tmpfldxz(i,1,bi,bj)
160       &                                      - vtop / delZ(1)       &                                      - vtop / delR(1)
161                    enddo                    enddo
162                  enddo                  enddo
163                enddo                enddo
# Line 185  cgg    surface level. This velocity is n Line 185  cgg    surface level. This velocity is n
185  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
186  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
187                        vtop = tmpfldxz(i,k,bi,bj)*                        vtop = tmpfldxz(i,k,bi,bj)*
188       &                maskW(i,j,k,bi,bj) * delZ(k) + vtop       &                maskW(i,j,k,bi,bj) * delR(k) + vtop
189  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
190                        if (maskW(i,j,k,bi,bj) .ne. 0.) then                        if (maskW(i,j,k,bi,bj) .ne. 0.) then
191                          tmpfldxz(i,k,bi,bj) = tmpfldxz(i,k,bi,bj)+ vbaro                          tmpfldxz(i,k,bi,bj) = tmpfldxz(i,k,bi,bj)+ vbaro
# Line 193  cgg    Add the barotropic velocity compo Line 193  cgg    Add the barotropic velocity compo
193                      enddo                      enddo
194  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
195                      tmpfldxz(i,1,bi,bj) = tmpfldxz(i,1,bi,bj)                      tmpfldxz(i,1,bi,bj) = tmpfldxz(i,1,bi,bj)
196       &                                      - vtop / delZ(1)       &                                      - vtop / delR(1)
197                    enddo                    enddo
198                  enddo                  enddo
199                enddo                enddo
# Line 264  cgg    surface level. This velocity is n Line 264  cgg    surface level. This velocity is n
264  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
265  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
266                        vtop = tmpfldxz(i,k,bi,bj)*                        vtop = tmpfldxz(i,k,bi,bj)*
267       &                maskS(i,j+jp1,k,bi,bj) * delZ(k) + vtop       &                maskS(i,j+jp1,k,bi,bj) * delR(k) + vtop
268  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
269                        if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then                        if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
270                          tmpfldxz(i,k,bi,bj) = tmpfldxz(i,k,bi,bj)+ vbaro                          tmpfldxz(i,k,bi,bj) = tmpfldxz(i,k,bi,bj)+ vbaro
# Line 272  cgg    Add the barotropic velocity compo Line 272  cgg    Add the barotropic velocity compo
272                      enddo                      enddo
273  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
274                      tmpfldxz(i,1,bi,bj) = tmpfldxz(i,1,bi,bj)                      tmpfldxz(i,1,bi,bj) = tmpfldxz(i,1,bi,bj)
275       &                                      - vtop / delZ(1)       &                                      - vtop / delR(1)
276                    enddo                    enddo
277                  enddo                  enddo
278                enddo                enddo
# Line 300  cgg    surface level. This velocity is n Line 300  cgg    surface level. This velocity is n
300  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
301  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
302                        vtop = tmpfldxz(i,k,bi,bj)*                        vtop = tmpfldxz(i,k,bi,bj)*
303       &                maskW(i,j,k,bi,bj) * delZ(k) + vtop       &                maskW(i,j,k,bi,bj) * delR(k) + vtop
304  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
305                        if (maskW(i,j,k,bi,bj) .ne. 0.) then                        if (maskW(i,j,k,bi,bj) .ne. 0.) then
306                          tmpfldxz(i,k,bi,bj) = tmpfldxz(i,k,bi,bj)+ vbaro                          tmpfldxz(i,k,bi,bj) = tmpfldxz(i,k,bi,bj)+ vbaro
# Line 308  cgg    Add the barotropic velocity compo Line 308  cgg    Add the barotropic velocity compo
308                      enddo                      enddo
309  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
310                      tmpfldxz(i,1,bi,bj) = tmpfldxz(i,1,bi,bj)                      tmpfldxz(i,1,bi,bj) = tmpfldxz(i,1,bi,bj)
311       &                                      - vtop / delZ(1)       &                                      - vtop / delR(1)
312                    enddo                    enddo
313                  enddo                  enddo
314                enddo                enddo
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22