/[MITgcm]/MITgcm/pkg/grdchk/grdchk_getxx.F
ViewVC logotype

Diff of /MITgcm/pkg/grdchk/grdchk_getxx.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.5 by heimbach, Mon Sep 16 18:11:58 2002 UTC revision 1.9 by heimbach, Fri Jul 18 21:10:16 2003 UTC
# Line 78  c--   == end of interface == Line 78  c--   == end of interface ==
78        ladinit      = .false.        ladinit      = .false.
79        dumiter      = 0        dumiter      = 0
80        dumtime      = 0. _d 0        dumtime      = 0. _d 0
81          write(fname(1:80),'(80a)') ' '
82    
83        if ( grdchkvarindex .eq. 0 ) then        if ( grdchkvarindex .eq. 0 ) then
84           STOP 'GRDCHK INDEX 0 NOT ALLOWED'           STOP 'GRDCHK INDEX 0 NOT ALLOWED'
# Line 85  c--   == end of interface == Line 86  c--   == end of interface ==
86  #ifdef ALLOW_THETA0_CONTROL  #ifdef ALLOW_THETA0_CONTROL
87        else if ( grdchkvarindex .eq. 1 ) then        else if ( grdchkvarindex .eq. 1 ) then
88           il=ilnblnk( xx_theta_file )           il=ilnblnk( xx_theta_file )
          write(fname(1:80),'(80a)') ' '  
89           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
90              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
91       &           yadmark, xx_theta_file(1:il),'.',optimcycle       &           yadmark, xx_theta_file(1:il),'.',optimcycle
# Line 94  c--   == end of interface == Line 94  c--   == end of interface ==
94       &           xx_theta_file(1:il),'.',optimcycle       &           xx_theta_file(1:il),'.',optimcycle
95           end if           end if
96    
97           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
98       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
99       &                         mythid, dummy)       &                         mythid, dummy)
100    
# Line 102  c--   == end of interface == Line 102  c--   == end of interface ==
102           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
103           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
104    
105           call active_write_xyz( fname, tmpfld3d, 1,           call active_write_xyz_loc( fname, tmpfld3d, 1,
106       &                          optimcycle,       &                          optimcycle,
107       &                          mythid, dummy)       &                          mythid, dummy)
108    
# Line 111  c--   == end of interface == Line 111  c--   == end of interface ==
111  #ifdef ALLOW_SALT0_CONTROL  #ifdef ALLOW_SALT0_CONTROL
112        else if ( grdchkvarindex .eq. 2 ) then        else if ( grdchkvarindex .eq. 2 ) then
113           il=ilnblnk( xx_salt_file )           il=ilnblnk( xx_salt_file )
          write(fname(1:80),'(80a)') ' '  
114           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
115              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
116       &           yadmark, xx_salt_file(1:il),'.',optimcycle       &           yadmark, xx_salt_file(1:il),'.',optimcycle
# Line 120  c--   == end of interface == Line 119  c--   == end of interface ==
119       &           xx_salt_file(1:il),'.',optimcycle       &           xx_salt_file(1:il),'.',optimcycle
120           end if           end if
121                    
122           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
123       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
124       &                         mythid, dummy)       &                         mythid, dummy)
125    
# Line 128  c--   == end of interface == Line 127  c--   == end of interface ==
127           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
128           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
129    
130           call active_write_xyz( fname, tmpfld3d, 1,           call active_write_xyz_loc( fname, tmpfld3d, 1,
131       &                          optimcycle,       &                          optimcycle,
132       &                          mythid, dummy)       &                          mythid, dummy)
133    
# Line 137  c--   == end of interface == Line 136  c--   == end of interface ==
136  #ifdef ALLOW_HFLUX_CONTROL  #ifdef ALLOW_HFLUX_CONTROL
137        else if ( grdchkvarindex .eq. 3 ) then        else if ( grdchkvarindex .eq. 3 ) then
138           il=ilnblnk( xx_hflux_file )           il=ilnblnk( xx_hflux_file )
          write(fname(1:80),'(80a)') ' '  
