Changeset 2715 for trunk/NEMOGCM/NEMO/LIM_SRC_3/iceini.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_3/iceini.F90 (modified) (5 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_3/iceini.F90

@@ -6 +6 @@
    !! History :  3.0  ! 2008-03  (M. Vancoppenolle) LIM-3 original code
    !!            3.3  ! 2010-12  (G. Madec) add call to lim_thd_init and lim_thd_sal_init
+   !!            4.0  ! 2011-02  (G. Madec) dynamical allocation
    !!----------------------------------------------------------------------
 #if defined key_lim3
@@ -13 +14 @@
    !!   ice_init       : sea-ice model initialization
    !!----------------------------------------------------------------------
-   USE phycst         ! physical constants
-   USE dom_oce        ! ocean domain
-   USE sbc_oce        ! Surface boundary condition: ocean fields
-   USE sbc_ice        ! Surface boundary condition: ice fields
-   USE ice            ! LIM: sea-ice variables
-   USE limmsh         ! LIM: mesh
-   USE limistate      ! LIM: initial state
-   USE limrst         ! LIM: restart
-   USE limthd         ! LIM: ice thermodynamics
-   USE limthd_sal     ! LIM: ice thermodynamics: salinity
-   USE par_ice        ! LIM: sea-ice parameters
-   USE limvar         ! LIM: variables
-   USE in_out_manager ! I/O manager
-   USE lib_mpp        ! MPP library
+   USE phycst           ! physical constants
+   USE dom_oce          ! ocean domain
+   USE sbc_oce          ! Surface boundary condition: ocean fields
+   USE sbc_ice          ! Surface boundary condition: ice   fields
+   USE ice              ! LIM variables
+   USE par_ice          ! LIM parameters
+   USE dom_ice          ! LIM domain
+   USE thd_ice          ! LIM thermodynamical variables
+   USE limitd_me        ! LIM ice thickness distribution
+   USE limrhg           ! LIM dynamics
+   USE limmsh           ! LIM mesh
+   USE limistate        ! LIM initial state
+   USE limrst           ! LIM restart
+   USE limthd           ! LIM ice thermodynamics
+   USE limthd_sal       ! LIM ice thermodynamics: salinity
+   USE limvar           ! LIM variables
+   USE limsbc           ! LIM surface boundary condition
+   USE in_out_manager   ! I/O manager
+   USE lib_mpp          ! MPP library
 
    IMPLICIT NONE
    PRIVATE
 
-   PUBLIC   ice_init   ! called by opa.F90
-
-   !!----------------------------------------------------------------------
-   !! NEMO/LIM3 3.3 , UCL - NEMO Consortium (2010)
+   PUBLIC   ice_init   ! called by sbcice_lim.F90
+
+   !!----------------------------------------------------------------------
+   !! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2011)
    !! $Id$
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -44 +50 @@
       !!                  ***  ROUTINE ice_init  ***
       !!
-      !! ** purpose :   
+      !! ** purpose :   Allocate all the dynamic arrays of the LIM-3 modules
       !!----------------------------------------------------------------------
-      !
+      INTEGER :: ierr
+      !!----------------------------------------------------------------------
+
+      !                                ! Allocate the ice arrays
+      ierr =        ice_alloc       ()       ! ice variables
+      ierr = ierr + dom_ice_alloc   ()       ! domain
+      ierr = ierr + sbc_ice_alloc   ()       ! surface forcing
+      ierr = ierr + thd_ice_alloc   ()       ! thermodynamics
+      ierr = ierr + lim_itd_me_alloc()       ! ice thickness distribution - mechanics
+      ierr = ierr + lim_rhg_alloc   ()       ! dynamics
+      !
+      IF( lk_mpp    )   CALL mpp_sum( ierr )
+      IF( ierr /= 0 )   CALL ctl_stop('STOP', 'ice_init : unable to allocate ice arrays')
+      !
+      !                                ! adequation jpk versus ice/snow layers/categories
+      IF( jpl   > jpk  .OR.  jpm    > jpk .OR.                                    &
+          jkmax > jpk  .OR.  nlay_s > jpk      )   CALL ctl_stop( 'STOP',         &
+         &     'ice_init: the 3rd dimension of workspace arrays is too small.',   &
+         &     'use more ocean levels or less ice/snow layers/categories.' )
+
       !                                ! Open the namelist file 
       CALL ctl_opn( numnam_ice, 'namelist_ice', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
       !
-      CALL ice_run                     ! namelist read some ice run parameters
-      !
-      CALL lim_thd_init                ! namelist read ice thermodynics parameters
-      !
-      CALL lim_thd_sal_init            ! namelist read ice salinity parameters
+      CALL ice_run                     ! set some ice run parameters
+      !
+      CALL lim_thd_init                ! set ice thermodynics parameters
+      !
+      CALL lim_thd_sal_init            ! set ice salinity parameters
       !
       rdt_ice = nn_fsbc * rdttra(1)    ! sea-ice timestep
@@ -60 +85 @@
       CALL lim_msh                     ! ice mesh initialization
       !
-      CALL lim_itd_ini                 ! initialize the ice thickness distribution
+      CALL lim_itd_ini                 ! ice thickness distribution initialization
+      !
+      CALL lim_sbc_init                ! ice surface boundary condition   
+
 
       !                                ! Initial sea-ice state
@@ -94 +122 @@
       !!
       !! ** Method  :   Read the namicerun namelist and check the parameter 
-      !!       values called at the first timestep (nit000)
+      !!              values called at the first timestep (nit000)
       !!
       !! ** input   :   Namelist namicerun