/[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.8 by heimbach, Tue Jan 21 19:18:35 2003 UTC revision 1.11 by heimbach, Fri Jul 18 21:10:16 2003 UTC
# Line 92  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 103  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 117  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 128  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 142  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 153  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 167  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 189  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 211  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 238  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 267  c--   bottom drag Line 263  c--   bottom drag
263        il=ilnblnk( xx_bottomdrag_file )        il=ilnblnk( xx_bottomdrag_file )
264        write(fnamebottomdrag(1:80),'(2a,i10.10)')        write(fnamebottomdrag(1:80),'(2a,i10.10)')
265       &     xx_bottomdrag_file(1:il),'.',optimcycle       &     xx_bottomdrag_file(1:il),'.',optimcycle
266        call active_read_xy ( fnamebottomdrag, tmpfld2d, 1,        call active_read_xy_loc ( fnamebottomdrag, tmpfld2d, 1,
267       &                      doglobalread, ladinit, optimcycle,       &                      doglobalread, ladinit, optimcycle,
268       &                      mythid, xx_bottomdrag_dummy )       &                      mythid, xx_bottomdrag_dummy )
269        do bj = jtlo,jthi        do bj = jtlo,jthi
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22