/[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.7 by heimbach, Fri Feb 28 02:34:56 2003 UTC revision 1.9 by heimbach, Fri Jul 18 21:10:16 2003 UTC
# 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    
334  #endif /* ALLOW_VWIND_CONTROL */  #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 )
# Line 344  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 352  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 369  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 377  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 394  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 402  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 421  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 429  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 443  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 462  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 470  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 488  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 496  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.7  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22