/[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.33 by jmc, Mon Sep 26 12:56:52 2011 UTC revision 1.34 by jmc, Wed Feb 8 03:14:13 2012 UTC
# Line 229  cph( Line 229  cph(
229  cph-print               print *, 'ph-grd ----- back from loc -----',  cph-print               print *, 'ph-grd ----- back from loc -----',
230  cph-print     &             icvrec, itilepos, jtilepos, layer, obcspos  cph-print     &             icvrec, itilepos, jtilepos, layer, obcspos
231  cph)  cph)
232                else
233                  itile    = 0
234                  jtile    = 0
235                  layer    = 0
236                  itilepos = 0
237                  jtilepos = 0
238              endif              endif
239              _BARRIER              _BARRIER
240    
# Line 263  c--   1. perturb control vector componen Line 269  c--   1. perturb control vector componen
269       &              itilepos, jtilepos,       &              itilepos, jtilepos,
270       &              xxmemo_ref, xxmemo_pert, localEps,       &              xxmemo_ref, xxmemo_pert, localEps,
271       &              mythid )       &              mythid )
272                else
273                  xxmemo_ref  = 0.
274                  xxmemo_pert = 0.
275              endif              endif
276              _BARRIER              _BARRIER
277    
# Line 319  c--   positive perturbation Line 328  c--   positive perturbation
328       &              itilepos, jtilepos,       &              itilepos, jtilepos,
329       &              xxmemo_ref, xxmemo_pert, localEps,       &              xxmemo_ref, xxmemo_pert, localEps,
330       &              mythid )       &              mythid )
331                else
332                  xxmemo_ref  = 0.
333                  xxmemo_pert = 0.
334              endif              endif
335              _BARRIER              _BARRIER
336    
# Line 360  c--   repeat the proceedure for a negati Line 372  c--   repeat the proceedure for a negati
372       &                 itilepos, jtilepos,       &                 itilepos, jtilepos,
373       &                 xxmemo_ref, xxmemo_pert, localEps,       &                 xxmemo_ref, xxmemo_pert, localEps,
374       &                 mythid )       &                 mythid )
375                   else
376                     xxmemo_ref  = 0.
377                     xxmemo_pert = 0.
378                 endif                 endif
379                 _BARRIER                 _BARRIER
380    
# Line 412  c*************************************** Line 427  c***************************************
427    
428              if ( myProcId .EQ. grdchkwhichproc .AND.              if ( myProcId .EQ. grdchkwhichproc .AND.
429       &           ierr .EQ. 0 ) then       &           ierr .EQ. 0 ) then
   
430                 tmpplot1(itilepos,jtilepos,layer,itile,jtile)                 tmpplot1(itilepos,jtilepos,layer,itile,jtile)
431       &              = gfd       &              = gfd
432                 tmpplot2(itilepos,jtilepos,layer,itile,jtile)                 tmpplot2(itilepos,jtilepos,layer,itile,jtile)
433       &              = ratio_ad       &              = ratio_ad
434                 tmpplot3(itilepos,jtilepos,layer,itile,jtile)                 tmpplot3(itilepos,jtilepos,layer,itile,jtile)
435       &              = ratio_ftl       &              = ratio_ftl
436                endif
437    
438                if ( ierr .EQ. 0 ) then
439                 fcrmem      ( ichknum ) = fcref                 fcrmem      ( ichknum ) = fcref
440                 fcppmem     ( ichknum ) = fcpertplus                 fcppmem     ( ichknum ) = fcpertplus
441                 fcpmmem     ( ichknum ) = fcpertminus                 fcpmmem     ( ichknum ) = fcpertminus
# Line 443  c*************************************** Line 459  c***************************************
459                 itestmem  ( ichknum ) = itest                 itestmem  ( ichknum ) = itest
460                 ierrmem   ( ichknum ) = ierr                 ierrmem   ( ichknum ) = ierr
461                 icglomem  ( ichknum ) = icglo                 icglomem  ( ichknum ) = icglo
462                endif
463    
464                if ( myProcId .EQ. grdchkwhichproc .AND.
465         &           ierr .EQ. 0 ) then
466    
467                 WRITE(standardmessageunit,'(A)')                 WRITE(standardmessageunit,'(A)')
468       &             'grad-res -------------------------------'       &             'grad-res -------------------------------'
Legend:
Removed from v.1.33  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22