- Timestamp:
- 2020-12-03T14:58:30+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep
- Property svn:externals
-
old new 8 8 9 9 # SETTE 10 ^/utils/CI/sette @13559sette10 ^/utils/CI/sette_wave@13990 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/TOP/PISCES/P4Z/p4zsms.F90
r13987 r14055 368 368 IF(lwp) WRITE(numout,*) '~~~~~~~' 369 369 ENDIF 370 CALL iom_rstput( kt, nitrst, numrtw, 'PH', hi(:,:,:) )371 CALL iom_rstput( kt, nitrst, numrtw, 'Silicalim', xksi(:,:) )370 CALL iom_rstput( kt, nitrst, numrtw, 'PH', hi(:,:,:) ) 371 CALL iom_rstput( kt, nitrst, numrtw, 'Silicalim', xksi(:,:) ) 372 372 CALL iom_rstput( kt, nitrst, numrtw, 'Silicamax', xksimax(:,:) ) 373 373 CALL iom_rstput( kt, nitrst, numrtw, 'tcflxcum', t_oce_co2_flx_cum ) -
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/TOP/PISCES/SED/sed.F90
r13987 r14055 43 43 REAL , PUBLIC :: sedmask 44 44 REAL(wp), PUBLIC :: denssol !: density of solid material 45 INTEGER , PUBLIC :: numrsr, numrsw !: logical unit for sed restart (read and write)46 45 LOGICAL , PUBLIC :: lrst_sed !: logical to control the trc restart write 47 46 LOGICAL , PUBLIC :: ln_rst_sed = .TRUE. !: initialisation from a restart file or not -
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/TOP/PISCES/SED/sedrst.F90
r13286 r14055 42 42 CHARACTER(LEN=50) :: clname ! trc output restart file name 43 43 CHARACTER(LEN=256) :: clpath ! full path to ocean output restart file 44 CHARACTER(LEN=52) :: clpname ! trc output restart file name including AGRIF 44 45 !!---------------------------------------------------------------------- 45 46 ! … … 80 81 IF(lwp) WRITE(numsed,*) & 81 82 ' open sed restart.output NetCDF file: ',TRIM(clpath)//clname 82 CALL iom_open( TRIM(clpath)//TRIM(clname), numrsw, ldwrt = .TRUE., kdlev = jpksed, cdcomp = 'SED' ) 83 IF(.NOT.lwxios) THEN 84 CALL iom_open( TRIM(clpath)//TRIM(clname), numrsw, ldwrt = .TRUE., kdlev = jpksed, cdcomp = 'SED' ) 85 ELSE 86 #if defined key_iomput 87 cw_sedrst_cxt = "rstws_"//TRIM(ADJUSTL(clkt)) 88 IF( TRIM(Agrif_CFixed()) == '0' ) THEN 89 clpname = clname 90 ELSE 91 clpname = TRIM(Agrif_CFixed())//"_"//clname 92 ENDIF 93 numrsw = iom_xios_setid(TRIM(clpath)//TRIM(clpname)) 94 CALL iom_init( cw_sedrst_cxt, kdid = numrsw, ld_closedef = .FALSE. ) 95 #else 96 CALL ctl_stop( 'Can not use XIOS in trc_rst_opn' ) 97 #endif 98 ENDIF 99 83 100 lrst_sed = .TRUE. 84 101 ENDIF … … 196 213 CALL pack_arr( jpoce, sedligand(1:jpoce,1:jpksed), & 197 214 & zdta2(1:jpi,1:jpj,1:jpksed), iarroce(1:jpoce) ) 198 199 215 IF( ln_timing ) CALL timing_stop('sed_rst_read') 200 216 … … 240 256 !! 1. WRITE in nutwrs 241 257 !! ------------------ 242 243 zinfo(1) = REAL( kt) 244 CALL iom_rstput( kt, nitrst, numrsw, 'kt', zinfo ) 258 ! zinfo(1) = REAL( kt) 259 CALL iom_rstput( kt, nitrst, numrsw, 'kt', REAL( kt , wp) ) 245 260 246 261 ! Back to 2D geometry … … 299 314 300 315 IF( kt == nitrst ) THEN 301 CALL iom_close( numrsw ) ! close the restart file (only at last time step) 316 IF(.NOT.lwxios) THEN 317 CALL iom_close( numrsw ) ! close the restart file (only at last time step) 318 ELSE 319 CALL iom_context_finalize( cw_sedrst_cxt ) 320 iom_file(numrsw)%nfid = 0 321 numrsw = 0 322 ENDIF 302 323 IF( l_offline .AND. ln_rst_list ) THEN 303 324 nrst_lst = nrst_lst + 1 … … 342 363 REAL(wp) :: zkt, zrdttrc1 343 364 REAL(wp) :: zndastp 365 CHARACTER(len = 82) :: clpname 344 366 345 367 ! Time domain : restart … … 353 375 354 376 IF( ln_rst_sed ) THEN 377 lxios_sini = .FALSE. 355 378 CALL iom_open( TRIM(cn_sedrst_indir)//'/'//cn_sedrst_in, numrsr ) 379 380 IF( lrxios) THEN 381 cr_sedrst_cxt = 'sed_rst' 382 IF(lwp) WRITE(numout,*) 'Enable restart reading by XIOS for SED' 383 ! IF( TRIM(Agrif_CFixed()) == '0' ) THEN 384 ! clpname = cn_sedrst_in 385 ! ELSE 386 ! clpname = TRIM(Agrif_CFixed())//"_"//cn_sedrst_in 387 ! ENDIF 388 CALL iom_init( cr_sedrst_cxt, kdid = numrsr, ld_closedef = .TRUE. ) 389 ENDIF 356 390 CALL iom_get ( numrsr, 'kt', zkt ) ! last time-step of previous run 357 358 391 IF(lwp) THEN 359 392 WRITE(numsed,*) ' *** Info read in restart : ' … … 402 435 IF(lwp) WRITE(numsed,*) 'trc_wri : write the TOP restart file (NetCDF) at it= ', kt, ' date= ', ndastp 403 436 IF(lwp) WRITE(numsed,*) '~~~~~~~' 437 IF( lwxios ) CALL iom_init_closedef(cw_sedrst_cxt) 404 438 ENDIF 405 439 CALL iom_rstput( kt, nitrst, numrsw, 'kt' , REAL( kt , wp) ) ! time-step 406 440 CALL iom_rstput( kt, nitrst, numrsw, 'ndastp' , REAL( ndastp, wp) ) ! date 407 CALL iom_rstput( kt, nitrst, numrsw, 'adatrj' , adatrj 408 ! ! the begining of the run [s]441 CALL iom_rstput( kt, nitrst, numrsw, 'adatrj' , adatrj ) ! number of elapsed days since 442 ! ! the begining of the run [s] 409 443 ENDIF 410 444 -
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/TOP/PISCES/SED/sedstp.F90
r12489 r14055 86 86 IF( kt == nitsed000 ) THEN 87 87 CALL iom_close( numrsr ) ! close input tracer restart file 88 ! IF(lwm) CALL FLUSH( numont ) ! flush namelist output 88 IF(lrxios) CALL iom_context_finalize( cr_sedrst_cxt ) 89 ! IF(lwm) CALL FLUSH( numont ) ! flush namelist output 89 90 ENDIF 90 91 IF( lrst_sed ) CALL sed_rst_wri( kt ) ! restart file output
Note: See TracChangeset
for help on using the changeset viewer.