139           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
140              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
141       &           yadmark, xx_hflux_file(1:il),'.',optimcycle       &           yadmark, xx_hflux_file(1:il),'.',optimcycle
# Line 146  c--   == end of interface == Line 144  c--   == end of interface ==
144       &           xx_hflux_file(1:il),'.',optimcycle       &           xx_hflux_file(1:il),'.',optimcycle
145           end if           end if
146    
147           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
148       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
149       &                        mythid, dummy)       &                        mythid, dummy)
150    
# Line 154  c--   == end of interface == Line 152  c--   == end of interface ==
152           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
153           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
154    
155           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
156       &                         optimcycle,       &                         optimcycle,
157       &                         mythid, dummy)       &                         mythid, dummy)
158    
# Line 163  c--   == end of interface == Line 161  c--   == end of interface ==
161  #ifdef ALLOW_SFLUX_CONTROL  #ifdef ALLOW_SFLUX_CONTROL
162        else if ( grdchkvarindex .eq. 4 ) then        else if ( grdchkvarindex .eq. 4 ) then
163           il=ilnblnk( xx_sflux_file )           il=ilnblnk( xx_sflux_file )
          write(fname(1:80),'(80a)') ' '  
164           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
165              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
166       &           yadmark, xx_sflux_file(1:il),'.',optimcycle       &           yadmark, xx_sflux_file(1:il),'.',optimcycle
# Line 172  c--   == end of interface == Line 169  c--   == end of interface ==
169       &           xx_sflux_file(1:il),'.',optimcycle       &           xx_sflux_file(1:il),'.',optimcycle
170           end if           end if
171    
172           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
173       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
174       &                        mythid, dummy)       &                        mythid, dummy)
175    
# Line 180  c--   == end of interface == Line 177  c--   == end of interface ==
177           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
178           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
179    
180           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
181       &                         optimcycle,       &                         optimcycle,
182       &                         mythid, dummy)       &                         mythid, dummy)
183    
# Line 189  c--   == end of interface == Line 186  c--   == end of interface ==
186  #ifdef ALLOW_USTRESS_CONTROL  #ifdef ALLOW_USTRESS_CONTROL
187        else if ( grdchkvarindex .eq. 5 ) then        else if ( grdchkvarindex .eq. 5 ) then
188           il=ilnblnk( xx_tauu_file )           il=ilnblnk( xx_tauu_file )
          write(fname(1:80),'(80a)') ' '  
189           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
190              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
191       &           yadmark, xx_tauu_file(1:il),'.',optimcycle       &           yadmark, xx_tauu_file(1:il),'.',optimcycle
# Line 198  c--   == end of interface == Line 194  c--   == end of interface ==
194       &           xx_tauu_file(1:il),'.',optimcycle       &           xx_tauu_file(1:il),'.',optimcycle
195           end if           end if
196    
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    
# Line 206  c--   == end of interface == Line 202  c--   == end of interface ==
202           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
203           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
204    
205           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
206       &                         optimcycle,       &                         optimcycle,
207       &                         mythid, dummy)       &                         mythid, dummy)
208    
# Line 215  c--   == end of interface == Line 211  c--   == end of interface ==
211  #ifdef ALLOW_VSTRESS_CONTROL  #ifdef ALLOW_VSTRESS_CONTROL
212        else if ( grdchkvarindex .eq. 6 ) then        else if ( grdchkvarindex .eq. 6 ) then
213           il=ilnblnk( xx_tauv_file )           il=ilnblnk( xx_tauv_file )
          write(fname(1:80),'(80a)') ' '  
