/[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.8 by heimbach, Tue Jun 24 16:08:45 2003 UTC revision 1.11 by heimbach, Thu Jul 28 13:54:36 2005 UTC
# Line 54  c     == routine arguments == Line 54  c     == routine arguments ==
54        _RL     localEps        _RL     localEps
55        integer mythid        integer mythid
56    
57  #ifdef ALLOW_GRADIENT_CHECK  #ifdef ALLOW_GRDCHK
58  c     == local variables ==  c     == local variables ==
59    
60        integer il        integer il
# 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 119  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 127  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 144  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 152  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 169  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 177  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 194  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 202  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 219  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 227  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    
# Line 244  c--   == end of interface == Line 244  c--   == end of interface ==
244       &           xx_atemp_file(1:il),'.',optimcycle       &           xx_atemp_file(1:il),'.',optimcycle
245           end if           end if
246    
247           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
248       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
249       &                        mythid, dummy)       &                        mythid, dummy)
250    
# Line 252  c--   == end of interface == Line 252  c--   == end of interface ==
252           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
253           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
254    
255           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
256       &                         optimcycle,       &                         optimcycle,
257       &                         mythid, dummy)       &                         mythid, dummy)
258    
# Line 269  c--   == end of interface == Line 269  c--   == end of interface ==
269       &           xx_aqh_file(1:il),'.',optimcycle       &           xx_aqh_file(1:il),'.',optimcycle
270           end if           end if
271    
272           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
273       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
274       &                        mythid, dummy)       &                        mythid, dummy)
275    
# Line 277  c--   == end of interface == Line 277  c--   == end of interface ==
277           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
278           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
279    
280           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
281       &                         optimcycle,       &                         optimcycle,
282       &                         mythid, dummy)       &                         mythid, dummy)
283    
# Line 294  c--   == end of interface == Line 294  c--   == end of interface ==
294       &           xx_uwind_file(1:il),'.',optimcycle       &           xx_uwind_file(1:il),'.',optimcycle
295           end if           end if
296    
297           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
298       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
299       &                        mythid, dummy)       &                        mythid, dummy)
300    
# Line 302  c--   == end of interface == Line 302  c--   == end of interface ==
302           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
303           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
304    
305           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
306       &                         optimcycle,       &                         optimcycle,
307       &                         mythid, dummy)       &                         mythid, dummy)
308    
# Line 319  c--   == end of interface == Line 319  c--   == end of interface ==
319       &           xx_vwind_file(1:il),'.',optimcycle       &           xx_vwind_file(1:il),'.',optimcycle
320           end if           end if
321    
322           call active_read_xy( fname, tmpfld2d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
323       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
324       &                        mythid, dummy)       &                        mythid, dummy)
325    
# Line 327  c--   == end of interface == Line 327  c--   == end of interface ==
327           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
328           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
329    
330           call active_write_xy( fname, tmpfld2d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
331       &                         optimcycle,       &                         optimcycle,
332       &                         mythid, dummy)       &                         mythid, dummy)
333    
# Line 344  c--   == end of interface == Line 344  c--   == end of interface ==
344       &           xx_obcsn_file(1:il),'.',optimcycle       &           xx_obcsn_file(1:il),'.',optimcycle
345           end if           end if
346    
347           call active_read_xz( fname, tmpfldxz, icvrec,           call active_read_xz_loc( fname, tmpfldxz, icvrec,
348       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
349       &                        mythid, dummy)       &                        mythid, dummy)
350    
# Line 352  c--   == end of interface == Line 352  c--   == end of interface ==
352           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
353           tmpfldxz( itilepos,layer,itile,jtile ) = xx_comp_pert           tmpfldxz( itilepos,layer,itile,jtile ) = xx_comp_pert
354    
355           call active_write_xz( fname, tmpfldxz, icvrec,           call active_write_xz_loc( fname, tmpfldxz, icvrec,
356       &                         optimcycle,       &                         optimcycle,
357       &                         mythid, dummy)       &                         mythid, dummy)
358    
# Line 369  c--   == end of interface == Line 369  c--   == end of interface ==
369       &           xx_obcss_file(1:il),'.',optimcycle       &           xx_obcss_file(1:il),'.',optimcycle
370           end if           end if
371    
372           call active_read_xz( fname, tmpfldxz, icvrec,           call active_read_xz_loc( fname, tmpfldxz, icvrec,
373       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
374       &                        mythid, dummy)       &                        mythid, dummy)
375    
# Line 377  c--   == end of interface == Line 377  c--   == end of interface ==
377           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
378           tmpfldxz( itilepos,layer,itile,jtile ) = xx_comp_pert           tmpfldxz( itilepos,layer,itile,jtile ) = xx_comp_pert
379    
380           call active_write_xz( fname, tmpfldxz, icvrec,           call active_write_xz_loc( fname, tmpfldxz, icvrec,
381       &                         optimcycle,       &                         optimcycle,
382       &                         mythid, dummy)       &                         mythid, dummy)
383    
# Line 394  c--   == end of interface == Line 394  c--   == end of interface ==
394       &           xx_obcsw_file(1:il),'.',optimcycle       &           xx_obcsw_file(1:il),'.',optimcycle
395           end if           end if
396    
397           call active_read_yz( fname, tmpfldyz, icvrec,           call active_read_yz_loc( fname, tmpfldyz, icvrec,
398       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
399       &                        mythid, dummy)       &                        mythid, dummy)
400    
# Line 402  c--   == end of interface == Line 402  c--   == end of interface ==
402           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
403           tmpfldyz( jtilepos,layer,itile,jtile ) = xx_comp_pert           tmpfldyz( jtilepos,layer,itile,jtile ) = xx_comp_pert
404    
405           call active_write_yz( fname, tmpfldyz, icvrec,           call active_write_yz_loc( fname, tmpfldyz, icvrec,
406       &                         optimcycle,       &                         optimcycle,
407       &                         mythid, dummy)       &                         mythid, dummy)
408    
# Line 419  c--   == end of interface == Line 419  c--   == end of interface ==
419       &           xx_obcse_file(1:il),'.',optimcycle       &           xx_obcse_file(1:il),'.',optimcycle
420           end if           end if
421    
422           call active_read_yz( fname, tmpfldyz, icvrec,           call active_read_yz_loc( fname, tmpfldyz, icvrec,
423       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
424       &                        mythid, dummy)       &                        mythid, dummy)
425    
# Line 427  c--   == end of interface == Line 427  c--   == end of interface ==
427           xx_comp_pert = xx_comp_ref + localEps           xx_comp_pert = xx_comp_ref + localEps
428           tmpfldyz( jtilepos,layer,itile,jtile ) = xx_comp_pert           tmpfldyz( jtilepos,layer,itile,jtile ) = xx_comp_pert
429    
430           call active_write_yz( fname, tmpfldyz, icvrec,           call active_write_yz_loc( fname, tmpfldyz, icvrec,
431       &                         optimcycle,       &                         optimcycle,
432       &                         mythid, dummy)       &                         mythid, dummy)
433    
434  #endif /* ALLOW_OBCSE_CONTROL */  #endif /* ALLOW_OBCSE_CONTROL */
435    
436  #ifdef ALLOW_TR10_CONTROL  #ifdef ALLOW_PRECIP_CONTROL
437        else if ( grdchkvarindex .eq. 17 ) then        else if ( grdchkvarindex .eq. 17 ) then
438           il=ilnblnk( xx_tr1_file )           il=ilnblnk( xx_precip_file )
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_precip_file(1:il),'.',optimcycle
442           else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then           else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
443              write(fname(1:80),'(2a,i10.10)')              write(fname(1:80),'(2a,i10.10)')
444       &           xx_tr1_file(1:il),'.',optimcycle       &           xx_precip_file(1:il),'.',optimcycle
445           end if           end if
446    
447           call active_read_xyz( fname, tmpfld3d, icvrec,           call active_read_xy_loc( fname, tmpfld2d, icvrec,
448       &                        doglobalread, ladinit, optimcycle,       &                        doglobalread, ladinit, optimcycle,
449       &                        mythid, dummy)       &                        mythid, dummy)
450    
451           xx_comp_ref = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )           xx_comp_ref = tmpfld2d( itilepos,jtilepos,itile,jtile )
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           tmpfld2d( itilepos,jtilepos,itile,jtile ) = xx_comp_pert
454    
455           call active_write_xyz( fname, tmpfld3d, icvrec,           call active_write_xy_loc( fname, tmpfld2d, icvrec,
456       &                         optimcycle,       &                         optimcycle,
457       &                         mythid, dummy)       &                         mythid, dummy)
458    
459  #endif /* ALLOW_TR10_CONTROL */  #endif /* ALLOW_PRECIP_CONTROL */
460    
461  #ifdef ALLOW_SST0_CONTROL  #ifdef ALLOW_SST0_CONTROL
462        else if ( grdchkvarindex .eq. 18 ) then        else if ( grdchkvarindex .eq. 18 ) then
# Line 469  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 477  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 494  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 502  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 521  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 529  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 543  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 562  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 570  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 588  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 596  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    
603  #endif /* ALLOW_EFLUXP0_CONTROL */  #endif /* ALLOW_EFLUXP0_CONTROL */
604    
605    #ifdef ALLOW_TR10_CONTROL
606          else if ( grdchkvarindex .eq. 32 ) then
607             il=ilnblnk( xx_tr1_file )
608             if ( theSimulationMode .EQ. TANGENT_SIMULATION ) then
609                write(fname(1:80),'(3a,i10.10)')
610         &           yadmark, xx_tr1_file(1:il),'.',optimcycle
611             else if ( theSimulationMode .EQ. FORWARD_SIMULATION ) then
612                write(fname(1:80),'(2a,i10.10)')
613         &           xx_tr1_file(1:il),'.',optimcycle
614             end if
615    
616             call active_read_xyz_loc( fname, tmpfld3d, icvrec,
617         &                        doglobalread, ladinit, optimcycle,
618         &                        mythid, dummy)
619    
620             xx_comp_ref = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
621             xx_comp_pert = xx_comp_ref + localEps
622             tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
623    
624             call active_write_xyz_loc( fname, tmpfld3d, icvrec,
625         &                         optimcycle,
626         &                         mythid, dummy)
627    
628    #endif /* ALLOW_TR10_CONTROL */
629    
630        else        else
631  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
632        endif        endif
633    
634  #endif /* ALLOW_GRADIENT_CHECK */  #endif /* ALLOW_GRDCHK */
635    
636        end        end
637    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22