/[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.21 by dfer, Thu Jan 17 20:49:27 2008 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CTRL_CPPOPTIONS.h"  #include "CTRL_CPPOPTIONS.h"
5    
# Line 48  c     == routine arguments == Line 49  c     == routine arguments ==
49        _RL     xx_comp        _RL     xx_comp
50        integer mythid        integer mythid
51    
52  #ifdef ALLOW_GRADIENT_CHECK  #ifdef ALLOW_GRDCHK
53  c     == local variables ==  c     == local variables ==
54    
55        integer il        integer il
# Line 61  c     == local variables == Line 62  c     == local variables ==
62    
63        character*(80) fname        character*(80) fname
64    
65          integer i,j,k
66    
67  c--   == external ==  c--   == external ==
68    
69        integer  ilnblnk        integer  ilnblnk
# Line 73  c--   == end of interface == Line 76  c--   == end of interface ==
76        dumiter      = 0        dumiter      = 0
77        dumtime      = 0. _d 0        dumtime      = 0. _d 0
78    
79        if      ( grdchkvarindex .eq. 1 ) then        if ( grdchkvarindex .eq. 0 ) then
80             STOP 'GRDCHK INDEX 0 NOT ALLOWED'
81    
82  #ifdef ALLOW_THETA0_CONTROL  #ifdef ALLOW_THETA0_CONTROL
83          else if ( grdchkvarindex .eq. 1 ) then
84           il=ilnblnk( xx_theta_file )           il=ilnblnk( xx_theta_file )
85           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
86           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
# Line 83  c--   == end of interface == Line 89  c--   == end of interface ==
89           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
90       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
91       &                         mythid, dummy)       &                         mythid, dummy)
   
92           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
   
93  #endif /* ALLOW_THETA0_CONTROL */  #endif /* ALLOW_THETA0_CONTROL */
94    
95  #ifdef ALLOW_SALT0_CONTROL  #ifdef ALLOW_SALT0_CONTROL
# Line 94  c--   == end of interface == Line 98  c--   == end of interface ==
98           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
99           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
100       &        yadmark, xx_salt_file(1:il),'.',optimcycle       &        yadmark, xx_salt_file(1:il),'.',optimcycle
           
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           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
105    
106  #endif /* ALLOW_SALT0_CONTROL */  #endif /* ALLOW_SALT0_CONTROL */
# Line 109  c--   == end of interface == Line 111  c--   == end of interface ==
111           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
112           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
113       &        yadmark, xx_hflux_file(1:il),'.',optimcycle       &        yadmark, xx_hflux_file(1:il),'.',optimcycle
   
114           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
115       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
116       &                        mythid, dummy)       &                        mythid, dummy)
   
117           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
118  #endif /* ALLOW_HFLUX_CONTROL */  #endif /* ALLOW_HFLUX_CONTROL */
119    
120  #ifdef ALLOW_SFLUX_CONTROL  #ifdef ALLOW_SFLUX_CONTROL
# Line 124  c--   == end of interface == Line 123  c--   == end of interface ==
123           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
124           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
125       &        yadmark, xx_sflux_file(1:il),'.',optimcycle       &        yadmark, xx_sflux_file(1:il),'.',optimcycle
   
126           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
127       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
128       &                        mythid, dummy)       &                        mythid, dummy)
   
129           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
130  #endif /* ALLOW_SFLUX_CONTROL */  #endif /* ALLOW_SFLUX_CONTROL */
131    
132  #ifdef ALLOW_USTRESS_CONTROL  #ifdef ALLOW_USTRESS_CONTROL
# Line 139  c--   == end of interface == Line 135  c--   == end of interface ==
135           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
136           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
137       &        yadmark, xx_tauu_file(1:il),'.',optimcycle       &        yadmark, xx_tauu_file(1:il),'.',optimcycle
   
138           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
139       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
140       &                        mythid, dummy)       &                        mythid, dummy)
   
