/[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.3 by heimbach, Thu May 30 22:47:26 2002 UTC revision 1.5 by heimbach, Mon Sep 16 18:11:58 2002 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                       xx_comp,       I                       adxx_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     xx_comp        _RL     adxx_comp
49        integer mythid        integer mythid
50    
51  #ifdef ALLOW_GRADIENT_CHECK  #ifdef ALLOW_GRADIENT_CHECK
# 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 84  c--   == end of interface == Line 87  c--   == end of interface ==
87       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
88       &                         mythid, dummy)       &                         mythid, dummy)
89    
90           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
91    
92  #endif /* ALLOW_THETA0_CONTROL */  #endif /* ALLOW_THETA0_CONTROL */
93    
# Line 99  c--   == end of interface == Line 102  c--   == end of interface ==
102       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
103       &                         mythid, dummy)       &                         mythid, dummy)
104    
105           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
106    
107  #endif /* ALLOW_SALT0_CONTROL */  #endif /* ALLOW_SALT0_CONTROL */
108    
# Line 114  c--   == end of interface == Line 117  c--   == end of interface ==
117       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
118       &                        mythid, dummy)       &                        mythid, dummy)
119    
120           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
121    
122  #endif /* ALLOW_HFLUX_CONTROL */  #endif /* ALLOW_HFLUX_CONTROL */
123    
# Line 129  c--   == end of interface == Line 132  c--   == end of interface ==
132       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
133       &                        mythid, dummy)       &                        mythid, dummy)
134    
135           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
136    
137  #endif /* ALLOW_SFLUX_CONTROL */  #endif /* ALLOW_SFLUX_CONTROL */
138    
# Line 144  c--   == end of interface == Line 147  c--   == end of interface ==
147       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
148       &                        mythid, dummy)       &                        mythid, dummy)
149    
150           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
151    
152  #endif /* ALLOW_USTRESS_CONTROL */  #endif /* ALLOW_USTRESS_CONTROL */
153    
# Line 159  c--   == end of interface == Line 162  c--   == end of interface ==
162       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
163       &                        mythid, dummy)       &                        mythid, dummy)
164    
165           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
166    
167  #endif /* ALLOW_VSTRESS_CONTROL */  #endif /* ALLOW_VSTRESS_CONTROL */
168    
# Line 174  c--   == end of interface == Line 177  c--   == end of interface ==
177       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
178       &                         mythid, dummy)       &                         mythid, dummy)
179    
180           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
181    
182  #endif /* ALLOW_TR10_CONTROL */  #endif /* ALLOW_TR10_CONTROL */
183    
# Line 189  c--   == end of interface == Line 192  c--   == end of interface ==
192       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
193       &                        mythid, dummy)       &                        mythid, dummy)
194    
195           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
196    
197  #endif /* ALLOW_SST0_CONTROL */  #endif /* ALLOW_SST0_CONTROL */
198    
# Line 204  c--   == end of interface == Line 207  c--   == end of interface ==
207       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
208       &                        mythid, dummy)       &                        mythid, dummy)
209    
210           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
211    
212  #endif /* ALLOW_SSS0_CONTROL */  #endif /* ALLOW_SSS0_CONTROL */
213    
# Line 221  c--   == end of interface == Line 224  c--   == end of interface ==
224       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
225       &                        mythid, dummy)       &                        mythid, dummy)
226    
227           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
228    
229  #else  #else
230    
# Line 229  c--   == end of interface == Line 232  c--   == end of interface ==
232       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
233       &                        mythid, dummy)       &                        mythid, dummy)
234    
235           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           adxx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
236    
237  #endif /* ALLOW_HFACC3D_CONTROL */  #endif /* ALLOW_HFACC3D_CONTROL */
238  #endif /* ALLOW_HFACC_CONTROL */  #endif /* ALLOW_HFACC_CONTROL */
239    
240    #ifdef ALLOW_EFLUXY0_CONTROL
241          else if ( grdchkvarindex .eq. 21 ) then
242             il=ilnblnk( xx_efluxy_file )
243             write(fname(1:80),'(80a)') ' '
244             write(fname(1:80),'(3a,i10.10)')
245         &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
246    
247             call active_read_xyz( fname, tmpfld3d, 1,
248         &                         doglobalread, ladinit, optimcycle,
249         &                         mythid, dummy)
250    
251             adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
252    
253    #endif /* ALLOW_THETA0_CONTROL */
254    
255    #ifdef ALLOW_EFLUXY0_CONTROL
256          else if ( grdchkvarindex .eq. 21 ) then
257             il=ilnblnk( xx_efluxy_file )
258             write(fname(1:80),'(80a)') ' '
259             write(fname(1:80),'(3a,i10.10)')
260         &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
261    
262             call active_read_xyz( fname, tmpfld3d, 1,
263         &                         doglobalread, ladinit, optimcycle,
264         &                         mythid, dummy)
265    
266             adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
267    
268    #endif /* ALLOW_EFLUXY0_CONTROL */
269    
270    #ifdef ALLOW_EFLUXP0_CONTROL
271          else if ( grdchkvarindex .eq. 22 ) then
272             il=ilnblnk( xx_efluxp_file )
273             write(fname(1:80),'(80a)') ' '
274             write(fname(1:80),'(3a,i10.10)')
275         &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
276    
277             call active_read_xyz( fname, tmpfld3d, 1,
278         &                         doglobalread, ladinit, optimcycle,
279         &                         mythid, dummy)
280    
281             adxx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
282    
283    #endif /* ALLOW_EFLUXP0_CONTROL */
284    
285        else        else
286  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
287        endif        endif
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22