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

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

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

revision 1.11 by heimbach, Mon Oct 27 22:32:55 2003 UTC revision 1.12 by heimbach, Tue Nov 4 20:47:42 2003 UTC
# Line 128  c     == local variables == Line 128  c     == local variables ==
128        _RL tmpplot2(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL tmpplot2(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
129        _RL tmpplot3(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL tmpplot3(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
130    
131          CHARACTER*(MAX_LEN_MBUF) msgBuf
132    
133  c     == end of interface ==  c     == end of interface ==
134  CEOP  CEOP
135    
# Line 178  cphadmtlm) Line 180  cphadmtlm)
180           grdchk_epsfac = 1. _d 0           grdchk_epsfac = 1. _d 0
181        end if        end if
182    
183        print *, 'ph-grd 3 -------------------------------'        print *, 'grad-res -------------------------------'
184        print ('(2a)'),        print ('(2a)'),
185       &     ' ph-grd 3  proc    #    i    j    k            fc ref',       &     ' grad-res  proc    #    i    j    k            fc ref',
186       &     '        fc + eps        fc - eps'       &     '        fc + eps        fc - eps'
187  #ifdef ALLOW_TANGENTLINEAR_RUN  #ifdef ALLOW_TANGENTLINEAR_RUN
188        print ('(2a)'),        print ('(2a)'),
189       &     ' ph-grd 3  proc    #    i    j    k          tlm grad',       &     ' grad-res  proc    #    i    j    k          tlm grad',
190       &     '         fd grad      1 - fd/tlm'       &     '         fd grad      1 - fd/tlm'
191  #else  #else
192        print ('(2a)'),        print ('(2a)'),
193       &     ' ph-grd 3  proc    #    i    j    k          adj grad',       &     ' grad-res  proc    #    i    j    k          adj grad',
194       &     '         fd grad      1 - fd/adj'       &     '         fd grad      1 - fd/adj'
195  #endif  #endif
196    
# Line 407  c*************************************** Line 409  c***************************************
409                 itestmem  ( ichknum ) = itest                 itestmem  ( ichknum ) = itest
410                 ierrmem   ( ichknum ) = ierr                 ierrmem   ( ichknum ) = ierr
411                                        
412                 print *, 'ph-grd 3 -------------------------------'                 print *, 'grad-res -------------------------------'
413                 print '(a,5I5,2x,3(1x,E15.9))', ' ph-grd 3 ',                 print '(a,5I5,2x,3(1x,E15.9))', ' grad-res ',
414       &              myprocid,ichknum,itilepos,jtilepos,layer,       &              myprocid,ichknum,itilepos,jtilepos,layer,
415       &              fcref, fcpertplus, fcpertminus       &              fcref, fcpertplus, fcpertminus
416  #ifdef ALLOW_TANGENTLINEAR_RUN  #ifdef ALLOW_TANGENTLINEAR_RUN
417                 print '(a,5I5,2x,3(1x,E15.9))', ' ph-grd 3 ',                 print '(a,5I5,2x,3(1x,E15.9))', ' grad-res ',
418       &              myprocid,ichknum,ichkmem(ichknum),       &              myprocid,ichknum,ichkmem(ichknum),
419       &              icompmem(ichknum),itestmem(ichknum),       &              icompmem(ichknum),itestmem(ichknum),
420       &              ftlxxmemo, gfd, ratio_ftl       &              ftlxxmemo, gfd, ratio_ftl
421                   WRITE(msgBuf,'(A34,2(1PE24.14,X))')
422         &              'precision_grdchk_result TLM ', fcref, ftlxxmemo
423                   CALL PRINT_MESSAGE
424         &              (msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
425  #else  #else
426                 print '(a,5I5,2x,3(1x,E15.9))', ' ph-grd 3 ',                 print '(a,5I5,2x,3(1x,E15.9))', ' grad-res ',
427       &              myprocid,ichknum,ichkmem(ichknum),       &              myprocid,ichknum,ichkmem(ichknum),
428       &              icompmem(ichknum),itestmem(ichknum),       &              icompmem(ichknum),itestmem(ichknum),
429       &              adxxmemo, gfd, ratio_ad       &              adxxmemo, gfd, ratio_ad
430                   WRITE(msgBuf,'(A34,2(1PE24.14,X))')
431         &              'precision_grdchk_result ADM ', fcref, adxxmemo
432                   CALL PRINT_MESSAGE
433         &              (msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
434  #endif  #endif
435    
436              endif              endif
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22