Line data Source code
1 : module restart_physics
2 :
3 : use shr_kind_mod, only: r8 => shr_kind_r8
4 : use spmd_utils, only: masterproc
5 : use co2_cycle, only: co2_transport
6 : use constituents, only: pcnst
7 :
8 : use radiation, only: radiation_define_restart, radiation_write_restart, &
9 : radiation_read_restart
10 :
11 : use ioFileMod
12 : use cam_abortutils, only: endrun
13 : use camsrfexch, only: cam_in_t, cam_out_t
14 : use cam_logfile, only: iulog
15 : use pio, only: file_desc_t, io_desc_t, var_desc_t, &
16 : pio_double, pio_int, pio_noerr, &
17 : pio_seterrorhandling, pio_bcast_error, &
18 : pio_inq_varid, &
19 : pio_def_var, pio_def_dim, &
20 : pio_put_var, pio_get_var
21 :
22 : implicit none
23 : private
24 : save
25 : !
26 : ! Public interfaces
27 : !
28 : public :: write_restart_physics ! Write the physics restart info out
29 : public :: read_restart_physics ! Read the physics restart info in
30 : public :: init_restart_physics
31 :
32 : !
33 : ! Private data
34 : !
35 :
36 : type(var_desc_t) :: flwds_desc, &
37 : solld_desc, co2prog_desc, co2diag_desc, sols_desc, soll_desc, &
38 : solsd_desc
39 :
40 : type(var_desc_t) :: bcphidry_desc, bcphodry_desc, ocphidry_desc, ocphodry_desc, &
41 : dstdry1_desc, dstdry2_desc, dstdry3_desc, dstdry4_desc
42 :
43 : type(var_desc_t) :: cflx_desc(pcnst)
44 :
45 : type(var_desc_t) :: wsx_desc
46 : type(var_desc_t) :: wsy_desc
47 : type(var_desc_t) :: shf_desc
48 :
49 : CONTAINS
50 3072 : subroutine init_restart_physics ( File, pbuf2d)
51 :
52 : use physics_buffer, only: pbuf_init_restart, physics_buffer_desc
53 : use ppgrid, only: pver, pverp
54 : use chemistry, only: chem_init_restart
55 : use prescribed_ozone, only: init_prescribed_ozone_restart
56 : use prescribed_ghg, only: init_prescribed_ghg_restart
57 : use prescribed_aero, only: init_prescribed_aero_restart
58 : use prescribed_volcaero, only: init_prescribed_volcaero_restart
59 : use cam_grid_support, only: cam_grid_write_attr, cam_grid_id
60 : use cam_grid_support, only: cam_grid_header_info_t
61 : use cam_pio_utils, only: cam_pio_def_dim
62 : use subcol_utils, only: is_subcol_on
63 : use subcol, only: subcol_init_restart
64 : use carma_intr, only: carma_restart_init
65 :
66 : type(file_desc_t), intent(inout) :: file
67 : type(physics_buffer_desc), pointer :: pbuf2d(:,:)
68 :
69 : integer :: grid_id
70 : integer :: hdimcnt, ierr, i
71 : integer :: dimids(4)
72 2304 : integer, allocatable :: hdimids(:)
73 2304 : type(cam_grid_header_info_t) :: info
74 : character(len=4) :: num
75 :
76 2304 : call pio_seterrorhandling(File, PIO_BCAST_ERROR)
77 : ! Probably should have the grid write this out.
78 2304 : grid_id = cam_grid_id('physgrid')
79 2304 : call cam_grid_write_attr(File, grid_id, info)
80 2304 : hdimcnt = info%num_hdims()
81 :
82 4608 : do i = 1, hdimcnt
83 4608 : dimids(i) = info%get_hdimid(i)
84 : end do
85 6912 : allocate(hdimids(hdimcnt))
86 4608 : hdimids(1:hdimcnt) = dimids(1:hdimcnt)
87 :
88 2304 : call pbuf_init_restart(File, pbuf2d)
89 :
90 2304 : call chem_init_restart(File)
91 :
92 2304 : call init_prescribed_ozone_restart(File)
93 2304 : call init_prescribed_ghg_restart(File)
94 2304 : call init_prescribed_aero_restart(File)
95 2304 : call init_prescribed_volcaero_restart(File)
96 :
97 2304 : ierr = pio_def_var(File, 'FLWDS', pio_double, hdimids, flwds_desc)
98 2304 : ierr = pio_def_var(File, 'SOLS', pio_double, hdimids, sols_desc)
99 2304 : ierr = pio_def_var(File, 'SOLL', pio_double, hdimids, soll_desc)
100 2304 : ierr = pio_def_var(File, 'SOLSD', pio_double, hdimids, solsd_desc)
101 2304 : ierr = pio_def_var(File, 'SOLLD', pio_double, hdimids, solld_desc)
102 :
103 2304 : ierr = pio_def_var(File, 'BCPHIDRY', pio_double, hdimids, bcphidry_desc)
104 2304 : ierr = pio_def_var(File, 'BCPHODRY', pio_double, hdimids, bcphodry_desc)
105 2304 : ierr = pio_def_var(File, 'OCPHIDRY', pio_double, hdimids, ocphidry_desc)
106 2304 : ierr = pio_def_var(File, 'OCPHODRY', pio_double, hdimids, ocphodry_desc)
107 2304 : ierr = pio_def_var(File, 'DSTDRY1', pio_double, hdimids, dstdry1_desc)
108 2304 : ierr = pio_def_var(File, 'DSTDRY2', pio_double, hdimids, dstdry2_desc)
109 2304 : ierr = pio_def_var(File, 'DSTDRY3', pio_double, hdimids, dstdry3_desc)
110 2304 : ierr = pio_def_var(File, 'DSTDRY4', pio_double, hdimids, dstdry4_desc)
111 :
112 2304 : if(co2_transport()) then
113 0 : ierr = pio_def_var(File, 'CO2PROG', pio_double, hdimids, co2prog_desc)
114 0 : ierr = pio_def_var(File, 'CO2DIAG', pio_double, hdimids, co2diag_desc)
115 : end if
116 :
117 : ! cam_import variables -- write the constituent surface fluxes as individual 2D arrays
118 : ! rather than as a single variable with a pcnst dimension. Note that the cflx components
119 : ! are only needed for those constituents that are not passed to the coupler. The restart
120 : ! for constituents passed through the coupler are handled by the .rs. restart file. But
121 : ! we don't currently have a mechanism to know whether the constituent is handled by the
122 : ! coupler or not, so we write all of cflx to the CAM restart file.
123 96768 : do i = 1, pcnst
124 94464 : write(num,'(i4.4)') i
125 96768 : ierr = pio_def_var(File, 'CFLX'//num, pio_double, hdimids, cflx_desc(i))
126 : end do
127 :
128 2304 : ierr = pio_def_var(File, 'wsx', pio_double, hdimids, wsx_desc)
129 2304 : ierr = pio_def_var(File, 'wsy', pio_double, hdimids, wsy_desc)
130 2304 : ierr = pio_def_var(File, 'shf', pio_double, hdimids, shf_desc)
131 :
132 2304 : call radiation_define_restart(file)
133 :
134 2304 : if (is_subcol_on()) then
135 0 : call subcol_init_restart(file, hdimids)
136 : end if
137 :
138 2304 : call carma_restart_init(file)
139 :
140 2304 : end subroutine init_restart_physics
141 :
142 4608 : subroutine write_restart_physics (File, cam_in, cam_out, pbuf2d)
143 :
144 : !-----------------------------------------------------------------------
145 2304 : use physics_buffer, only: physics_buffer_desc, pbuf_write_restart
146 : use phys_grid, only: phys_decomp
147 :
148 : use ppgrid, only: begchunk, endchunk, pcols, pverp
149 : use chemistry, only: chem_write_restart
150 : use prescribed_ozone, only: write_prescribed_ozone_restart
151 : use prescribed_ghg, only: write_prescribed_ghg_restart
152 : use prescribed_aero, only: write_prescribed_aero_restart
153 : use prescribed_volcaero, only: write_prescribed_volcaero_restart
154 :
155 : use cam_history_support, only: fillvalue
156 : use spmd_utils, only: iam
157 : use cam_grid_support, only: cam_grid_write_dist_array, cam_grid_id
158 : use cam_grid_support, only: cam_grid_get_decomp, cam_grid_dimensions
159 : use cam_grid_support, only: cam_grid_write_var
160 : use pio, only: pio_write_darray
161 : use subcol_utils, only: is_subcol_on
162 : use subcol, only: subcol_write_restart
163 : use carma_intr, only: carma_restart_write
164 : !
165 : ! Input arguments
166 : !
167 : type(file_desc_t), intent(inout) :: File
168 : type(cam_in_t), intent(in) :: cam_in(begchunk:endchunk)
169 : type(cam_out_t), intent(in) :: cam_out(begchunk:endchunk)
170 : type(physics_buffer_desc), pointer :: pbuf2d(:,:)
171 : !
172 : ! Local workspace
173 : !
174 : type(io_desc_t), pointer :: iodesc
175 4608 : real(r8):: tmpfield(pcols, begchunk:endchunk)
176 : integer :: i, m ! loop index
177 : integer :: ncol ! number of vertical columns
178 : integer :: ierr
179 : integer :: physgrid
180 : integer :: dims(3), gdims(3)
181 : integer :: nhdims
182 : !-----------------------------------------------------------------------
183 :
184 : ! Write grid vars
185 2304 : call cam_grid_write_var(File, phys_decomp)
186 :
187 : ! Physics buffer
188 2304 : if (is_subcol_on()) then
189 0 : call subcol_write_restart(File)
190 : end if
191 :
192 2304 : call pbuf_write_restart(File, pbuf2d)
193 :
194 2304 : physgrid = cam_grid_id('physgrid')
195 2304 : call cam_grid_dimensions(physgrid, gdims(1:2), nhdims)
196 :
197 : ! data for chemistry
198 2304 : call chem_write_restart(File)
199 :
200 2304 : call write_prescribed_ozone_restart(File)
201 2304 : call write_prescribed_ghg_restart(File)
202 2304 : call write_prescribed_aero_restart(File)
203 2304 : call write_prescribed_volcaero_restart(File)
204 :
205 : ! cam_in/out variables
206 : ! This is a group of surface variables so can reuse dims
207 2304 : dims(1) = pcols
208 2304 : dims(2) = endchunk - begchunk + 1
209 0 : call cam_grid_get_decomp(physgrid, dims(1:2), gdims(1:nhdims), &
210 2304 : pio_double, iodesc)
211 :
212 11592 : do i = begchunk, endchunk
213 9288 : ncol = cam_out(i)%ncol
214 155088 : tmpfield(:ncol, i) = cam_out(i)%flwds(:ncol)
215 : ! Only have to do this once (cam_in/out vars all same shape)
216 11592 : if (ncol < pcols) then
217 5112 : tmpfield(ncol+1:, i) = fillvalue
218 : end if
219 : end do
220 2304 : call pio_write_darray(File, flwds_desc, iodesc, tmpfield, ierr)
221 :
222 11592 : do i = begchunk, endchunk
223 9288 : ncol = cam_out(i)%ncol
224 157392 : tmpfield(:ncol, i) = cam_out(i)%sols(:ncol)
225 : end do
226 2304 : call pio_write_darray(File, sols_desc, iodesc, tmpfield, ierr)
227 :
228 11592 : do i = begchunk, endchunk
229 9288 : ncol = cam_out(i)%ncol
230 157392 : tmpfield(:ncol, i) = cam_out(i)%soll(:ncol)
231 : end do
232 2304 : call pio_write_darray(File, soll_desc, iodesc, tmpfield, ierr)
233 :
234 11592 : do i = begchunk, endchunk
235 9288 : ncol = cam_out(i)%ncol
236 157392 : tmpfield(:ncol, i) = cam_out(i)%solsd(:ncol)
237 : end do
238 2304 : call pio_write_darray(File, solsd_desc, iodesc, tmpfield, ierr)
239 :
240 11592 : do i = begchunk, endchunk
241 9288 : ncol = cam_out(i)%ncol
242 157392 : tmpfield(:ncol, i) = cam_out(i)%solld(:ncol)
243 : end do
244 2304 : call pio_write_darray(File, solld_desc, iodesc, tmpfield, ierr)
245 :
246 11592 : do i = begchunk, endchunk
247 9288 : ncol = cam_out(i)%ncol
248 157392 : tmpfield(:ncol, i) = cam_out(i)%bcphidry(:ncol)
249 : end do
250 2304 : call pio_write_darray(File, bcphidry_desc, iodesc, tmpfield, ierr)
251 :
252 11592 : do i = begchunk, endchunk
253 9288 : ncol = cam_out(i)%ncol
254 157392 : tmpfield(:ncol, i) = cam_out(i)%bcphodry(:ncol)
255 : end do
256 2304 : call pio_write_darray(File, bcphodry_desc, iodesc, tmpfield, ierr)
257 :
258 11592 : do i = begchunk, endchunk
259 9288 : ncol = cam_out(i)%ncol
260 157392 : tmpfield(:ncol, i) = cam_out(i)%ocphidry(:ncol)
261 : end do
262 2304 : call pio_write_darray(File, ocphidry_desc, iodesc, tmpfield, ierr)
263 :
264 11592 : do i = begchunk, endchunk
265 9288 : ncol = cam_out(i)%ncol
266 157392 : tmpfield(:ncol, i) = cam_out(i)%ocphodry(:ncol)
267 : end do
268 2304 : call pio_write_darray(File, ocphodry_desc, iodesc, tmpfield, ierr)
269 :
270 11592 : do i = begchunk, endchunk
271 9288 : ncol = cam_out(i)%ncol
272 157392 : tmpfield(:ncol, i) = cam_out(i)%dstdry1(:ncol)
273 : end do
274 2304 : call pio_write_darray(File, dstdry1_desc, iodesc, tmpfield, ierr)
275 :
276 11592 : do i = begchunk, endchunk
277 9288 : ncol = cam_out(i)%ncol
278 157392 : tmpfield(:ncol, i) = cam_out(i)%dstdry2(:ncol)
279 : end do
280 2304 : call pio_write_darray(File, dstdry2_desc, iodesc, tmpfield, ierr)
281 :
282 11592 : do i = begchunk, endchunk
283 9288 : ncol = cam_out(i)%ncol
284 157392 : tmpfield(:ncol, i) = cam_out(i)%dstdry3(:ncol)
285 : end do
286 2304 : call pio_write_darray(File, dstdry3_desc, iodesc, tmpfield, ierr)
287 :
288 11592 : do i = begchunk, endchunk
289 9288 : ncol = cam_out(i)%ncol
290 157392 : tmpfield(:ncol, i) = cam_out(i)%dstdry4(:ncol)
291 : end do
292 2304 : call pio_write_darray(File, dstdry4_desc, iodesc, tmpfield, ierr)
293 :
294 2304 : if (co2_transport()) then
295 0 : do i = begchunk, endchunk
296 0 : ncol = cam_out(i)%ncol
297 0 : tmpfield(:ncol, i) = cam_out(i)%co2prog(:ncol)
298 : end do
299 0 : call pio_write_darray(File, co2prog_desc, iodesc, tmpfield, ierr)
300 :
301 0 : do i = begchunk, endchunk
302 0 : ncol = cam_out(i)%ncol
303 0 : tmpfield(:ncol, i) = cam_out(i)%co2diag(:ncol)
304 : end do
305 0 : call pio_write_darray(File, co2diag_desc, iodesc, tmpfield, ierr)
306 : end if
307 :
308 : ! cam_in components
309 96768 : do m = 1, pcnst
310 475272 : do i = begchunk, endchunk
311 380808 : ncol = cam_in(i)%ncol
312 6453072 : tmpfield(:ncol, i) = cam_in(i)%cflx(:ncol, m)
313 : end do
314 96768 : call pio_write_darray(File, cflx_desc(m), iodesc, tmpfield, ierr)
315 : end do
316 :
317 11592 : do i = begchunk, endchunk
318 9288 : ncol = cam_in(i)%ncol
319 157392 : tmpfield(:ncol,i) = cam_in(i)%wsx(:ncol)
320 : end do
321 2304 : call pio_write_darray(File, wsx_desc, iodesc, tmpfield, ierr)
322 :
323 11592 : do i = begchunk, endchunk
324 9288 : ncol = cam_in(i)%ncol
325 157392 : tmpfield(:ncol,i) = cam_in(i)%wsy(:ncol)
326 : end do
327 2304 : call pio_write_darray(File, wsy_desc, iodesc, tmpfield, ierr)
328 :
329 11592 : do i = begchunk, endchunk
330 9288 : ncol = cam_in(i)%ncol
331 157392 : tmpfield(:ncol,i) = cam_in(i)%shf(:ncol)
332 : end do
333 2304 : call pio_write_darray(File, shf_desc, iodesc, tmpfield, ierr)
334 :
335 2304 : call radiation_write_restart(file)
336 2304 : call carma_restart_write(file)
337 :
338 2304 : end subroutine write_restart_physics
339 :
340 : !#######################################################################
341 :
342 768 : subroutine read_restart_physics(File, cam_in, cam_out, pbuf2d)
343 :
344 : !-----------------------------------------------------------------------
345 2304 : use physics_buffer, only: physics_buffer_desc, pbuf_read_restart
346 :
347 : use ppgrid, only: begchunk, endchunk, pcols, pver, pverp
348 : use chemistry, only: chem_read_restart
349 : use cam_grid_support, only: cam_grid_read_dist_array, cam_grid_id
350 : use cam_grid_support, only: cam_grid_get_decomp, cam_grid_dimensions
351 : use cam_history_support, only: fillvalue
352 :
353 : use prescribed_ozone, only: read_prescribed_ozone_restart
354 : use prescribed_ghg, only: read_prescribed_ghg_restart
355 : use prescribed_aero, only: read_prescribed_aero_restart
356 : use prescribed_volcaero, only: read_prescribed_volcaero_restart
357 : use subcol_utils, only: is_subcol_on
358 : use subcol, only: subcol_read_restart
359 : use pio, only: pio_read_darray
360 : use carma_intr, only: carma_restart_read
361 : !
362 : ! Arguments
363 : !
364 : type(file_desc_t), intent(inout) :: File
365 : type(cam_in_t), pointer :: cam_in(:)
366 : type(cam_out_t), pointer :: cam_out(:)
367 : type(physics_buffer_desc), pointer :: pbuf2d(:,:)
368 : !
369 : ! Local workspace
370 : !
371 768 : real(r8), allocatable :: tmpfield2(:,:)
372 : integer :: i, c, m ! loop index
373 : integer :: ierr ! I/O status
374 : type(io_desc_t), pointer :: iodesc
375 : type(var_desc_t) :: vardesc
376 : integer :: csize, vsize
377 : character(len=4) :: num
378 : integer :: dims(3), gdims(3), nhdims
379 : integer :: err_handling
380 : integer :: physgrid
381 : !-----------------------------------------------------------------------
382 :
383 : ! subcol_read_restart must be called before pbuf_read_restart
384 768 : if (is_subcol_on()) then
385 0 : call subcol_read_restart(File)
386 : end if
387 :
388 768 : call pbuf_read_restart(File, pbuf2d)
389 :
390 768 : csize=endchunk-begchunk+1
391 768 : dims(1) = pcols
392 768 : dims(2) = csize
393 :
394 768 : physgrid = cam_grid_id('physgrid')
395 :
396 768 : call cam_grid_dimensions(physgrid, gdims(1:2))
397 :
398 768 : if (gdims(2) == 1) then
399 : nhdims = 1
400 : else
401 0 : nhdims = 2
402 : end if
403 0 : call cam_grid_get_decomp(physgrid, dims(1:2), gdims(1:nhdims), pio_double, &
404 768 : iodesc)
405 :
406 : ! data for chemistry
407 768 : call chem_read_restart(File)
408 :
409 768 : call read_prescribed_ozone_restart(File)
410 768 : call read_prescribed_ghg_restart(File)
411 768 : call read_prescribed_aero_restart(File)
412 768 : call read_prescribed_volcaero_restart(File)
413 :
414 2304 : allocate(tmpfield2(pcols, begchunk:endchunk))
415 53400 : tmpfield2 = fillvalue
416 :
417 768 : ierr = pio_inq_varid(File, 'FLWDS', vardesc)
418 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
419 3864 : do c=begchunk,endchunk
420 53400 : do i=1,pcols
421 52632 : cam_out(c)%flwds(i) = tmpfield2(i, c)
422 : end do
423 : end do
424 :
425 768 : ierr = pio_inq_varid(File, 'SOLS', vardesc)
426 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
427 3864 : do c=begchunk,endchunk
428 53400 : do i=1,pcols
429 52632 : cam_out(c)%sols(i) = tmpfield2(i, c)
430 : end do
431 : end do
432 :
433 768 : ierr = pio_inq_varid(File, 'SOLL', vardesc)
434 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
435 3864 : do c=begchunk,endchunk
436 53400 : do i=1,pcols
437 52632 : cam_out(c)%soll(i) = tmpfield2(i, c)
438 : end do
439 : end do
440 :
441 768 : ierr = pio_inq_varid(File, 'SOLSD', vardesc)
442 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
443 3864 : do c=begchunk,endchunk
444 53400 : do i=1,pcols
445 52632 : cam_out(c)%solsd(i) = tmpfield2(i, c)
446 : end do
447 : end do
448 :
449 768 : ierr = pio_inq_varid(File, 'SOLLD', vardesc)
450 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
451 3864 : do c=begchunk,endchunk
452 53400 : do i=1,pcols
453 52632 : cam_out(c)%solld(i) = tmpfield2(i, c)
454 : end do
455 : end do
456 :
457 768 : ierr = pio_inq_varid(File, 'BCPHIDRY', vardesc)
458 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
459 3864 : do c=begchunk,endchunk
460 53400 : do i=1,pcols
461 52632 : cam_out(c)%bcphidry(i) = tmpfield2(i, c)
462 : end do
463 : end do
464 :
465 768 : ierr = pio_inq_varid(File, 'BCPHODRY', vardesc)
466 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
467 3864 : do c=begchunk,endchunk
468 53400 : do i=1,pcols
469 52632 : cam_out(c)%bcphodry(i) = tmpfield2(i, c)
470 : end do
471 : end do
472 :
473 768 : ierr = pio_inq_varid(File, 'OCPHIDRY', vardesc)
474 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
475 3864 : do c=begchunk,endchunk
476 53400 : do i=1,pcols
477 52632 : cam_out(c)%ocphidry(i) = tmpfield2(i, c)
478 : end do
479 : end do
480 :
481 768 : ierr = pio_inq_varid(File, 'OCPHODRY', vardesc)
482 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
483 3864 : do c=begchunk,endchunk
484 53400 : do i=1,pcols
485 52632 : cam_out(c)%ocphodry(i) = tmpfield2(i, c)
486 : end do
487 : end do
488 :
489 768 : ierr = pio_inq_varid(File, 'DSTDRY1', vardesc)
490 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
491 3864 : do c=begchunk,endchunk
492 53400 : do i=1,pcols
493 52632 : cam_out(c)%dstdry1(i) = tmpfield2(i, c)
494 : end do
495 : end do
496 :
497 768 : ierr = pio_inq_varid(File, 'DSTDRY2', vardesc)
498 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
499 3864 : do c=begchunk,endchunk
500 53400 : do i=1,pcols
501 52632 : cam_out(c)%dstdry2(i) = tmpfield2(i, c)
502 : end do
503 : end do
504 :
505 768 : ierr = pio_inq_varid(File, 'DSTDRY3', vardesc)
506 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
507 3864 : do c=begchunk,endchunk
508 53400 : do i=1,pcols
509 52632 : cam_out(c)%dstdry3(i) = tmpfield2(i, c)
510 : end do
511 : end do
512 :
513 768 : ierr = pio_inq_varid(File, 'DSTDRY4', vardesc)
514 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
515 3864 : do c=begchunk,endchunk
516 53400 : do i=1,pcols
517 52632 : cam_out(c)%dstdry4(i) = tmpfield2(i, c)
518 : end do
519 : end do
520 :
521 768 : if (co2_transport()) then
522 0 : ierr = pio_inq_varid(File, 'CO2PROG', vardesc)
523 0 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
524 0 : do c=begchunk,endchunk
525 0 : do i=1,pcols
526 0 : cam_out(c)%co2prog(i) = tmpfield2(i, c)
527 : end do
528 : end do
529 :
530 0 : ierr = pio_inq_varid(File, 'CO2DIAG', vardesc)
531 0 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
532 0 : do c=begchunk,endchunk
533 0 : do i=1,pcols
534 0 : cam_out(c)%co2diag(i) = tmpfield2(i, c)
535 : end do
536 : end do
537 : end if
538 :
539 : ! Reading the CFLX* components from the restart is optional for
540 : ! backwards compatibility. These fields were not needed for an
541 : ! exact restart until the UNICON scheme was added. More generally,
542 : ! these components are only needed if they are not handled by the
543 : ! coupling layer restart (the ".rs." file), and if the values are
544 : ! used in the tphysbc physics before the tphysac code has a chance
545 : ! to update the values that are coming from boundary datasets.
546 32256 : do m = 1, pcnst
547 :
548 31488 : write(num,'(i4.4)') m
549 :
550 31488 : call pio_seterrorhandling(File, PIO_BCAST_ERROR, err_handling)
551 31488 : ierr = pio_inq_varid(File, 'CFLX'//num, vardesc)
552 31488 : call pio_seterrorhandling(File, err_handling)
553 :
554 32256 : if (ierr == PIO_NOERR) then ! CFLX variable found on restart file
555 31488 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
556 158424 : do c= begchunk, endchunk
557 2189400 : do i = 1, pcols
558 2157912 : cam_in(c)%cflx(i,m) = tmpfield2(i, c)
559 : end do
560 : end do
561 : end if
562 :
563 : end do
564 :
565 768 : call pio_seterrorhandling(File, PIO_BCAST_ERROR, err_handling)
566 768 : ierr = pio_inq_varid(File, 'wsx', vardesc)
567 768 : if (ierr == PIO_NOERR) then ! variable found on restart file
568 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
569 3864 : do c= begchunk, endchunk
570 53400 : do i = 1, pcols
571 52632 : cam_in(c)%wsx(i) = tmpfield2(i, c)
572 : end do
573 : end do
574 : end if
575 768 : ierr = pio_inq_varid(File, 'wsy', vardesc)
576 768 : if (ierr == PIO_NOERR) then ! variable found on restart file
577 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
578 3864 : do c= begchunk, endchunk
579 53400 : do i = 1, pcols
580 52632 : cam_in(c)%wsy(i) = tmpfield2(i, c)
581 : end do
582 : end do
583 : end if
584 768 : ierr = pio_inq_varid(File, 'shf', vardesc)
585 768 : if (ierr == PIO_NOERR) then ! variable found on restart file
586 768 : call pio_read_darray(File, vardesc, iodesc, tmpfield2, ierr)
587 3864 : do c= begchunk, endchunk
588 53400 : do i = 1, pcols
589 52632 : cam_in(c)%shf(i) = tmpfield2(i, c)
590 : end do
591 : end do
592 : endif
593 768 : call pio_seterrorhandling(File, err_handling)
594 :
595 768 : deallocate(tmpfield2)
596 :
597 768 : call radiation_read_restart(file)
598 768 : call carma_restart_read(file)
599 :
600 10752 : end subroutine read_restart_physics
601 :
602 : end module restart_physics
|
