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

Diff of /MITgcm/pkg/ctrl/ctrl_getobcsw.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.5 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 212  cgg     & Line 212  cgg     &
212    
213          if ( (obcswfirst) .or. (obcswchanged)) then          if ( (obcswfirst) .or. (obcswchanged)) then
214                        
215  cgg(    This is a terribly long way to do it. However, the dimensions don't exactly  cgg(    This is a terribly long way to do it. However, the dimensions do not exactly
216  cgg     match up. I will blame Fortran for the ugliness.  cgg     match up. I will blame Fortran for the ugliness.
217    
218            do bj = jtlo,jthi            do bj = jtlo,jthi
# Line 265  cgg    surface level. This velocity is n Line 265  cgg    surface level. This velocity is n
265  cgg    exactly balance the volume flux, since we are dealing with  cgg    exactly balance the volume flux, since we are dealing with
266  cgg    the baroclinic velocity structure..  cgg    the baroclinic velocity structure..
267                        utop = tmpfldyz(j,k,bi,bj)*                        utop = tmpfldyz(j,k,bi,bj)*
268       &                maskW(i+ip1,j,k,bi,bj) * delZ(k) + utop       &                maskW(i+ip1,j,k,bi,bj) * delR(k) + utop
269  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
270                        if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then                        if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then
271                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro
# Line 273  cgg    Add the barotropic velocity compo Line 273  cgg    Add the barotropic velocity compo
273                      enddo                      enddo
274  cgg    Compute the baroclinic velocity at level 1. Should balance flux.  cgg    Compute the baroclinic velocity at level 1. Should balance flux.
275                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)
276       &                                      - utop / delZ(1)       &                                      - utop / delR(1)
277                    enddo                    enddo
278                  enddo                  enddo
279                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                        utop = tmpfldyz(j,k,bi,bj)*                        utop = tmpfldyz(j,k,bi,bj)*
303       &                maskS(i,j,k,bi,bj) * delZ(k) + utop       &                maskS(i,j,k,bi,bj) * delR(k) + utop
304  cgg    Add the barotropic velocity component.  cgg    Add the barotropic velocity component.
305                        if (maskS(i,j,k,bi,bj) .ne. 0.) then                        if (maskS(i,j,k,bi,bj) .ne. 0.) then
306                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro                          tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro
# 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                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)                      tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj)
311       &                                      - utop / delZ(1)       &                                      - utop / delR(1)
312                    enddo                    enddo
313                  enddo                  enddo
314                enddo                enddo
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22