2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
5 |
|
CBOP |
6 |
|
C !ROUTINE: UPDATE_CG2D |
7 |
|
C !INTERFACE: |
8 |
SUBROUTINE UPDATE_CG2D( myTime, myIter, myThid ) |
SUBROUTINE UPDATE_CG2D( myTime, myIter, myThid ) |
9 |
C /==========================================================\ |
C !DESCRIPTION: \bv |
10 |
C | SUBROUTINE UPDATE_CG2D | |
C *==========================================================* |
11 |
C | o Update 2d conjugate gradient solver operators | |
C | SUBROUTINE UPDATE_CG2D |
12 |
C | account for Free-Surf effect on total column thickness | |
C | o Update 2d conjugate gradient solver operators |
13 |
C |==========================================================| |
C | account for Free-Surf effect on total column thickness |
14 |
C | This routine is based on INI_CG2D, and simplified | |
C *==========================================================* |
15 |
C \==========================================================/ |
C | This routine is based on INI_CG2D, and simplified. It is |
16 |
IMPLICIT NONE |
C | only active when the non-linear free surface mode of |
17 |
|
C | equations is active. |
18 |
|
C *==========================================================* |
19 |
|
C \ev |
20 |
|
|
21 |
|
C !USES: |
22 |
|
IMPLICIT NONE |
23 |
C === Global variables === |
C === Global variables === |
24 |
#include "SIZE.h" |
#include "SIZE.h" |
25 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
31 |
#include "OBCS.h" |
#include "OBCS.h" |
32 |
#endif |
#endif |
33 |
|
|
34 |
|
C !INPUT/OUTPUT PARAMETERS: |
35 |
C === Routine arguments === |
C === Routine arguments === |
36 |
C myTime - Current time in simulation |
C myTime - Current time in simulation |
37 |
C myIter - Current iteration number in simulation |
C myIter - Current iteration number in simulation |
39 |
_RL myTime |
_RL myTime |
40 |
INTEGER myIter |
INTEGER myIter |
41 |
INTEGER myThid |
INTEGER myThid |
42 |
|
|
43 |
|
C !LOCAL VARIABLES: |
44 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
45 |
C-- Note : compared to "INI_CG2D", no needs to compute again |
C-- Note : compared to "INI_CG2D", no needs to compute again |
46 |
C the solver norn=malisation factor of the solver tolerance |
C the solver norn=malisation factor of the solver tolerance |
|
|
|
47 |
C === Local variables === |
C === Local variables === |
48 |
C bi,bj - Loop counters |
C bi,bj - Loop counters |
49 |
C I,J,K |
C I,J,K |
52 |
INTEGER I, J, K |
INTEGER I, J, K |
53 |
_RL faceArea |
_RL faceArea |
54 |
_RL aC, aCw, aCs |
_RL aC, aCw, aCs |
55 |
|
CEOP |
56 |
|
|
57 |
C-- Initialise laplace operator |
C-- Initialise laplace operator |
58 |
C aW2d: integral in Z Ax/dX |
C aW2d: integral in Z Ax/dX |