/[MITgcm]/MITgcm_contrib/SOSE/BoxAdj/code_ad/ctrl_getobcss.F
ViewVC logotype

Diff of /MITgcm_contrib/SOSE/BoxAdj/code_ad/ctrl_getobcss.F

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

revision 1.1 by mmazloff, Tue Jan 18 19:33:08 2011 UTC revision 1.2 by mmazloff, Wed Jan 19 12:39:13 2011 UTC
# Line 120  cih Line 120  cih
120       &                         (obcsscount0-1)*nobcs+iobcs,       &                         (obcsscount0-1)*nobcs+iobcs,
121       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
122       &                         mythid, xx_obcss_dummy )       &                         mythid, xx_obcss_dummy )
123    CNMA
124    #ifdef CTRL_OBCS_MODES
125    cih
126            if ( optimcycle .ge. 0) then            if ( optimcycle .ge. 0) then
127  cih    If normal velocity        cih    If normal velocity      
128              if (iobcs .eq. 3) then              if (iobcs .eq. 3) then
# Line 150  cih    Compute absolute velocities from Line 152  cih    Compute absolute velocities from
152                                end if                                end if
153                             end do                               end do  
154                             do k = 1,Nr                             do k = 1,Nr
155                                tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)                          tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)
156         &                           *recip_hFacS(i,j+jp1,k,bi,bj)
157                             end do                                                   end do                      
158  cih    End if open boundary.  cih    End if open boundary.
159                          end if                                                  end if                        
# Line 190  cih    Compute absolute velocities from Line 193  cih    Compute absolute velocities from
193                                end if                                end if
194                             end do                             end do
195                             do k = 1,Nr                             do k = 1,Nr
196                                tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)                          tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)
197         &                           *recip_hFacW(i,j,k,bi,bj)
198                             end do                                                     end do                        
199  cih    End if open boundary.  cih    End if open boundary.
200                          end if                                                  end if                        
# Line 202  cih    End if iobcs = 4. Line 206  cih    End if iobcs = 4.
206              end if              end if
207  cih    End if optimcycle > 0 .    cih    End if optimcycle > 0 .  
208           end if           end if
209    CNMA
210    #endif
211  cih  cih
212            do bj = jtlo,jthi            do bj = jtlo,jthi
213              do bi = itlo,ithi              do bi = itlo,ithi
# Line 243  cih Line 249  cih
249       &                         (obcsscount1-1)*nobcs+iobcs,       &                         (obcsscount1-1)*nobcs+iobcs,
250       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
251       &                         mythid, xx_obcss_dummy )       &                         mythid, xx_obcss_dummy )
252    CNMA
253    #ifdef CTRL_OBCS_MODES
254  cih          cih        
255            if ( optimcycle .ge. 0) then                      if ( optimcycle .ge. 0) then          
256  cih    If normal velocity        cih    If normal velocity      
# Line 273  cih    Compute absolute velocities from Line 281  cih    Compute absolute velocities from
281                               end if                               end if
282                             end do                               end do  
283                             do k = 1,Nr                             do k = 1,Nr
284                                tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)                          tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)
285         &                           *recip_hFacS(i,j+jp1,k,bi,bj)
286                             end do                                                   end do                      
287  cih    End if open boundary.  cih    End if open boundary.
288                          end if                                                  end if                        
# Line 313  cih    Compute absolute velocities from Line 322  cih    Compute absolute velocities from
322                               end if                               end if
323                             end do                               end do  
324                             do k = 1,Nr                             do k = 1,Nr
325                                tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)  c                              tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)
326                            tmpfldxz(i,k,bi,bj) = tmpz(k,bi,bj)
327         &                           *recip_hFacW(i,j,k,bi,bj)
328                             end do                                                   end do                      
329  cih    End if open boundary.  cih    End if open boundary.
330                          end if                                                  end if                        
# Line 325  cih    End if iobcs = 4. Line 336  cih    End if iobcs = 4.
336              end if              end if
337  cih    End if optimcycle > 0 .    cih    End if optimcycle > 0 .  
338           end if           end if
339    CNMA
340    #endif
341  cih  cih
342            do bj = jtlo,jthi            do bj = jtlo,jthi
343              do bi = itlo,ithi              do bi = itlo,ithi
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22