/[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.6 by heimbach, Thu Dec 19 14:00:08 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)')
84       &        yadmark, xx_theta_file(1:il),'.',optimcycle       &        yadprefix, 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           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 93  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       &        yadmark, xx_salt_file(1:il),'.',optimcycle       &        yadprefix, 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           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 108  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       &        yadmark, xx_hflux_file(1:il),'.',optimcycle       &        yadprefix, 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           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 123  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       &        yadmark, xx_sflux_file(1:il),'.',optimcycle       &        yadprefix, 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           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 138  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       &        yadmark, xx_tauu_file(1:il),'.',optimcycle       &        yadprefix, 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           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 153  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       &        yadmark, xx_tauv_file(1:il),'.',optimcycle       &        yadprefix, 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           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 168  c--   == end of interface == Line 171  c--   == end of interface ==
171           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_tr1_file )
172           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
173           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
174       &        yadmark, xx_tr1_file(1:il),'.',optimcycle       &        yadprefix, xx_tr1_file(1:il),'.',optimcycle
175                    
176           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
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 183  c--   == end of interface == Line 186  c--   == end of interface ==
186           il=ilnblnk( xx_sst_file )           il=ilnblnk( xx_sst_file )
187           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
188           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
189       &        yadmark, xx_sst_file(1:il),'.',optimcycle       &        yadprefix, xx_sst_file(1:il),'.',optimcycle
190    
191           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
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 198  c--   == end of interface == Line 201  c--   == end of interface ==
201           il=ilnblnk( xx_sss_file )           il=ilnblnk( xx_sss_file )
202           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
203           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
204       &        yadmark, xx_sss_file(1:il),'.',optimcycle       &        yadprefix, xx_sss_file(1:il),'.',optimcycle
205    
206           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
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 213  c--   == end of interface == Line 216  c--   == end of interface ==
216           il=ilnblnk( xx_hfacc_file )           il=ilnblnk( xx_hfacc_file )
217           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
218           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
219       &        yadmark, xx_hfacc_file(1:il),'.',optimcycle       &        yadprefix, xx_hfacc_file(1:il),'.',optimcycle
220    
221  #ifdef ALLOW_HFACC3D_CONTROL  #ifdef ALLOW_HFACC3D_CONTROL
222    
# 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         &        yadprefix, 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         &        yadprefix, 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         &        yadprefix, 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.6
  ViewVC Help
Powered by ViewVC 1.1.22