- Timestamp:
- 2013-12-11T18:34:00+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90
r4161 r4333 29 29 30 30 LOGICAL :: linit = .TRUE. ! initialization flag (set to flase after the 1st call) 31 REAL(wp) :: epsi04 = 1 e-04 ! constant31 REAL(wp) :: epsi04 = 1.e-04 ! constant 32 32 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: efact ! metric coefficient 33 33 … … 54 54 REAL(wp), DIMENSION(jpi,jpj), INTENT( inout ) :: ptab ! Field on which the diffusion is applied 55 55 ! 56 INTEGER :: ji, jj ! dummy loop indices57 INTEGER :: its, iter, ierr ! local integers58 REAL(wp) :: zalfa, zrlxint, zconv , zeps! local scalars56 INTEGER :: ji, jj ! dummy loop indices 57 INTEGER :: its, iter, ierr ! local integers 58 REAL(wp) :: zalfa, zrlxint, zconv ! local scalars 59 59 REAL(wp), POINTER, DIMENSION(:,:) :: zrlx, zflu, zflv, zdiv0, zdiv, ztab0 60 60 CHARACTER(lc) :: charout ! local character … … 79 79 zalfa = 0.5_wp ! =1.0/0.5/0.0 = implicit/Cranck-Nicholson/explicit 80 80 its = 100 ! Maximum number of iteration 81 zeps = 2._wp * epsi0482 81 ! 83 82 ztab0(:, : ) = ptab(:,:) ! Arrays initialization … … 94 93 iter = 0 95 94 ! 96 DO WHILE( zconv > zeps.AND. iter <= its ) ! Sub-time step loop95 DO WHILE( zconv > ( 2._wp * epsi04 ) .AND. iter <= its ) ! Sub-time step loop 97 96 ! 98 97 iter = iter + 1 ! incrementation of the sub-time step number
Note: See TracChangeset
for help on using the changeset viewer.