/[MITgcm]/MITgcm/pkg/streamice/streamice_vel_solve.F
ViewVC logotype

Diff of /MITgcm/pkg/streamice/streamice_vel_solve.F

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

revision 1.3 by jmc, Thu Jun 13 19:08:43 2013 UTC revision 1.4 by dgoldberg, Sat Aug 24 20:34:27 2013 UTC
# Line 21  C     === Global variables === Line 21  C     === Global variables ===
21  #include "PARAMS.h"  #include "PARAMS.h"
22  #include "STREAMICE.h"  #include "STREAMICE.h"
23  #include "STREAMICE_CG.h"  #include "STREAMICE_CG.h"
24  #ifdef ALLOW_PETSC  !#ifdef ALLOW_PETSC
25  #include "finclude/petsc.h"  !#include "finclude/petsc.h"
26  #endif  !#endif
27    
28  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
29  # include "tamc.h"  # include "tamc.h"
# Line 46  C     LOCAL VARIABLES Line 46  C     LOCAL VARIABLES
46    
47        INTEGER conv_flag, i, j, k, l, iter, cg_iters, bi, bj        INTEGER conv_flag, i, j, k, l, iter, cg_iters, bi, bj
48        INTEGER iter_numconv        INTEGER iter_numconv
49        INTEGER ikey_nl        INTEGER ikey_nl, myThidTemp
50        _RL err_max, err_tempu, err_tempv, err_init, area, err_max_fp        _RL err_max, err_tempu, err_tempv, err_init, area, err_max_fp
51        _RL max_vel, tempu, tempv, err_lastchange, cgtol        _RL max_vel, tempu, tempv, err_lastchange, cgtol
52        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
53  #ifdef ALLOW_PETSC  #ifdef ALLOW_PETSC
54        PetscErrorCode ierr  !      myThidTemp = myThid
55    !      call streamice_initialize_petsc (myThidTemp)
56  #endif  #endif
57  !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
58  !       _RL taudy (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  !       _RL taudy (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 191  C     LOCAL VARIABLES Line 192  C     LOCAL VARIABLES
192    
193  C START NL ITER. LOOP  C START NL ITER. LOOP
194  C-------------------------------------------------------------------  C-------------------------------------------------------------------
 #ifdef ALLOW_PETSC  
       call petscInitialize(PETSC_NULL_CHARACTER,ierr)  
 #endif  
195    
196        DO iter=1,streamice_max_nl_iter        DO iter=1,streamice_max_nl_iter
197    
# Line 453  C       but after convergence loop does Line 451  C       but after convergence loop does
451        ENDDO        ENDDO
452    
453  #ifdef ALLOW_PETSC  #ifdef ALLOW_PETSC
454        call PetscFinalize(ierr)  !      call streamice_finalize_petsc (myThidTemp)
455    !      call streamice_finalize_petsc (myThid)
456  #endif  #endif
457    
458  C END NL ITER. LOOP  C END NL ITER. LOOP
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22