Changeset 508 for trunk/NEMO/LIM_SRC/iceini.F90

Timestamp:

2006-10-03T17:58:55+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_071:RB: add iom for restart and reorganization of restart

File:

• trunk/NEMO/LIM_SRC/iceini.F90 (modified) (7 diffs)

trunk/NEMO/LIM_SRC/iceini.F90

@@ -4 +4 @@
    !!   Sea-ice model : LIM 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
+   !!----------------------------------------------------------------------
 #if defined key_ice_lim
    !!----------------------------------------------------------------------
    !!   'key_ice_lim' :                                   LIM sea-ice model
    !!----------------------------------------------------------------------
+   !!----------------------------------------------------------------------
    !!   ice_init       : sea-ice model initialization
+   !!   ice_run        : Definition some run parameter for ice model
    !!----------------------------------------------------------------------
    USE dom_oce
@@ -19 +24 @@
    USE limmsh
    USE limistate
-   USE limrst
+   USE limrst   
    USE ini1d           ! initialization of the 1D configuration
-
+      
    IMPLICIT NONE
    PRIVATE
 
-   !! * Routine accessibility
-   PUBLIC ice_init                 ! called by opa.F90
+   PUBLIC   ice_init                 ! called by opa.F90
 
-   !! * Share Module variables
-   LOGICAL , PUBLIC  ::   & !!! ** init namelist (namicerun) **
-      ln_limdyn   = .TRUE.   !: flag for ice dynamics (T) or not (F)
-   INTEGER , PUBLIC  ::   &  !:
-      nstart ,            &  !: iteration number of the begining of the run 
-      nlast  ,            &  !: iteration number of the end of the run 
-      nitrun ,            &  !: number of iteration
-      numit                  !: iteration number
-   REAL(wp), PUBLIC  ::   &  !:
-      hsndif = 0.e0 ,     &  !: computation of temp. in snow (0) or not (9999)
-      hicdif = 0.e0 ,     &  !: computation of temp. in ice (0) or not (9999)
-      tpstot                 !: time of the run in seconds
-   REAL(wp), PUBLIC, DIMENSION(2)  ::  &  !:
-      acrit  = (/ 1.e-06 , 1.e-06 /)    !: minimum fraction for leads in 
-      !                                   !  north and south hemisphere
+   LOGICAL , PUBLIC               ::   ln_limdyn = .TRUE.   !: flag for ice dynamics (T) or not (F)
+   REAL(wp), PUBLIC               ::   hsndif = 0.e0        !: computation of temp. in snow (0) or not (9999)
+   REAL(wp), PUBLIC               ::   hicdif = 0.e0        !: computation of temp. in ice (0) or not (9999)
+   REAL(wp), PUBLIC, DIMENSION(2) ::   acrit  = (/ 1.e-06 , 1.e-06 /)    !: minimum fraction for leads in 
+      !                                                                  !  north and south hemisphere
    !!----------------------------------------------------------------------
    !!   LIM 2.0,  UCL-LOCEAN-IPSL (2005) 
@@ -56 +50 @@
       !!
       !! ** purpose :   
-      !!
-      !! History :
-      !!   8.5  !  02-08  (G. Madec)  F90: Free form and modules
       !!----------------------------------------------------------------------
-       CHARACTER(len=80) :: namelist_icename
-       
+      CHARACTER(len=80) :: namelist_icename
+      !!----------------------------------------------------------------------
+      !
       ! Open the namelist file 
       namelist_icename = 'namelist_ice'
-           
       CALL ctlopn(numnam_ice,namelist_icename,'OLD', 'FORMATTED', 'SEQUENTIAL',   &
                      1,numout,.FALSE.,1)      
-
       CALL ice_run                    !  read in namelist some run parameters
                  
@@ -83 +73 @@
       ! Initial sea-ice state
       IF( .NOT.ln_rstart ) THEN
-         numit = 0
          CALL lim_istate              ! start from rest: sea-ice deduced from sst
       ELSE
-         CALL lim_rst_read( numit )   ! start from a restart file
+         CALL lim_rst_read            ! start from a restart file
       ENDIF
       
@@ -94 +83 @@
       alb_ice(:,:) = albege(:,:)      ! sea-ice albedo
 # endif
-      
-      nstart = numit  + nfice      
-      nitrun = nitend - nit000 + 1 
-      nlast  = numit  + nitrun 
-
-      IF( nstock == 0  )  nstock = nlast + 1
-
+      !
    END SUBROUTINE ice_init
 
@@ -114 +97 @@
       !!
       !! ** input   :   Namelist namicerun
-      !!
-      !! history :
-      !!   2.0  !  03-08 (C. Ethe)  Original code
       !!-------------------------------------------------------------------
-
       NAMELIST/namicerun/ ln_limdyn, acrit, hsndif, hicdif
       !!-------------------------------------------------------------------
-
-      !                                           ! Read Namelist namicerun 
-      REWIND ( numnam_ice )
+      !                    
+      REWIND ( numnam_ice )                       ! Read Namelist namicerun 
       READ   ( numnam_ice , namicerun )
 
-      IF( lk_cfg_1d  )  ln_limdyn = .FALSE.       ! No ice transport in 1D configuration
+      IF( lk_cfg_1d  )   ln_limdyn = .FALSE.      ! No ice transport in 1D configuration
 
       IF(lwp) THEN
@@ -137 +115 @@
          WRITE(numout,*) '   computation of temp. in ice  (=0) or not (=9999) hicdif = ', hicdif
       ENDIF
+      !
    END SUBROUTINE ice_run