/[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.14 by heimbach, Wed Aug 31 00:03:44 2005 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 85  c--   == end of interface == Line 86  c--   == end of interface ==
86           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
87       &        yadmark, xx_theta_file(1:il),'.',optimcycle       &        yadmark, xx_theta_file(1:il),'.',optimcycle
88    
89           call active_read_xyz_loc( 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 99  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_loc( 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 114  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_loc( 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 129  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_loc( 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 144  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_loc( 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 159  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,
          call active_read_xy_loc( fname, tmpfld2d, icvrec,  
151       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
152       &                        mythid, dummy)       &                        mythid, dummy)
   
153           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
154  #endif /* ALLOW_VSTRESS_CONTROL */  #endif /* ALLOW_VSTRESS_CONTROL */
155    
156  #ifdef ALLOW_ATEMP_CONTROL  #ifdef ALLOW_ATEMP_CONTROL
# Line 174  c--   == end of interface == Line 159  c--   == end of interface ==
159           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
160           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
161       &        yadmark, xx_atemp_file(1:il),'.',optimcycle       &        yadmark, xx_atemp_file(1:il),'.',optimcycle
162             call active_read_xy( fname, tmpfld2d, icvrec,
          call active_read_xy_loc( fname, tmpfld2d, icvrec,  
163       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
164       &                        mythid, dummy)       &                        mythid, dummy)
   
165           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
166  #endif /* ALLOW_ATEMP_CONTROL */  #endif /* ALLOW_ATEMP_CONTROL */
167    
168  #ifdef ALLOW_AQH_CONTROL  #ifdef ALLOW_AQH_CONTROL
# Line 189  c--   == end of interface == Line 171  c--   == end of interface ==
171           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
172           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
173       &        yadmark, xx_aqh_file(1:il),'.',optimcycle       &        yadmark, xx_aqh_file(1:il),'.',optimcycle
174             call active_read_xy( fname, tmpfld2d, icvrec,
          call active_read_xy_loc( fname, tmpfld2d, icvrec,  
175       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
176       &                        mythid, dummy)       &                        mythid, dummy)
   
177           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
178  #endif /* ALLOW_AQH_CONTROL */  #endif /* ALLOW_AQH_CONTROL */
179    
180  #ifdef ALLOW_UWIND_CONTROL  #ifdef ALLOW_UWIND_CONTROL
# Line 204  c--   == end of interface == Line 183  c--   == end of interface ==
183           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
184           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
185       &        yadmark, xx_uwind_file(1:il),'.',optimcycle       &        yadmark, xx_uwind_file(1:il),'.',optimcycle
186             call active_read_xy( fname, tmpfld2d, icvrec,
          call active_read_xy_loc( fname, tmpfld2d, icvrec,  
187       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
188       &                        mythid, dummy)       &                        mythid, dummy)
   
189           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
190  #endif /* ALLOW_UWIND_CONTROL */  #endif /* ALLOW_UWIND_CONTROL */
191    
192  #ifdef ALLOW_VWIND_CONTROL  #ifdef ALLOW_VWIND_CONTROL
# Line 219  c--   == end of interface == Line 195  c--   == end of interface ==
195           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
196           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
197       &        yadmark, xx_vwind_file(1:il),'.',optimcycle       &        yadmark, xx_vwind_file(1:il),'.',optimcycle
198             call active_read_xy( fname, tmpfld2d, icvrec,
          call active_read_xy_loc( fname, tmpfld2d, icvrec,  
199       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
200       &                        mythid, dummy)       &                        mythid, dummy)
   
201           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
202  #endif /* ALLOW_VWIND_CONTROL */  #endif /* ALLOW_VWIND_CONTROL */
203    
204  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
# Line 235  c--   == end of interface == Line 208  c--   == end of interface ==
208           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
209       &        yadmark, xx_obcsn_file(1:il),'.',optimcycle       &        yadmark, xx_obcsn_file(1:il),'.',optimcycle
210    
211           call active_read_xz_loc( fname, tmpfldxz, icvrec,           call active_read_xz( fname, tmpfldxz, icvrec,
212       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
213       &                        mythid, dummy)       &                        mythid, dummy)
214    
# Line 250  c--   == end of interface == Line 223  c--   == end of interface ==
223           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
224       &        yadmark, xx_obcss_file(1:il),'.',optimcycle       &        yadmark, xx_obcss_file(1:il),'.',optimcycle
225    
226           call active_read_xz_loc( fname, tmpfldxz, icvrec,           call active_read_xz( fname, tmpfldxz, icvrec,
227       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
228       &                        mythid, dummy)       &                        mythid, dummy)
229    
# Line 265  c--   == end of interface == Line 238  c--   == end of interface ==
238           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
239       &        yadmark, xx_obcsw_file(1:il),'.',optimcycle       &        yadmark, xx_obcsw_file(1:il),'.',optimcycle
240    
241           call active_read_yz_loc( fname, tmpfldyz, icvrec,           call active_read_yz( fname, tmpfldyz, icvrec,
242       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
243       &                        mythid, dummy)       &                        mythid, dummy)
244    
# Line 280  c--   == end of interface == Line 253  c--   == end of interface ==
253           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
254       &        yadmark, xx_obcse_file(1:il),'.',optimcycle       &        yadmark, xx_obcse_file(1:il),'.',optimcycle
255    
256           call active_read_yz_loc( fname, tmpfldyz, icvrec,           call active_read_yz( fname, tmpfldyz, icvrec,
257       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
258       &                        mythid, dummy)       &                        mythid, dummy)
259    
# Line 288  c--   == end of interface == Line 261  c--   == end of interface ==
261    
262  #endif /* ALLOW_OBCSE_CONTROL */  #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
296           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_tr1_file )
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_loc( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
302       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
303       &                         mythid, dummy)       &                         mythid, dummy)
304    
# Line 309  c--   == end of interface == Line 312  c--   == end of interface ==
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_loc( 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  #if (defined (ALLOW_SSS_CONTROL) || defined (ALLOW_SSS0_CONTROL))  #if (defined (ALLOW_SSS_CONTROL) || defined (ALLOW_SSS0_CONTROL))
# Line 324  c--   == end of interface == Line 324  c--   == end of interface ==
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,
          call active_read_xy_loc( fname, tmpfld2d, icvrec,  
328       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
329       &                        mythid, dummy)       &                        mythid, dummy)
   
330           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
331  #endif /* ALLOW_SSS0_CONTROL */  #endif /* ALLOW_SSS0_CONTROL */
332    
333  #ifdef ALLOW_HFACC_CONTROL  #ifdef ALLOW_DEPTH_CONTROL
334        else if ( grdchkvarindex .eq. 20 ) then        else if ( grdchkvarindex .eq. 20 ) then
335           il=ilnblnk( xx_hfacc_file )           il=ilnblnk( xx_depth_file )
336           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
337           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
338       &        yadmark, xx_hfacc_file(1:il),'.',optimcycle       &        yadmark, xx_depth_file(1:il),'.',optimcycle
339             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,  
340       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
341       &                        mythid, dummy)       &                        mythid, dummy)
   
342           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
343    #endif /* ALLOW_DEPTH_CONTROL */
 #endif /* ALLOW_HFACC3D_CONTROL */  
 #endif /* ALLOW_HFACC_CONTROL */  
344    
345  #ifdef ALLOW_EFLUXY0_CONTROL  #ifdef ALLOW_EFLUXY0_CONTROL
346        else if ( grdchkvarindex .eq. 21 ) then        else if ( grdchkvarindex .eq. 21 ) then
# Line 366  c--   == end of interface == Line 349  c--   == end of interface ==
349           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
350       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle       &        yadmark, xx_efluxy_file(1:il),'.',optimcycle
351    
352           call active_read_xyz_loc( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
353       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
354       &                         mythid, dummy)       &                         mythid, dummy)
355    
# Line 381  c--   == end of interface == Line 364  c--   == end of interface ==
364           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
365       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle       &        yadmark, xx_efluxp_file(1:il),'.',optimcycle
366    
367           call active_read_xyz_loc( fname, tmpfld3d, 1,           call active_read_xyz( fname, tmpfld3d, 1,
368       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
369       &                         mythid, dummy)       &                         mythid, dummy)
370    
# Line 389  c--   == end of interface == Line 372  c--   == end of interface ==
372    
373  #endif /* ALLOW_EFLUXP0_CONTROL */  #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 )
385    #endif /* ALLOW_HFLUXM_CONTROL */
386    
387  #ifdef ALLOW_PRECIP_CONTROL  #ifdef ALLOW_PRECIP_CONTROL
388        else if ( grdchkvarindex .eq. 32 ) then        else if ( grdchkvarindex .eq. 32 ) then
389           il=ilnblnk( xx_precip_file )           il=ilnblnk( xx_precip_file )
390           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
391           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
392       &        yadmark, xx_precip_file(1:il),'.',optimcycle       &        yadmark, xx_precip_file(1:il),'.',optimcycle
393             call active_read_xy( fname, tmpfld2d, icvrec,
          call active_read_xy_loc( fname, tmpfld2d, icvrec,  
394       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
395       &                        mythid, dummy)       &                        mythid, dummy)
   
