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/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/LIM_SRC_3 – NEMO

Context Navigation

source: branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90 @ 7772

Last change on this file since 7772 was 7772, checked in by clem, 7 years ago
fix ice restart for categories>9, see ticket #1819
Property svn:keywords set to `Id`
File size: 25.1 KB

Line
1	MODULE 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 dom_oce ! ocean domain
21	USE sbc_oce ! Surface boundary condition: ocean fields
22	USE sbc_ice ! Surface boundary condition: ice fields
23	USE in_out_manager ! I/O manager
24	USE iom ! I/O library
25	USE lib_mpp ! MPP library
26	USE wrk_nemo ! work arrays
27	USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)
28	USE limctl
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 sbc_lim_init
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$
43	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
44	!!----------------------------------------------------------------------
45	CONTAINS
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=256) :: 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 - 2nn_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	IF( ln_icectl ) CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - Beginning the time step - ' ) ! control print
96	END SUBROUTINE lim_rst_opn
97
98
99	SUBROUTINE lim_rst_write( kt )
100	!!----------------------------------------------------------------------
101	!! * lim_rst_write *
102	!!
103	!! ** purpose : output of sea-ice variable in a netcdf file
104	!!----------------------------------------------------------------------
105	INTEGER, INTENT(in) :: kt ! number of iteration
106	!!
107	INTEGER :: ji, jj, jk ,jl ! dummy loop indices
108	INTEGER :: iter
109	CHARACTER(len=15) :: znam
110	CHARACTER(len=2) :: zchar, zchar1
111	REAL(wp), POINTER, DIMENSION(:,:) :: z2d
112	!!----------------------------------------------------------------------
113
114	CALL wrk_alloc( jpi, jpj, z2d )
115
116	iter = kt + nn_fsbc - 1 ! ice restarts are written at kt == nitrst - nn_fsbc + 1
117
118	IF( iter == nitrst ) THEN
119	IF(lwp) WRITE(numout,*)
120	IF(lwp) WRITE(numout,*) 'lim_rst_write : write ice restart file kt =', kt
121	IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~'
122	ENDIF
123
124	! Write in numriw (if iter == nitrst)
125	! ------------------
126	! ! calendar control
127	CALL iom_rstput( iter, nitrst, numriw, 'nn_fsbc', REAL( nn_fsbc, wp ) ) ! time-step
128	CALL iom_rstput( iter, nitrst, numriw, 'kt_ice' , REAL( iter , wp ) ) ! date
129
130	! Prognostic variables
131	DO jl = 1, jpl
132	WRITE(zchar,'(I2)') jl
133	znam = 'v_i'//'_htc'//TRIM(ADJUSTL(zchar))
134	z2d(:,:) = v_i(:,:,jl)
135	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
136	znam = 'v_s'//'_htc'//TRIM(ADJUSTL(zchar))
137	z2d(:,:) = v_s(:,:,jl)
138	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
139	znam = 'smv_i'//'_htc'//TRIM(ADJUSTL(zchar))
140	z2d(:,:) = smv_i(:,:,jl)
141	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
142	znam = 'oa_i'//'_htc'//TRIM(ADJUSTL(zchar))
143	z2d(:,:) = oa_i(:,:,jl)
144	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
145	znam = 'a_i'//'_htc'//TRIM(ADJUSTL(zchar))
146	z2d(:,:) = a_i(:,:,jl)
147	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
148	znam = 't_su'//'_htc'//TRIM(ADJUSTL(zchar))
149	z2d(:,:) = t_su(:,:,jl)
150	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
151	znam = 'tempt_sl1'//'_htc'//TRIM(ADJUSTL(zchar))
152	z2d(:,:) = e_s(:,:,1,jl)
153	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
154	DO jk = 1, nlay_i
155	WRITE(zchar1,'(I2)') jk
156	znam = 'tempt'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
157	z2d(:,:) = e_i(:,:,jk,jl)
158	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
159	END DO
160	END DO
161
162	CALL iom_rstput( iter, nitrst, numriw, 'u_ice' , u_ice )
163	CALL iom_rstput( iter, nitrst, numriw, 'v_ice' , v_ice )
164	CALL iom_rstput( iter, nitrst, numriw, 'stress1_i' , stress1_i )
165	CALL iom_rstput( iter, nitrst, numriw, 'stress2_i' , stress2_i )
166	CALL iom_rstput( iter, nitrst, numriw, 'stress12_i' , stress12_i )
167	CALL iom_rstput( iter, nitrst, numriw, 'snwice_mass' , snwice_mass )
168	CALL iom_rstput( iter, nitrst, numriw, 'snwice_mass_b', snwice_mass_b )
169
170	DO jl = 1, jpl
171	WRITE(zchar,'(I2)') jl
172	znam = 'sxice'//'_htc'//TRIM(ADJUSTL(zchar))
173	z2d(:,:) = sxice(:,:,jl)
174	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
175	znam = 'syice'//'_htc'//TRIM(ADJUSTL(zchar))
176	z2d(:,:) = syice(:,:,jl)
177	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
178	znam = 'sxxice'//'_htc'//TRIM(ADJUSTL(zchar))
179	z2d(:,:) = sxxice(:,:,jl)
180	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
181	znam = 'syyice'//'_htc'//TRIM(ADJUSTL(zchar))
182	z2d(:,:) = syyice(:,:,jl)
183	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
184	znam = 'sxyice'//'_htc'//TRIM(ADJUSTL(zchar))
185	z2d(:,:) = sxyice(:,:,jl)
186	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
187	znam = 'sxsn'//'_htc'//TRIM(ADJUSTL(zchar))
188	z2d(:,:) = sxsn(:,:,jl)
189	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
190	znam = 'sysn'//'_htc'//TRIM(ADJUSTL(zchar))
191	z2d(:,:) = sysn(:,:,jl)
192	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
193	znam = 'sxxsn'//'_htc'//TRIM(ADJUSTL(zchar))
194	z2d(:,:) = sxxsn(:,:,jl)
195	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
196	znam = 'syysn'//'_htc'//TRIM(ADJUSTL(zchar))
197	z2d(:,:) = syysn(:,:,jl)
198	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
199	znam = 'sxysn'//'_htc'//TRIM(ADJUSTL(zchar))
200	z2d(:,:) = sxysn(:,:,jl)
201	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
202	znam = 'sxa'//'_htc'//TRIM(ADJUSTL(zchar))
203	z2d(:,:) = sxa(:,:,jl)
204	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
205	znam = 'sya'//'_htc'//TRIM(ADJUSTL(zchar))
206	z2d(:,:) = sya(:,:,jl)
207	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
208	znam = 'sxxa'//'_htc'//TRIM(ADJUSTL(zchar))
209	z2d(:,:) = sxxa(:,:,jl)
210	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
211	znam = 'syya'//'_htc'//TRIM(ADJUSTL(zchar))
212	z2d(:,:) = syya(:,:,jl)
213	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
214	znam = 'sxya'//'_htc'//TRIM(ADJUSTL(zchar))
215	z2d(:,:) = sxya(:,:,jl)
216	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
217	znam = 'sxc0'//'_htc'//TRIM(ADJUSTL(zchar))
218	z2d(:,:) = sxc0(:,:,jl)
219	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
220	znam = 'syc0'//'_htc'//TRIM(ADJUSTL(zchar))
221	z2d(:,:) = syc0(:,:,jl)
222	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
223	znam = 'sxxc0'//'_htc'//TRIM(ADJUSTL(zchar))
224	z2d(:,:) = sxxc0(:,:,jl)
225	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
226	znam = 'syyc0'//'_htc'//TRIM(ADJUSTL(zchar))
227	z2d(:,:) = syyc0(:,:,jl)
228	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
229	znam = 'sxyc0'//'_htc'//TRIM(ADJUSTL(zchar))
230	z2d(:,:) = sxyc0(:,:,jl)
231	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
232	znam = 'sxsal'//'_htc'//TRIM(ADJUSTL(zchar))
233	z2d(:,:) = sxsal(:,:,jl)
234	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
235	znam = 'sysal'//'_htc'//TRIM(ADJUSTL(zchar))
236	z2d(:,:) = sysal(:,:,jl)
237	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
238	znam = 'sxxsal'//'_htc'//TRIM(ADJUSTL(zchar))
239	z2d(:,:) = sxxsal(:,:,jl)
240	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
241	znam = 'syysal'//'_htc'//TRIM(ADJUSTL(zchar))
242	z2d(:,:) = syysal(:,:,jl)
243	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
244	znam = 'sxysal'//'_htc'//TRIM(ADJUSTL(zchar))
245	z2d(:,:) = sxysal(:,:,jl)
246	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
247	znam = 'sxage'//'_htc'//TRIM(ADJUSTL(zchar))
248	z2d(:,:) = sxage(:,:,jl)
249	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
250	znam = 'syage'//'_htc'//TRIM(ADJUSTL(zchar))
251	z2d(:,:) = syage(:,:,jl)
252	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
253	znam = 'sxxage'//'_htc'//TRIM(ADJUSTL(zchar))
254	z2d(:,:) = sxxage(:,:,jl)
255	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
256	znam = 'syyage'//'_htc'//TRIM(ADJUSTL(zchar))
257	z2d(:,:) = syyage(:,:,jl)
258	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
259	znam = 'sxyage'//'_htc'//TRIM(ADJUSTL(zchar))
260	z2d(:,:) = sxyage(:,:,jl)
261	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
262	END DO
263
264	CALL iom_rstput( iter, nitrst, numriw, 'sxopw ' , sxopw )
265	CALL iom_rstput( iter, nitrst, numriw, 'syopw ' , syopw )
266	CALL iom_rstput( iter, nitrst, numriw, 'sxxopw' , sxxopw )
267	CALL iom_rstput( iter, nitrst, numriw, 'syyopw' , syyopw )
268	CALL iom_rstput( iter, nitrst, numriw, 'sxyopw' , sxyopw )
269
270	DO jl = 1, jpl
271	WRITE(zchar,'(I2)') jl
272	DO jk = 1, nlay_i
273	WRITE(zchar1,'(I2)') jk
274	znam = 'sxe'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
275	z2d(:,:) = sxe(:,:,jk,jl)
276	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
277	znam = 'sye'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
278	z2d(:,:) = sye(:,:,jk,jl)
279	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
280	znam = 'sxxe'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
281	z2d(:,:) = sxxe(:,:,jk,jl)
282	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
283	znam = 'syye'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
284	z2d(:,:) = syye(:,:,jk,jl)
285	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
286	znam = 'sxye'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
287	z2d(:,:) = sxye(:,:,jk,jl)
288	CALL iom_rstput( iter, nitrst, numriw, znam , z2d )
289	END DO
290	END DO
291
292	IF( iter == nitrst ) THEN
293	CALL iom_close( numriw ) ! close the restart file
294	lrst_ice = .FALSE.
295	ENDIF
296	!
297	CALL wrk_dealloc( jpi, jpj, z2d )
298	!
299	END SUBROUTINE lim_rst_write
300
301
302	SUBROUTINE lim_rst_read
303	!!----------------------------------------------------------------------
304	!! * lim_rst_read *
305	!!
306	!! ** purpose : read of sea-ice variable restart in a netcdf file
307	!!----------------------------------------------------------------------
308	INTEGER :: ji, jj, jk, jl
309	REAL(wp) :: zfice, ziter
310	REAL(wp), POINTER, DIMENSION(:,:) :: z2d
311	CHARACTER(len=15) :: znam
312	CHARACTER(len=2) :: zchar, zchar1
313	INTEGER :: jlibalt = jprstlib
314	LOGICAL :: llok
315	!!----------------------------------------------------------------------
316
317	CALL wrk_alloc( jpi, jpj, z2d )
318
319	IF(lwp) THEN
320	WRITE(numout,*)
321	WRITE(numout,*) 'lim_rst_read : read ice NetCDF restart file'
322	WRITE(numout,*) '~~~~~~~~~~~~~'
323	ENDIF
324
325	IF ( jprstlib == jprstdimg ) THEN
326	! eventually read netcdf file (monobloc) for restarting on different number of processors
327	! if {cn_icerst_in}.nc exists, then set jlibalt to jpnf90
328	INQUIRE( FILE = TRIM(cn_icerst_indir)//'/'//TRIM(cn_icerst_in)//'.nc', EXIST = llok )
329	IF ( llok ) THEN ; jlibalt = jpnf90 ; ELSE ; jlibalt = jprstlib ; ENDIF
330	ENDIF
331
332	CALL iom_open ( TRIM(cn_icerst_indir)//'/'//cn_icerst_in, numrir, kiolib = jprstlib )
333
334	CALL iom_get( numrir, 'nn_fsbc', zfice )
335	CALL iom_get( numrir, 'kt_ice' , ziter )
336	IF(lwp) WRITE(numout,*) ' read ice restart file at time step : ', ziter
337	IF(lwp) WRITE(numout,*) ' in any case we force it to nit000 - 1 : ', nit000 - 1
338
339	!Control of date
340
341	IF( ( nit000 - NINT(ziter) ) /= 1 .AND. ABS( nrstdt ) == 1 ) &
342	& CALL ctl_stop( 'lim_rst_read ===>>>> : problem with nit000 in ice restart', &
343	& ' verify the file or rerun with the value 0 for the', &
344	& ' control of time parameter nrstdt' )
345	IF( NINT(zfice) /= nn_fsbc .AND. ABS( nrstdt ) == 1 ) &
346	& CALL ctl_stop( 'lim_rst_read ===>>>> : problem with nn_fsbc in ice restart', &
347	& ' verify the file or rerun with the value 0 for the', &
348	& ' control of time parameter nrstdt' )
349
350	DO jl = 1, jpl
351	WRITE(zchar,'(I2)') jl
352	znam = 'v_i'//'_htc'//TRIM(ADJUSTL(zchar))
353	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
354	v_i(:,:,jl) = z2d(:,:)
355	znam = 'v_s'//'_htc'//TRIM(ADJUSTL(zchar))
356	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
357	v_s(:,:,jl) = z2d(:,:)
358	znam = 'smv_i'//'_htc'//TRIM(ADJUSTL(zchar))
359	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
360	smv_i(:,:,jl) = z2d(:,:)
361	znam = 'oa_i'//'_htc'//TRIM(ADJUSTL(zchar))
362	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
363	oa_i(:,:,jl) = z2d(:,:)
364	znam = 'a_i'//'_htc'//TRIM(ADJUSTL(zchar))
365	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
366	a_i(:,:,jl) = z2d(:,:)
367	znam = 't_su'//'_htc'//TRIM(ADJUSTL(zchar))
368	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
369	t_su(:,:,jl) = z2d(:,:)
370	znam = 'tempt_sl1'//'_htc'//TRIM(ADJUSTL(zchar))
371	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
372	e_s(:,:,1,jl) = z2d(:,:)
373	DO jk = 1, nlay_i
374	WRITE(zchar1,'(I2)') jk
375	znam = 'tempt'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
376	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
377	e_i(:,:,jk,jl) = z2d(:,:)
378	END DO
379	END DO
380
381	CALL iom_get( numrir, jpdom_autoglo, 'u_ice' , u_ice )
382	CALL iom_get( numrir, jpdom_autoglo, 'v_ice' , v_ice )
383	CALL iom_get( numrir, jpdom_autoglo, 'stress1_i' , stress1_i )
384	CALL iom_get( numrir, jpdom_autoglo, 'stress2_i' , stress2_i )
385	CALL iom_get( numrir, jpdom_autoglo, 'stress12_i', stress12_i )
386	CALL iom_get( numrir, jpdom_autoglo, 'snwice_mass' , snwice_mass )
387	CALL iom_get( numrir, jpdom_autoglo, 'snwice_mass_b', snwice_mass_b )
388
389	DO jl = 1, jpl
390	WRITE(zchar,'(I2)') jl
391	znam = 'sxice'//'_htc'//TRIM(ADJUSTL(zchar))
392	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
393	sxice(:,:,jl) = z2d(:,:)
394	znam = 'syice'//'_htc'//TRIM(ADJUSTL(zchar))
395	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
396	syice(:,:,jl) = z2d(:,:)
397	znam = 'sxxice'//'_htc'//TRIM(ADJUSTL(zchar))
398	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
399	sxxice(:,:,jl) = z2d(:,:)
400	znam = 'syyice'//'_htc'//TRIM(ADJUSTL(zchar))
401	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
402	syyice(:,:,jl) = z2d(:,:)
403	znam = 'sxyice'//'_htc'//TRIM(ADJUSTL(zchar))
404	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
405	sxyice(:,:,jl) = z2d(:,:)
406	znam = 'sxsn'//'_htc'//TRIM(ADJUSTL(zchar))
407	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
408	sxsn(:,:,jl) = z2d(:,:)
409	znam = 'sysn'//'_htc'//TRIM(ADJUSTL(zchar))
410	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
411	sysn(:,:,jl) = z2d(:,:)
412	znam = 'sxxsn'//'_htc'//TRIM(ADJUSTL(zchar))
413	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
414	sxxsn(:,:,jl) = z2d(:,:)
415	znam = 'syysn'//'_htc'//TRIM(ADJUSTL(zchar))
416	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
417	syysn(:,:,jl) = z2d(:,:)
418	znam = 'sxysn'//'_htc'//TRIM(ADJUSTL(zchar))
419	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
420	sxysn(:,:,jl) = z2d(:,:)
421	znam = 'sxa'//'_htc'//TRIM(ADJUSTL(zchar))
422	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
423	sxa(:,:,jl) = z2d(:,:)
424	znam = 'sya'//'_htc'//TRIM(ADJUSTL(zchar))
425	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
426	sya(:,:,jl) = z2d(:,:)
427	znam = 'sxxa'//'_htc'//TRIM(ADJUSTL(zchar))
428	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
429	sxxa(:,:,jl) = z2d(:,:)
430	znam = 'syya'//'_htc'//TRIM(ADJUSTL(zchar))
431	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
432	syya(:,:,jl) = z2d(:,:)
433	znam = 'sxya'//'_htc'//TRIM(ADJUSTL(zchar))
434	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
435	sxya(:,:,jl) = z2d(:,:)
436	znam = 'sxc0'//'_htc'//TRIM(ADJUSTL(zchar))
437	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
438	sxc0(:,:,jl) = z2d(:,:)
439	znam = 'syc0'//'_htc'//TRIM(ADJUSTL(zchar))
440	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
441	syc0(:,:,jl) = z2d(:,:)
442	znam = 'sxxc0'//'_htc'//TRIM(ADJUSTL(zchar))
443	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
444	sxxc0(:,:,jl) = z2d(:,:)
445	znam = 'syyc0'//'_htc'//TRIM(ADJUSTL(zchar))
446	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
447	syyc0(:,:,jl) = z2d(:,:)
448	znam = 'sxyc0'//'_htc'//TRIM(ADJUSTL(zchar))
449	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
450	sxyc0(:,:,jl) = z2d(:,:)
451	znam = 'sxsal'//'_htc'//TRIM(ADJUSTL(zchar))
452	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
453	sxsal(:,:,jl) = z2d(:,:)
454	znam = 'sysal'//'_htc'//TRIM(ADJUSTL(zchar))
455	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
456	sysal(:,:,jl) = z2d(:,:)
457	znam = 'sxxsal'//'_htc'//TRIM(ADJUSTL(zchar))
458	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
459	sxxsal(:,:,jl) = z2d(:,:)
460	znam = 'syysal'//'_htc'//TRIM(ADJUSTL(zchar))
461	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
462	syysal(:,:,jl) = z2d(:,:)
463	znam = 'sxysal'//'_htc'//TRIM(ADJUSTL(zchar))
464	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
465	sxysal(:,:,jl) = z2d(:,:)
466	znam = 'sxage'//'_htc'//TRIM(ADJUSTL(zchar))
467	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
468	sxage(:,:,jl) = z2d(:,:)
469	znam = 'syage'//'_htc'//TRIM(ADJUSTL(zchar))
470	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
471	syage(:,:,jl) = z2d(:,:)
472	znam = 'sxxage'//'_htc'//TRIM(ADJUSTL(zchar))
473	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
474	sxxage(:,:,jl) = z2d(:,:)
475	znam = 'syyage'//'_htc'//TRIM(ADJUSTL(zchar))
476	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
477	syyage(:,:,jl) = z2d(:,:)
478	znam = 'sxyage'//'_htc'//TRIM(ADJUSTL(zchar))
479	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
480	sxyage(:,:,jl)= z2d(:,:)
481	END DO
482
483	CALL iom_get( numrir, jpdom_autoglo, 'sxopw ' , sxopw )
484	CALL iom_get( numrir, jpdom_autoglo, 'syopw ' , syopw )
485	CALL iom_get( numrir, jpdom_autoglo, 'sxxopw' , sxxopw )
486	CALL iom_get( numrir, jpdom_autoglo, 'syyopw' , syyopw )
487	CALL iom_get( numrir, jpdom_autoglo, 'sxyopw' , sxyopw )
488
489	DO jl = 1, jpl
490	WRITE(zchar,'(I2)') jl
491	DO jk = 1, nlay_i
492	WRITE(zchar1,'(I2)') jk
493	znam = 'sxe'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
494	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
495	sxe(:,:,jk,jl) = z2d(:,:)
496	znam = 'sye'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
497	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
498	sye(:,:,jk,jl) = z2d(:,:)
499	znam = 'sxxe'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
500	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
501	sxxe(:,:,jk,jl) = z2d(:,:)
502	znam = 'syye'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
503	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
504	syye(:,:,jk,jl) = z2d(:,:)
505	znam = 'sxye'//'_il'//TRIM(ADJUSTL(zchar1))//'_htc'//TRIM(ADJUSTL(zchar))
506	CALL iom_get( numrir, jpdom_autoglo, znam , z2d )
507	sxye(:,:,jk,jl) = z2d(:,:)
508	END DO
509	END DO
510	!
511	! clem: I do not understand why the following IF is needed
512	! I suspect something inconsistent in the main code with option nn_icesal=1
513	IF( nn_icesal == 1 ) THEN
514	DO jl = 1, jpl
515	sm_i(:,:,jl) = rn_icesal
516	DO jk = 1, nlay_i
517	s_i(:,:,jk,jl) = rn_icesal
518	END DO
519	END DO
520	ENDIF
521	!
522	!CALL iom_close( numrir ) !clem: closed in sbcice_lim.F90
523	!
524	CALL wrk_dealloc( jpi, jpj, z2d )
525	!
526	END SUBROUTINE lim_rst_read
527
528	#else
529	!!----------------------------------------------------------------------
530	!! Default option : Empty module NO LIM sea-ice model
531	!!----------------------------------------------------------------------
532	CONTAINS
533	SUBROUTINE lim_rst_read ! Empty routine
534	END SUBROUTINE lim_rst_read
535	SUBROUTINE lim_rst_write ! Empty routine
536	END SUBROUTINE lim_rst_write
537	#endif
538
539	!!======================================================================
540	END MODULE limrst

Note: See TracBrowser for help on using the repository browser.

Download in other formats: