70 |
CALL DUMMY_IN_STEPPING( myTime, myIter, myThid ) |
CALL DUMMY_IN_STEPPING( myTime, myIter, myThid ) |
71 |
#endif |
#endif |
72 |
|
|
73 |
|
#ifdef NONLIN_FRSURF |
74 |
|
C-- compute the future surface level thickness |
75 |
|
C according to the current Eta field |
76 |
|
IF ( nonlinFreeSurf.GT.0) THEN |
77 |
|
CALL CALC_SURF_DR(etaN, myTime, myIter, myThid ) |
78 |
|
ENDIF |
79 |
|
#endif |
80 |
|
|
81 |
C-- Load forcing/external data fields. |
C-- Load forcing/external data fields. |
82 |
#ifdef INCLUDE_EXTERNAL_FORCING_PACKAGE |
#ifdef INCLUDE_EXTERNAL_FORCING_PACKAGE |
83 |
C NOTE, that although the exf package is part of the |
C NOTE, that although the exf package is part of the |
111 |
#endif |
#endif |
112 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
113 |
|
|
114 |
|
#ifdef NONLIN_FRSURF |
115 |
|
C-- update hfacC,W,S and recip_hFac according to the current Eta field |
116 |
|
IF ( nonlinFreeSurf.GT.0) THEN |
117 |
|
CALL UPDATE_SURF_DR( myTime, myIter, myThid ) |
118 |
|
ENDIF |
119 |
|
C- update also CG2D matrix (and preconditioner) |
120 |
|
IF ( nonlinFreeSurf.GT.2) THEN |
121 |
|
CALL UPDATE_CG2D( myTime, myIter, myThid ) |
122 |
|
ENDIF |
123 |
|
#endif |
124 |
|
|
125 |
C-- This block has been moved to the_correction_step(). We have |
C-- This block has been moved to the_correction_step(). We have |
126 |
C left this code here to indicate where the filters used to be |
C left this code here to indicate where the filters used to be |
127 |
C in the algorithm before JMC moved them to after the pressure |
C in the algorithm before JMC moved them to after the pressure |