Changeset 11467 for NEMO/branches/UKMO/dev_r10037_GPU/src/ICE/icevar.F90

Timestamp:

2019-08-22T11:49:08+02:00 (5 years ago)

Author:

andmirek

Message:

Ticket #2197 allocate arrays at the beggining of the run

File:

• NEMO/branches/UKMO/dev_r10037_GPU/src/ICE/icevar.F90 (modified) (12 diffs)

NEMO/branches/UKMO/dev_r10037_GPU/src/ICE/icevar.F90

@@ -89 +89 @@
       !!              all-ice variables, i.e. it turns VGLO into VAGG
       !!-------------------------------------------------------------------
+      USE scoce, ONLY : z1_at_i => scr2D4, z1_vt_i => scr2D2, z1_vt_s => scr2D3
       INTEGER, INTENT( in ) ::   kn     ! =1 state variables only
       !                                 ! >1 state variables + others
       !
       INTEGER ::   ji, jj, jk, jl   ! dummy loop indices
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   z1_at_i, z1_vt_i, z1_vt_s
       !!-------------------------------------------------------------------
       !
@@ -112 +112 @@
       IF( kn > 1 ) THEN
          !
-         ALLOCATE( z1_at_i(jpi,jpj) , z1_vt_i(jpi,jpj) , z1_vt_s(jpi,jpj) )
          WHERE( at_i(:,:) > epsi20 )   ;   z1_at_i(:,:) = 1._wp / at_i(:,:)
          ELSEWHERE                     ;   z1_at_i(:,:) = 0._wp
@@ -152 +151 @@
          END WHERE
 
-         DEALLOCATE( z1_at_i , z1_vt_i , z1_vt_s )
       ENDIF
       !
@@ -165 +163 @@
       !!              global variables, i.e. it turns VGLO into VEQV
       !!-------------------------------------------------------------------
+      USE scice, ONLY : z1_a_i => scr3i, z1_v_i => scr4i
       INTEGER  ::   ji, jj, jk, jl   ! dummy loop indices
       REAL(wp) ::   ze_i             ! local scalars
@@ -170 +169 @@
       REAL(wp) ::   zhmax, z1_zhmax                 !   -      -
       REAL(wp) ::   zlay_i, zlay_s                  !   -      -
-      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   z1_a_i, z1_v_i
       !!-------------------------------------------------------------------
 
@@ -294 +292 @@
       !! ** References : Vancoppenolle et al., 2007
       !!-------------------------------------------------------------------
+      USE scice , ONLY : z_slope_s => scr5i, zalpha =>scr6i
       INTEGER  ::   ji, jj, jk, jl   ! dummy loop index
       REAL(wp) ::   zsal, z1_dS
       REAL(wp) ::   zargtemp , zs0, zs
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   z_slope_s, zalpha    ! case 2 only
       REAL(wp), PARAMETER :: zsi0 = 3.5_wp
       REAL(wp), PARAMETER :: zsi1 = 4.5_wp
@@ -315 +313 @@
       CASE( 2 )       !  time varying salinity with linear profile  !
          !            !---------------------------------------------!
-         !
-         ALLOCATE( z_slope_s(jpi,jpj,jpl) , zalpha(jpi,jpj,jpl) )
          !
          DO jl = 1, jpl
@@ -353 +349 @@
          END DO
          !
-         DEALLOCATE( z_slope_s , zalpha )
-         !
          !            !-------------------------------------------!
       CASE( 3 )       ! constant salinity with a fix profile      ! (Schwarzacher (1959) multiyear salinity profile
@@ -386 +380 @@
       !!                Works with 1d vectors and is used by thermodynamic modules
       !!-------------------------------------------------------------------
+      USE scice, ONLY : z_slope_s => scr1D1i, zalpha => scr1D2i
       INTEGER  ::   ji, jk    ! dummy loop indices
       REAL(wp) ::   zargtemp, zsal, z1_dS   ! local scalars
       REAL(wp) ::   zs, zs0              !   -      -
       !
-      REAL(wp), ALLOCATABLE, DIMENSION(:) ::   z_slope_s, zalpha   !
       REAL(wp), PARAMETER :: zsi0 = 3.5_wp
       REAL(wp), PARAMETER :: zsi1 = 4.5_wp
@@ -405 +399 @@
       CASE( 2 )       !  time varying salinity with linear profile  !
          !            !---------------------------------------------!
-         !
-         ALLOCATE( z_slope_s(jpij), zalpha(jpij) )
          !
          !                                      ! Slope of the linear profile 
@@ -430 +422 @@
          END DO
          !
-         DEALLOCATE( z_slope_s, zalpha )
-
          !            !-------------------------------------------!
       CASE( 3 )       ! constant salinity with a fix profile      ! (Schwarzacher (1959) multiyear salinity profile
@@ -458 +448 @@
       !! ** Purpose :   Remove too small sea ice areas and correct fluxes
       !!-------------------------------------------------------------------
+      USE scoce, ONLY : zswitch => scr2D2
       INTEGER  ::   ji, jj, jl, jk   ! dummy loop indices
-      REAL(wp), DIMENSION(jpi,jpj) ::   zswitch
       !!-------------------------------------------------------------------
       !