New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
limistate_2.F90 in trunk/NEMOGCM/NEMO/LIM_SRC_2 – NEMO

source: trunk/NEMOGCM/NEMO/LIM_SRC_2/limistate_2.F90 @ 4696

Last change on this file since 4696 was 4624, checked in by acc, 10 years ago

#1305. Fix slow start-up problems on some systems by introducing and using lwm logical to restrict output of merged namelists to the first (or only) processor. lwm is true only on the first processor regardless of ln_ctl. Small changes to all flavours of nemogcm.F90 are also required to write namctl and namcfg after the call to mynode which now opens output.namelist.dyn and writes nammpp.
  • Property svn:keywords set to Id
File size: 10.4 KB
RevLine 
[821]1MODULE limistate_2
[3]2   !!======================================================================
[821]3   !!                     ***  MODULE  limistate_2  ***
[3]4   !!              Initialisation of diagnostics ice variables
5   !!======================================================================
[888]6   !! History :   1.0  !  01-04  (C. Ethe, G. Madec)  Original code
7   !!             2.0  !  03-08  (G. Madec)  add lim_istate_init
[508]8   !!                  !  04-04  (S. Theetten) initialization from a file
9   !!                  !  06-07  (S. Masson)  IOM to read the restart
[888]10   !!                  !  07-10  (G. Madec)  surface module
[508]11   !!--------------------------------------------------------------------
[821]12#if defined key_lim2
[3]13   !!----------------------------------------------------------------------
[821]14   !!   'key_lim2' :                                  LIM 2.0 sea-ice model
[3]15   !!----------------------------------------------------------------------
[508]16   !!----------------------------------------------------------------------
[821]17   !!   lim_istate_2      :  Initialisation of diagnostics ice variables
18   !!   lim_istate_init_2 :  initialization of ice state and namelist read
[3]19   !!----------------------------------------------------------------------
20   USE phycst
[821]21   USE par_ice_2       ! ice parameters
22   USE dom_ice_2
[1037]23   USE eosbn2          ! equation of state
[508]24   USE lbclnk
[888]25   USE oce
[821]26   USE ice_2
[508]27   USE iom
28   USE in_out_manager
[3625]29   USE lib_fortran     ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 
[3]30
31   IMPLICIT NONE
32   PRIVATE
33
[821]34   PUBLIC lim_istate_2      ! routine called by lim_init_2.F90
[3]35
[4147]36   !!! **  namelist (namiceini) **
37   LOGICAL  ::   ln_limini   !: Ice initialization state
38   REAL(wp) ::   ttest       ! threshold water temperature for initial sea ice
39   REAL(wp) ::   hninn       ! initial snow thickness in the north
40   REAL(wp) ::   hginn       ! initial ice thickness in the north
41   REAL(wp) ::   alinn       ! initial leads area in the north
42   REAL(wp) ::   hnins       ! initial snow thickness in the south
43   REAL(wp) ::   hgins       ! initial ice thickness in the south
44   REAL(wp) ::   alins       ! initial leads area in the south
[1228]45   
46   REAL(wp) ::   zero      = 0.e0     ! constant value = 0
47   REAL(wp) ::   zone      = 1.e0     ! constant value = 1
[3]48   !!----------------------------------------------------------------------
[2528]49   !! NEMO/LIM2 3.3 , UCL - NEMO Consortium (2010)
[1156]50   !! $Id$
[2528]51   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[3]52   !!----------------------------------------------------------------------
53
54CONTAINS
55
[821]56   SUBROUTINE lim_istate_2
[3]57      !!-------------------------------------------------------------------
[821]58      !!                    ***  ROUTINE lim_istate_2  ***
[3]59      !!
[77]60      !! ** Purpose :   defined the sea-ice initial state
61      !!
62      !! ** Method  :   restart from a state defined in a binary file
63      !!                or from arbitrary sea-ice conditions
[3]64      !!--------------------------------------------------------------------
[508]65      INTEGER  ::   ji, jj, jk                ! dummy loop indices
[1045]66      REAL(wp) ::   zidto                     ! temporary scalar
[3]67      !--------------------------------------------------------------------
[888]68 
69      CALL lim_istate_init_2     !  reading the initials parameters of the ice
[3]70
[888]71      IF( .NOT. ln_limini ) THEN 
[419]72         
[3294]73         tfu(:,:) = tfreez( tsn(:,:,1,jp_sal) ) * tmask(:,:,1)       ! freezing/melting point of sea water [Celcius]
[1037]74
[419]75         DO jj = 1, jpj
76            DO ji = 1, jpi
[1037]77               !                     ! ice if sst <= t-freez + ttest
[3294]78               IF( tsn(ji,jj,1,jp_tem)  - tfu(ji,jj) >= ttest ) THEN   ;   zidto = 0.e0      ! no ice
79               ELSE                                                    ;   zidto = 1.e0      !    ice
[1037]80               ENDIF
81               !
[419]82               IF( fcor(ji,jj) >= 0.e0 ) THEN     !--  Northern hemisphere.
83                  hicif(ji,jj)   = zidto * hginn
84                  frld(ji,jj)    = zidto * alinn + ( 1.0 - zidto ) * 1.0
85                  hsnif(ji,jj)   = zidto * hninn
86               ELSE                               !---  Southern hemisphere.
87                  hicif(ji,jj)   = zidto * hgins
88                  frld(ji,jj)    = zidto * alins + ( 1.0 - zidto ) * 1.0
89                  hsnif(ji,jj)   = zidto * hnins
90               ENDIF
91            END DO
92         END DO
[1037]93
94         tfu(:,:) = tfu(:,:) + rt0       ! ftu converted from Celsius to Kelvin (rt0 over land)
[419]95         
96         sist  (:,:)   = tfu(:,:)
97         tbif  (:,:,1) = tfu(:,:)
98         tbif  (:,:,2) = tfu(:,:)
99         tbif  (:,:,3) = tfu(:,:)
[888]100
[419]101      ENDIF
[888]102     
[3]103      fsbbq (:,:)   = 0.e0
104      qstoif(:,:)   = 0.e0
[1470]105      u_ice (:,:)   = 0.e0
106      v_ice (:,:)   = 0.e0
[3]107
[77]108      !---  Moments for advection.             
[3]109
[77]110      sxice (:,:)  = 0.e0   ;   sxsn (:,:)  = 0.e0   ;   sxa  (:,:)  = 0.e0
111      syice (:,:)  = 0.e0   ;   sysn (:,:)  = 0.e0   ;   sya  (:,:)  = 0.e0
112      sxxice(:,:)  = 0.e0   ;   sxxsn(:,:)  = 0.e0   ;   sxxa (:,:)  = 0.e0
113      syyice(:,:)  = 0.e0   ;   syysn(:,:)  = 0.e0   ;   syya (:,:)  = 0.e0
114      sxyice(:,:)  = 0.e0   ;   sxysn(:,:)  = 0.e0   ;   sxya (:,:)  = 0.e0
[3]115
[77]116      sxc0  (:,:)  = 0.e0   ;   sxc1 (:,:)  = 0.e0   ;   sxc2 (:,:)  = 0.e0
117      syc0  (:,:)  = 0.e0   ;   syc1 (:,:)  = 0.e0   ;   syc2 (:,:)  = 0.e0
118      sxxc0 (:,:)  = 0.e0   ;   sxxc1(:,:)  = 0.e0   ;   sxxc2(:,:)  = 0.e0
119      syyc0 (:,:)  = 0.e0   ;   syyc1(:,:)  = 0.e0   ;   syyc2(:,:)  = 0.e0
120      sxyc0 (:,:)  = 0.e0   ;   sxyc1(:,:)  = 0.e0   ;   sxyc2(:,:)  = 0.e0
[3]121
[77]122      sxst  (:,:)  = 0.e0
123      syst  (:,:)  = 0.e0
124      sxxst (:,:)  = 0.e0
125      syyst (:,:)  = 0.e0
126      sxyst (:,:)  = 0.e0
[2855]127#if ! defined key_lim2_vp
128      stress1_i (:,:) = 0._wp                          ! EVP rheology
129      stress2_i (:,:) = 0._wp
130      stress12_i(:,:) = 0._wp
131#endif
[3]132
[77]133      !-- lateral boundary conditions
134      CALL lbc_lnk( hicif, 'T', 1. )
135      CALL lbc_lnk( frld , 'T', 1. )
[299]136
137      ! C A U T I O N  frld = 1 over land and lbc_lnk put zero along
138      ! *************  closed boundaries herefore we force to one over land
139      frld(:,:) = tms(:,:) * frld(:,:) + ( 1. - tms(:,:) )   
140
[77]141      CALL lbc_lnk( hsnif, 'T', 1. )
[888]142      CALL lbc_lnk( sist , 'T', 1. , pval = rt0 )      ! set rt0 on closed boundary (required by bulk formulation)
[77]143      DO jk = 1, jplayersp1
144         CALL lbc_lnk(tbif(:,:,jk), 'T', 1. )
145      END DO
146      CALL lbc_lnk( fsbbq  , 'T', 1. )
147      CALL lbc_lnk( qstoif , 'T', 1. )
[888]148
[821]149   END SUBROUTINE lim_istate_2
[3]150
151   
[821]152   SUBROUTINE lim_istate_init_2
[3]153      !!-------------------------------------------------------------------
[821]154      !!                   ***  ROUTINE lim_istate_init_2  ***
[3]155      !!       
[77]156      !! ** Purpose :   Definition of initial state of the ice
[3]157      !!
[77]158      !! ** Method  :   Read the namiceini namelist and check the parameter
[888]159      !!       values called at the first timestep (nit000)
[3]160      !!
[77]161      !! ** input   :   Namelist namiceini
[3]162      !!-------------------------------------------------------------------
[508]163      INTEGER :: inum_ice
[673]164      INTEGER :: ji,jj
[4147]165      INTEGER :: ios                 ! Local integer output status for namelist read
[419]166
[675]167      NAMELIST/namiceini/ ln_limini, ttest, hninn, hginn, alinn, &
[419]168         &                hnins, hgins, alins
[3]169      !!-------------------------------------------------------------------
[4147]170                   
171      REWIND( numnam_ice_ref )              ! Namelist namiceini in reference namelist : Ice initial state
172      READ  ( numnam_ice_ref, namiceini, IOSTAT = ios, ERR = 901)
173901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceini in reference namelist', lwp )
174
175      REWIND( numnam_ice_cfg )              ! Namelist namiceini in configuration namelist : Ice initial state
176      READ  ( numnam_ice_cfg, namiceini, IOSTAT = ios, ERR = 902 )
177902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceini in configuration namelist', lwp )
[4624]178      IF(lwm) WRITE ( numoni, namiceini )
[888]179      !
180      IF(lwp) THEN
181         WRITE(numout,*)
182         WRITE(numout,*) 'lim_istate_init_2 : ice parameters inititialisation '
183         WRITE(numout,*) '~~~~~~~~~~~~~~~~~'
184         WRITE(numout,*) '         threshold water temp. for initial sea-ice    ttest      = ', ttest
185         WRITE(numout,*) '         initial snow thickness in the north          hninn      = ', hninn
186         WRITE(numout,*) '         initial ice thickness in the north           hginn      = ', hginn 
187         WRITE(numout,*) '         initial leads area in the north              alinn      = ', alinn           
188         WRITE(numout,*) '         initial snow thickness in the south          hnins      = ', hnins 
189         WRITE(numout,*) '         initial ice thickness in the south           hgins      = ', hgins
190         WRITE(numout,*) '         initial leads area in the south              alins      = ', alins
191         WRITE(numout,*) '         Ice state initialization using input file    ln_limini  = ', ln_limini
[3]192      ENDIF
[419]193
194      IF( ln_limini ) THEN                      ! Ice initialization using input file
[888]195         !
[508]196         CALL iom_open( 'Ice_initialization.nc', inum_ice )
[888]197         !
[508]198         IF( inum_ice > 0 ) THEN
[888]199            IF(lwp) WRITE(numout,*)
200            IF(lwp) WRITE(numout,*) '                  ice state initialization with : Ice_initialization.nc'
[3]201           
[888]202            CALL iom_get( inum_ice, jpdom_data, 'hicif', hicif )     
203            CALL iom_get( inum_ice, jpdom_data, 'hsnif', hsnif )     
204            CALL iom_get( inum_ice, jpdom_data, 'frld' , frld  )     
205            CALL iom_get( inum_ice, jpdom_data, 'ts'   , sist  )
[673]206            CALL iom_get( inum_ice, jpdom_unknown, 'tbif', tbif(1:nlci,1:nlcj,:),   &
[508]207                 &        kstart = (/ mig(1),mjg(1),1 /), kcount = (/ nlci,nlcj,jplayersp1 /) )
[673]208            ! put some values in the extra-halo...
209            DO jj = nlcj+1, jpj   ;   tbif(1:nlci,jj,:) = tbif(1:nlci,nlej,:)   ;   END DO
210            DO ji = nlci+1, jpi   ;   tbif(ji    ,: ,:) = tbif(nlei  ,:   ,:)   ;   END DO
[508]211
212            CALL iom_close( inum_ice)
[888]213            !
[419]214         ENDIF
215      ENDIF
[888]216      !     
[821]217   END SUBROUTINE lim_istate_init_2
[3]218
219#else
220   !!----------------------------------------------------------------------
[821]221   !!   Default option :         Empty module      NO LIM 2.0 sea-ice model
[3]222   !!----------------------------------------------------------------------
223CONTAINS
[821]224   SUBROUTINE lim_istate_2        ! Empty routine
225   END SUBROUTINE lim_istate_2
[3]226#endif
227
228   !!======================================================================
[821]229END MODULE limistate_2
Note: See TracBrowser for help on using the repository browser.