Changeset 2599

Timestamp:

2011-02-20T16:02:36+01:00 (13 years ago)

Author:

gm

Message:

dynamic mem: #785 ; LIM-3 case: DOCTOR norm alb_ice.. changed in zalb_ice...

File:

• branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90 (modified) (10 diffs)

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

@@ -4 +4 @@
    !! Surface module :  update the ocean surface boundary condition over ice
    !!       &           covered area using LIM sea-ice model
-   !! Sea-Ice model  :  LIM 3.0 Sea ice model time-stepping
+   !! Sea-Ice model  :  LIM-3 Sea ice model time-stepping
    !!=====================================================================
    !! History :  2.0  ! 2006-12  (M. Vancoppenolle) Original code
@@ -10 +10 @@
    !!             -   ! 2008-04  (G. Madec)  sltyle and lim_ctl routine
    !!            3.3  ! 2010-11  (G. Madec) ice-ocean stress always computed at each ocean time-step
+   !!            4.0  ! 2011-01  (A Porter)  dynamical allocation
    !!----------------------------------------------------------------------
 #if defined key_lim3
@@ -62 +63 @@
 #  include "vectopt_loop_substitute.h90"
    !!----------------------------------------------------------------------
-   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
+   !! NEMO/OPA 4.0 , UCL NEMO Consortium (2010)
    !! $Id$
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -89 +90 @@
       !!---------------------------------------------------------------------
       USE wrk_nemo, ONLY: wrk_use, wrk_release
-      USE wrk_nemo, ONLY: alb_ice_os => wrk_3d_1 ! albedo of the ice under overcast sky
-      USE wrk_nemo, ONLY: alb_ice_os => wrk_3d_2 ! albedo of ice under clear sky
+      USE wrk_nemo, ONLY: zalb_ice_os => wrk_3d_1 ! albedo of the ice under overcast sky
+      USE wrk_nemo, ONLY: zalb_ice_cs => wrk_3d_2 ! albedo of ice under clear sky
       !!
       INTEGER, INTENT(in) ::   kt      ! ocean time step
@@ -99 +100 @@
       !!----------------------------------------------------------------------
 
-      IF(.NOT. wrk_use(3, 1,2))THEN
-         CALL ctl_stop('sbc_ice_lim: requested workspace arrays are unavailable.')
-         RETURN
-      ELSE IF(jpl > jpk)THEN
-         CALL ctl_stop('sbc_ice_lim: extent of 3rd dimension of workspace arrays needs to exceed jpk.')
-         RETURN
-      END IF
+      IF( .NOT. wrk_use(3, 1,2) ) THEN
+         CALL ctl_stop( 'sbc_ice_lim: requested workspace arrays are unavailable.' )   ;   RETURN
+      ELSEIF( jpl > jpk ) THEN
+         CALL ctl_stop( 'sbc_ice_lim: extent of 3rd dimension of workspace arrays needs to exceed jpk.' )   ;   RETURN
+      ENDIF
 
       IF( kt == nit000 ) THEN
@@ -131 +130 @@
          t_bo(:,:) = tfreez( sss_m ) +  rt0          ! masked sea surface freezing temperature [Kelvin]
          !                                           ! (set to rt0 over land)
-         CALL albedo_ice( t_su, ht_i, ht_s, alb_ice_cs, alb_ice_os )  ! ... ice albedo
+         CALL albedo_ice( t_su, ht_i, ht_s, zalb_ice_cs, zalb_ice_os )  ! ... ice albedo
 
          DO jl = 1, jpl
@@ -146 +145 @@
          SELECT CASE( kblk )
          CASE( 3 )                                       ! CLIO bulk formulation
-            CALL blk_ice_clio( t_su , alb_ice_cs, alb_ice_os,                           &
+            CALL blk_ice_clio( t_su , zalb_ice_cs, zalb_ice_os,                           &
                &                      utau_ice  , vtau_ice  , qns_ice   , qsr_ice   ,   &
                &                      qla_ice   , dqns_ice  , dqla_ice  ,               &
@@ -153 +152 @@
             !         
          CASE( 4 )                                       ! CORE bulk formulation
-            CALL blk_ice_core( t_su , u_ice     , v_ice     , alb_ice_cs,               &
+            CALL blk_ice_core( t_su , u_ice     , v_ice     , zalb_ice_cs,               &
                &                      utau_ice  , vtau_ice  , qns_ice   , qsr_ice   ,   &
                &                      qla_ice   , dqns_ice  , dqla_ice  ,               &
@@ -255 +254 @@
 !!gm   remark, the ocean-ice stress is not saved in ice diag call above .....  find a solution!!!
       !
-      IF(.NOT. wrk_release(3, 1,2))THEN
-         CALL ctl_stop('sbc_ice_lim: failed to release workspace arrays.')
+      IF( .NOT. wrk_release(3, 1,2) ) THEN
+         CALL ctl_stop( 'sbc_ice_lim: failed to release workspace arrays.' )
       END IF
       !
@@ -277 +276 @@
 
       inb_altests = 10
-      inb_alp(:) = 0
+      inb_alp(:)  = 0
 
       ! Alert if incompatible volume and concentration