/[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.11 by heimbach, Fri Jul 18 21:10:16 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 91  c--   Temperature field. Line 92  c--   Temperature field.
92        il=ilnblnk( xx_theta_file )        il=ilnblnk( xx_theta_file )
93        write(fnametheta(1:80),'(2a,i10.10)')        write(fnametheta(1:80),'(2a,i10.10)')
94       &     xx_theta_file(1:il),'.',optimcycle       &     xx_theta_file(1:il),'.',optimcycle
95        call active_read_xyz( fnametheta, tmpfld3d, 1,        call active_read_xyz_loc( fnametheta, tmpfld3d, 1,
96       &                      doglobalread, ladinit, optimcycle,       &                      doglobalread, ladinit, optimcycle,
97       &                      mythid, xx_theta_dummy )       &                      mythid, xx_theta_dummy )
98    
# 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 116  c--   Temperature field. Line 117  c--   Temperature field.
117        il=ilnblnk( xx_salt_file )        il=ilnblnk( xx_salt_file )
118        write(fnamesalt(1:80),'(2a,i10.10)')        write(fnamesalt(1:80),'(2a,i10.10)')
119       &     xx_salt_file(1:il),'.',optimcycle       &     xx_salt_file(1:il),'.',optimcycle
120        call active_read_xyz( fnamesalt, tmpfld3d, 1,        call active_read_xyz_loc( fnamesalt, tmpfld3d, 1,
121       &                      doglobalread, ladinit, optimcycle,       &                      doglobalread, ladinit, optimcycle,
122       &                      mythid, xx_salt_dummy )       &                      mythid, xx_salt_dummy )
123    
# 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 141  c--   Temperature field. Line 140  c--   Temperature field.
140        il=ilnblnk( xx_tr1_file )        il=ilnblnk( xx_tr1_file )
141        write(fnametr1(1:80),'(2a,i10.10)')        write(fnametr1(1:80),'(2a,i10.10)')
142       &     xx_tr1_file(1:il),'.',optimcycle       &     xx_tr1_file(1:il),'.',optimcycle
143        call active_read_xyz( fnametr1, tmpfld3d, 1,        call active_read_xyz_loc( fnametr1, tmpfld3d, 1,
144       &                      doglobalread, ladinit, optimcycle,       &                      doglobalread, ladinit, optimcycle,
145       &                      mythid, xx_tr1_dummy )       &                      mythid, xx_tr1_dummy )
146    
# 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 166  c--   diffkr. Line 163  c--   diffkr.
163        il=ilnblnk( xx_diffkr_file )        il=ilnblnk( xx_diffkr_file )
164        write(fnamediffkr(1:80),'(2a,i10.10)')        write(fnamediffkr(1:80),'(2a,i10.10)')
165       &     xx_diffkr_file(1:il),'.',optimcycle       &     xx_diffkr_file(1:il),'.',optimcycle
166        call active_read_xyz( fnamediffkr, tmpfld3d, 1,        call active_read_xyz_loc( fnamediffkr, tmpfld3d, 1,
167       &                      doglobalread, ladinit, optimcycle,       &                      doglobalread, ladinit, optimcycle,
168       &                      mythid, xx_diffkr_dummy )       &                      mythid, xx_diffkr_dummy )
169        do bj = jtlo,jthi        do bj = jtlo,jthi
# Line 188  c--   kapgm. Line 185  c--   kapgm.
185        il=ilnblnk( xx_kapgm_file )        il=ilnblnk( xx_kapgm_file )
186        write(fnamekapgm(1:80),'(2a,i10.10)')        write(fnamekapgm(1:80),'(2a,i10.10)')
187       &     xx_kapgm_file(1:il),'.',optimcycle       &     xx_kapgm_file(1:il),'.',optimcycle
188        call active_read_xyz( fnamekapgm, tmpfld3d, 1,        call active_read_xyz_loc( fnamekapgm, tmpfld3d, 1,
189       &                      doglobalread, ladinit, optimcycle,       &                      doglobalread, ladinit, optimcycle,
190       &                      mythid, xx_kapgm_dummy )       &                      mythid, xx_kapgm_dummy )
191        do bj = jtlo,jthi        do bj = jtlo,jthi
# Line 210  c--   y-component EP-flux field. Line 207  c--   y-component EP-flux field.
207        il=ilnblnk( xx_efluxy_file )        il=ilnblnk( xx_efluxy_file )
208        write(fnameefluxy(1:80),'(2a,i10.10)')        write(fnameefluxy(1:80),'(2a,i10.10)')
209       &     xx_efluxy_file(1:il),'.',optimcycle       &     xx_efluxy_file(1:il),'.',optimcycle
210        call active_read_xyz( fnameefluxy, tmpfld3d, 1,        call active_read_xyz_loc( fnameefluxy, tmpfld3d, 1,
211       &                      doglobalread, ladinit, optimcycle,       &                      doglobalread, ladinit, optimcycle,
212       &                      mythid, xx_efluxy_dummy )       &                      mythid, xx_efluxy_dummy )
213    
# Line 237  c--   p-component EP-flux field. Line 234  c--   p-component EP-flux field.
234        il=ilnblnk( xx_efluxp_file )        il=ilnblnk( xx_efluxp_file )
235        write(fnameefluxp(1:80),'(2a,i10.10)')        write(fnameefluxp(1:80),'(2a,i10.10)')
236       &     xx_efluxp_file(1:il),'.',optimcycle       &     xx_efluxp_file(1:il),'.',optimcycle
237        call active_read_xyz( fnameefluxp, tmpfld3d, 1,        call active_read_xyz_loc( fnameefluxp, tmpfld3d, 1,
238       &                      doglobalread, ladinit, optimcycle,       &                      doglobalread, ladinit, optimcycle,
239       &                      mythid, xx_efluxp_dummy )       &                      mythid, xx_efluxp_dummy )
240    
# 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_loc ( 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.11
  ViewVC Help
Powered by ViewVC 1.1.22