214           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
215              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
216       &           yadmark, xx_tauv_file(1:il),'.',optimcycle       &           yadmark, xx_tauv_file(1:il),'.',optimcycle
# Line 224  c--   == end of interface == Line 219  c--   == end of interface ==
219       &           xx_tauv_file(1:il),'.',optimcycle       &           xx_tauv_file(1:il),'.',optimcycle
220           end if           end if
221    
222           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
223       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
224       &                        mythid, dummy)       &                        mythid, dummy)
225    
# Line 232  c--   == end of interface == Line 227  c--   == end of interface ==
227           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
228           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
229    
230           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
231       &                         optimcycle,       &                         optimcycle,
232       &                         mythid, dummy)       &                         mythid, dummy)
233    
234  #endif /* ALLOW_VSTRESS_CONTROL */  #endif /* ALLOW_VSTRESS_CONTROL */
235    
236    #ifdef ALLOW_ATEMP_CONTROL
237          else if ( grdchkvarindex .eq. 7 ) then
238             il=ilnblnk( xx_atemp_file )
239             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
240                write(fname(1:80),'(3a,i10.10)')
241         &           yadmark, xx_atemp_file(1:il),'.',optimcycle
242             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
243                write(fname(1:80),'(2a,i10.10)')
244         &           xx_atemp_file(1:il),'.',optimcycle
245             end if
246    
247             call active_read_xy_loc( fname, tmpfld2d, icvrec,
248         &                        doglobalread, ladinit, optimcycle,
249         &                        mythid, dummy)
250    
251             xx_comp_ref = tmpfld2d( itilepos,jtilepos,itile,jtile )
252             xx_comp_pert = xx_comp_ref + localEps
253             tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
254    
255             call active_write_xy_loc( fname, tmpfld2d, icvrec,
256         &                         optimcycle,
257         &                         mythid, dummy)
258    
259    #endif /* ALLOW_ATEMP_CONTROL */
260    
261    #ifdef ALLOW_AQH_CONTROL
262          else if ( grdchkvarindex .eq. 8 ) then
263             il=ilnblnk( xx_aqh_file )
264             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
265                write(fname(1:80),'(3a,i10.10)')
266         &           yadmark, xx_aqh_file(1:il),'.',optimcycle
267             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
268                write(fname(1:80),'(2a,i10.10)')
269         &           xx_aqh_file(1:il),'.',optimcycle
270             end if
271    
272             call active_read_xy_loc( fname, tmpfld2d, icvrec,
273         &                        doglobalread, ladinit, optimcycle,
274         &                        mythid, dummy)
275    
276             xx_comp_ref = tmpfld2d( itilepos,jtilepos,itile,jtile )
277             xx_comp_pert = xx_comp_ref + localEps
278             tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
279    
280             call active_write_xy_loc( fname, tmpfld2d, icvrec,
281         &                         optimcycle,
282         &                         mythid, dummy)
283    
284    #endif /* ALLOW_AQH_CONTROL */
285    
286    #ifdef ALLOW_UWIND_CONTROL
287          else if ( grdchkvarindex .eq. 9 ) then
288             il=ilnblnk( xx_uwind_file )
289             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
290                write(fname(1:80),'(3a,i10.10)')
291         &           yadmark, xx_uwind_file(1:il),'.',optimcycle
292             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
293                write(fname(1:80),'(2a,i10.10)')
294         &           xx_uwind_file(1:il),'.',optimcycle
295             end if
296    
297             call active_read_xy_loc( fname, tmpfld2d, icvrec,
298         &                        doglobalread, ladinit, optimcycle,
299         &                        mythid, dummy)
300    
301             xx_comp_ref = tmpfld2d( itilepos,jtilepos,itile,jtile )
302             xx_comp_pert = xx_comp_ref + localEps
303             tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
304    
305             call active_write_xy_loc( fname, tmpfld2d, icvrec,
306         &                         optimcycle,
307         &                         mythid, dummy)
308    
309    #endif /* ALLOW_UWIND_CONTROL */
310    
311    #ifdef ALLOW_VWIND_CONTROL
312          else if ( grdchkvarindex .eq. 10 ) then
313             il=ilnblnk( xx_vwind_file )
314             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
315                write(fname(1:80),'(3a,i10.10)')
316         &           yadmark, xx_vwind_file(1:il),'.',optimcycle
317             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
318                write(fname(1:80),'(2a,i10.10)')
319         &           xx_vwind_file(1:il),'.',optimcycle
320             end if
321    
322             call active_read_xy_loc( fname, tmpfld2d, icvrec,
323         &                        doglobalread, ladinit, optimcycle,
324         &                        mythid, dummy)
325    
326             xx_comp_ref = tmpfld2d( itilepos,jtilepos,itile,jtile )
327             xx_comp_pert = xx_comp_ref + localEps
328             tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
329    
330             call active_write_xy_loc( fname, tmpfld2d, icvrec,
331         &                         optimcycle,
332         &                         mythid, dummy)
333    
334    #endif /* ALLOW_VWIND_CONTROL */
335    
336    #ifdef ALLOW_OBCSN_CONTROL
337          else if ( grdchkvarindex .eq. 11 ) then
338             il=ilnblnk( xx_obcsn_file )
339             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
340                write(fname(1:80),'(3a,i10.10)')
341         &           yadmark, xx_obcsn_file(1:il),'.',optimcycle
342             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
343                write(fname(1:80),'(2a,i10.10)')
344         &           xx_obcsn_file(1:il),'.',optimcycle
345             end if
346    
347             call active_read_xz_loc( fname, tmpfldxz, icvrec,
348         &                        doglobalread, ladinit, optimcycle,
349         &                        mythid, dummy)
350    
351             xx_comp_ref = tmpfldxz( itilepos,layer,itile,jtile )
352             xx_comp_pert = xx_comp_ref + localEps
353             tmpfldxz( itilepos,layer,itile,jtile ) = xx_comp_pert
354    
355             call active_write_xz_loc( fname, tmpfldxz, icvrec,
356         &                         optimcycle,
357         &                         mythid, dummy)
358    
359    #endif /* ALLOW_OBCSN_CONTROL */
360    
361    #ifdef ALLOW_OBCSS_CONTROL
362          else if ( grdchkvarindex .eq. 12 ) then
363             il=ilnblnk( xx_obcss_file )
364             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
365                write(fname(1:80),'(3a,i10.10)')
366         &           yadmark, xx_obcss_file(1:il),'.',optimcycle
367             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
368                write(fname(1:80),'(2a,i10.10)')
369         &           xx_obcss_file(1:il),'.',optimcycle
370             end if
371    
372             call active_read_xz_loc( fname, tmpfldxz, icvrec,
373         &                        doglobalread, ladinit, optimcycle,
374         &                        mythid, dummy)
375    
376             xx_comp_ref = tmpfldxz( itilepos,layer,itile,jtile )
377             xx_comp_pert = xx_comp_ref + localEps
378             tmpfldxz( itilepos,layer,itile,jtile ) = xx_comp_pert
379    
380             call active_write_xz_loc( fname, tmpfldxz, icvrec,
381         &                         optimcycle,
382         &                         mythid, dummy)
383    
384    #endif /* ALLOW_OBCSS_CONTROL */
385    
386    #ifdef ALLOW_OBCSW_CONTROL
387          else if ( grdchkvarindex .eq. 13 ) then
388             il=ilnblnk( xx_obcsw_file )
389             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
390                write(fname(1:80),'(3a,i10.10)')
391         &           yadmark, xx_obcsw_file(1:il),'.',optimcycle
392             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
393                write(fname(1:80),'(2a,i10.10)')
394         &           xx_obcsw_file(1:il),'.',optimcycle
395             end if
396    
397             call active_read_yz_loc( fname, tmpfldyz, icvrec,
398         &                        doglobalread, ladinit, optimcycle,
399         &                        mythid, dummy)
400    
401             xx_comp_ref = tmpfldyz( jtilepos,layer,itile,jtile )
402             xx_comp_pert = xx_comp_ref + localEps
403             tmpfldyz( jtilepos,layer,itile,jtile ) = xx_comp_pert
404    
405             call active_write_yz_loc( fname, tmpfldyz, icvrec,
406         &                         optimcycle,
407         &                         mythid, dummy)
408    
409    #endif /* ALLOW_OBCSW_CONTROL */
410    
411    #ifdef ALLOW_OBCSE_CONTROL
412          else if ( grdchkvarindex .eq. 14 ) then
413             il=ilnblnk( xx_obcse_file )
414             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
415                write(fname(1:80),'(3a,i10.10)')
416         &           yadmark, xx_obcse_file(1:il),'.',optimcycle
417             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
418                write(fname(1:80),'(2a,i10.10)')
419         &           xx_obcse_file(1:il),'.',optimcycle
420             end if
421    
422             call active_read_yz_loc( fname, tmpfldyz, icvrec,
423         &                        doglobalread, ladinit, optimcycle,
424         &                        mythid, dummy)
425    
426             xx_comp_ref = tmpfldyz( jtilepos,layer,itile,jtile )
427             xx_comp_pert = xx_comp_ref + localEps
428             tmpfldyz( jtilepos,layer,itile,jtile ) = xx_comp_pert
429    
430             call active_write_yz_loc( fname, tmpfldyz, icvrec,
431         &                         optimcycle,
432         &                         mythid, dummy)
433    
434    #endif /* ALLOW_OBCSE_CONTROL */
435    
436  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_TR10_CONTROL
437        else if ( grdchkvarindex .eq. 17 ) then        else if ( grdchkvarindex .eq. 17 ) then
438           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_tr1_file )
          write(fname(1:80),'(80a)') ' '  