141           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
142  #endif /* ALLOW_USTRESS_CONTROL */  #endif /* ALLOW_USTRESS_CONTROL */
143    
144  #ifdef ALLOW_VSTRESS_CONTROL  #ifdef ALLOW_VSTRESS_CONTROL
# Line 154  c--   == end of interface == Line 147  c--   == end of interface ==
147           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
148           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
149       &        yadmark, xx_tauv_file(1:il),'.',optimcycle       &        yadmark, xx_tauv_file(1:il),'.',optimcycle
150             call active_read_xy( fname, tmpfld2d, icvrec,
151         &                        doglobalread, ladinit, optimcycle,
152         &                        mythid, dummy)
153             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
154    #endif /* ALLOW_VSTRESS_CONTROL */
155    
156    #ifdef ALLOW_ATEMP_CONTROL
157          else if ( grdchkvarindex .eq. 7 ) then
158             il=ilnblnk( xx_atemp_file )
159             write(fname(1:80),'(80a)') ' '
160             write(fname(1:80),'(3a,i10.10)')
161         &        yadmark, xx_atemp_file(1:il),'.',optimcycle
162           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
163       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
164       &                        mythid, dummy)       &                        mythid, dummy)
165             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
166    #endif /* ALLOW_ATEMP_CONTROL */
167    
168    #ifdef ALLOW_AQH_CONTROL
169          else if ( grdchkvarindex .eq. 8 ) then
170             il=ilnblnk( xx_aqh_file )
171             write(fname(1:80),'(80a)') ' '
172             write(fname(1:80),'(3a,i10.10)')
173         &        yadmark, xx_aqh_file(1:il),'.',optimcycle
174             call active_read_xy( fname, tmpfld2d, icvrec,
175         &                        doglobalread, ladinit, optimcycle,
176         &                        mythid, dummy)
177           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
178    #endif /* ALLOW_AQH_CONTROL */
179    
180  #endif /* ALLOW_VSTRESS_CONTROL */  #ifdef ALLOW_UWIND_CONTROL
181          else if ( grdchkvarindex .eq. 9 ) then
182             il=ilnblnk( xx_uwind_file )
183             write(fname(1:80),'(80a)') ' '
184             write(fname(1:80),'(3a,i10.10)')
185         &        yadmark, xx_uwind_file(1:il),'.',optimcycle
186             call active_read_xy( fname, tmpfld2d, icvrec,
187         &                        doglobalread, ladinit, optimcycle,
188         &                        mythid, dummy)
189             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
190    #endif /* ALLOW_UWIND_CONTROL */
191    
192    #ifdef ALLOW_VWIND_CONTROL
193          else if ( grdchkvarindex .eq. 10 ) then
194             il=ilnblnk( xx_vwind_file )
195             write(fname(1:80),'(80a)') ' '
196             write(fname(1:80),'(3a,i10.10)')
197         &        yadmark, xx_vwind_file(1:il),'.',optimcycle
198             call active_read_xy( fname, tmpfld2d, icvrec,
199         &                        doglobalread, ladinit, optimcycle,
200         &                        mythid, dummy)
201             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
202    #endif /* ALLOW_VWIND_CONTROL */
203    
204    #ifdef ALLOW_OBCSN_CONTROL
205          else if ( grdchkvarindex .eq. 11 ) then
206             il=ilnblnk( xx_obcsn_file )
207             write(fname(1:80),'(80a)') ' '
208             write(fname(1:80),'(3a,i10.10)')
209         &        yadmark, xx_obcsn_file(1:il),'.',optimcycle
210    
211             call active_read_xz( fname, tmpfldxz, icvrec,
212         &                        doglobalread, ladinit, optimcycle,
213         &                        mythid, dummy)
214    
215             xx_comp = tmpfldxz( itilepos,layer,itile,jtile )
216    
217    #endif /* ALLOW_OBCSN_CONTROL */
218    
219    #ifdef ALLOW_OBCSS_CONTROL
220          else if ( grdchkvarindex .eq. 12 ) then
221             il=ilnblnk( xx_obcss_file )
222             write(fname(1:80),'(80a)') ' '
223             write(fname(1:80),'(3a,i10.10)')
224         &        yadmark, xx_obcss_file(1:il),'.',optimcycle
225    
226             call active_read_xz( fname, tmpfldxz, icvrec,
227         &                        doglobalread, ladinit, optimcycle,
228         &                        mythid, dummy)
229    
230             xx_comp = tmpfldxz( itilepos,layer,itile,jtile )
231    
232    #endif /* ALLOW_OBCSS_CONTROL */
233    
234    #ifdef ALLOW_OBCSW_CONTROL
235          else if ( grdchkvarindex .eq. 13 ) then
236             il=ilnblnk( xx_obcsw_file )
237             write(fname(1:80),'(80a)') ' '
238             write(fname(1:80),'(3a,i10.10)')
239         &        yadmark, xx_obcsw_file(1:il),'.',optimcycle
240    
241             call active_read_yz( fname, tmpfldyz, icvrec,
242         &                        doglobalread, ladinit, optimcycle,
243         &                        mythid, dummy)
244    
245             xx_comp = tmpfldyz( jtilepos,layer,itile,jtile )
246    
247    #endif /* ALLOW_OBCSW_CONTROL */
248    
249    #ifdef ALLOW_OBCSE_CONTROL
250          else if ( grdchkvarindex .eq. 14 ) then
251             il=ilnblnk( xx_obcse_file )
252             write(fname(1:80),'(80a)') ' '
253             write(fname(1:80),'(3a,i10.10)')
254         &        yadmark, xx_obcse_file(1:il),'.',optimcycle
255    
256             call active_read_yz( fname, tmpfldyz, icvrec,
257         &                        doglobalread, ladinit, optimcycle,
258         &                        mythid, dummy)
259    
260             xx_comp = tmpfldyz( jtilepos,layer,itile,jtile )
261    
262    #endif /* ALLOW_OBCSE_CONTROL */
263    
264    #ifdef ALLOW_DIFFKR_CONTROL
265          else if ( grdchkvarindex .eq. 15 ) then
266             il=ilnblnk( xx_diffkr_file )
267             write(fname(1:80),'(80a)') ' '
268             write(fname(1:80),'(3a,i10.10)')
269         &        yadmark, xx_diffkr_file(1:il),'.',optimcycle
270    
271             call active_read_xyz( fname, tmpfld3d, 1,
272         &                         doglobalread, ladinit, optimcycle,
273         &                         mythid, dummy)
274    
275             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
276    
277    #endif /* ALLOW_DIFFKR_CONTROL */
278    
279    #ifdef ALLOW_KAPGM_CONTROL
280          else if ( grdchkvarindex .eq. 16 ) then
281             il=ilnblnk( xx_kapgm_file )
282             write(fname(1:80),'(80a)') ' '
283             write(fname(1:80),'(3a,i10.10)')
284         &        yadmark, xx_kapgm_file(1:il),'.',optimcycle
285    
286             call active_read_xyz( fname, tmpfld3d, 1,
287         &                         doglobalread, ladinit, optimcycle,
288         &                         mythid, dummy)
289    
290             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
291    
292    #endif /* ALLOW_KAPGM_CONTROL */
293    
294  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
295        else if ( grdchkvarindex .eq. 17 ) then        else if ( grdchkvarindex .eq. 17 ) then
# Line 169  c--   == end of interface == Line 297  c--   == end of interface ==
297           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
298           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
299       &        yadmark, xx_tr1_file(1:il),'.',optimcycle       &        yadmark, xx_tr1_file(1:il),'.',optimcycle
300            
301           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
302       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
303       &                         mythid, dummy)       &                         mythid, dummy)
# Line 178  c--   == end of interface == Line 306  c--   == end of interface ==
306    
307  #endif /* ALLOW_TR10_CONTROL */  #endif /* ALLOW_TR10_CONTROL */
308    
309  #ifdef ALLOW_SST0_CONTROL  #if (defined (ALLOW_SST_CONTROL) || defined (ALLOW_SST0_CONTROL))
310        else if ( grdchkvarindex .eq. 18 ) then        else if ( grdchkvarindex .eq. 18 ) then
311           il=ilnblnk( xx_sst_file )           il=ilnblnk( xx_sst_file )
312           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
313           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
314       &        yadmark, xx_sst_file(1:il),'.',optimcycle       &        yadmark, xx_sst_file(1:il),'.',optimcycle
   
