/[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.16 by heimbach, Wed Jun 7 01:55:14 2006 UTC revision 1.17 by heimbach, Fri Oct 27 05:16:54 2006 UTC
# Line 88  c--   == end of interface == Line 88  c--   == end of interface ==
88           call active_read_xyz_loc( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
89       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
90       &                         mythid, dummy)       &                         mythid, dummy)
   
91           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
   
92  #endif /* ALLOW_THETA0_CONTROL */  #endif /* ALLOW_THETA0_CONTROL */
93    
94  #ifdef ALLOW_SALT0_CONTROL  #ifdef ALLOW_SALT0_CONTROL
# Line 99  c--   == end of interface == Line 97  c--   == end of interface ==
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       &        yadmark, xx_salt_file(1:il),'.',optimcycle       &        yadmark, xx_salt_file(1:il),'.',optimcycle
           
100           call active_read_xyz_loc( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
101       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
102       &                         mythid, dummy)       &                         mythid, dummy)
   
103           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
104    
105  #endif /* ALLOW_SALT0_CONTROL */  #endif /* ALLOW_SALT0_CONTROL */
# Line 114  c--   == end of interface == Line 110  c--   == end of interface ==
110           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
111           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
112       &        yadmark, xx_hflux_file(1:il),'.',optimcycle       &        yadmark, xx_hflux_file(1:il),'.',optimcycle
   
113           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
114       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
115       &                        mythid, dummy)       &                        mythid, dummy)
   
116           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
117  #endif /* ALLOW_HFLUX_CONTROL */  #endif /* ALLOW_HFLUX_CONTROL */
118    
119  #ifdef ALLOW_SFLUX_CONTROL  #ifdef ALLOW_SFLUX_CONTROL
# Line 129  c--   == end of interface == Line 122  c--   == end of interface ==
122           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
123           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
124       &        yadmark, xx_sflux_file(1:il),'.',optimcycle       &        yadmark, xx_sflux_file(1:il),'.',optimcycle
   
125           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
126       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
127       &                        mythid, dummy)       &                        mythid, dummy)
   
128           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
129  #endif /* ALLOW_SFLUX_CONTROL */  #endif /* ALLOW_SFLUX_CONTROL */
130    
131  #ifdef ALLOW_USTRESS_CONTROL  #ifdef ALLOW_USTRESS_CONTROL
# Line 144  c--   == end of interface == Line 134  c--   == end of interface ==
134           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
135           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
136       &        yadmark, xx_tauu_file(1:il),'.',optimcycle       &        yadmark, xx_tauu_file(1:il),'.',optimcycle
   
137           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
138       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
139       &                        mythid, dummy)       &                        mythid, dummy)
   
140           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
141  #endif /* ALLOW_USTRESS_CONTROL */  #endif /* ALLOW_USTRESS_CONTROL */
142    
143  #ifdef ALLOW_VSTRESS_CONTROL  #ifdef ALLOW_VSTRESS_CONTROL
# Line 159  c--   == end of interface == Line 146  c--   == end of interface ==
146           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
147           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
148       &        yadmark, xx_tauv_file(1:il),'.',optimcycle       &        yadmark, xx_tauv_file(1:il),'.',optimcycle
   
149           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
150       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
151       &                        mythid, dummy)       &                        mythid, dummy)
   
152           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
153  #endif /* ALLOW_VSTRESS_CONTROL */  #endif /* ALLOW_VSTRESS_CONTROL */
154    
155  #ifdef ALLOW_ATEMP_CONTROL  #ifdef ALLOW_ATEMP_CONTROL
# Line 174  c--   == end of interface == Line 158  c--   == end of interface ==
158           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
159           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
160       &        yadmark, xx_atemp_file(1:il),'.',optimcycle       &        yadmark, xx_atemp_file(1:il),'.',optimcycle
   
161           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
162       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
163       &                        mythid, dummy)       &                        mythid, dummy)
   
164           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
165  #endif /* ALLOW_ATEMP_CONTROL */  #endif /* ALLOW_ATEMP_CONTROL */
166    
167  #ifdef ALLOW_AQH_CONTROL  #ifdef ALLOW_AQH_CONTROL
# Line 189  c--   == end of interface == Line 170  c--   == end of interface ==
170           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
171           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
172       &        yadmark, xx_aqh_file(1:il),'.',optimcycle       &        yadmark, xx_aqh_file(1:il),'.',optimcycle
   
