Changeset 4332 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90

Timestamp:

2013-12-11T15:38:42+01:00 (10 years ago)

Author:

clem

Message:

update LIM3 to fix remaining bugs. Now working in global and regional config.

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90 (modified) (4 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90

@@ -29 +29 @@
 
    LOGICAL  ::   linit = .TRUE.              ! initialization flag (set to flase after the 1st call)
-   REAL(wp) ::   epsi04 = 1e-04              ! constant
+   REAL(wp) ::   epsi04 = 1.e-04              ! constant
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   efact   ! metric coefficient
 
@@ -54 +54 @@
       REAL(wp), DIMENSION(jpi,jpj), INTENT( inout ) ::   ptab    ! Field on which the diffusion is applied
       !
-      INTEGER  ::  ji, jj            ! dummy loop indices
-      INTEGER  ::  its, iter, ierr   ! local integers
-      REAL(wp) ::   zalfa, zrlxint, zconv, zeps   ! local scalars
+      INTEGER  ::  ji, jj                   ! dummy loop indices
+      INTEGER  ::  its, iter, ierr          ! local integers
+      REAL(wp) ::   zalfa, zrlxint, zconv   ! local scalars
       REAL(wp), POINTER, DIMENSION(:,:) ::   zrlx, zflu, zflv, zdiv0, zdiv, ztab0
       CHARACTER(lc) ::   charout   ! local character
@@ -79 +79 @@
       zalfa = 0.5_wp                      ! =1.0/0.5/0.0 = implicit/Cranck-Nicholson/explicit
       its   = 100                         ! Maximum number of iteration
-      zeps  =  2._wp * epsi04
       !
       ztab0(:, : ) = ptab(:,:)      ! Arrays initialization
@@ -94 +93 @@
       iter  = 0
       !
-      DO WHILE( zconv > zeps .AND. iter <= its )   ! Sub-time step loop
+      DO WHILE( zconv > ( 2._wp * epsi04 ) .AND. iter <= its )   ! Sub-time step loop
          !
          iter = iter + 1                                 ! incrementation of the sub-time step number