Changeset 8866 for branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/NST_SRC/agrif_opa_update.F90

Timestamp:

2017-12-01T07:22:09+01:00 (6 years ago)

Author:

gm

Message:

#1911 (ENHANCE-09): Agrif: remove Agrif_update_tke + style cleaning

File:

• branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/NST_SRC/agrif_opa_update.F90 (modified) (5 diffs)

branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/NST_SRC/agrif_opa_update.F90

@@ -7 +7 @@
    !! AGRIF: update package for the ocean dynamics (OPA)
    !!======================================================================
-   !! History :  2.0  !  2002-06  (XXX)  Original cade
-   !!             -   !  2005-11  (XXX) 
+   !! History :  2.0  !  2002-06  (L. Debreu)  Original code
    !!            3.2  !  2009-04  (R. Benshila) 
    !!            3.6  !  2014-09  (R. Benshila) 
@@ -29 +28 @@
 
    PUBLIC   Agrif_Update_Tra, Agrif_Update_Dyn, Update_Scales
-   PUBLIC   Agrif_Update_Tke
 
    !!----------------------------------------------------------------------
@@ -160 +158 @@
    END SUBROUTINE Agrif_Update_Dyn
 
-!!gm Missing GLS case !!!!!
-
-   SUBROUTINE Agrif_Update_Tke( kt )
-      !!----------------------------------------------------------------------
-      !!                   *** ROUTINE Agrif_Update_Tke ***
-      !!----------------------------------------------------------------------
-      INTEGER, INTENT(in) ::   kt   ! timestep index
-      !!----------------------------------------------------------------------
-      !
-!!gm test on kt/=0  ????  why not nit000-1  ?  doesn't seem logic
-      IF( ( Agrif_NbStepint() /= 0 ) .AND. kt /= 0 )   RETURN
-#  if defined TWO_WAY
-      !
-      Agrif_UseSpecialValueInUpdate = .TRUE.
-      Agrif_SpecialValueFineGrid    = 0._wp
-      !
-      CALL Agrif_Update_Variable(  en_id, locupdate=(/0,0/), procname=updateEN  )
-      CALL Agrif_Update_Variable( avt_id, locupdate=(/0,0/), procname=updateAVT )
-      CALL Agrif_Update_Variable( avm_id, locupdate=(/0,0/), procname=updateAVM )
-      !
-      Agrif_UseSpecialValueInUpdate = .FALSE.
-      !
-#  endif
-      !
-   END SUBROUTINE Agrif_Update_Tke
-   
 
    SUBROUTINE updateTS( tabres, i1, i2, j1, j2, k1, k2, n1, n2, before )
@@ -215 +187 @@
                   DO jj = j1, j2
                      DO ji = i1, i2
-                        IF( tabres(ji,jj,jk,jn) .NE. 0. ) THEN
+                        IF( tabres(ji,jj,jk,jn) /= 0._wp ) THEN
                            tsb(ji,jj,jk,jn) = tsb(ji,jj,jk,jn) & 
                               &             + atfp * ( tabres(ji,jj,jk,jn) - tsn(ji,jj,jk,jn) ) * tmask(ji,jj,jk)
@@ -228 +200 @@
                DO jj=j1,j2
                   DO ji=i1,i2
-                     IF( tabres(ji,jj,jk,jn) .NE. 0. ) THEN 
+                     IF( tabres(ji,jj,jk,jn) /= 0._wp ) THEN 
                         tsn(ji,jj,jk,jn) = tabres(ji,jj,jk,jn) * tmask(ji,jj,jk)
                      END IF