/[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.14 by heimbach, Tue Mar 23 19:42:53 2004 UTC revision 1.15 by heimbach, Tue Oct 26 20:10:25 2004 UTC
# Line 155  cph   assuming all xx_ fields are initia Line 155  cph   assuming all xx_ fields are initia
155        ierr_grdchk = 0        ierr_grdchk = 0
156  cphadmtlm(  cphadmtlm(
157        fcref = fc        fcref = fc
158  cphadmtlm      fcref = objf_state_final(45,4,1,1)  cphadmtlm      fcref = objf_state_final(45,4,1,1,1)
159  cphadmtlm)  cphadmtlm)
160    
161        print *, 'ph-check fcref = ', fcref        print *, 'ph-check fcref = ', fcref
# Line 255  cphadmtlm( Line 255  cphadmtlm(
255              g_fc = 0.              g_fc = 0.
256  cphadmtlm            do j=1,sny  cphadmtlm            do j=1,sny
257  cphadmtlm               do i=1,snx  cphadmtlm               do i=1,snx
258  cphadmtlm                  g_objf_state_final(i,j,1,1) = 0.  cphadmtlm                  g_objf_state_final(i,j,1,1,1) = 0.
259    cphadmtlm                  g_objf_state_final(i,j,1,1,2) = 0.
260  cphadmtlm               enddo  cphadmtlm               enddo
261  cphadmtlm            enddo  cphadmtlm            enddo
262  cphadmtlm)  cphadmtlm)
263              call g_the_main_loop( mytime, myiter, mythid )              call g_the_main_loop( mytime, myiter, mythid )
264  cphadmtlm(  cphadmtlm(
265              ftlxxmemo = g_fc              ftlxxmemo = g_fc
266  cphadmtlm            ftlxxmemo = g_objf_state_final(45,4,1,1)  cphadmtlm            ftlxxmemo = g_objf_state_final(45,4,1,1,1)
267  cphadmtlm)  cphadmtlm)
268              _BARRIER              _BARRIER
269  c--  c--
# Line 303  c--   forward run with perturbed control Line 304  c--   forward run with perturbed control
304              call the_main_loop( mytime, myiter, mythid )              call the_main_loop( mytime, myiter, mythid )
305  cphadmtlm(  cphadmtlm(
306              fcpertplus = fc              fcpertplus = fc
307  cphadmtlm            fcpertplus = objf_state_final(45,4,1,1)  cphadmtlm            fcpertplus = objf_state_final(45,4,1,1,1)
308  cphadmtlm)  cphadmtlm)
309              print *, 'ph-check fcpertplus = ', fcpertplus              print *, 'ph-check fcpertplus = ', fcpertplus
310              _BARRIER              _BARRIER
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22