Changeset 13149 for NEMO/branches/UKMO/NEMO_4.0_mirror_SI3_GPU/src/OCE/DYN/dynldf_lap_blp.F90

Timestamp:

2020-06-24T09:03:45+02:00 (4 years ago)

Author:

andmirek

Message:

Ticket #2482: Dissable restart and use allocatable arrays

File:

• NEMO/branches/UKMO/NEMO_4.0_mirror_SI3_GPU/src/OCE/DYN/dynldf_lap_blp.F90 (modified) (3 diffs)

NEMO/branches/UKMO/NEMO_4.0_mirror_SI3_GPU/src/OCE/DYN/dynldf_lap_blp.F90

@@ -56 +56 @@
       REAL(wp) ::   zua, zva     ! local scalars
       REAL(wp), DIMENSION(jpi,jpj) ::   zcur, zdiv
+!     REAL(wp), ALLOCATABLE, DIMENSION(:, :) :: zcur, zdiv
       !!----------------------------------------------------------------------
       !
-      IF( kt == nit000 .AND. lwp ) THEN
-         WRITE(numout,*)
-         WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacian) operator, pass=', kpass
-         WRITE(numout,*) '~~~~~~~ '
+      IF( kt == nit000 ) THEN
+         IF(lwp) THEN
+            WRITE(numout,*)
+            WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacian) operator, pass=', kpass
+            WRITE(numout,*) '~~~~~~~ '
+         ENDIF
+!        ALLOCATE(zcur(jpi,jpj), zdiv(jpi,jpj))
       ENDIF
       !
@@ -122 +126 @@
       REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) ::   pua, pva   ! momentum trend
       !
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zulap, zvlap   ! laplacian at u- and v-point
+!     REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zulap, zvlap   ! laplacian at u- and v-point
+      REAL(wp), ALLOCATABLE, DIMENSION(:, :, :) :: zulap, zvlap
       !!----------------------------------------------------------------------
       !
@@ -129 +134 @@
          IF(lwp) WRITE(numout,*) 'dyn_ldf_blp : bilaplacian operator momentum '
          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~'
+         ALLOCATE(zulap(jpi,jpj,jpk), zvlap(jpi,jpj,jpk))
       ENDIF
       !