315           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
316       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
317       &                        mythid, dummy)       &                        mythid, dummy)
   
318           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
319  #endif /* ALLOW_SST0_CONTROL */  #endif /* ALLOW_SST0_CONTROL */
320    
321  #ifdef ALLOW_SSS0_CONTROL  #if (defined (ALLOW_SSS_CONTROL) || defined (ALLOW_SSS0_CONTROL))
322        else if ( grdchkvarindex .eq. 19 ) then        else if ( grdchkvarindex .eq. 19 ) then
323           il=ilnblnk( xx_sss_file )           il=ilnblnk( xx_sss_file )
324           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
325           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
326       &        yadmark, xx_sss_file(1:il),'.',optimcycle       &        yadmark, xx_sss_file(1:il),'.',optimcycle
327             call active_read_xy( fname, tmpfld2d, icvrec,
328         &                        doglobalread, ladinit, optimcycle,
329         &                        mythid, dummy)
330             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
331    #endif /* ALLOW_SSS0_CONTROL */
332    
333    #ifdef ALLOW_DEPTH_CONTROL
334          else if ( grdchkvarindex .eq. 20 ) then
335             il=ilnblnk( xx_depth_file )
336             write(fname(1:80),'(80a)') ' '
337             write(fname(1:80),'(3a,i10.10)')
338         &        yadmark, xx_depth_file(1:il),'.',optimcycle
339           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
340       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
341       &                        mythid, dummy)       &                        mythid, dummy)
342             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
343    #endif /* ALLOW_DEPTH_CONTROL */
344    
345    #ifdef ALLOW_EFLUXY0_CONTROL
346          else if ( grdchkvarindex .eq. 21 ) then
347             il=ilnblnk( xx_efluxy_file )
348             write(fname(1:80),'(80a)') ' '
349             write(fname(1:80),'(3a,i10.10)')
350         &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
351    
352             call active_read_xyz( fname, tmpfld3d, 1,
353         &                         doglobalread, ladinit, optimcycle,
354         &                         mythid, dummy)
355    
356             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
357    
358    #endif /* ALLOW_EFLUXY0_CONTROL */
359    
360    #ifdef ALLOW_EFLUXP0_CONTROL
361          else if ( grdchkvarindex .eq. 22 ) then
362             il=ilnblnk( xx_efluxp_file )
363             write(fname(1:80),'(80a)') ' '
364             write(fname(1:80),'(3a,i10.10)')
365         &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
366    
367             call active_read_xyz( fname, tmpfld3d, 1,
368         &                         doglobalread, ladinit, optimcycle,
369         &                         mythid, dummy)
370    
371             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
372    
373    #endif /* ALLOW_EFLUXP0_CONTROL */
374    
375    #ifdef ALLOW_HFLUXM_CONTROL
376          else if ( grdchkvarindex .eq. 24 ) then
377             il=ilnblnk( xx_hfluxm_file )
378             write(fname(1:80),'(80a)') ' '
379             write(fname(1:80),'(3a,i10.10)')
380         &        yadmark, xx_hfluxm_file(1:il),'.',optimcycle
381             call active_read_xy( fname, tmpfld2d, icvrec,
382         &                        doglobalread, ladinit, optimcycle,
383         &                        mythid, dummy)
384           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
385    #endif /* ALLOW_HFLUXM_CONTROL */
386    
387  #endif /* ALLOW_SSS0_CONTROL */  #ifdef ALLOW_PRECIP_CONTROL
388          else if ( grdchkvarindex .eq. 32 ) then
389             il=ilnblnk( xx_precip_file )
390             write(fname(1:80),'(80a)') ' '
391             write(fname(1:80),'(3a,i10.10)')
392         &        yadmark, xx_precip_file(1:il),'.',optimcycle
393             call active_read_xy( fname, tmpfld2d, icvrec,
394         &                        doglobalread, ladinit, optimcycle,
395         &                        mythid, dummy)
396             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
397    #endif /* ALLOW_PRECIP_CONTROL */
398    
399  #ifdef ALLOW_HFACC_CONTROL  #ifdef ALLOW_SWFLUX_CONTROL
400        else if ( grdchkvarindex .eq. 20 ) then        else if ( grdchkvarindex .eq. 33 ) then
401           il=ilnblnk( xx_hfacc_file )           il=ilnblnk( xx_swflux_file )
402           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
403           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
404       &        yadmark, xx_hfacc_file(1:il),'.',optimcycle       &        yadmark, xx_swflux_file(1:il),'.',optimcycle
405             call active_read_xy( fname, tmpfld2d, icvrec,
406         &                        doglobalread, ladinit, optimcycle,
407         &                        mythid, dummy)
408             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
409    #endif /* ALLOW_SWFLUX_CONTROL */
410    
411  #ifdef ALLOW_HFACC3D_CONTROL  #ifdef ALLOW_SWDOWN_CONTROL
412          else if ( grdchkvarindex .eq. 34 ) then
413             il=ilnblnk( xx_swdown_file )
414             write(fname(1:80),'(80a)') ' '
415             write(fname(1:80),'(3a,i10.10)')
416         &        yadmark, xx_swdown_file(1:il),'.',optimcycle
417             call active_read_xy( fname, tmpfld2d, icvrec,
418         &                        doglobalread, ladinit, optimcycle,
419         &                        mythid, dummy)
420             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
421    #endif /* ALLOW_SWDOWN_CONTROL */
422    
423           call active_read_xyz( fname, tmpfld3d, icvrec,  #ifdef ALLOW_LWFLUX_CONTROL
424          else if ( grdchkvarindex .eq. 35 ) then
425             il=ilnblnk( xx_lwflux_file )
426             write(fname(1:80),'(80a)') ' '
427             write(fname(1:80),'(3a,i10.10)')
428         &        yadmark, xx_lwflux_file(1:il),'.',optimcycle
429             call active_read_xy( fname, tmpfld2d, icvrec,
430       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
431       &                        mythid, dummy)       &                        mythid, dummy)
432             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
433    #endif /* ALLOW_LWFLUX_CONTROL */
434    
435           xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )  #ifdef ALLOW_LWDOWN_CONTROL
436          else if ( grdchkvarindex .eq. 36 ) then
437             il=ilnblnk( xx_lwdown_file )
438             write(fname(1:80),'(80a)') ' '
439             write(fname(1:80),'(3a,i10.10)')
440         &        yadmark, xx_lwdown_file(1:il),'.',optimcycle
441             call active_read_xy( fname, tmpfld2d, icvrec,
442         &                        doglobalread, ladinit, optimcycle,
443         &                        mythid, dummy)
444             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
445    #endif /* ALLOW_LWDOWN_CONTROL */
446    
447    #ifdef ALLOW_EVAP_CONTROL
448          else if ( grdchkvarindex .eq. 37 ) then
449             il=ilnblnk( xx_evap_file )
450             write(fname(1:80),'(80a)') ' '
451             write(fname(1:80),'(3a,i10.10)')
452         &        yadmark, xx_evap_file(1:il),'.',optimcycle
453             call active_read_xy( fname, tmpfld2d, icvrec,
454         &                        doglobalread, ladinit, optimcycle,
455         &                        mythid, dummy)
456             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
457    #endif /* ALLOW_EVAP_CONTROL */
458    
459    #ifdef ALLOW_SNOWPRECIP_CONTROL
460          else if ( grdchkvarindex .eq. 38 ) then
461             il=ilnblnk( xx_snowprecip_file )
462             write(fname(1:80),'(80a)') ' '
463             write(fname(1:80),'(3a,i10.10)')
464         &        yadmark, xx_snowprecip_file(1:il),'.',optimcycle
465             call active_read_xy( fname, tmpfld2d, icvrec,
466         &                        doglobalread, ladinit, optimcycle,
467         &                        mythid, dummy)
468             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
469    #endif /* ALLOW_SNOWPRECIP_CONTROL */
470    
471    #ifdef ALLOW_APRESSURE_CONTROL
472          else if ( grdchkvarindex .eq. 39 ) then
473             il=ilnblnk( xx_apressure_file )
474             write(fname(1:80),'(80a)') ' '
475             write(fname(1:80),'(3a,i10.10)')
476         &        yadmark, xx_apressure_file(1:il),'.',optimcycle
477    
478  #else           call active_read_xy( fname, tmpfld2d, icvrec,
479         &                        doglobalread, ladinit, optimcycle,
480         &                        mythid, dummy)
481             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
482    #endif /* ALLOW_APRESSURE_CONTROL */
483    
484    #ifdef ALLOW_RUNOFF_CONTROL
485          else if ( grdchkvarindex .eq. 40 ) then
486             il=ilnblnk( xx_runoff_file )
487             write(fname(1:80),'(80a)') ' '
488             write(fname(1:80),'(3a,i10.10)')
489         &        yadmark, xx_runoff_file(1:il),'.',optimcycle
490           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy( fname, tmpfld2d, icvrec,
491       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
492       &                        mythid, dummy)       &                        mythid, dummy)
493             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
494    #endif /* ALLOW_RUNOFF_CONTROL */
495    
496    #ifdef ALLOW_SIAREA_CONTROL
497          else if ( grdchkvarindex .eq. 41 ) then
498             il=ilnblnk( xx_siarea_file )
499             write(fname(1:80),'(80a)') ' '
500             write(fname(1:80),'(3a,i10.10)')
501         &        yadmark, xx_siarea_file(1:il),'.',optimcycle
502             call active_read_xy( fname, tmpfld2d, icvrec,
503         &                        doglobalread, ladinit, optimcycle,
504         &                        mythid, dummy)
505           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
506    #endif /* ALLOW_SIAREA_CONTROL */
507    
508  #endif /* ALLOW_HFACC3D_CONTROL */  #ifdef ALLOW_SIHEFF_CONTROL
509  #endif /* ALLOW_HFACC_CONTROL */        else if ( grdchkvarindex .eq. 42 ) then
510             il=ilnblnk( xx_siheff_file )
511             write(fname(1:80),'(80a)') ' '
512             write(fname(1:80),'(3a,i10.10)')
513         &        yadmark, xx_siheff_file(1:il),'.',optimcycle
514             call active_read_xy( fname, tmpfld2d, icvrec,
515         &                        doglobalread, ladinit, optimcycle,
516         &                        mythid, dummy)
517             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
518    #endif /* ALLOW_SIHEFF_CONTROL */
519    
520    #ifdef ALLOW_SIHSNOW_CONTROL
521          else if ( grdchkvarindex .eq. 43 ) then
522             il=ilnblnk( xx_sihsnow_file )
523             write(fname(1:80),'(80a)') ' '
524             write(fname(1:80),'(3a,i10.10)')
525         &        yadmark, xx_sihsnow_file(1:il),'.',optimcycle
526             call active_read_xy( fname, tmpfld2d, icvrec,
527         &                        doglobalread, ladinit, optimcycle,
528         &                        mythid, dummy)
529             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
530    #endif /* ALLOW_SIHSNOW_CONTROL */
531    
532        else        else
533  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
534        endif        endif
535    
536  #endif /* ALLOW_GRADIENT_CHECK */  #endif /* ALLOW_GRDCHK */
537    
538        end        end
539    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22