/[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.18 by heimbach, Fri Nov 10 04:57:10 2006 UTC revision 1.19 by heimbach, Sun Nov 19 23:17:17 2006 UTC
# Line 109  c     == local variables == Line 109  c     == local variables ==
109        integer obcspos        integer obcspos
110        integer itilepos        integer itilepos
111        integer jtilepos        integer jtilepos
112          integer icglo
113        integer itest        integer itest
114        integer ierr        integer ierr
115        integer ierr_grdchk        integer ierr_grdchk
# Line 184  cphadmtlm) Line 185  cphadmtlm)
185    
186        print *, 'grad-res -------------------------------'        print *, 'grad-res -------------------------------'
187        print ('(2a)'),        print ('(2a)'),
188       &    ' grad-res  proc    #    i    j    k iobc',       &    ' grad-res  proc    #    i    j    k   bi   bj iobc',
189       &    '               fc ref           fc + eps           fc - eps'       &    '               fc ref           fc + eps           fc - eps'
190  #ifdef ALLOW_TANGENTLINEAR_RUN  #ifdef ALLOW_TANGENTLINEAR_RUN
191        print ('(2a)'),        print ('(2a)'),
192       &    ' grad-res  proc    #    i    j    k iobc',       &    ' grad-res  proc    #    i    j    k   bi   bj iobc',
193       &    '             tlm grad            fd grad         1 - fd/tlm'       &    '             tlm grad            fd grad         1 - fd/tlm'
194  #else  #else
195        print ('(2a)'),        print ('(2a)'),
196       &    ' grad-res  proc    #    i    j    k iobc',       &    ' grad-res  proc    #    i    j    k   bi   bj iobc',
197       &    '             adj grad            fd grad         1 - fd/adj'       &    '             adj grad            fd grad         1 - fd/adj'
198  #endif  #endif
199    
# Line 207  c--   gradient checks. Line 208  c--   gradient checks.
208    
209           ichknum = (icomp - nbeg)/nstep + 1           ichknum = (icomp - nbeg)/nstep + 1
210    
211    cph(
212    cph-print         print *, 'ph-grd _main: nbeg, icomp, ichknum ',
213    cph-print     &        nbeg, icomp, ichknum
214    cph)
215           if (ichknum .le. maxgrdchecks ) then           if (ichknum .le. maxgrdchecks ) then
216    
217  c--         Determine the location of icomp on the grid.  c--         Determine the location of icomp on the grid.
218              if ( myProcId .EQ. grdchkwhichproc ) then              if ( myProcId .EQ. grdchkwhichproc ) then
219                 call grdchk_loc( icomp, ichknum,                 call grdchk_loc( icomp, ichknum,
220       &              icvrec, itile, jtile, layer, obcspos,       &              icvrec, itile, jtile, layer, obcspos,
221       &              itilepos, jtilepos, itest, ierr,       &              itilepos, jtilepos, icglo, itest, ierr,
222       &              mythid )       &              mythid )
223  cph(  cph(
224                 print *, 'ph-grd ----- back from loc -----',  cph-print               print *, 'ph-grd ----- back from loc -----',
225       &             icvrec, itilepos, jtilepos, layer, obcspos  cph-print     &             icvrec, itilepos, jtilepos, layer, obcspos
226  cph)  cph)
227              endif              endif
228              _BARRIER              _BARRIER
# Line 420  c*************************************** Line 425  c***************************************
425                 ichkmem   ( ichknum ) = ichknum                 ichkmem   ( ichknum ) = ichknum
426                 itestmem  ( ichknum ) = itest                 itestmem  ( ichknum ) = itest
427                 ierrmem   ( ichknum ) = ierr                 ierrmem   ( ichknum ) = ierr
428                                     icglomem  ( ichknum ) = icglo
429    
430                 print *, 'grad-res -------------------------------'                 print *, 'grad-res -------------------------------'
431                 print '(a,6I5,2x,3(1x,E18.12))', ' grad-res ',                 print '(a,8I5,2x,3(1x,E18.12))', ' grad-res ',
432       &              myprocid,ichknum,itilepos,jtilepos,layer,obcspos,       &              myprocid,ichknum,itilepos,jtilepos,
433         &              layer,itile,jtile,obcspos,
434       &              fcref, fcpertplus, fcpertminus       &              fcref, fcpertplus, fcpertminus
435  #ifdef ALLOW_TANGENTLINEAR_RUN  #ifdef ALLOW_TANGENTLINEAR_RUN
436                 print '(a,6I5,2x,3(1x,E18.12))', ' grad-res ',                 print '(a,8I5,2x,3(1x,E18.12))', ' grad-res ',
437       &              myprocid,ichknum,ichkmem(ichknum),       &              myprocid,ichknum,ichkmem(ichknum),
438       &              icompmem(ichknum),itestmem(ichknum),obcspos,       &              icompmem(ichknum),itestmem(ichknum),
439         &              bimem(ichknum),bjmem(ichknum),iobcsmem(ichknum),
440       &              ftlxxmemo, gfd, ratio_ftl       &              ftlxxmemo, gfd, ratio_ftl
441                 WRITE(msgBuf,'(A34,2(1PE24.14,X))')                 WRITE(msgBuf,'(A34,2(1PE24.14,X))')
442       &              'precision_grdchk_result TLM ', fcref, ftlxxmemo       &              'precision_grdchk_result TLM ', fcref, ftlxxmemo
443                 CALL PRINT_MESSAGE                 CALL PRINT_MESSAGE
444       &              (msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)       &              (msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
445  #else  #else
446                 print '(a,6I5,2x,3(1x,E18.12))', ' grad-res ',                 print '(a,8I5,2x,3(1x,E18.12))', ' grad-res ',
447       &              myprocid,ichknum,ichkmem(ichknum),       &              myprocid,ichknum,ichkmem(ichknum),
448       &              icompmem(ichknum),itestmem(ichknum),obcspos,       &              icompmem(ichknum),itestmem(ichknum),
449         &              bimem(ichknum),bjmem(ichknum),obcspos,
450       &              adxxmemo, gfd, ratio_ad       &              adxxmemo, gfd, ratio_ad
451                 WRITE(msgBuf,'(A34,2(1PE24.14,X))')                 WRITE(msgBuf,'(A34,2(1PE24.14,X))')
452       &              'precision_grdchk_result ADM ', fcref, adxxmemo       &              'precision_grdchk_result ADM ', fcref, adxxmemo
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22