/[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.2 by heimbach, Fri Jul 13 14:50:46 2001 UTC revision 1.8 by heimbach, Tue Jun 24 16:08:45 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_OBCSN_CONTROL
230          else if ( grdchkvarindex .eq. 11 ) then
231             il=ilnblnk( xx_obcsn_file )
232             write(fname(1:80),'(80a)') ' '
233             write(fname(1:80),'(3a,i10.10)')
234         &        yadmark, xx_obcsn_file(1:il),'.',optimcycle
235    
236             call active_read_xz( fname, tmpfldxz, icvrec,
237         &                        doglobalread, ladinit, optimcycle,
238         &                        mythid, dummy)
239    
240             xx_comp = tmpfldxz( itilepos,layer,itile,jtile )
241    
242    #endif /* ALLOW_OBCSN_CONTROL */
243    
244    #ifdef ALLOW_OBCSS_CONTROL
245          else if ( grdchkvarindex .eq. 12 ) then
246             il=ilnblnk( xx_obcss_file )
247             write(fname(1:80),'(80a)') ' '
248             write(fname(1:80),'(3a,i10.10)')
249         &        yadmark, xx_obcss_file(1:il),'.',optimcycle
250    
251             call active_read_xz( fname, tmpfldxz, icvrec,
252         &                        doglobalread, ladinit, optimcycle,
253         &                        mythid, dummy)
254    
255             xx_comp = tmpfldxz( itilepos,layer,itile,jtile )
256    
257    #endif /* ALLOW_OBCSS_CONTROL */
258    
259    #ifdef ALLOW_OBCSW_CONTROL
260          else if ( grdchkvarindex .eq. 13 ) then
261             il=ilnblnk( xx_obcsw_file )
262             write(fname(1:80),'(80a)') ' '
263             write(fname(1:80),'(3a,i10.10)')
264         &        yadmark, xx_obcsw_file(1:il),'.',optimcycle
265    
266             call active_read_yz( fname, tmpfldyz, icvrec,
267         &                        doglobalread, ladinit, optimcycle,
268         &                        mythid, dummy)
269    
270             xx_comp = tmpfldyz( jtilepos,layer,itile,jtile )
271    
272    #endif /* ALLOW_OBCSW_CONTROL */
273    
274    #ifdef ALLOW_OBCSE_CONTROL
275          else if ( grdchkvarindex .eq. 14 ) then
276             il=ilnblnk( xx_obcse_file )
277             write(fname(1:80),'(80a)') ' '
278             write(fname(1:80),'(3a,i10.10)')
279         &        yadmark, xx_obcse_file(1:il),'.',optimcycle
280    
281             call active_read_yz( fname, tmpfldyz, icvrec,
282         &                        doglobalread, ladinit, optimcycle,
283         &                        mythid, dummy)
284    
285             xx_comp = tmpfldyz( jtilepos,layer,itile,jtile )
286    
287    #endif /* ALLOW_OBCSE_CONTROL */
288    
289  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
290        else if ( grdchkvarindex .eq. 17 ) then        else if ( grdchkvarindex .eq. 17 ) then
291           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_tr1_file )
# Line 208  c--   == end of interface == Line 331  c--   == end of interface ==
331    
332  #endif /* ALLOW_SSS0_CONTROL */  #endif /* ALLOW_SSS0_CONTROL */
333    
334    #ifdef ALLOW_HFACC_CONTROL
335          else if ( grdchkvarindex .eq. 20 ) then
336             il=ilnblnk( xx_hfacc_file )
337             write(fname(1:80),'(80a)') ' '
338             write(fname(1:80),'(3a,i10.10)')
339         &        yadmark, xx_hfacc_file(1:il),'.',optimcycle
340    
341    #ifdef ALLOW_HFACC3D_CONTROL
342    
343             call active_read_xyz( fname, tmpfld3d, icvrec,
344         &                        doglobalread, ladinit, optimcycle,
345         &                        mythid, dummy)
346    
347             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
348    
349    #else
350    
351             call active_read_xy( fname, tmpfld2d, icvrec,
352         &                        doglobalread, ladinit, optimcycle,
353         &                        mythid, dummy)
354    
355             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
356    
357    #endif /* ALLOW_HFACC3D_CONTROL */
358    #endif /* ALLOW_HFACC_CONTROL */
359    
360    #ifdef ALLOW_EFLUXY0_CONTROL
361          else if ( grdchkvarindex .eq. 21 ) then
362             il=ilnblnk( xx_efluxy_file )
363             write(fname(1:80),'(80a)') ' '
364             write(fname(1:80),'(3a,i10.10)')
365         &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
366    
367             call active_read_xyz( fname, tmpfld3d, 1,
368         &                         doglobalread, ladinit, optimcycle,
369         &                         mythid, dummy)
370    
371             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
372    
373    #endif /* ALLOW_EFLUXY0_CONTROL */
374    
375    #ifdef ALLOW_EFLUXP0_CONTROL
376          else if ( grdchkvarindex .eq. 22 ) then
377             il=ilnblnk( xx_efluxp_file )
378             write(fname(1:80),'(80a)') ' '
379             write(fname(1:80),'(3a,i10.10)')
380         &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
381    
382             call active_read_xyz( fname, tmpfld3d, 1,
383         &                         doglobalread, ladinit, optimcycle,
384         &                         mythid, dummy)
385    
386             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
387    
388    #endif /* ALLOW_EFLUXP0_CONTROL */
389    
390        else        else
391  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
392        endif        endif
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22