/[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.5 by heimbach, Mon Sep 16 18:11:58 2002 UTC revision 1.9 by heimbach, Fri Jul 18 21:10:16 2003 UTC
# Line 10  C $Header$ Line 10  C $Header$
10       I                       layer,       I                       layer,
11       I                       itilepos,       I                       itilepos,
12       I                       jtilepos,       I                       jtilepos,
13       I                       adxx_comp,       I                       xx_comp,
14       I                       mythid       I                       mythid
15       &                     )       &                     )
16    
# Line 45  c     == routine arguments == Line 45  c     == routine arguments ==
45        integer layer        integer layer
46        integer itilepos        integer itilepos
47        integer jtilepos        integer jtilepos
48        _RL     adxx_comp        _RL     xx_comp
49        integer mythid        integer mythid
50    
51  #ifdef ALLOW_GRADIENT_CHECK  #ifdef ALLOW_GRADIENT_CHECK
# Line 61  c     == local variables == Line 61  c     == local variables ==
61    
62        character*(80) fname        character*(80) fname
63    
64          integer i,j,k
65    
66  c--   == external ==  c--   == external ==
67    
68        integer  ilnblnk        integer  ilnblnk
# Line 83  c--   == end of interface == Line 85  c--   == end of interface ==
85           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
86       &        yadmark, xx_theta_file(1:il),'.',optimcycle       &        yadmark, xx_theta_file(1:il),'.',optimcycle
87    
88           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
89       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
90       &                         mythid, dummy)       &                         mythid, dummy)
91    
92           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
93    
94  #endif /* ALLOW_THETA0_CONTROL */  #endif /* ALLOW_THETA0_CONTROL */
95    
# Line 98  c--   == end of interface == Line 100  c--   == end of interface ==
100           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
101       &        yadmark, xx_salt_file(1:il),'.',optimcycle       &        yadmark, xx_salt_file(1:il),'.',optimcycle
102                    
103           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
104       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
105       &                         mythid, dummy)       &                         mythid, dummy)
106    
107           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
108    
109  #endif /* ALLOW_SALT0_CONTROL */  #endif /* ALLOW_SALT0_CONTROL */
110    
# Line 113  c--   == end of interface == Line 115  c--   == end of interface ==
115           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
116       &        yadmark, xx_hflux_file(1:il),'.',optimcycle       &        yadmark, xx_hflux_file(1:il),'.',optimcycle
117    
118           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
119       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
120       &                        mythid, dummy)       &                        mythid, dummy)
121    
122           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
123    
124  #endif /* ALLOW_HFLUX_CONTROL */  #endif /* ALLOW_HFLUX_CONTROL */
125    
# Line 128  c--   == end of interface == Line 130  c--   == end of interface ==
130           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
131       &        yadmark, xx_sflux_file(1:il),'.',optimcycle       &        yadmark, xx_sflux_file(1:il),'.',optimcycle
132    
133           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
134       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
135       &                        mythid, dummy)       &                        mythid, dummy)
136    
137           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
138    
139  #endif /* ALLOW_SFLUX_CONTROL */  #endif /* ALLOW_SFLUX_CONTROL */
140    
# Line 143  c--   == end of interface == Line 145  c--   == end of interface ==
145           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
146       &        yadmark, xx_tauu_file(1:il),'.',optimcycle       &        yadmark, xx_tauu_file(1:il),'.',optimcycle
147    
148           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
149       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
150       &                        mythid, dummy)       &                        mythid, dummy)
151    
152           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
153    
154  #endif /* ALLOW_USTRESS_CONTROL */  #endif /* ALLOW_USTRESS_CONTROL */
155    
# Line 158  c--   == end of interface == Line 160  c--   == end of interface ==
160           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
161       &        yadmark, xx_tauv_file(1:il),'.',optimcycle       &        yadmark, xx_tauv_file(1:il),'.',optimcycle
162    
163           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
164       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
165       &                        mythid, dummy)       &                        mythid, dummy)
166    
167           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
168    
169  #endif /* ALLOW_VSTRESS_CONTROL */  #endif /* ALLOW_VSTRESS_CONTROL */
170    
171    #ifdef ALLOW_ATEMP_CONTROL
172          else if ( grdchkvarindex .eq. 7 ) then
173             il=ilnblnk( xx_atemp_file )
174             write(fname(1:80),'(80a)') ' '
175             write(fname(1:80),'(3a,i10.10)')
176         &        yadmark, xx_atemp_file(1:il),'.',optimcycle
177    
178             call active_read_xy_loc( fname, tmpfld2d, icvrec,
179         &                        doglobalread, ladinit, optimcycle,
180         &                        mythid, dummy)
181    
182             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
183    
184    #endif /* ALLOW_ATEMP_CONTROL */
185    
186    #ifdef ALLOW_AQH_CONTROL
187          else if ( grdchkvarindex .eq. 8 ) then
188             il=ilnblnk( xx_aqh_file )
189             write(fname(1:80),'(80a)') ' '
190             write(fname(1:80),'(3a,i10.10)')
191         &        yadmark, xx_aqh_file(1:il),'.',optimcycle
192    
193             call active_read_xy_loc( fname, tmpfld2d, icvrec,
194         &                        doglobalread, ladinit, optimcycle,
195         &                        mythid, dummy)
196    
197             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
198    
199    #endif /* ALLOW_AQH_CONTROL */
200    
201    #ifdef ALLOW_UWIND_CONTROL
202          else if ( grdchkvarindex .eq. 9 ) then
203             il=ilnblnk( xx_uwind_file )
204             write(fname(1:80),'(80a)') ' '
205             write(fname(1:80),'(3a,i10.10)')
206         &        yadmark, xx_uwind_file(1:il),'.',optimcycle
207    
208             call active_read_xy_loc( fname, tmpfld2d, icvrec,
209         &                        doglobalread, ladinit, optimcycle,
210         &                        mythid, dummy)
211    
212             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
213    
214    #endif /* ALLOW_UWIND_CONTROL */
215    
216    #ifdef ALLOW_VWIND_CONTROL
217          else if ( grdchkvarindex .eq. 10 ) then
218             il=ilnblnk( xx_vwind_file )
219             write(fname(1:80),'(80a)') ' '
220             write(fname(1:80),'(3a,i10.10)')
221         &        yadmark, xx_vwind_file(1:il),'.',optimcycle
222    
223             call active_read_xy_loc( fname, tmpfld2d, icvrec,
224         &                        doglobalread, ladinit, optimcycle,
225         &                        mythid, dummy)
226    
227             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
228    
229    #endif /* ALLOW_VWIND_CONTROL */
230    
231    #ifdef ALLOW_OBCSN_CONTROL
232          else if ( grdchkvarindex .eq. 11 ) then
233             il=ilnblnk( xx_obcsn_file )
234             write(fname(1:80),'(80a)') ' '
235             write(fname(1:80),'(3a,i10.10)')
236         &        yadmark, xx_obcsn_file(1:il),'.',optimcycle
237    
238             call active_read_xz_loc( fname, tmpfldxz, icvrec,
239         &                        doglobalread, ladinit, optimcycle,
240         &                        mythid, dummy)
241    
242             xx_comp = tmpfldxz( itilepos,layer,itile,jtile )
243    
244    #endif /* ALLOW_OBCSN_CONTROL */
245    
246    #ifdef ALLOW_OBCSS_CONTROL
247          else if ( grdchkvarindex .eq. 12 ) then
248             il=ilnblnk( xx_obcss_file )
249             write(fname(1:80),'(80a)') ' '
250             write(fname(1:80),'(3a,i10.10)')
251         &        yadmark, xx_obcss_file(1:il),'.',optimcycle
252    
253             call active_read_xz_loc( fname, tmpfldxz, icvrec,
254         &                        doglobalread, ladinit, optimcycle,
255         &                        mythid, dummy)
256    
257             xx_comp = tmpfldxz( itilepos,layer,itile,jtile )
258    
259    #endif /* ALLOW_OBCSS_CONTROL */
260    
261    #ifdef ALLOW_OBCSW_CONTROL
262          else if ( grdchkvarindex .eq. 13 ) then
263             il=ilnblnk( xx_obcsw_file )
264             write(fname(1:80),'(80a)') ' '
265             write(fname(1:80),'(3a,i10.10)')
266         &        yadmark, xx_obcsw_file(1:il),'.',optimcycle
267    
268             call active_read_yz_loc( fname, tmpfldyz, icvrec,
269         &                        doglobalread, ladinit, optimcycle,
270         &                        mythid, dummy)
271    
272             xx_comp = tmpfldyz( jtilepos,layer,itile,jtile )
273    
274    #endif /* ALLOW_OBCSW_CONTROL */
275    
276    #ifdef ALLOW_OBCSE_CONTROL
277          else if ( grdchkvarindex .eq. 14 ) then
278             il=ilnblnk( xx_obcse_file )
279             write(fname(1:80),'(80a)') ' '
280             write(fname(1:80),'(3a,i10.10)')
281         &        yadmark, xx_obcse_file(1:il),'.',optimcycle
282    
283             call active_read_yz_loc( fname, tmpfldyz, icvrec,
284         &                        doglobalread, ladinit, optimcycle,
285         &                        mythid, dummy)
286    
287             xx_comp = tmpfldyz( jtilepos,layer,itile,jtile )
288    
289    #endif /* ALLOW_OBCSE_CONTROL */
290    
291  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
292        else if ( grdchkvarindex .eq. 17 ) then        else if ( grdchkvarindex .eq. 17 ) then
293           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_tr1_file )
# Line 173  c--   == end of interface == Line 295  c--   == end of interface ==
295           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
296       &        yadmark, xx_tr1_file(1:il),'.',optimcycle       &        yadmark, xx_tr1_file(1:il),'.',optimcycle
297                    
298           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
299       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
300       &                         mythid, dummy)       &                         mythid, dummy)
301    
302           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
303    
304  #endif /* ALLOW_TR10_CONTROL */  #endif /* ALLOW_TR10_CONTROL */
305    
# Line 188  c--   == end of interface == Line 310  c--   == end of interface ==
310           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
311       &        yadmark, xx_sst_file(1:il),'.',optimcycle       &        yadmark, xx_sst_file(1:il),'.',optimcycle
312    
313           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
314       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
315       &                        mythid, dummy)       &                        mythid, dummy)
316    
317           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
318    
319  #endif /* ALLOW_SST0_CONTROL */  #endif /* ALLOW_SST0_CONTROL */
320    
# Line 203  c--   == end of interface == Line 325  c--   == end of interface ==
325           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
326       &        yadmark, xx_sss_file(1:il),'.',optimcycle       &        yadmark, xx_sss_file(1:il),'.',optimcycle
327    
328           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
329       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
330       &                        mythid, dummy)       &                        mythid, dummy)
331    
332           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
333    
334  #endif /* ALLOW_SSS0_CONTROL */  #endif /* ALLOW_SSS0_CONTROL */
335    
# Line 220  c--   == end of interface == Line 342  c--   == end of interface ==
342    
343  #ifdef ALLOW_HFACC3D_CONTROL  #ifdef ALLOW_HFACC3D_CONTROL
344    
345           call active_read_xyz( fname, tmpfld3d, icvrec,           call active_read_xyz_loc( fname, tmpfld3d, icvrec,
346       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
347       &                        mythid, dummy)       &                        mythid, dummy)
348    
349           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
350    
351  #else  #else
352    
353           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
354       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
355       &                        mythid, dummy)       &                        mythid, dummy)
356    
357           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
358    
359  #endif /* ALLOW_HFACC3D_CONTROL */  #endif /* ALLOW_HFACC3D_CONTROL */
360  #endif /* ALLOW_HFACC_CONTROL */  #endif /* ALLOW_HFACC_CONTROL */
# Line 244  c--   == end of interface == Line 366  c--   == end of interface ==
366           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
367       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
368    
369           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
      &                         doglobalread, ladinit, optimcycle,  
      &                         mythid, dummy)  
   
          adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )  
   
 #endif /* ALLOW_THETA0_CONTROL */  
   
 #ifdef ALLOW_EFLUXY0_CONTROL  
       else if ( grdchkvarindex .eq. 21 ) then  
          il=ilnblnk( xx_efluxy_file )  
          write(fname(1:80),'(80a)') ' '  
          write(fname(1:80),'(3a,i10.10)')  
      &        yadmark, xx_efluxy_file(1:il),'.',optimcycle  
   
          call active_read_xyz( fname, tmpfld3d, 1,  
370       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
371       &                         mythid, dummy)       &                         mythid, dummy)
372    
373           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
374    
375  #endif /* ALLOW_EFLUXY0_CONTROL */  #endif /* ALLOW_EFLUXY0_CONTROL */
376    
# Line 274  c--   == end of interface == Line 381  c--   == end of interface ==
381           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
382       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
383    
384           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
385       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
386       &                         mythid, dummy)       &                         mythid, dummy)
387    
388           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
389    
390  #endif /* ALLOW_EFLUXP0_CONTROL */  #endif /* ALLOW_EFLUXP0_CONTROL */
391    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22