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

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

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

revision 1.6 by heimbach, Sat Jul 13 02:47:32 2002 UTC revision 1.9 by heimbach, Tue Jun 24 16:07:06 2003 UTC
# Line 56  c     == local variables == Line 56  c     == local variables ==
56        character*( 80)   fnamekapgm        character*( 80)   fnamekapgm
57        character*( 80)   fnameefluxy        character*( 80)   fnameefluxy
58        character*( 80)   fnameefluxp        character*( 80)   fnameefluxp
59          character*( 80)   fnamebottomdrag
60    
61        _RL     fac        _RL     fac
62    
# Line 70  CEOP Line 71  CEOP
71        jthi = mybyhi(mythid)        jthi = mybyhi(mythid)
72        itlo = mybxlo(mythid)        itlo = mybxlo(mythid)
73        ithi = mybxhi(mythid)        ithi = mybxhi(mythid)
74        jmin = 1-oly        jmin = 1
75        jmax = sny+oly        jmax = sny
76        imin = 1-olx        imin = 1
77        imax = snx+olx        imax = snx
78    
79        doglobalread = .false.        doglobalread = .false.
80        ladinit      = .false.        ladinit      = .false.
# Line 102  c--   Temperature field. Line 103  c--   Temperature field.
103                do i = imin,imax                do i = imin,imax
104                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
105       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
106  cph                gtNm1(i,j,k,bi,bj) = gtNm1(i,j,k,bi,bj) +                  if(theta(i,j,k,bi,bj).lt.-2.0)
107  cph     &                               fac*tmpfld3d(i,j,k,bi,bj)       &               theta(i,j,k,bi,bj)= -2.0  
108                enddo                enddo
109              enddo              enddo
110            enddo            enddo
# Line 127  c--   Temperature field. Line 128  c--   Temperature field.
128                do i = imin,imax                do i = imin,imax
129                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +
130       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
 cph                gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) +  
 cph     &                               fac*tmpfld3d(i,j,k,bi,bj)  
131                enddo                enddo
132              enddo              enddo
133            enddo            enddo
# Line 152  c--   Temperature field. Line 151  c--   Temperature field.
151                do i = imin,imax                do i = imin,imax
152                  tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) +                  tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) +
153       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
 cph                gtr1Nm1(i,j,k,bi,bj) = gtr1Nm1(i,j,k,bi,bj) +  
 cph     &                               fac*tmpfld3d(i,j,k,bi,bj)  
154                enddo                enddo
155              enddo              enddo
156            enddo            enddo
# Line 261  cph     & Line 258  cph     &
258        enddo        enddo
259  #endif  #endif
260    
261    #ifdef ALLOW_BOTTOMDRAG_CONTROL
262    c--   bottom drag
263          il=ilnblnk( xx_bottomdrag_file )
264          write(fnamebottomdrag(1:80),'(2a,i10.10)')
265         &     xx_bottomdrag_file(1:il),'.',optimcycle
266          call active_read_xy ( fnamebottomdrag, tmpfld2d, 1,
267         &                      doglobalread, ladinit, optimcycle,
268         &                      mythid, xx_bottomdrag_dummy )
269          do bj = jtlo,jthi
270            do bi = itlo,ithi
271              do j = jmin,jmax
272                do i = imin,imax
273                  bottomdragfld(i,j,bi,bj) = bottomdragfld(i,j,bi,bj)
274         &                                   + tmpfld2d(i,j,bi,bj)
275                enddo
276              enddo
277            enddo
278          enddo
279    #endif
280    
281    
282  c--   Update the tile edges.  c--   Update the tile edges.
283    
# Line 288  cph      _EXCH_XYZ_R8( gTr1Nm1, mythid ) Line 305  cph      _EXCH_XYZ_R8( gTr1Nm1, mythid )
305  #ifdef ALLOW_EFLUXP0_CONTROL  #ifdef ALLOW_EFLUXP0_CONTROL
306        _EXCH_XYZ_R8( EfluxP, mythid )        _EXCH_XYZ_R8( EfluxP, mythid )
307  #endif  #endif
308    #ifdef ALLOW_BOTTOMDRAG_CONTROL
309          _EXCH_XY_R8( bottomdragfld, mythid )
310    #endif
311    
312    
313        return        return
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22