/[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.3 by heimbach, Thu May 30 22:47:26 2002 UTC revision 1.4 by heimbach, Sat Jul 13 02:55:58 2002 UTC
# Line 75  c--   == end of interface == Line 75  c--   == end of interface ==
75        dumiter      = 0        dumiter      = 0
76        dumtime      = 0. _d 0        dumtime      = 0. _d 0
77    
78        if      ( grdchkvarindex .eq. 1 ) then        if ( grdchkvarindex .eq. 0 ) then
79             STOP 'GRDCHK INDEX 0 NOT ALLOWED'
80    
81  #ifdef ALLOW_THETA0_CONTROL  #ifdef ALLOW_THETA0_CONTROL
82          else if ( grdchkvarindex .eq. 1 ) then
83           il=ilnblnk( xx_theta_file )           il=ilnblnk( xx_theta_file )
84           write(fname(1:80),'(80a)') ' '           write(fname(1:80),'(80a)') ' '
85           write(fname(1:80),'(2a,i10.10)')           write(fname(1:80),'(2a,i10.10)')
# Line 301  c--   == end of interface == Line 304  c--   == end of interface ==
304    
305  #endif /* ALLOW_HFACC3D_CONTROL */  #endif /* ALLOW_HFACC3D_CONTROL */
306  #endif /* ALLOW_HFACC_CONTROL */  #endif /* ALLOW_HFACC_CONTROL */
307    
308    #ifdef ALLOW_EFLUXY0_CONTROL
309          else if ( grdchkvarindex .eq. 21 ) then
310             il=ilnblnk( xx_efluxy_file )
311             write(fname(1:80),'(80a)') ' '
312             write(fname(1:80),'(2a,i10.10)')
313         &        xx_efluxy_file(1:il),'.',optimcycle
314    
315             call active_read_xyz( fname, tmpfld3d, 1,
316         &                         doglobalread, ladinit, optimcycle,
317         &                         mythid, dummy)
318    
319             xx_comp_ref = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
320             xx_comp_pert = xx_comp_ref + grdchk_eps
321             tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
322    
323             call active_write_xyz( fname, tmpfld3d, 1,
324         &                          optimcycle,
325         &                          mythid, dummy)
326    
327    #endif /* ALLOW_EFLUXY0_CONTROL */
328    
329    #ifdef ALLOW_EFLUXP0_CONTROL
330          else if ( grdchkvarindex .eq. 22 ) then
331             il=ilnblnk( xx_efluxp_file )
332             write(fname(1:80),'(80a)') ' '
333             write(fname(1:80),'(2a,i10.10)')
334         &        xx_efluxp_file(1:il),'.',optimcycle
335    
336             call active_read_xyz( fname, tmpfld3d, 1,
337         &                         doglobalread, ladinit, optimcycle,
338         &                         mythid, dummy)
339    
340             xx_comp_ref = tmpfld3d( itilepos,jtilepos,layer,itile,jtile )
341             xx_comp_pert = xx_comp_ref + grdchk_eps
342             tmpfld3d( itilepos,jtilepos,layer,itile,jtile ) = xx_comp_pert
343    
344             call active_write_xyz( fname, tmpfld3d, 1,
345         &                          optimcycle,
346         &                          mythid, dummy)
347    
348    #endif /* ALLOW_EFLUXP0_CONTROL */
349    
350        else        else
351  ce      --> this index does not exist yet.  ce      --> this index does not exist yet.
352        endif        endif
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22