/[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.8 by heimbach, Tue Jun 24 16:08:45 2003 UTC revision 1.15 by heimbach, Fri May 12 02:17:03 2006 UTC
# Line 48  c     == routine arguments == Line 48  c     == routine arguments ==
48        _RL     xx_comp        _RL     xx_comp
49        integer mythid        integer mythid
50    
51  #ifdef ALLOW_GRADIENT_CHECK  #ifdef ALLOW_GRDCHK
52  c     == local variables ==  c     == local variables ==
53    
54        integer il        integer il
# 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    
# 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    
# 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    
# 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    
# 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    
# 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    
# Line 173  c--   == end of interface == Line 175  c--   == end of interface ==
175           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
176       &        yadmark, xx_atemp_file(1:il),'.',optimcycle       &        yadmark, xx_atemp_file(1:il),'.',optimcycle
177    
178           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
179       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
180       &                        mythid, dummy)       &                        mythid, dummy)
181    
# Line 188  c--   == end of interface == Line 190  c--   == end of interface ==
190           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
191       &        yadmark, xx_aqh_file(1:il),'.',optimcycle       &        yadmark, xx_aqh_file(1:il),'.',optimcycle
192    
193           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
194       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
195       &                        mythid, dummy)       &                        mythid, dummy)
196    
# Line 203  c--   == end of interface == Line 205  c--   == end of interface ==
205           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
206       &        yadmark, xx_uwind_file(1:il),'.',optimcycle       &        yadmark, xx_uwind_file(1:il),'.',optimcycle
207    
208           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
209       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
210       &                        mythid, dummy)       &                        mythid, dummy)
211    
# Line 218  c--   == end of interface == Line 220  c--   == end of interface ==
220           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
221       &        yadmark, xx_vwind_file(1:il),'.',optimcycle       &        yadmark, xx_vwind_file(1:il),'.',optimcycle
222    
223           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
224       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
225       &                        mythid, dummy)       &                        mythid, dummy)
226    
# Line 233  c--   == end of interface == Line 235  c--   == end of interface ==
235           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
236       &        yadmark, xx_obcsn_file(1:il),'.',optimcycle       &        yadmark, xx_obcsn_file(1:il),'.',optimcycle
237    
238           call active_read_xz( fname, tmpfldxz, icvrec,           call active_read_xz_loc( fname, tmpfldxz, icvrec,
239       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
240       &                        mythid, dummy)       &                        mythid, dummy)
241    
# Line 248  c--   == end of interface == Line 250  c--   == end of interface ==
250           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
251       &        yadmark, xx_obcss_file(1:il),'.',optimcycle       &        yadmark, xx_obcss_file(1:il),'.',optimcycle
252    
253           call active_read_xz( fname, tmpfldxz, icvrec,           call active_read_xz_loc( fname, tmpfldxz, icvrec,
254       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
255       &                        mythid, dummy)       &                        mythid, dummy)
256    
# Line 263  c--   == end of interface == Line 265  c--   == end of interface ==
265           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
266       &        yadmark, xx_obcsw_file(1:il),'.',optimcycle       &        yadmark, xx_obcsw_file(1:il),'.',optimcycle
267    
268           call active_read_yz( fname, tmpfldyz, icvrec,           call active_read_yz_loc( fname, tmpfldyz, icvrec,
269       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
270       &                        mythid, dummy)       &                        mythid, dummy)
271    
# Line 278  c--   == end of interface == Line 280  c--   == end of interface ==
280           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
281       &        yadmark, xx_obcse_file(1:il),'.',optimcycle       &        yadmark, xx_obcse_file(1:il),'.',optimcycle
282    
283           call active_read_yz( fname, tmpfldyz, icvrec,           call active_read_yz_loc( fname, tmpfldyz, icvrec,
284       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
285       &                        mythid, dummy)       &                        mythid, dummy)
286    
# Line 286  c--   == end of interface == Line 288  c--   == end of interface ==
288    
289  #endif /* ALLOW_OBCSE_CONTROL */  #endif /* ALLOW_OBCSE_CONTROL */
290    
291    #ifdef ALLOW_DIFFKR_CONTROL
292          else if ( grdchkvarindex .eq. 15 ) then
293             il=ilnblnk( xx_diffkr_file )
294             write(fname(1:80),'(80a)') ' '
295             write(fname(1:80),'(3a,i10.10)')
296         &        yadmark, xx_diffkr_file(1:il),'.',optimcycle
297    
298             call active_read_xyz_loc( fname, tmpfld3d, 1,
299         &                         doglobalread, ladinit, optimcycle,
300         &                         mythid, dummy)
301    
302             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
303    
304    #endif /* ALLOW_DIFFKR_CONTROL */
305    
306    #ifdef ALLOW_KAPGM_CONTROL
307          else if ( grdchkvarindex .eq. 16 ) then
308             il=ilnblnk( xx_kapgm_file )
309             write(fname(1:80),'(80a)') ' '
310             write(fname(1:80),'(3a,i10.10)')
311         &        yadmark, xx_kapgm_file(1:il),'.',optimcycle
312    
313             call active_read_xyz_loc( fname, tmpfld3d, 1,
314         &                         doglobalread, ladinit, optimcycle,
315         &                         mythid, dummy)
316    
317             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
318    
319    #endif /* ALLOW_KAPGM_CONTROL */
320    
321  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
322        else if ( grdchkvarindex .eq. 17 ) then        else if ( grdchkvarindex .eq. 17 ) then
323           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_tr1_file )
# Line 293  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_tr1_file(1:il),'.',optimcycle       &        yadmark, xx_tr1_file(1:il),'.',optimcycle
327                    
328           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
329       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
330       &                         mythid, dummy)       &                         mythid, dummy)
331    
# Line 301  c--   == end of interface == Line 333  c--   == end of interface ==
333    
334  #endif /* ALLOW_TR10_CONTROL */  #endif /* ALLOW_TR10_CONTROL */
335    
336  #ifdef ALLOW_SST0_CONTROL  #if (defined (ALLOW_SST_CONTROL) || defined (ALLOW_SST0_CONTROL))
337        else if ( grdchkvarindex .eq. 18 ) then        else if ( grdchkvarindex .eq. 18 ) then
338           il=ilnblnk( xx_sst_file )           il=ilnblnk( xx_sst_file )
339           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
340           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
341       &        yadmark, xx_sst_file(1:il),'.',optimcycle       &        yadmark, xx_sst_file(1:il),'.',optimcycle
342    
343           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
344       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
345       &                        mythid, dummy)       &                        mythid, dummy)
346    
# Line 316  c--   == end of interface == Line 348  c--   == end of interface ==
348    
349  #endif /* ALLOW_SST0_CONTROL */  #endif /* ALLOW_SST0_CONTROL */
350    
351  #ifdef ALLOW_SSS0_CONTROL  #if (defined (ALLOW_SSS_CONTROL) || defined (ALLOW_SSS0_CONTROL))
352        else if ( grdchkvarindex .eq. 19 ) then        else if ( grdchkvarindex .eq. 19 ) then
353           il=ilnblnk( xx_sss_file )           il=ilnblnk( xx_sss_file )
354           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
355           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
356       &        yadmark, xx_sss_file(1:il),'.',optimcycle       &        yadmark, xx_sss_file(1:il),'.',optimcycle
357    
358           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
359       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
360       &                        mythid, dummy)       &                        mythid, dummy)
361    
# Line 340  c--   == end of interface == Line 372  c--   == end of interface ==
372    
373  #ifdef ALLOW_HFACC3D_CONTROL  #ifdef ALLOW_HFACC3D_CONTROL
374    
375           call active_read_xyz( fname, tmpfld3d, icvrec,           call active_read_xyz_loc( fname, tmpfld3d, icvrec,
376       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
377       &                        mythid, dummy)       &                        mythid, dummy)
378    
# Line 348  c--   == end of interface == Line 380  c--   == end of interface ==
380    
381  #else  #else
382    
383           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
384       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
385       &                        mythid, dummy)       &                        mythid, dummy)
386    
# Line 364  c--   == end of interface == Line 396  c--   == end of interface ==
396           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
397       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
398    
399           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
400       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
401       &                         mythid, dummy)       &                         mythid, dummy)
402    
# Line 379  c--   == end of interface == Line 411  c--   == end of interface ==
411           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
412       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
413    
414           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
415       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
416       &                         mythid, dummy)       &                         mythid, dummy)
417    
# Line 387  c--   == end of interface == Line 419  c--   == end of interface ==
419    
420  #endif /* ALLOW_EFLUXP0_CONTROL */  #endif /* ALLOW_EFLUXP0_CONTROL */
421    
422    #ifdef ALLOW_PRECIP_CONTROL
423          else if ( grdchkvarindex .eq. 32 ) then
424             il=ilnblnk( xx_precip_file )
425             write(fname(1:80),'(80a)') ' '
426             write(fname(1:80),'(3a,i10.10)')
427         &        yadmark, xx_precip_file(1:il),'.',optimcycle
428    
429             call active_read_xy_loc( fname, tmpfld2d, icvrec,
430         &                        doglobalread, ladinit, optimcycle,
431         &                        mythid, dummy)
432    
433             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
434    
435    #endif /* ALLOW_PRECIP_CONTROL */
436    
437    #ifdef ALLOW_SWFLUX_CONTROL
438          else if ( grdchkvarindex .eq. 33 ) then
439             il=ilnblnk( xx_swflux_file )
440             write(fname(1:80),'(80a)') ' '
441             write(fname(1:80),'(3a,i10.10)')
442         &        yadmark, xx_swflux_file(1:il),'.',optimcycle
443    
444             call active_read_xy_loc( fname, tmpfld2d, icvrec,
445         &                        doglobalread, ladinit, optimcycle,
446         &                        mythid, dummy)
447    
448             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
449    
450    #endif /* ALLOW_SWFLUX_CONTROL */
451    
452    #ifdef ALLOW_SWDOWN_CONTROL
453          else if ( grdchkvarindex .eq. 34 ) then
454             il=ilnblnk( xx_swdown_file )
455             write(fname(1:80),'(80a)') ' '
456             write(fname(1:80),'(3a,i10.10)')
457         &        yadmark, xx_swdown_file(1:il),'.',optimcycle
458    
459             call active_read_xy_loc( fname, tmpfld2d, icvrec,
460         &                        doglobalread, ladinit, optimcycle,
461         &                        mythid, dummy)
462    
463             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
464    
465    #endif /* ALLOW_SWDOWN_CONTROL */
466    
467        else        else
468  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
469        endif        endif
470    
471  #endif /* ALLOW_GRADIENT_CHECK */  #endif /* ALLOW_GRDCHK */
472    
473        end        end
474    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22