/[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.9 by heimbach, Fri Jul 18 21:10:16 2003 UTC revision 1.19 by heimbach, Thu Jun 21 04:06:21 2007 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 85  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_loc( fname, tmpfld3d, 1,           call active_read_xyz( 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( 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( 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( 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( 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( 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( 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( 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( 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( 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 235  c--   == end of interface == Line 207  c--   == end of interface ==
207           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
208       &        yadmark, xx_obcsn_file(1:il),'.',optimcycle       &        yadmark, xx_obcsn_file(1:il),'.',optimcycle
209    
210           call active_read_xz_loc( fname, tmpfldxz, icvrec,           call active_read_xz( fname, tmpfldxz, icvrec,
211       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
212       &                        mythid, dummy)       &                        mythid, dummy)
213    
# Line 250  c--   == end of interface == Line 222  c--   == end of interface ==
222           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
223       &        yadmark, xx_obcss_file(1:il),'.',optimcycle       &        yadmark, xx_obcss_file(1:il),'.',optimcycle
224    
225           call active_read_xz_loc( fname, tmpfldxz, icvrec,           call active_read_xz( fname, tmpfldxz, icvrec,
226       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
227       &                        mythid, dummy)       &                        mythid, dummy)
228    
# Line 265  c--   == end of interface == Line 237  c--   == end of interface ==
237           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
238       &        yadmark, xx_obcsw_file(1:il),'.',optimcycle       &        yadmark, xx_obcsw_file(1:il),'.',optimcycle
239    
240           call active_read_yz_loc( fname, tmpfldyz, icvrec,           call active_read_yz( fname, tmpfldyz, icvrec,
241       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
242       &                        mythid, dummy)       &                        mythid, dummy)
243    
# Line 280  c--   == end of interface == Line 252  c--   == end of interface ==
252           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
253       &        yadmark, xx_obcse_file(1:il),'.',optimcycle       &        yadmark, xx_obcse_file(1:il),'.',optimcycle
254    
255           call active_read_yz_loc( fname, tmpfldyz, icvrec,           call active_read_yz( fname, tmpfldyz, icvrec,
256       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
257       &                        mythid, dummy)       &                        mythid, dummy)
258    
# Line 288  c--   == end of interface == Line 260  c--   == end of interface ==
260    
261  #endif /* ALLOW_OBCSE_CONTROL */  #endif /* ALLOW_OBCSE_CONTROL */
262    
263    #ifdef ALLOW_DIFFKR_CONTROL
264          else if ( grdchkvarindex .eq. 15 ) then
265             il=ilnblnk( xx_diffkr_file )
266             write(fname(1:80),'(80a)') ' '
267             write(fname(1:80),'(3a,i10.10)')
268         &        yadmark, xx_diffkr_file(1:il),'.',optimcycle
269    
270             call active_read_xyz( fname, tmpfld3d, 1,
271         &                         doglobalread, ladinit, optimcycle,
272         &                         mythid, dummy)
273    
274             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
275    
276    #endif /* ALLOW_DIFFKR_CONTROL */
277    
278    #ifdef ALLOW_KAPGM_CONTROL
279          else if ( grdchkvarindex .eq. 16 ) then
280             il=ilnblnk( xx_kapgm_file )
281             write(fname(1:80),'(80a)') ' '
282             write(fname(1:80),'(3a,i10.10)')
283         &        yadmark, xx_kapgm_file(1:il),'.',optimcycle
284    
285             call active_read_xyz( fname, tmpfld3d, 1,
286         &                         doglobalread, ladinit, optimcycle,
287         &                         mythid, dummy)
288    
289             xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
290    
291    #endif /* ALLOW_KAPGM_CONTROL */
292    
293  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
294        else if ( grdchkvarindex .eq. 17 ) then        else if ( grdchkvarindex .eq. 17 ) then
295           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_tr1_file )
# Line 295  c--   == end of interface == Line 297  c--   == end of interface ==
297           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
298       &        yadmark, xx_tr1_file(1:il),'.',optimcycle       &        yadmark, xx_tr1_file(1:il),'.',optimcycle
299                    
300           call active_read_xyz_loc( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
301       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
302       &                         mythid, dummy)       &                         mythid, dummy)
303    
# Line 303  c--   == end of interface == Line 305  c--   == end of interface ==
305    
306  #endif /* ALLOW_TR10_CONTROL */  #endif /* ALLOW_TR10_CONTROL */
307    
308  #ifdef ALLOW_SST0_CONTROL  #if (defined (ALLOW_SST_CONTROL) || defined (ALLOW_SST0_CONTROL))
309        else if ( grdchkvarindex .eq. 18 ) then        else if ( grdchkvarindex .eq. 18 ) then
310           il=ilnblnk( xx_sst_file )           il=ilnblnk( xx_sst_file )
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( 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  #ifdef ALLOW_SSS0_CONTROL  #if (defined (ALLOW_SSS_CONTROL) || defined (ALLOW_SSS0_CONTROL))
321        else if ( grdchkvarindex .eq. 19 ) then        else if ( grdchkvarindex .eq. 19 ) then
322           il=ilnblnk( xx_sss_file )           il=ilnblnk( xx_sss_file )
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( 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_HFACC_CONTROL  #ifdef ALLOW_DEPTH_CONTROL
333        else if ( grdchkvarindex .eq. 20 ) then        else if ( grdchkvarindex .eq. 20 ) then
334           il=ilnblnk( xx_hfacc_file )           il=ilnblnk( xx_depth_file )
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_hfacc_file(1:il),'.',optimcycle       &        yadmark, xx_depth_file(1:il),'.',optimcycle
338             call active_read_xy( fname, tmpfld2d, icvrec,
 #ifdef ALLOW_HFACC3D_CONTROL  
   
          call active_read_xyz_loc( fname, tmpfld3d, icvrec,  
      &                        doglobalread, ladinit, optimcycle,  
      &                        mythid, dummy)  
   
          xx_comp = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )  
   
 #else  
   
          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_HFACC3D_CONTROL */  
 #endif /* ALLOW_HFACC_CONTROL */  
343    
344  #ifdef ALLOW_EFLUXY0_CONTROL  #ifdef ALLOW_EFLUXY0_CONTROL
345        else if ( grdchkvarindex .eq. 21 ) then        else if ( grdchkvarindex .eq. 21 ) then
# Line 366  c--   == end of interface == Line 348  c--   == end of interface ==
348           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
349       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
350    
351           call active_read_xyz_loc( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
352       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
353       &                         mythid, dummy)       &                         mythid, dummy)
354    
# Line 381  c--   == end of interface == Line 363  c--   == end of interface ==
363           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
364       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
365    
366           call active_read_xyz_loc( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
367       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
368       &                         mythid, dummy)       &                         mythid, dummy)
369    
# Line 389  c--   == end of interface == Line 371  c--   == end of interface ==
371    
372  #endif /* ALLOW_EFLUXP0_CONTROL */  #endif /* ALLOW_EFLUXP0_CONTROL */
373    
374    #ifdef ALLOW_PRECIP_CONTROL
375          else if ( grdchkvarindex .eq. 32 ) then
376             il=ilnblnk( xx_precip_file )
377             write(fname(1:80),'(80a)') ' '
378             write(fname(1:80),'(3a,i10.10)')
379         &        yadmark, xx_precip_file(1:il),'.',optimcycle
380             call active_read_xy( fname, tmpfld2d, icvrec,
381         &                        doglobalread, ladinit, optimcycle,
382         &                        mythid, dummy)
383             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
384    #endif /* ALLOW_PRECIP_CONTROL */
385    
386    #ifdef ALLOW_SWFLUX_CONTROL
387          else if ( grdchkvarindex .eq. 33 ) then
388             il=ilnblnk( xx_swflux_file )
389             write(fname(1:80),'(80a)') ' '
390             write(fname(1:80),'(3a,i10.10)')
391         &        yadmark, xx_swflux_file(1:il),'.',optimcycle
392             call active_read_xy( fname, tmpfld2d, icvrec,
393         &                        doglobalread, ladinit, optimcycle,
394         &                        mythid, dummy)
395             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
396    #endif /* ALLOW_SWFLUX_CONTROL */
397    
398    #ifdef ALLOW_SWDOWN_CONTROL
399          else if ( grdchkvarindex .eq. 34 ) then
400             il=ilnblnk( xx_swdown_file )
401             write(fname(1:80),'(80a)') ' '
402             write(fname(1:80),'(3a,i10.10)')
403         &        yadmark, xx_swdown_file(1:il),'.',optimcycle
404             call active_read_xy( 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( fname, tmpfld2d, icvrec,
417         &                        doglobalread, ladinit, optimcycle,
418         &                        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( fname, tmpfld2d, icvrec,
429         &                        doglobalread, ladinit, optimcycle,
430         &                        mythid, dummy)
431             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
432    #endif /* ALLOW_LWDOWN_CONTROL */
433    
434    #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( 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( 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( 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( 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    #ifdef ALLOW_SIAREA_CONTROL
484          else if ( grdchkvarindex .eq. 41 ) then
485             il=ilnblnk( xx_siarea_file )
486             write(fname(1:80),'(80a)') ' '
487             write(fname(1:80),'(3a,i10.10)')
488         &        yadmark, xx_siarea_file(1:il),'.',optimcycle
489             call active_read_xy( fname, tmpfld2d, icvrec,
490         &                        doglobalread, ladinit, optimcycle,
491         &                        mythid, dummy)
492             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
493    #endif /* ALLOW_SIAREA_CONTROL */
494    
495    #ifdef ALLOW_SIHEFF_CONTROL
496          else if ( grdchkvarindex .eq. 42 ) then
497             il=ilnblnk( xx_siheff_file )
498             write(fname(1:80),'(80a)') ' '
499             write(fname(1:80),'(3a,i10.10)')
500         &        yadmark, xx_siheff_file(1:il),'.',optimcycle
501             call active_read_xy( fname, tmpfld2d, icvrec,
502         &                        doglobalread, ladinit, optimcycle,
503         &                        mythid, dummy)
504             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
505    #endif /* ALLOW_SIHEFF_CONTROL */
506    
507    #ifdef ALLOW_SIHSNOW_CONTROL
508          else if ( grdchkvarindex .eq. 43 ) then
509             il=ilnblnk( xx_sihsnow_file )
510             write(fname(1:80),'(80a)') ' '
511             write(fname(1:80),'(3a,i10.10)')
512         &        yadmark, xx_sihsnow_file(1:il),'.',optimcycle
513             call active_read_xy( fname, tmpfld2d, icvrec,
514         &                        doglobalread, ladinit, optimcycle,
515         &                        mythid, dummy)
516             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
517    #endif /* ALLOW_SIHSNOW_CONTROL */
518    
519        else        else
520  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
521        endif        endif
522    
523  #endif /* ALLOW_GRADIENT_CHECK */  #endif /* ALLOW_GRDCHK */
524    
525        end        end
526    
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22