Changeset 2777 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

Timestamp:

2011-06-07T11:55:02+02:00 (13 years ago)

Author:

smasson

Message:

LIM3 compiling and (partly?) running in v3_3_1, see ticket#817

File:

• trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

@@ -508 +508 @@
       !! ** Inputs / Ouputs : kstrngth (what kind of ice strength we are using)
       !!----------------------------------------------------------------------
-      USE wrk_nemo, ONLY: zworka => wrk_2d_1    ! 2D workspace
+      USE wrk_nemo, ONLY: zworka => wrk_2d_3    ! 2D workspace
       !
       INTEGER, INTENT(in) ::   kstrngth    ! = 1 for Rothrock formulation, 0 for Hibler (1979)
@@ -519 +519 @@
       !!----------------------------------------------------------------------
 
-      IF( wrk_in_use(2, 1) ) THEN
+      IF( wrk_in_use(2, 3) ) THEN
          CALL ctl_stop('lim_itd_me_icestrength : requested workspace array unavailable')   ;   RETURN
       ENDIF
@@ -675 +675 @@
       CALL lbc_lnk( strength, 'T', 1. )      ! Boundary conditions
 
-      IF( wrk_not_released(2, 1) )   CALL ctl_stop('lim_itd_me_icestrength: failed to release workspace array')
+      IF( wrk_not_released(2, 3) )   CALL ctl_stop('lim_itd_me_icestrength: failed to release workspace array')
       !
    END SUBROUTINE lim_itd_me_icestrength
@@ -894 +894 @@
 
       ! Normalization factor : aksum, ensures mass conservation
-      aksum(:,:) = athorn(ji,jj,0)
+      aksum(:,:) = athorn(:,:,0)
       DO jl = 1, jpl 
          aksum(:,:)    = aksum(:,:) + aridge(:,:,jl) * ( 1._wp - 1._wp / krdg(:,:,jl) )    &