439           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
440              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
441       &           yadmark, xx_tr1_file(1:il),'.',optimcycle       &           yadmark, xx_tr1_file(1:il),'.',optimcycle
# Line 250  c--   == end of interface == Line 444  c--   == end of interface ==
444       &           xx_tr1_file(1:il),'.',optimcycle       &           xx_tr1_file(1:il),'.',optimcycle
445           end if           end if
446    
447           call active_read_xyz( fname, tmpfld3d, icvrec,           call active_read_xyz_loc( fname, tmpfld3d, icvrec,
448       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
449       &                        mythid, dummy)       &                        mythid, dummy)
450    
# Line 258  c--   == end of interface == Line 452  c--   == end of interface ==
452           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
453           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
454    
455           call active_write_xyz( fname, tmpfld3d, icvrec,           call active_write_xyz_loc( fname, tmpfld3d, icvrec,
456       &                         optimcycle,       &                         optimcycle,
457       &                         mythid, dummy)       &                         mythid, dummy)
458    
# Line 267  c--   == end of interface == Line 461  c--   == end of interface ==
461  #ifdef ALLOW_SST0_CONTROL  #ifdef ALLOW_SST0_CONTROL
462        else if ( grdchkvarindex .eq. 18 ) then        else if ( grdchkvarindex .eq. 18 ) then
463           il=ilnblnk( xx_sst_file )           il=ilnblnk( xx_sst_file )
          write(fname(1:80),'(80a)') ' '  
