/[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.25 by jmc, Sat Sep 15 17:32:41 2007 UTC revision 1.28 by heimbach, Fri May 29 06:14:47 2009 UTC
# Line 156  cph   Gradient via adjoint has already b Line 156  cph   Gradient via adjoint has already b
156  cph   and so has unperturbed cost function,  cph   and so has unperturbed cost function,
157  cph   assuming all xx_ fields are initialised to zero.  cph   assuming all xx_ fields are initialised to zero.
158    
159          ierr      = 0
160        ierr_grdchk = 0        ierr_grdchk = 0
161          adxxmemo  = 0.
162          ftlxxmemo = 0.
163  #ifdef ALLOW_ADMTLM  #ifdef ALLOW_ADMTLM
164        fcref = objf_state_final(idep,jdep,1,1,1)        fcref = objf_state_final(idep,jdep,1,1,1)
165  #else  #else
# Line 250  c--   (B): Get gradient component g_fc f Line 253  c--   (B): Get gradient component g_fc f
253  c******************************************************  c******************************************************
254  c--  c--
255  c--   1. perturb control vector component: xx(i)=1.  c--   1. perturb control vector component: xx(i)=1.
             ftlxxmemo = 0.  
256    
257              if ( myProcId .EQ. grdchkwhichproc .AND.              if ( myProcId .EQ. grdchkwhichproc .AND.
258       &           ierr .EQ. 0 ) then       &           ierr .EQ. 0 ) then
# Line 472  c Line 474  c
474       &              icompmem(ichknum),itestmem(ichknum),       &              icompmem(ichknum),itestmem(ichknum),
475       &              bimem(ichknum),bjmem(ichknum),obcspos,       &              bimem(ichknum),bjmem(ichknum),obcspos,
476       &              adxxmemo, gfd, ratio_ad       &              adxxmemo, gfd, ratio_ad
477                 WRITE(msgBuf,'(A34,2(1PE24.14,X))')  c              WRITE(msgBuf,'(A34,2(1PE24.14,X))')
478       &              'precision_grdchk_result ADM ', fcref, adxxmemo  c    &              'precision_grdchk_result ADM ', fcref, adxxmemo
479                 CALL PRINT_MESSAGE  c              CALL PRINT_MESSAGE
480       &              (msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)  c    &              (msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
 c  
481                 WRITE(msgBuf,'(A34,1PE24.14)')                 WRITE(msgBuf,'(A34,1PE24.14)')
482       &              ' ADM  precision_derivative_cost =', fcref       &              ' ADM  precision_derivative_cost =', fcref
483                 CALL PRINT_MESSAGE                 CALL PRINT_MESSAGE
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22