/[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.2 by heimbach, Fri Jul 13 13:40:17 2001 UTC revision 1.5 by heimbach, Fri Sep 28 15:15:55 2001 UTC
# Line 2  C $Header$ Line 2  C $Header$
2    
3  #include "CTRL_CPPOPTIONS.h"  #include "CTRL_CPPOPTIONS.h"
4    
5    CBOP
6    C     !ROUTINE: ctrl_map_ini
7    C     !INTERFACE:
8          subroutine ctrl_map_ini( mythid )
9    
10    C     !DESCRIPTION: \bv
11    c     *=================================================================
12    c     | SUBROUTINE ctrl_map_ini
13    c     | Add the temperature, salinity, and diffusivity parts of the
14    c     | control vector to the model state and update the tile halos.
15    c     | The control vector is defined in the header file "ctrl.h".
16    c     *=================================================================
17    C     \ev
18    
19        subroutine ctrl_map_ini(  C     !USES:
      I                     mythid  
      &                   )  
   
 c     ==================================================================  
 c     SUBROUTINE ctrl_map_ini  
 c     ==================================================================  
 c  
 c     o Add the temperature and salinity parts of the control vector to  
 c       the model state and update the tile edges. The control vector is  
 c       defined in the header file "ctrl.h".  
 c  
 c     started: Christian Eckert eckert@mit.edu 30-Jun-1999  
 c  
 c     changed: Christian Eckert eckert@mit.edu 23-Feb-2000  
 c  
 c              - Restructured the code in order to create a package  
 c                for the MITgcmUV.  
 c  
 c     ==================================================================  
 c     SUBROUTINE ctrl_map_ini  
 c     ==================================================================  
   
20        implicit none        implicit none
21    
22  c     == global variables ==  c     == global variables ==
   
23  #include "EEPARAMS.h"  #include "EEPARAMS.h"
24  #include "SIZE.h"  #include "SIZE.h"
25  #include "DYNVARS.h"  #include "DYNVARS.h"
26  #include "TR1.h"  #include "TR1.h"
   
27  #include "ctrl.h"  #include "ctrl.h"
28  #include "ctrl_dummy.h"  #include "ctrl_dummy.h"
29  #include "optim.h"  #include "optim.h"
30    
31    C     !INPUT/OUTPUT PARAMETERS:
32  c     == routine arguments ==  c     == routine arguments ==
   
33        integer mythid        integer mythid
34    
35    C     !LOCAL VARIABLES:
36  c     == local variables ==  c     == local variables ==
37    
       _RL     fac  
38        integer bi,bj        integer bi,bj
39        integer i,j,k        integer i,j,k
40        integer itlo,ithi        integer itlo,ithi
# Line 61  c     == local variables == Line 50  c     == local variables ==
50        character*( 80)   fnametheta        character*( 80)   fnametheta
51        character*( 80)   fnamesalt        character*( 80)   fnamesalt
52        character*( 80)   fnametr1        character*( 80)   fnametr1
53          character*( 80)   fnamediffkr
54          character*( 80)   fnamekapgm
55    
56  c     == external ==        _RL     fac
57    
58    c     == external ==
59        integer  ilnblnk        integer  ilnblnk
60        external ilnblnk        external ilnblnk
61    
62  c     == end of interface ==  c     == end of interface ==
63    CEOP
64    
65        jtlo = mybylo(mythid)        jtlo = mybylo(mythid)
66        jthi = mybyhi(mythid)        jthi = mybyhi(mythid)
# Line 105  c--   Temperature field. Line 98  c--   Temperature field.
98                do i = imin,imax                do i = imin,imax
99                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
100       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
101                  gtNm1(i,j,k,bi,bj) = gtNm1(i,j,k,bi,bj) +  cph                gtNm1(i,j,k,bi,bj) = gtNm1(i,j,k,bi,bj) +
102       &                               fac*tmpfld3d(i,j,k,bi,bj)  cph     &                               fac*tmpfld3d(i,j,k,bi,bj)
103                enddo                enddo
104              enddo              enddo
105            enddo            enddo
# Line 130  c--   Temperature field. Line 123  c--   Temperature field.
123                do i = imin,imax                do i = imin,imax
124                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +
125       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
126                  gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) +  cph                gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) +
127       &                               fac*tmpfld3d(i,j,k,bi,bj)  cph     &                               fac*tmpfld3d(i,j,k,bi,bj)
128                enddo                enddo
129              enddo              enddo
130            enddo            enddo
# Line 155  c--   Temperature field. Line 148  c--   Temperature field.
148                do i = imin,imax                do i = imin,imax
149                  tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) +                  tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) +
150       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
151                  gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) +  cph                gtr1Nm1(i,j,k,bi,bj) = gtr1Nm1(i,j,k,bi,bj) +
152       &                               fac*tmpfld3d(i,j,k,bi,bj)  cph     &                               fac*tmpfld3d(i,j,k,bi,bj)
153                  enddo
154                enddo
155              enddo
156           enddo
157          enddo
158    #endif
159    
160    #ifdef ALLOW_DIFFKR_CONTROL
161    c--   diffkr.
162          il=ilnblnk( xx_diffkr_file )
163          write(fnamediffkr(1:80),'(2a,i10.10)')
164         &     xx_diffkr_file(1:il),'.',optimcycle
165          call active_read_xyz( fnamediffkr, tmpfld3d, 1,
166         &                      doglobalread, ladinit, optimcycle,
167         &                      mythid, xx_diffkr_dummy )
168          do bj = jtlo,jthi
169            do bi = itlo,ithi
170              do k = 1,nr
171                do j = jmin,jmax
172                  do i = imin,imax
173                    diffkr(i,j,k,bi,bj) = diffkr(i,j,k,bi,bj) +
174         &                                tmpfld3d(i,j,k,bi,bj)
175                  enddo
176                enddo
177              enddo
178           enddo
179          enddo
180    #endif
181    
182    #ifdef ALLOW_KAPGM_CONTROL
183    c--   kapgm.
184          il=ilnblnk( xx_kapgm_file )
185          write(fnamekapgm(1:80),'(2a,i10.10)')
186         &     xx_kapgm_file(1:il),'.',optimcycle
187          call active_read_xyz( fnamekapgm, tmpfld3d, 1,
188         &                      doglobalread, ladinit, optimcycle,
189         &                      mythid, xx_kapgm_dummy )
190          do bj = jtlo,jthi
191            do bi = itlo,ithi
192              do k = 1,nr
193                do j = jmin,jmax
194                  do i = imin,imax
195                    kapgm(i,j,k,bi,bj) = kapgm(i,j,k,bi,bj) +
196         &                               tmpfld3d(i,j,k,bi,bj)
197                enddo                enddo
198              enddo              enddo
199            enddo            enddo
# Line 176  c--   Update the tile edges. Line 213  c--   Update the tile edges.
213        _EXCH_XYZ_R8( gsNm1, mythid )        _EXCH_XYZ_R8( gsNm1, mythid )
214  #endif  #endif
215  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
216        _EXCH_XYZ_R8(  tr1, mythid )        _EXCH_XYZ_R8(     tr1, mythid )
217        _EXCH_XYZ_R8( gTr1Nm1, mythid )        _EXCH_XYZ_R8( gTr1Nm1, mythid )
218  #endif  #endif
219    #ifdef ALLOW_DIFFKR_CONTROL
220          _EXCH_XYZ_R8( diffkr, mythid)
221    #endif
222    #ifdef ALLOW_KAPGM_CONTROL
223          _EXCH_XYZ_R8( kapgm, mythid)
224    #endif
225    
226    
227        return        return
228        end        end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22