464           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
465              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
466       &           yadmark, xx_sst_file(1:il),'.',optimcycle       &           yadmark, xx_sst_file(1:il),'.',optimcycle
# Line 276  c--   == end of interface == Line 469  c--   == end of interface ==
469       &           xx_sst_file(1:il),'.',optimcycle       &           xx_sst_file(1:il),'.',optimcycle
470           end if           end if
471    
472           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
473       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
474       &                        mythid, dummy)       &                        mythid, dummy)
475    
# Line 284  c--   == end of interface == Line 477  c--   == end of interface ==
477           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
478           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
479    
480           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
481       &                         optimcycle,       &                         optimcycle,
482       &                         mythid, dummy)       &                         mythid, dummy)
483    
# Line 293  c--   == end of interface == Line 486  c--   == end of interface ==
486  #ifdef ALLOW_SSS0_CONTROL  #ifdef ALLOW_SSS0_CONTROL
487        else if ( grdchkvarindex .eq. 19 ) then        else if ( grdchkvarindex .eq. 19 ) then
488           il=ilnblnk( xx_sss_file )           il=ilnblnk( xx_sss_file )
          write(fname(1:80),'(80a)') ' '  
489           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
490              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
491       &           yadmark, xx_sss_file(1:il),'.',optimcycle       &           yadmark, xx_sss_file(1:il),'.',optimcycle
# Line 302  c--   == end of interface == Line 494  c--   == end of interface ==
494       &           xx_sss_file(1:il),'.',optimcycle       &           xx_sss_file(1:il),'.',optimcycle
495           end if           end if
496    
497           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
498       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
499       &                        mythid, dummy)       &                        mythid, dummy)
500    
# Line 310  c--   == end of interface == Line 502  c--   == end of interface ==
502           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
503           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
504    
505           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
506       &                         optimcycle,       &                         optimcycle,
507       &                         mythid, dummy)       &                         mythid, dummy)
508    
# Line 319  c--   == end of interface == Line 511  c--   == end of interface ==
511  #ifdef ALLOW_HFACC_CONTROL  #ifdef ALLOW_HFACC_CONTROL
512        else if ( grdchkvarindex .eq. 20 ) then        else if ( grdchkvarindex .eq. 20 ) then
513           il=ilnblnk( xx_hfacc_file )           il=ilnblnk( xx_hfacc_file )
          write(fname(1:80),'(80a)') ' '  