396           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
397  #endif /* ALLOW_PRECIP_CONTROL */  #endif /* ALLOW_PRECIP_CONTROL */
398    
399  #ifdef ALLOW_SWFLUX_CONTROL  #ifdef ALLOW_SWFLUX_CONTROL
# Line 410  c--   == end of interface == Line 402  c--   == end of interface ==
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_swflux_file(1:il),'.',optimcycle       &        yadmark, xx_swflux_file(1:il),'.',optimcycle
405             call active_read_xy( fname, tmpfld2d, icvrec,
          call active_read_xy_loc( fname, tmpfld2d, icvrec,  
406       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
407       &                        mythid, dummy)       &                        mythid, dummy)
   
408           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
   
409  #endif /* ALLOW_SWFLUX_CONTROL */  #endif /* ALLOW_SWFLUX_CONTROL */
410    
411  #ifdef ALLOW_SWDOWN_CONTROL  #ifdef ALLOW_SWDOWN_CONTROL
# Line 425  c--   == end of interface == Line 414  c--   == end of interface ==
414           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
415           write(fname(1:80),'(3a,i10.10)')           write(fname(1:80),'(3a,i10.10)')
416       &        yadmark, xx_swdown_file(1:il),'.',optimcycle       &        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    #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,
431         &                        mythid, dummy)
432             xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
433    #endif /* ALLOW_LWFLUX_CONTROL */
434    
435           call active_read_xy_loc( fname, tmpfld2d, icvrec,  #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,       &                        doglobalread, ladinit, optimcycle,
443       &                        mythid, dummy)       &                        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 )           xx_comp = tmpfld2d( itilepos,jtilepos,itile,jtile )
457    #endif /* ALLOW_EVAP_CONTROL */
458    
459  #endif /* ALLOW_SWDOWN_CONTROL */  #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             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,
491         &                        doglobalread, ladinit, optimcycle,
492         &                        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 )
506    #endif /* ALLOW_SIAREA_CONTROL */
507    
508    #ifdef ALLOW_SIHEFF_CONTROL
509          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.
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22