173           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
174       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
175       &                        mythid, dummy)       &                        mythid, dummy)
   
176           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
177  #endif /* ALLOW_AQH_CONTROL */  #endif /* ALLOW_AQH_CONTROL */
178    
179  #ifdef ALLOW_UWIND_CONTROL  #ifdef ALLOW_UWIND_CONTROL
# Line 204  c--   == end of interface == Line 182  c--   == end of interface ==
182           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
183           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
184       &        yadmark, xx_uwind_file(1:il),'.',optimcycle       &        yadmark, xx_uwind_file(1:il),'.',optimcycle
   
185           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
186       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
187       &                        mythid, dummy)       &                        mythid, dummy)
   
188           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
189  #endif /* ALLOW_UWIND_CONTROL */  #endif /* ALLOW_UWIND_CONTROL */
190    
191  #ifdef ALLOW_VWIND_CONTROL  #ifdef ALLOW_VWIND_CONTROL
# Line 219  c--   == end of interface == Line 194  c--   == end of interface ==
194           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
195           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
196       &        yadmark, xx_vwind_file(1:il),'.',optimcycle       &        yadmark, xx_vwind_file(1:il),'.',optimcycle
   
197           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
198       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
199       &                        mythid, dummy)       &                        mythid, dummy)
   
200           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
201  #endif /* ALLOW_VWIND_CONTROL */  #endif /* ALLOW_VWIND_CONTROL */
202    
203  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
# Line 339  c--   == end of interface == Line 311  c--   == end of interface ==
311           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
312           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
313       &        yadmark, xx_sst_file(1:il),'.',optimcycle       &        yadmark, xx_sst_file(1:il),'.',optimcycle
   
314           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
315       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
316       &                        mythid, dummy)       &                        mythid, dummy)
   
317           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
318  #endif /* ALLOW_SST0_CONTROL */  #endif /* ALLOW_SST0_CONTROL */
319    
320  #if (defined (ALLOW_SSS_CONTROL) || defined (ALLOW_SSS0_CONTROL))  #if (defined (ALLOW_SSS_CONTROL) || defined (ALLOW_SSS0_CONTROL))
# Line 354  c--   == end of interface == Line 323  c--   == end of interface ==
323           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
324           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
325       &        yadmark, xx_sss_file(1:il),'.',optimcycle       &        yadmark, xx_sss_file(1:il),'.',optimcycle
   
326           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
327       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
328       &                        mythid, dummy)       &                        mythid, dummy)
   
329           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
330  #endif /* ALLOW_SSS0_CONTROL */  #endif /* ALLOW_SSS0_CONTROL */
331    
332  #ifdef ALLOW_DEPTH_CONTROL  #ifdef ALLOW_DEPTH_CONTROL
# Line 369  c--   == end of interface == Line 335  c--   == end of interface ==
335           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
336           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
337       &        yadmark, xx_depth_file(1:il),'.',optimcycle       &        yadmark, xx_depth_file(1:il),'.',optimcycle
   
338           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
339       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
340       &                        mythid, dummy)       &                        mythid, dummy)
   
341           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
342  #endif /* ALLOW_DEPTH_CONTROL */  #endif /* ALLOW_DEPTH_CONTROL */
343    
344  #ifdef ALLOW_EFLUXY0_CONTROL  #ifdef ALLOW_EFLUXY0_CONTROL
# Line 414  c--   == end of interface == Line 377  c--   == end of interface ==
377           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
378           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
379       &        yadmark, xx_precip_file(1:il),'.',optimcycle       &        yadmark, xx_precip_file(1:il),'.',optimcycle
   
380           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
381       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
382       &                        mythid, dummy)       &                        mythid, dummy)
   
383           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
384  #endif /* ALLOW_PRECIP_CONTROL */  #endif /* ALLOW_PRECIP_CONTROL */
385    
386  #ifdef ALLOW_SWFLUX_CONTROL  #ifdef ALLOW_SWFLUX_CONTROL
# Line 429  c--   == end of interface == Line 389  c--   == end of interface ==
389           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
390           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
391       &        yadmark, xx_swflux_file(1:il),'.',optimcycle       &        yadmark, xx_swflux_file(1:il),'.',optimcycle
   
392           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
393       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
394       &                        mythid, dummy)       &                        mythid, dummy)
   
