/[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.6 by heimbach, Thu Dec 19 14:00:08 2002 UTC revision 1.8 by heimbach, Tue Jun 24 16:08:45 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 81  c--   == end of interface == Line 81  c--   == end of interface ==
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)')
84       &        yadprefix, xx_theta_file(1:il),'.',optimcycle       &        yadmark, xx_theta_file(1:il),'.',optimcycle
85    
86           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
87       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
88       &                         mythid, dummy)       &                         mythid, dummy)
89    
90           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
91    
92  #endif /* ALLOW_THETA0_CONTROL */  #endif /* ALLOW_THETA0_CONTROL */
93    
# Line 96  c--   == end of interface == Line 96  c--   == end of interface ==
96           il=ilnblnk( xx_salt_file )           il=ilnblnk( xx_salt_file )
97           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
98           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
99       &        yadprefix, xx_salt_file(1:il),'.',optimcycle       &        yadmark, xx_salt_file(1:il),'.',optimcycle
100                    
101           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
102       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
103       &                         mythid, dummy)       &                         mythid, dummy)
104    
105           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
106    
107  #endif /* ALLOW_SALT0_CONTROL */  #endif /* ALLOW_SALT0_CONTROL */
108    
# Line 111  c--   == end of interface == Line 111  c--   == end of interface ==
111           il=ilnblnk( xx_hflux_file )           il=ilnblnk( xx_hflux_file )
112           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
113           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
114       &        yadprefix, xx_hflux_file(1:il),'.',optimcycle       &        yadmark, xx_hflux_file(1:il),'.',optimcycle
115    
116           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
117       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
118       &                        mythid, dummy)       &                        mythid, dummy)
119    
120           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
121    
122  #endif /* ALLOW_HFLUX_CONTROL */  #endif /* ALLOW_HFLUX_CONTROL */
123    
# Line 126  c--   == end of interface == Line 126  c--   == end of interface ==
126           il=ilnblnk( xx_sflux_file )           il=ilnblnk( xx_sflux_file )
127           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
128           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
129       &        yadprefix, xx_sflux_file(1:il),'.',optimcycle       &        yadmark, xx_sflux_file(1:il),'.',optimcycle
130    
131           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
132       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
133       &                        mythid, dummy)       &                        mythid, dummy)
134    
135           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
136    
137  #endif /* ALLOW_SFLUX_CONTROL */  #endif /* ALLOW_SFLUX_CONTROL */
138    
# Line 141  c--   == end of interface == Line 141  c--   == end of interface ==
141           il=ilnblnk( xx_tauu_file )           il=ilnblnk( xx_tauu_file )
142           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
143           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
144       &        yadprefix, xx_tauu_file(1:il),'.',optimcycle       &        yadmark, xx_tauu_file(1:il),'.',optimcycle
145    
146           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
147       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
148       &                        mythid, dummy)       &                        mythid, dummy)
149    
150           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
151    
152  #endif /* ALLOW_USTRESS_CONTROL */  #endif /* ALLOW_USTRESS_CONTROL */
153    
# Line 156  c--   == end of interface == Line 156  c--   == end of interface ==
156           il=ilnblnk( xx_tauv_file )           il=ilnblnk( xx_tauv_file )
157           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
158           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
159       &        yadprefix, xx_tauv_file(1:il),'.',optimcycle       &        yadmark, xx_tauv_file(1:il),'.',optimcycle
160    
161           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
162       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
163       &                        mythid, dummy)       &                        mythid, dummy)
164    
165           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
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 )
292           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
293           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
294       &        yadprefix, xx_tr1_file(1:il),'.',optimcycle       &        yadmark, xx_tr1_file(1:il),'.',optimcycle
295                    
296           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
297       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
298       &                         mythid, dummy)       &                         mythid, dummy)
299    
300           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
301    
302  #endif /* ALLOW_TR10_CONTROL */  #endif /* ALLOW_TR10_CONTROL */
303    
# Line 186  c--   == end of interface == Line 306  c--   == end of interface ==
306           il=ilnblnk( xx_sst_file )           il=ilnblnk( xx_sst_file )
307           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
308           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
309       &        yadprefix, xx_sst_file(1:il),'.',optimcycle       &        yadmark, xx_sst_file(1:il),'.',optimcycle
310    
311           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
312       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
313       &                        mythid, dummy)       &                        mythid, dummy)
314    
315           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
316    
317  #endif /* ALLOW_SST0_CONTROL */  #endif /* ALLOW_SST0_CONTROL */
318    
# Line 201  c--   == end of interface == Line 321  c--   == end of interface ==
321           il=ilnblnk( xx_sss_file )           il=ilnblnk( xx_sss_file )
322           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
323           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
324       &        yadprefix, xx_sss_file(1:il),'.',optimcycle       &        yadmark, xx_sss_file(1:il),'.',optimcycle
325    
326           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
327       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
328       &                        mythid, dummy)       &                        mythid, dummy)
329    
330           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
331    
332  #endif /* ALLOW_SSS0_CONTROL */  #endif /* ALLOW_SSS0_CONTROL */
333    
# Line 216  c--   == end of interface == Line 336  c--   == end of interface ==
336           il=ilnblnk( xx_hfacc_file )           il=ilnblnk( xx_hfacc_file )
337           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
338           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
339       &        yadprefix, xx_hfacc_file(1:il),'.',optimcycle       &        yadmark, xx_hfacc_file(1:il),'.',optimcycle
340    
341  #ifdef ALLOW_HFACC3D_CONTROL  #ifdef ALLOW_HFACC3D_CONTROL
342    
# Line 224  c--   == end of interface == Line 344  c--   == end of interface ==
344       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
345       &                        mythid, dummy)       &                        mythid, dummy)
346    
347           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
348    
349  #else  #else
350    
# Line 232  c--   == end of interface == Line 352  c--   == end of interface ==
352       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
353       &                        mythid, dummy)       &                        mythid, dummy)
354    
355           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
356    
357  #endif /* ALLOW_HFACC3D_CONTROL */  #endif /* ALLOW_HFACC3D_CONTROL */
358  #endif /* ALLOW_HFACC_CONTROL */  #endif /* ALLOW_HFACC_CONTROL */
# Line 242  c--   == end of interface == Line 362  c--   == end of interface ==
362           il=ilnblnk( xx_efluxy_file )           il=ilnblnk( xx_efluxy_file )
363           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
364           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
365       &        yadprefix, xx_efluxy_file(1:il),'.',optimcycle       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
   
          call active_read_xyz( 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)')  
      &        yadprefix, xx_efluxy_file(1:il),'.',optimcycle  
366    
367           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
368       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
369       &                         mythid, dummy)       &                         mythid, dummy)
370    
371           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
372    
373  #endif /* ALLOW_EFLUXY0_CONTROL */  #endif /* ALLOW_EFLUXY0_CONTROL */
374    
# Line 272  c--   == end of interface == Line 377  c--   == end of interface ==
377           il=ilnblnk( xx_efluxp_file )           il=ilnblnk( xx_efluxp_file )
378           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
379           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
380       &        yadprefix, xx_efluxp_file(1:il),'.',optimcycle       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
381    
382           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
383       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
384       &                         mythid, dummy)       &                         mythid, dummy)
385    
386           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
387    
388  #endif /* ALLOW_EFLUXP0_CONTROL */  #endif /* ALLOW_EFLUXP0_CONTROL */
389    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22