Changeset 2715 for trunk/NEMOGCM/NEMO/LIM_SRC_2/iceini_2.F90

Timestamp:

2011-03-30T17:58:35+02:00 (13 years ago)

Author:

rblod

Message:

First attempt to put dynamic allocation on the trunk

File:

• trunk/NEMOGCM/NEMO/LIM_SRC_2/iceini_2.F90 (modified) (5 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_2/iceini_2.F90

@@ -4 +4 @@
    !!   Sea-ice model : LIM 2.0 Sea ice model Initialization
    !!======================================================================
-   !! History :   1.0  !  02-08  (G. Madec)  F90: Free form and modules
-   !!             2.0  !  03-08  (C. Ethe)  add ice_run
-   !!             3.3  !  09-05  (G.Garric, C. Bricaud) addition of the lim2_evp case
+   !! History :  1.0  ! 2002-08  (G. Madec)  F90: Free form and modules
+   !!            2.0  ! 2003-08  (C. Ethe)  add ice_run
+   !!            3.3  ! 2009-05  (G. Garric, C. Bricaud) addition of the lim2_evp case
+   !!            4.0  ! 2011-02  (G. Madec) dynamical allocation
    !!----------------------------------------------------------------------
 #if defined key_lim2
@@ -12 +13 @@
    !!   'key_lim2' :                                  LIM 2.0 sea-ice model
    !!----------------------------------------------------------------------
-   !!----------------------------------------------------------------------
    !!   ice_init_2       : sea-ice model initialization
    !!   ice_run_2        : Definition some run parameter for ice model
    !!----------------------------------------------------------------------
+   USE phycst           ! physical constants
    USE dom_oce          ! ocean domain
-   USE dom_ice_2        ! LIM2: ice domain
    USE sbc_oce          ! surface boundary condition: ocean
-   USE sbc_ice          ! surface boundary condition: ice
-   USE phycst           ! Define parameters for the routines
-   USE ice_2            ! LIM2: ice variable
-   USE limmsh_2         ! LIM2: mesh
-   USE limistate_2      ! LIM2: initial state
-   USE limrst_2         ! LIM2: restart
+   USE sbc_ice          ! LIM2 surface boundary condition
+   USE dom_ice_2        ! LIM2 ice domain
+   USE par_ice_2        ! LIM2 parameters
+   USE thd_ice_2        ! LIM2 thermodynamical variables
+   USE limrhg           ! LIM2 rheology
+   USE ice_2            ! LIM2 ice variable
+   USE limmsh_2         ! LIM2 mesh
+   USE limistate_2      ! LIM2 initial state
+   USE limrst_2         ! LIM2 restart
+   USE limsbc_2         ! LIM2 surface boundary condition
    USE in_out_manager   ! I/O manager
-      
+   USE lib_mpp          ! MPP library
+
    IMPLICIT NONE
    PRIVATE
@@ -33 +38 @@
 
    !!----------------------------------------------------------------------
-   !! NEMO/LIM2 3.3 , UCL - NEMO Consortium (2010)
+   !! NEMO/LIM2 4.0 , UCL - NEMO Consortium (2011)
    !! $Id$ 
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -45 +50 @@
       !! ** purpose :   initialisation of LIM-2 domain and variables  
       !!----------------------------------------------------------------------
+      INTEGER :: ierr
+      !!----------------------------------------------------------------------
       !
+      IF(lwp) THEN
+         WRITE(numout,*)
+         WRITE(numout,*) 'ice_init_2 : LIM-2 sea-ice - initialization'
+         WRITE(numout,*) '~~~~~~~~~~~   '
+      ENDIF
+      !                                ! Allocate the ice arrays
+      ierr =        ice_alloc_2    ()       ! ice variables
+      ierr = ierr + dom_ice_alloc_2()       ! domain
+      ierr = ierr + sbc_ice_alloc  ()       ! surface forcing
+      ierr = ierr + thd_ice_alloc_2()       ! thermodynamics
+#if ! defined key_lim2_vp
+      ierr = ierr + lim_rhg_alloc  ()
+#endif
+      IF( lk_mpp    )   CALL mpp_sum( ierr )
+      IF( ierr /= 0 )   CALL ctl_stop( 'STOP', 'ice_init_2 : unable to allocate ice arrays' )
+
+      !                                ! adequation jpk versus ice/snow layers
+      IF( jpl > jpk  .OR.  jplayersp1 > jpk  )   CALL ctl_stop( 'STOP',           &
+         &     'ice_init: the 3rd dimension of workspace arrays is too small.',   &
+         &     'use more ocean levels or less ice layers/categories.' )
+
       !                                ! Open the namelist file 
       CALL ctl_opn( numnam_ice, 'namelist_ice', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 
@@ -61 +89 @@
       ENDIF
       !
-      tn_ice(:,:,1) = sist(:,:)        ! initialisation of ice temperature   
-      fr_i  (:,:) = 1.0 - frld(:,:)    ! initialisation of sea-ice fraction    
+      tn_ice(:,:,1) = sist(:,:)        ! ice temperature  known by the ocean
+      fr_i  (:,:)   = 1.0 - frld(:,:)  ! sea-ice fraction known by the ocean
+      !
+      CALL lim_sbc_init_2              ! ice surface boundary condition   
+      !
+      IF( lk_lim2_vp )   THEN   ;   WRITE(numout,*) '                VP  rheology - B-grid case'
+      ELSE                      ;   WRITE(numout,*) '                EVP rheology - C-grid case'
+      ENDIF
       !
    END SUBROUTINE ice_init_2