395           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
396  #endif /* ALLOW_SWFLUX_CONTROL */  #endif /* ALLOW_SWFLUX_CONTROL */
397    
398  #ifdef ALLOW_SWDOWN_CONTROL  #ifdef ALLOW_SWDOWN_CONTROL
# Line 444  c--   == end of interface == Line 401  c--   == end of interface ==
401           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
402           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
403       &        yadmark, xx_swdown_file(1:il),'.',optimcycle       &        yadmark, xx_swdown_file(1:il),'.',optimcycle
404             call active_read_xy_loc( fname, tmpfld2d, icvrec,
405         &                        doglobalread, ladinit, optimcycle,
406         &                        mythid, dummy)
407             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
408    #endif /* ALLOW_SWDOWN_CONTROL */
409    
410    #ifdef ALLOW_LWFLUX_CONTROL
411          else if ( grdchkvarindex .eq. 35 ) then
412             il=ilnblnk( xx_lwflux_file )
413             write(fname(1:80),'(80a)') ' '
414             write(fname(1:80),'(3a,i10.10)')
415         &        yadmark, xx_lwflux_file(1:il),'.',optimcycle
416           call active_read_xy_loc( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
417       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
418       &                        mythid, dummy)       &                        mythid, dummy)
419             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
420    #endif /* ALLOW_LWFLUX_CONTROL */
421    
422    #ifdef ALLOW_LWDOWN_CONTROL
423          else if ( grdchkvarindex .eq. 36 ) then
424             il=ilnblnk( xx_lwdown_file )
425             write(fname(1:80),'(80a)') ' '
426             write(fname(1:80),'(3a,i10.10)')
427         &        yadmark, xx_lwdown_file(1:il),'.',optimcycle
428             call active_read_xy_loc( fname, tmpfld2d, icvrec,
429         &                        doglobalread, ladinit, optimcycle,
430         &                        mythid, dummy)
431           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
432    #endif /* ALLOW_LWDOWN_CONTROL */
433    
434  #endif /* ALLOW_SWDOWN_CONTROL */  #ifdef ALLOW_EVAP_CONTROL
435          else if ( grdchkvarindex .eq. 37 ) then
436             il=ilnblnk( xx_evap_file )
437             write(fname(1:80),'(80a)') ' '
438             write(fname(1:80),'(3a,i10.10)')
439         &        yadmark, xx_evap_file(1:il),'.',optimcycle
440             call active_read_xy_loc( fname, tmpfld2d, icvrec,
441         &                        doglobalread, ladinit, optimcycle,
442         &                        mythid, dummy)
443             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
444    #endif /* ALLOW_EVAP_CONTROL */
445    
446    #ifdef ALLOW_SNOWPRECIP_CONTROL
447          else if ( grdchkvarindex .eq. 38 ) then
448             il=ilnblnk( xx_snowprecip_file )
449             write(fname(1:80),'(80a)') ' '
450             write(fname(1:80),'(3a,i10.10)')
451         &        yadmark, xx_snowprecip_file(1:il),'.',optimcycle
452             call active_read_xy_loc( fname, tmpfld2d, icvrec,
453         &                        doglobalread, ladinit, optimcycle,
454         &                        mythid, dummy)
455             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
456    #endif /* ALLOW_SNOWPRECIP_CONTROL */
457    
458    #ifdef ALLOW_APRESSURE_CONTROL
459          else if ( grdchkvarindex .eq. 39 ) then
460             il=ilnblnk( xx_apressure_file )
461             write(fname(1:80),'(80a)') ' '
462             write(fname(1:80),'(3a,i10.10)')
463         &        yadmark, xx_apressure_file(1:il),'.',optimcycle
464    
465             call active_read_xy_loc( fname, tmpfld2d, icvrec,
466         &                        doglobalread, ladinit, optimcycle,
467         &                        mythid, dummy)
468             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
469    #endif /* ALLOW_APRESSURE_CONTROL */
470    
471    #ifdef ALLOW_RUNOFF_CONTROL
472          else if ( grdchkvarindex .eq. 40 ) then
473             il=ilnblnk( xx_runoff_file )
474             write(fname(1:80),'(80a)') ' '
475             write(fname(1:80),'(3a,i10.10)')
476         &        yadmark, xx_runoff_file(1:il),'.',optimcycle
477             call active_read_xy_loc( fname, tmpfld2d, icvrec,
478         &                        doglobalread, ladinit, optimcycle,
479         &                        mythid, dummy)
480             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
481    #endif /* ALLOW_RUNOFF_CONTROL */
482    
483        else        else
484  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22