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.
limrst.F90 in branches/2014/dev_r4650_UKMO11_restart_functionality/NEMOGCM/NEMO/LIM_SRC_3 – NEMO

source: branches/2014/dev_r4650_UKMO11_restart_functionality/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90 @ 5231

Last change on this file since 5231 was 5231, checked in by davestorkey, 9 years ago

Svn keywords deactivated using "svn propdel" in
branch 2014/dev_r4650_UKMO11_restart_functionality.
File size: 24.2 KB
Line 
1MODULE limrst
2   !!======================================================================
3   !!                     ***  MODULE  limrst  ***
4   !! Ice restart :  write the ice restart file
5   !!======================================================================
6   !! History:   -   ! 2005-04 (M. Vancoppenolle) Original code
7   !!           3.0  ! 2008-03 (C. Ethe) restart files in using IOM interface
8   !!           4.0  ! 2011-02 (G. Madec) dynamical allocation
9   !!----------------------------------------------------------------------
10#if defined key_lim3
11   !!----------------------------------------------------------------------
12   !!   'key_lim3' :                                   LIM sea-ice model
13   !!----------------------------------------------------------------------
14   !!   lim_rst_opn   : open ice restart file
15   !!   lim_rst_write : write of the restart file
16   !!   lim_rst_read  : read  the restart file
17   !!----------------------------------------------------------------------
18   USE ice            ! sea-ice variables
19   USE oce     , ONLY :  snwice_mass, snwice_mass_b
20   USE par_ice        ! sea-ice parameters
21   USE dom_oce        ! ocean domain
22   USE sbc_oce        ! Surface boundary condition: ocean fields
23   USE sbc_ice        ! Surface boundary condition: ice fields
24   USE in_out_manager ! I/O manager
25   USE iom            ! I/O library
26   USE lib_mpp        ! MPP library
27   USE wrk_nemo       ! work arrays
28   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 
29
30   IMPLICIT NONE
31   PRIVATE
32
33   PUBLIC   lim_rst_opn    ! routine called by icestep.F90
34   PUBLIC   lim_rst_write  ! routine called by icestep.F90
35   PUBLIC   lim_rst_read   ! routine called by iceini.F90
36
37   LOGICAL, PUBLIC ::   lrst_ice         !: logical to control the ice restart write
38   INTEGER, PUBLIC ::   numrir, numriw   !: logical unit for ice restart (read and write)
39
40   !!----------------------------------------------------------------------
41   !! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2011)
42   !! $Id: limrst.F90 5208 2015-04-13 13:08:59Z davestorkey $
43   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
44   !!----------------------------------------------------------------------
45CONTAINS
46
47   SUBROUTINE lim_rst_opn( kt )
48      !!----------------------------------------------------------------------
49      !!                    ***  lim_rst_opn  ***
50      !!
51      !! ** purpose  :   output of sea-ice variable in a netcdf file
52      !!----------------------------------------------------------------------
53      INTEGER, INTENT(in) ::   kt       ! number of iteration
54      !
55      CHARACTER(LEN=20)   ::   clkt     ! ocean time-step define as a character
56      CHARACTER(LEN=50)   ::   clname   ! ice output restart file name
57      CHARACTER(len=150)  ::   clpath   ! full path to ice output restart file
58      !!----------------------------------------------------------------------
59      !
60      IF( kt == nit000 )   lrst_ice = .FALSE.   ! default definition
61
62      ! in order to get better performances with NetCDF format, we open and define the ice restart file
63      ! one ice time step before writing the data (-> at nitrst - 2*nn_fsbc + 1), except if we write ice
64      ! restart files every ice time step or if an ice restart file was writen at nitend - 2*nn_fsbc + 1
65      IF( kt == nitrst - 2*nn_fsbc + 1 .OR. nstock == nn_fsbc    &
66         &                             .OR. ( kt == nitend - nn_fsbc + 1 .AND. .NOT. lrst_ice ) ) THEN
67         IF( nitrst <= nitend .AND. nitrst > 0 ) THEN
68            ! beware of the format used to write kt (default is i8.8, that should be large enough...)
69            IF( nitrst > 99999999 ) THEN   ;   WRITE(clkt, *       ) nitrst
70            ELSE                           ;   WRITE(clkt, '(i8.8)') nitrst
71            ENDIF
72            ! create the file
73            clname = TRIM(cexper)//"_"//TRIM(ADJUSTL(clkt))//"_"//TRIM(cn_icerst_out)
74            clpath = TRIM(cn_icerst_outdir) 
75            IF( clpath(LEN_TRIM(clpath):) /= '/' ) clpath = TRIM(clpath)//'/'
76            IF(lwp) THEN
77               WRITE(numout,*)
78               SELECT CASE ( jprstlib )
79               CASE ( jprstdimg )
80                  WRITE(numout,*) '             open ice restart binary file: ',TRIM(clpath)//clname
81               CASE DEFAULT
82                  WRITE(numout,*) '             open ice restart NetCDF file: ',TRIM(clpath)//clname
83               END SELECT
84               IF( kt == nitrst - 2*nn_fsbc + 1 ) THEN   
85                  WRITE(numout,*)         '             kt = nitrst - 2*nn_fsbc + 1 = ', kt,' date= ', ndastp
86               ELSE   ;   WRITE(numout,*) '             kt = '                         , kt,' date= ', ndastp
87               ENDIF
88            ENDIF
89            !
90            CALL iom_open( TRIM(clpath)//TRIM(clname), numriw, ldwrt = .TRUE., kiolib = jprstlib )
91            lrst_ice = .TRUE.
92         ENDIF
93      ENDIF
94      !
95   END SUBROUTINE lim_rst_opn
96
97
98   SUBROUTINE lim_rst_write( kt )
99      !!----------------------------------------------------------------------
100      !!                    ***  lim_rst_write  ***
101      !!
102      !! ** purpose  :   output of sea-ice variable in a netcdf file
103      !!----------------------------------------------------------------------
104      INTEGER, INTENT(in) ::   kt     ! number of iteration
105      !!
106      INTEGER ::   ji, jj, jk ,jl   ! dummy loop indices
107      INTEGER ::   iter
108      CHARACTER(len=15) ::   znam
109      CHARACTER(len=1)  ::   zchar, zchar1
110      REAL(wp), POINTER, DIMENSION(:,:) :: z2d
111      !!----------------------------------------------------------------------
112
113      CALL wrk_alloc( jpi, jpj, z2d )
114
115      iter = kt + nn_fsbc - 1   ! ice restarts are written at kt == nitrst - nn_fsbc + 1
116
117      IF( iter == nitrst ) THEN
118         IF(lwp) WRITE(numout,*)
119         IF(lwp) WRITE(numout,*) 'lim_rst_write : write ice restart file  kt =', kt
120         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~'         
121      ENDIF
122
123      ! Write in numriw (if iter == nitrst)
124      ! ------------------
125      !                                                                        ! calendar control
126      CALL iom_rstput( iter, nitrst, numriw, 'nn_fsbc', REAL( nn_fsbc, wp ) )      ! time-step
127      CALL iom_rstput( iter, nitrst, numriw, 'kt_ice' , REAL( iter   , wp ) )      ! date
128
129      ! Prognostic variables
130      DO jl = 1, jpl 
131         WRITE(zchar,'(I1)') jl
132         znam = 'v_i'//'_htc'//zchar
133         z2d(:,:) = v_i(:,:,jl)
134         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
135         znam = 'v_s'//'_htc'//zchar
136         z2d(:,:) = v_s(:,:,jl)
137         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
138         znam = 'smv_i'//'_htc'//zchar
139         z2d(:,:) = smv_i(:,:,jl)
140         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
141         znam = 'oa_i'//'_htc'//zchar
142         z2d(:,:) = oa_i(:,:,jl)
143         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
144         znam = 'a_i'//'_htc'//zchar
145         z2d(:,:) = a_i(:,:,jl)
146         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
147         znam = 't_su'//'_htc'//zchar
148         z2d(:,:) = t_su(:,:,jl)
149         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
150      END DO
151
152      DO jl = 1, jpl 
153         WRITE(zchar,'(I1)') jl
154         znam = 'tempt_sl1'//'_htc'//zchar
155         z2d(:,:) = e_s(:,:,1,jl)
156         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
157      END DO
158
159      DO jl = 1, jpl 
160         WRITE(zchar,'(I1)') jl
161         DO jk = 1, nlay_i 
162            WRITE(zchar1,'(I1)') jk
163            znam = 'tempt'//'_il'//zchar1//'_htc'//zchar
164            z2d(:,:) = e_i(:,:,jk,jl)
165            CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
166         END DO
167      END DO
168
169      CALL iom_rstput( iter, nitrst, numriw, 'u_ice'        , u_ice      )
170      CALL iom_rstput( iter, nitrst, numriw, 'v_ice'        , v_ice      )
171      CALL iom_rstput( iter, nitrst, numriw, 'stress1_i'    , stress1_i  )
172      CALL iom_rstput( iter, nitrst, numriw, 'stress2_i'    , stress2_i  )
173      CALL iom_rstput( iter, nitrst, numriw, 'stress12_i'   , stress12_i )
174      CALL iom_rstput( iter, nitrst, numriw, 'snwice_mass'  , snwice_mass )   !clem modif
175      CALL iom_rstput( iter, nitrst, numriw, 'snwice_mass_b', snwice_mass_b ) !clem modif
176
177      DO jl = 1, jpl 
178         WRITE(zchar,'(I1)') jl
179         znam = 'sxice'//'_htc'//zchar
180         z2d(:,:) = sxice(:,:,jl)
181         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
182         znam = 'syice'//'_htc'//zchar
183         z2d(:,:) = syice(:,:,jl)
184         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
185         znam = 'sxxice'//'_htc'//zchar
186         z2d(:,:) = sxxice(:,:,jl)
187         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
188         znam = 'syyice'//'_htc'//zchar
189         z2d(:,:) = syyice(:,:,jl)
190         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
191         znam = 'sxyice'//'_htc'//zchar
192         z2d(:,:) = sxyice(:,:,jl)
193         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
194         znam = 'sxsn'//'_htc'//zchar
195         z2d(:,:) = sxsn(:,:,jl)
196         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
197         znam = 'sysn'//'_htc'//zchar
198         z2d(:,:) = sysn(:,:,jl)
199         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
200         znam = 'sxxsn'//'_htc'//zchar
201         z2d(:,:) = sxxsn(:,:,jl)
202         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
203         znam = 'syysn'//'_htc'//zchar
204         z2d(:,:) = syysn(:,:,jl)
205         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
206         znam = 'sxysn'//'_htc'//zchar
207         z2d(:,:) = sxysn(:,:,jl)
208         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
209         znam = 'sxa'//'_htc'//zchar
210         z2d(:,:) = sxa(:,:,jl)
211         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
212         znam = 'sya'//'_htc'//zchar
213         z2d(:,:) = sya(:,:,jl)
214         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
215         znam = 'sxxa'//'_htc'//zchar
216         z2d(:,:) = sxxa(:,:,jl)
217         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
218         znam = 'syya'//'_htc'//zchar
219         z2d(:,:) = syya(:,:,jl)
220         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
221         znam = 'sxya'//'_htc'//zchar
222         z2d(:,:) = sxya(:,:,jl)
223         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
224         znam = 'sxc0'//'_htc'//zchar
225         z2d(:,:) = sxc0(:,:,jl)
226         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
227         znam = 'syc0'//'_htc'//zchar
228         z2d(:,:) = syc0(:,:,jl)
229         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
230         znam = 'sxxc0'//'_htc'//zchar
231         z2d(:,:) = sxxc0(:,:,jl)
232         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
233         znam = 'syyc0'//'_htc'//zchar
234         z2d(:,:) = syyc0(:,:,jl)
235         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
236         znam = 'sxyc0'//'_htc'//zchar
237         z2d(:,:) = sxyc0(:,:,jl)
238         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
239         znam = 'sxsal'//'_htc'//zchar
240         z2d(:,:) = sxsal(:,:,jl)
241         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
242         znam = 'sysal'//'_htc'//zchar
243         z2d(:,:) = sysal(:,:,jl)
244         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
245         znam = 'sxxsal'//'_htc'//zchar
246         z2d(:,:) = sxxsal(:,:,jl)
247         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
248         znam = 'syysal'//'_htc'//zchar
249         z2d(:,:) = syysal(:,:,jl)
250         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
251         znam = 'sxysal'//'_htc'//zchar
252         z2d(:,:) = sxysal(:,:,jl)
253         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
254         znam = 'sxage'//'_htc'//zchar
255         z2d(:,:) = sxage(:,:,jl)
256         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
257         znam = 'syage'//'_htc'//zchar
258         z2d(:,:) = syage(:,:,jl)
259         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
260         znam = 'sxxage'//'_htc'//zchar
261         z2d(:,:) = sxxage(:,:,jl)
262         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
263         znam = 'syyage'//'_htc'//zchar
264         z2d(:,:) = syyage(:,:,jl)
265         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
266         znam = 'sxyage'//'_htc'//zchar
267         z2d(:,:) = sxyage(:,:,jl)
268         CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
269      END DO
270
271      CALL iom_rstput( iter, nitrst, numriw, 'sxopw ' ,  sxopw  )
272      CALL iom_rstput( iter, nitrst, numriw, 'syopw ' ,  syopw  )
273      CALL iom_rstput( iter, nitrst, numriw, 'sxxopw' ,  sxxopw )
274      CALL iom_rstput( iter, nitrst, numriw, 'syyopw' ,  syyopw )
275      CALL iom_rstput( iter, nitrst, numriw, 'sxyopw' ,  sxyopw )
276
277      DO jl = 1, jpl 
278         WRITE(zchar,'(I1)') jl
279         DO jk = 1, nlay_i 
280            WRITE(zchar1,'(I1)') jk
281            znam = 'sxe'//'_il'//zchar1//'_htc'//zchar
282            z2d(:,:) = sxe(:,:,jk,jl)
283            CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
284            znam = 'sye'//'_il'//zchar1//'_htc'//zchar
285            z2d(:,:) = sye(:,:,jk,jl)
286            CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
287            znam = 'sxxe'//'_il'//zchar1//'_htc'//zchar
288            z2d(:,:) = sxxe(:,:,jk,jl)
289            CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
290            znam = 'syye'//'_il'//zchar1//'_htc'//zchar
291            z2d(:,:) = syye(:,:,jk,jl)
292            CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
293            znam = 'sxye'//'_il'//zchar1//'_htc'//zchar
294            z2d(:,:) = sxye(:,:,jk,jl)
295            CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
296         END DO
297      END DO
298
299      IF( iter == nitrst ) THEN
300         CALL iom_close( numriw )                         ! close the restart file
301         lrst_ice = .FALSE.
302      ENDIF
303      !
304      CALL wrk_dealloc( jpi, jpj, z2d )
305      !
306   END SUBROUTINE lim_rst_write
307
308
309   SUBROUTINE lim_rst_read
310      !!----------------------------------------------------------------------
311      !!                    ***  lim_rst_read  ***
312      !!
313      !! ** purpose  :   read of sea-ice variable restart in a netcdf file
314      !!----------------------------------------------------------------------
315      INTEGER :: ji, jj, jk, jl, indx
316      REAL(wp) ::   zfice, ziter
317      REAL(wp) ::   zs_inf, z_slope_s, zsmax, zsmin, zalpha   ! local scalars used for the salinity profile
318      REAL(wp), POINTER, DIMENSION(:)   ::   zs_zero 
319      REAL(wp), POINTER, DIMENSION(:,:) ::   z2d
320      CHARACTER(len=15) ::   znam
321      CHARACTER(len=1)  ::   zchar, zchar1
322      INTEGER           ::   jlibalt = jprstlib
323      LOGICAL           ::   llok
324      !!----------------------------------------------------------------------
325
326      CALL wrk_alloc( nlay_i, zs_zero )
327      CALL wrk_alloc( jpi, jpj, z2d )
328
329      IF(lwp) THEN
330         WRITE(numout,*)
331         WRITE(numout,*) 'lim_rst_read : read ice NetCDF restart file'
332         WRITE(numout,*) '~~~~~~~~~~~~~'
333      ENDIF
334
335      IF ( jprstlib == jprstdimg ) THEN
336        ! eventually read netcdf file (monobloc)  for restarting on different number of processors
337        ! if {cn_icerst_in}.nc exists, then set jlibalt to jpnf90
338        INQUIRE( FILE = TRIM(cn_icerst_indir)//'/'//TRIM(cn_icerst_in)//'.nc', EXIST = llok )
339        IF ( llok ) THEN ; jlibalt = jpnf90  ; ELSE ; jlibalt = jprstlib ; ENDIF
340      ENDIF
341
342      CALL iom_open ( TRIM(cn_icerst_indir)//'/'//cn_icerst_in, numrir, kiolib = jprstlib )
343
344      CALL iom_get( numrir, 'nn_fsbc', zfice )
345      CALL iom_get( numrir, 'kt_ice' , ziter )   
346      IF(lwp) WRITE(numout,*) '   read ice restart file at time step    : ', ziter
347      IF(lwp) WRITE(numout,*) '   in any case we force it to nit000 - 1 : ', nit000 - 1
348
349      !Control of date
350
351      IF( ( nit000 - NINT(ziter) ) /= 1 .AND. ABS( nrstdt ) == 1 )   &
352         &     CALL ctl_stop( 'lim_rst_read ===>>>> : problem with nit000 in ice restart',  &
353         &                   '   verify the file or rerun with the value 0 for the',        &
354         &                   '   control of time parameter  nrstdt' )
355      IF( NINT(zfice) /= nn_fsbc          .AND. ABS( nrstdt ) == 1 )   &
356         &     CALL ctl_stop( 'lim_rst_read ===>>>> : problem with nn_fsbc in ice restart',  &
357         &                   '   verify the file or rerun with the value 0 for the',         &
358         &                   '   control of time parameter  nrstdt' )
359
360      DO jl = 1, jpl 
361         WRITE(zchar,'(I1)') jl
362         znam = 'v_i'//'_htc'//zchar
363         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
364         v_i(:,:,jl) = z2d(:,:)
365         znam = 'v_s'//'_htc'//zchar
366         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
367         v_s(:,:,jl) = z2d(:,:) 
368         znam = 'smv_i'//'_htc'//zchar
369         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
370         smv_i(:,:,jl) = z2d(:,:)
371         znam = 'oa_i'//'_htc'//zchar
372         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
373         oa_i(:,:,jl) = z2d(:,:)
374         znam = 'a_i'//'_htc'//zchar
375         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
376         a_i(:,:,jl) = z2d(:,:)
377         znam = 't_su'//'_htc'//zchar
378         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
379         t_su(:,:,jl) = z2d(:,:)
380      END DO
381
382      DO jl = 1, jpl 
383         WRITE(zchar,'(I1)') jl
384         znam = 'tempt_sl1'//'_htc'//zchar
385         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
386         e_s(:,:,1,jl) = z2d(:,:)
387      END DO
388
389      DO jl = 1, jpl 
390         WRITE(zchar,'(I1)') jl
391         DO jk = 1, nlay_i 
392            WRITE(zchar1,'(I1)') jk
393            znam = 'tempt'//'_il'//zchar1//'_htc'//zchar
394            CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
395            e_i(:,:,jk,jl) = z2d(:,:)
396         END DO
397      END DO
398
399      CALL iom_get( numrir, jpdom_autoglo, 'u_ice'     , u_ice      )
400      CALL iom_get( numrir, jpdom_autoglo, 'v_ice'     , v_ice      )
401      CALL iom_get( numrir, jpdom_autoglo, 'stress1_i' , stress1_i  )
402      CALL iom_get( numrir, jpdom_autoglo, 'stress2_i' , stress2_i  )
403      CALL iom_get( numrir, jpdom_autoglo, 'stress12_i', stress12_i )
404      CALL iom_get( numrir, jpdom_autoglo, 'snwice_mass'  , snwice_mass )   !clem modif
405      CALL iom_get( numrir, jpdom_autoglo, 'snwice_mass_b', snwice_mass_b ) !clem modif
406
407      DO jl = 1, jpl 
408         WRITE(zchar,'(I1)') jl
409         znam = 'sxice'//'_htc'//zchar
410         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
411         sxice(:,:,jl) = z2d(:,:)
412         znam = 'syice'//'_htc'//zchar
413         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
414         syice(:,:,jl) = z2d(:,:)
415         znam = 'sxxice'//'_htc'//zchar
416         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
417         sxxice(:,:,jl) = z2d(:,:)
418         znam = 'syyice'//'_htc'//zchar
419         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
420         syyice(:,:,jl) = z2d(:,:)
421         znam = 'sxyice'//'_htc'//zchar
422         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
423         sxyice(:,:,jl) = z2d(:,:)
424         znam = 'sxsn'//'_htc'//zchar
425         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
426         sxsn(:,:,jl) = z2d(:,:)
427         znam = 'sysn'//'_htc'//zchar
428         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
429         sysn(:,:,jl) = z2d(:,:)
430         znam = 'sxxsn'//'_htc'//zchar
431         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
432         sxxsn(:,:,jl) = z2d(:,:)
433         znam = 'syysn'//'_htc'//zchar
434         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
435         syysn(:,:,jl) = z2d(:,:)
436         znam = 'sxysn'//'_htc'//zchar
437         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
438         sxysn(:,:,jl) = z2d(:,:)
439         znam = 'sxa'//'_htc'//zchar
440         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
441         sxa(:,:,jl) = z2d(:,:)
442         znam = 'sya'//'_htc'//zchar
443         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
444         sya(:,:,jl) = z2d(:,:)
445         znam = 'sxxa'//'_htc'//zchar
446         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
447         sxxa(:,:,jl) = z2d(:,:)
448         znam = 'syya'//'_htc'//zchar
449         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
450         syya(:,:,jl) = z2d(:,:)
451         znam = 'sxya'//'_htc'//zchar
452         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
453         sxya(:,:,jl) = z2d(:,:)
454         znam = 'sxc0'//'_htc'//zchar
455         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
456         sxc0(:,:,jl) = z2d(:,:)
457         znam = 'syc0'//'_htc'//zchar
458         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
459         syc0(:,:,jl) = z2d(:,:)
460         znam = 'sxxc0'//'_htc'//zchar
461         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
462         sxxc0(:,:,jl) = z2d(:,:)
463         znam = 'syyc0'//'_htc'//zchar
464         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
465         syyc0(:,:,jl) = z2d(:,:)
466         znam = 'sxyc0'//'_htc'//zchar
467         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
468         sxyc0(:,:,jl) = z2d(:,:)
469         znam = 'sxsal'//'_htc'//zchar
470         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
471         sxsal(:,:,jl) = z2d(:,:)
472         znam = 'sysal'//'_htc'//zchar
473         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
474         sysal(:,:,jl) = z2d(:,:)
475         znam = 'sxxsal'//'_htc'//zchar
476         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
477         sxxsal(:,:,jl) = z2d(:,:)
478         znam = 'syysal'//'_htc'//zchar
479         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
480         syysal(:,:,jl) = z2d(:,:)
481         znam = 'sxysal'//'_htc'//zchar
482         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
483         sxysal(:,:,jl) = z2d(:,:)
484         znam = 'sxage'//'_htc'//zchar
485         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
486         sxage(:,:,jl) = z2d(:,:)
487         znam = 'syage'//'_htc'//zchar
488         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
489         syage(:,:,jl) = z2d(:,:)
490         znam = 'sxxage'//'_htc'//zchar
491         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
492         sxxage(:,:,jl) = z2d(:,:)
493         znam = 'syyage'//'_htc'//zchar
494         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
495         syyage(:,:,jl) = z2d(:,:)
496         znam = 'sxyage'//'_htc'//zchar
497         CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
498         sxyage(:,:,jl)= z2d(:,:)
499      END DO
500
501      CALL iom_get( numrir, jpdom_autoglo, 'sxopw ' ,  sxopw  )
502      CALL iom_get( numrir, jpdom_autoglo, 'syopw ' ,  syopw  )
503      CALL iom_get( numrir, jpdom_autoglo, 'sxxopw' ,  sxxopw )
504      CALL iom_get( numrir, jpdom_autoglo, 'syyopw' ,  syyopw )
505      CALL iom_get( numrir, jpdom_autoglo, 'sxyopw' ,  sxyopw )
506
507      DO jl = 1, jpl 
508         WRITE(zchar,'(I1)') jl
509         DO jk = 1, nlay_i 
510            WRITE(zchar1,'(I1)') jk
511            znam = 'sxe'//'_il'//zchar1//'_htc'//zchar
512            CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
513            sxe(:,:,jk,jl) = z2d(:,:)
514            znam = 'sye'//'_il'//zchar1//'_htc'//zchar
515            CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
516            sye(:,:,jk,jl) = z2d(:,:)
517            znam = 'sxxe'//'_il'//zchar1//'_htc'//zchar
518            CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
519            sxxe(:,:,jk,jl) = z2d(:,:)
520            znam = 'syye'//'_il'//zchar1//'_htc'//zchar
521            CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
522            syye(:,:,jk,jl) = z2d(:,:)
523            znam = 'sxye'//'_il'//zchar1//'_htc'//zchar
524            CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
525            sxye(:,:,jk,jl) = z2d(:,:)
526         END DO
527      END DO
528      !
529      ! clem: I do not understand why the following IF is needed
530      !       I suspect something inconsistent in the main code with option num_sal=1
531      IF( num_sal == 1 ) THEN
532         DO jl = 1, jpl 
533            sm_i(:,:,jl) = bulk_sal
534            DO jk = 1, nlay_i 
535               s_i(:,:,jk,jl) = bulk_sal
536            END DO
537         END DO
538      ENDIF
539      !
540      !CALL iom_close( numrir ) !clem: closed in sbcice_lim.F90
541      !
542      CALL wrk_dealloc( nlay_i, zs_zero )
543      CALL wrk_dealloc( jpi, jpj, z2d )
544      !
545   END SUBROUTINE lim_rst_read
546
547#else
548   !!----------------------------------------------------------------------
549   !!   Default option :       Empty module            NO LIM sea-ice model
550   !!----------------------------------------------------------------------
551CONTAINS
552   SUBROUTINE lim_rst_read             ! Empty routine
553   END SUBROUTINE lim_rst_read
554   SUBROUTINE lim_rst_write            ! Empty routine
555   END SUBROUTINE lim_rst_write
556#endif
557
558   !!======================================================================
559END MODULE limrst
Note: See TracBrowser for help on using the repository browser.