Changeset 5180 for trunk/NEMOGCM

Timestamp:

2015-03-27T11:16:58+01:00 (9 years ago)

Author:

vancop

Message:

Remove out-of-bound exceptions

Location:

trunk/NEMOGCM/NEMO/LIM_SRC_3

Files:

• limthd_dif.F90 (modified) (1 diff)
• limtrp.F90 (modified) (2 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

@@ -679 +679 @@
          END DO
 
-         DO numeq = nlay_i + nlay_s + 1, nlay_s + 2, -1
+         DO numeq = nlay_i + nlay_s, nlay_s + 2, -1
             DO ji = kideb , kiut
                jk    =  numeq - nlay_s - 1

trunk/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90

@@ -80 +80 @@
       IF( nn_timing == 1 )  CALL timing_start('limtrp')
 
-      CALL wrk_alloc( jpi,jpj,           zsm, zatold, zeiold, zesold )
-      CALL wrk_alloc( jpi,jpj,jpl,       z0ice, z0snw, z0ai, z0es , z0smi , z0oi )
-      CALL wrk_alloc( jpi,jpj,1,         z0opw )
-      CALL wrk_alloc( jpi,jpj,nlay_i+1,jpl, z0ei )
-      CALL wrk_alloc( jpi,jpj,jpl,       zhimax, zviold, zvsold, zsmvold )
+      CALL wrk_alloc( jpi,jpj,            zsm, zatold, zeiold, zesold )
+      CALL wrk_alloc( jpi,jpj,jpl,        z0ice, z0snw, z0ai, z0es , z0smi , z0oi )
+      CALL wrk_alloc( jpi,jpj,1,          z0opw )
+      CALL wrk_alloc( jpi,jpj,nlay_i,jpl, z0ei )
+      CALL wrk_alloc( jpi,jpj,jpl,        zhimax, zviold, zvsold, zsmvold )
 
       IF( numit == nstart .AND. lwp ) THEN
@@ -459 +459 @@
       IF( ln_icectl )   CALL lim_prt( kt, iiceprt, jiceprt,-1, ' - ice dyn & trp - ' )
       !
-      CALL wrk_dealloc( jpi,jpj,           zsm, zatold, zeiold, zesold )
-      CALL wrk_dealloc( jpi,jpj,jpl,       z0ice, z0snw, z0ai, z0es , z0smi , z0oi )
-      CALL wrk_dealloc( jpi,jpj,1,         z0opw )
-      CALL wrk_dealloc( jpi,jpj,nlay_i+1,jpl, z0ei )
-      CALL wrk_dealloc( jpi,jpj,jpl,       zviold, zvsold, zhimax, zsmvold )
+      CALL wrk_dealloc( jpi,jpj,            zsm, zatold, zeiold, zesold )
+      CALL wrk_dealloc( jpi,jpj,jpl,        z0ice, z0snw, z0ai, z0es , z0smi , z0oi )
+      CALL wrk_dealloc( jpi,jpj,1,          z0opw )
+      CALL wrk_dealloc( jpi,jpj,nlay_i,jpl, z0ei )
+      CALL wrk_dealloc( jpi,jpj,jpl,        zviold, zvsold, zhimax, zsmvold )
       !
       IF( nn_timing == 1 )  CALL timing_stop('limtrp')