/[MITgcm]/MITgcm/pkg/grdchk/grdchk_getadxx.F
ViewVC logotype

Diff of /MITgcm/pkg/grdchk/grdchk_getadxx.F

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

revision 1.3 by heimbach, Thu May 30 22:47:26 2002 UTC revision 1.7 by heimbach, Fri Feb 28 02:34:56 2003 UTC
# Line 73  c--   == end of interface == Line 73  c--   == end of interface ==
73        dumiter      = 0        dumiter      = 0
74        dumtime      = 0. _d 0        dumtime      = 0. _d 0
75    
76        if      ( grdchkvarindex .eq. 1 ) then        if ( grdchkvarindex .eq. 0 ) then
77             STOP 'GRDCHK INDEX 0 NOT ALLOWED'
78    
79  #ifdef ALLOW_THETA0_CONTROL  #ifdef ALLOW_THETA0_CONTROL
80          else if ( grdchkvarindex .eq. 1 ) then
81           il=ilnblnk( xx_theta_file )           il=ilnblnk( xx_theta_file )
82           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
83           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
# Line 163  c--   == end of interface == Line 166  c--   == end of interface ==
166    
167  #endif /* ALLOW_VSTRESS_CONTROL */  #endif /* ALLOW_VSTRESS_CONTROL */
168    
169    #ifdef ALLOW_ATEMP_CONTROL
170          else if ( grdchkvarindex .eq. 7 ) then
171             il=ilnblnk( xx_atemp_file )
172             write(fname(1:80),'(80a)') ' '
173             write(fname(1:80),'(3a,i10.10)')
174         &        yadmark, xx_atemp_file(1:il),'.',optimcycle
175    
176             call active_read_xy( fname, tmpfld2d, icvrec,
177         &                        doglobalread, ladinit, optimcycle,
178         &                        mythid, dummy)
179    
180             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
181    
182    #endif /* ALLOW_ATEMP_CONTROL */
183    
184    #ifdef ALLOW_AQH_CONTROL
185          else if ( grdchkvarindex .eq. 8 ) then
186             il=ilnblnk( xx_aqh_file )
187             write(fname(1:80),'(80a)') ' '
188             write(fname(1:80),'(3a,i10.10)')
189         &        yadmark, xx_aqh_file(1:il),'.',optimcycle
190    
191             call active_read_xy( fname, tmpfld2d, icvrec,
192         &                        doglobalread, ladinit, optimcycle,
193         &                        mythid, dummy)
194    
195             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
196    
197    #endif /* ALLOW_AQH_CONTROL */
198    
199    #ifdef ALLOW_UWIND_CONTROL
200          else if ( grdchkvarindex .eq. 9 ) then
201             il=ilnblnk( xx_uwind_file )
202             write(fname(1:80),'(80a)') ' '
203             write(fname(1:80),'(3a,i10.10)')
204         &        yadmark, xx_uwind_file(1:il),'.',optimcycle
205    
206             call active_read_xy( fname, tmpfld2d, icvrec,
207         &                        doglobalread, ladinit, optimcycle,
208         &                        mythid, dummy)
209    
210             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
211    
212    #endif /* ALLOW_UWIND_CONTROL */
213    
214    #ifdef ALLOW_VWIND_CONTROL
215          else if ( grdchkvarindex .eq. 10 ) then
216             il=ilnblnk( xx_vwind_file )
217             write(fname(1:80),'(80a)') ' '
218             write(fname(1:80),'(3a,i10.10)')
219         &        yadmark, xx_vwind_file(1:il),'.',optimcycle
220    
221             call active_read_xy( fname, tmpfld2d, icvrec,
222         &                        doglobalread, ladinit, optimcycle,
223         &                        mythid, dummy)
224    
225             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
226    
227    #endif /* ALLOW_VWIND_CONTROL */
228    
229  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
230        else if ( grdchkvarindex .eq. 17 ) then        else if ( grdchkvarindex .eq. 17 ) then
231           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_tr1_file )
# Line 234  c--   == end of interface == Line 297  c--   == end of interface ==
297  #endif /* ALLOW_HFACC3D_CONTROL */  #endif /* ALLOW_HFACC3D_CONTROL */
298  #endif /* ALLOW_HFACC_CONTROL */  #endif /* ALLOW_HFACC_CONTROL */
299    
300    #ifdef ALLOW_EFLUXY0_CONTROL
301          else if ( grdchkvarindex .eq. 21 ) then
302             il=ilnblnk( xx_efluxy_file )
303             write(fname(1:80),'(80a)') ' '
304             write(fname(1:80),'(3a,i10.10)')
305         &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
306    
307             call active_read_xyz( fname, tmpfld3d, 1,
308         &                         doglobalread, ladinit, optimcycle,
309         &                         mythid, dummy)
310    
311             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
312    
313    #endif /* ALLOW_EFLUXY0_CONTROL */
314    
315    #ifdef ALLOW_EFLUXP0_CONTROL
316          else if ( grdchkvarindex .eq. 22 ) then
317             il=ilnblnk( xx_efluxp_file )
318             write(fname(1:80),'(80a)') ' '
319             write(fname(1:80),'(3a,i10.10)')
320         &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
321    
322             call active_read_xyz( fname, tmpfld3d, 1,
323         &                         doglobalread, ladinit, optimcycle,
324         &                         mythid, dummy)
325    
326             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
327    
328    #endif /* ALLOW_EFLUXP0_CONTROL */
329    
330        else        else
331  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
332        endif        endif
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22