514           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then           if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
515              write(fname(1:80),'(3a,i10.10)')              write(fname(1:80),'(3a,i10.10)')
516       &           yadmark, xx_hfacc_file(1:il),'.',optimcycle       &           yadmark, xx_hfacc_file(1:il),'.',optimcycle
# Line 330  c--   == end of interface == Line 521  c--   == end of interface ==
521    
522  #ifdef ALLOW_HFACC3D_CONTROL  #ifdef ALLOW_HFACC3D_CONTROL
523    
524           call active_read_xyz( fname, tmpfld3d, icvrec,           call active_read_xyz_loc( fname, tmpfld3d, icvrec,
525       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
526       &                        mythid, dummy)       &                        mythid, dummy)
527    
# Line 338  c--   == end of interface == Line 529  c--   == end of interface ==
529           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
530           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
531    
532           call active_write_xyz( fname, tmpfld3d, icvrec,           call active_write_xyz_loc( fname, tmpfld3d, icvrec,
533       &                         optimcycle,       &                         optimcycle,
534       &                         mythid, dummy)       &                         mythid, dummy)
535    
536  #else  #else
537    
538           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
539       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
540       &                        mythid, dummy)       &                        mythid, dummy)
541    
# Line 352  c--   == end of interface == Line 543  c--   == end of interface ==
543           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
544           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
545    
546           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
547       &                         optimcycle,       &                         optimcycle,
548       &                         mythid, dummy)       &                         mythid, dummy)
549    
# Line 371  c--   == end of interface == Line 562  c--   == end of interface ==
562       &           xx_efluxy_file(1:il),'.',optimcycle       &           xx_efluxy_file(1:il),'.',optimcycle
563           end if           end if
564    
565           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
566       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
567       &                         mythid, dummy)       &                         mythid, dummy)
568    
# Line 379  c--   == end of interface == Line 570  c--   == end of interface ==
570           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
571           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
572    
573           call active_write_xyz( fname, tmpfld3d, 1,           call active_write_xyz_loc( fname, tmpfld3d, 1,
574       &                          optimcycle,       &                          optimcycle,
575       &                          mythid, dummy)       &                          mythid, dummy)
576    
# Line 397  c--   == end of interface == Line 588  c--   == end of interface ==
588       &           xx_efluxp_file(1:il),'.',optimcycle       &           xx_efluxp_file(1:il),'.',optimcycle
589           end if           end if
590    
591           call active_read_xyz( fname, tmpfld3d, 1,           call active_read_xyz_loc( fname, tmpfld3d, 1,
592       &                         doglobalread, ladinit, optimcycle,       &                         doglobalread, ladinit, optimcycle,
593       &                         mythid, dummy)       &                         mythid, dummy)
594    
# Line 405  c--   == end of interface == Line 596  c--   == end of interface ==
596           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
597           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert           tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
598    
599           call active_write_xyz( fname, tmpfld3d, 1,           call active_write_xyz_loc( fname, tmpfld3d, 1,
600       &                          optimcycle,       &                          optimcycle,
601       &                          mythid, dummy)       &                          mythid, dummy)
602    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22