/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_vel_solve.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_vel_solve.F

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

revision 1.12 by dgoldberg, Thu Mar 7 15:23:19 2013 UTC revision 1.13 by dgoldberg, Sat Apr 6 17:43:41 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
25    #include "finclude/petsc.h"
26    #endif
27    
28  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
29  # include "tamc.h"  # include "tamc.h"
30  #endif  #endif
# Line 46  C     LOCAL VARIABLES Line 50  C     LOCAL VARIABLES
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
54          PetscErrorCode ierr
55    #endif
56  !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+Oly,nSx,nSy)  !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+Oly,nSx,nSy)
57  !       _RL taudy (1-OLx:sNx+OLx,1-OLy:sNy+Oly,nSx,nSy)  !       _RL taudy (1-OLx:sNx+OLx,1-OLy:sNy+Oly,nSx,nSy)
58    
# Line 184  C     LOCAL VARIABLES Line 191  C     LOCAL VARIABLES
191    
192  C START NL ITER. LOOP  C START NL ITER. LOOP
193  C-------------------------------------------------------------------  C-------------------------------------------------------------------
194    #ifdef ALLOW_PETSC
195          call petscInitialize(PETSC_NULL_CHARACTER,ierr)
196    #endif
197    
198        DO iter=1,streamice_max_nl_iter        DO iter=1,streamice_max_nl_iter
199    
# Line 442  C       but after convergence loop does Line 452  C       but after convergence loop does
452         ENDIF ! (err_max .GT. streamice_nonlin_tol * err_init)         ENDIF ! (err_max .GT. streamice_nonlin_tol * err_init)
453        ENDDO        ENDDO
454    
455    #ifdef ALLOW_PETSC
456          call PetscFinalize(ierr)
457    #endif
458    
459  C END NL ITER. LOOP  C END NL ITER. LOOP
460  C-------------------------------------------------------------------  C-------------------------------------------------------------------
461    
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22