MODULE trcrst !!====================================================================== !! *** MODULE trcrst *** !! TOP : Manage the passive tracer restart !!====================================================================== !! History : - ! 1991-03 () original code !! 1.0 ! 2005-03 (O. Aumont, A. El Moussaoui) F90 !! - ! 2005-10 (C. Ethe) print control !! 2.0 ! 2005-10 (C. Ethe, G. Madec) revised architecture !!---------------------------------------------------------------------- #if defined key_top !!---------------------------------------------------------------------- !! 'key_top' TOP models !!---------------------------------------------------------------------- !!---------------------------------------------------------------------- !! trc_rst : Restart for passive tracer !!---------------------------------------------------------------------- !!---------------------------------------------------------------------- !! 'key_top' TOP models !!---------------------------------------------------------------------- !! trc_rst_opn : open restart file !! trc_rst_read : read restart file !! trc_rst_wri : write restart file !!---------------------------------------------------------------------- USE oce_trc USE trc USE trcnam_trp USE iom USE trcrst_cfc ! CFC USE trcrst_lobster ! LOBSTER restart USE trcrst_pisces ! PISCES restart USE trcrst_c14b ! C14 bomb restart USE trcrst_my_trc ! MY_TRC restart USE daymod IMPLICIT NONE PRIVATE PUBLIC trc_rst_opn ! called by ??? PUBLIC trc_rst_read ! called by ??? PUBLIC trc_rst_wri ! called by ??? INTEGER, PUBLIC :: numrtr, numrtw !: logical unit for trc restart (read and write) !! * Substitutions # include "top_substitute.h90" CONTAINS SUBROUTINE trc_rst_opn( kt ) !!---------------------------------------------------------------------- !! *** trc_rst_opn *** !! !! ** purpose : output of sea-trc variable in a netcdf file !!---------------------------------------------------------------------- INTEGER, INTENT(in) :: kt ! number of iteration ! CHARACTER(LEN=20) :: clkt ! ocean time-step define as a character CHARACTER(LEN=50) :: clname ! trc output restart file name !!---------------------------------------------------------------------- ! IF( lk_offline ) THEN IF( kt == nit000 ) THEN lrst_trc = .FALSE. nitrst = nitend ENDIF IF( MOD( kt - 1, nstock ) == 0 ) THEN ! we use kt - 1 and not kt - nit000 to keep the same periodicity from the beginning of the experiment nitrst = kt + nstock - 1 ! define the next value of nitrst for restart writing IF( nitrst > nitend ) nitrst = nitend ! make sure we write a restart at the end of the run ENDIF ELSE IF( kt == nit000 ) lrst_trc = .FALSE. ENDIF ! to get better performances with NetCDF format: ! we open and define the tracer restart file one tracer time step before writing the data (-> at nitrst - 2*nn_dttrc + 1) ! except if we write tracer restart files every tracer time step or if a tracer restart file was writen at nitend - 2*nn_dttrc + 1 IF( kt == nitrst - 2*nn_dttrc + 1 .OR. nstock == nn_dttrc .OR. ( kt == nitend - nn_dttrc + 1 .AND. .NOT. lrst_trc ) ) THEN ! beware of the format used to write kt (default is i8.8, that should be large enough) IF( nitrst > 1.0e9 ) THEN ; WRITE(clkt,* ) nitrst ELSE ; WRITE(clkt,'(i8.8)') nitrst ENDIF ! create the file IF(lwp) WRITE(numout,*) clname = TRIM(cexper)//"_"//TRIM(ADJUSTL(clkt))//"_"//TRIM(cn_trcrst_out) IF(lwp) WRITE(numout,*) ' open trc restart.output NetCDF file: '//clname CALL iom_open( clname, numrtw, ldwrt = .TRUE., kiolib = jprstlib ) lrst_trc = .TRUE. ENDIF ! END SUBROUTINE trc_rst_opn SUBROUTINE trc_rst_read !!---------------------------------------------------------------------- !! *** trc_rst_opn *** !! !! ** purpose : read passive tracer fields in restart files !!---------------------------------------------------------------------- INTEGER :: jn INTEGER :: jlibalt = jprstlib LOGICAL :: llok !!---------------------------------------------------------------------- IF(lwp) WRITE(numout,*) IF(lwp) WRITE(numout,*) 'trc_rst_read : read the TOP restart file' IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' IF ( jprstlib == jprstdimg ) THEN ! eventually read netcdf file (monobloc) for restarting on different number of processors ! if {cn_trcrst_in}.nc exists, then set jlibalt to jpnf90 INQUIRE( FILE = TRIM(cn_trcrst_in)//'.nc', EXIST = llok ) IF ( llok ) THEN ; jlibalt = jpnf90 ; ELSE ; jlibalt = jprstlib ; ENDIF ENDIF CALL iom_open( cn_trcrst_in, numrtr, kiolib = jlibalt ) ! Time domain : restart ! --------------------- CALL trc_rst_cal( nit000, 'READ' ) ! calendar ! READ prognostic variables and computes diagnostic variable DO jn = 1, jptra CALL iom_get( numrtr, jpdom_autoglo, 'TRN'//ctrcnm(jn), trn(:,:,:,jn) ) END DO DO jn = 1, jptra CALL iom_get( numrtr, jpdom_autoglo, 'TRB'//ctrcnm(jn), trb(:,:,:,jn) ) END DO IF( lk_lobster ) CALL trc_rst_read_lobster( numrtr ) ! LOBSTER bio-model IF( lk_pisces ) CALL trc_rst_read_pisces ( numrtr ) ! PISCES bio-model IF( lk_cfc ) CALL trc_rst_read_cfc ( numrtr ) ! CFC tracers IF( lk_c14b ) CALL trc_rst_read_c14b ( numrtr ) ! C14 bomb tracer IF( lk_my_trc ) CALL trc_rst_read_my_trc ( numrtr ) ! MY_TRC tracers CALL iom_close( numrtr ) ! END SUBROUTINE trc_rst_read SUBROUTINE trc_rst_wri( kt ) !!---------------------------------------------------------------------- !! *** trc_rst_wri *** !! !! ** purpose : write passive tracer fields in restart files !!---------------------------------------------------------------------- INTEGER, INTENT( in ) :: kt ! ocean time-step index !! INTEGER :: jn REAL(wp) :: zarak0 !!---------------------------------------------------------------------- CALL trc_rst_cal( kt, 'WRITE' ) ! calendar CALL iom_rstput( kt, nitrst, numrtw, 'rdttrc1', rdttrc(1) ) ! surface passive tracer time step ! prognostic variables ! -------------------- DO jn = 1, jptra CALL iom_rstput( kt, nitrst, numrtw, 'TRN'//ctrcnm(jn), trn(:,:,:,jn) ) END DO DO jn = 1, jptra CALL iom_rstput( kt, nitrst, numrtw, 'TRB'//ctrcnm(jn), trb(:,:,:,jn) ) END DO IF( lk_lobster ) CALL trc_rst_wri_lobster( kt, nitrst, numrtw ) ! LOBSTER bio-model IF( lk_pisces ) CALL trc_rst_wri_pisces ( kt, nitrst, numrtw ) ! PISCES bio-model IF( lk_cfc ) CALL trc_rst_wri_cfc ( kt, nitrst, numrtw ) ! CFC tracers IF( lk_c14b ) CALL trc_rst_wri_c14b ( kt, nitrst, numrtw ) ! C14 bomb tracer IF( lk_my_trc ) CALL trc_rst_wri_my_trc ( kt, nitrst, numrtw ) ! MY_TRC tracers IF( kt == nitrst ) THEN CALL trc_rst_stat ! statistics CALL iom_close( numrtw ) ! close the restart file (only at last time step) #if ! defined key_trdmld_trc lrst_trc = .FALSE. #endif ENDIF ! END SUBROUTINE trc_rst_wri SUBROUTINE trc_rst_cal( kt, cdrw ) !!--------------------------------------------------------------------- !! *** ROUTINE trc_rst_cal *** !! !! ** Purpose : Read or write calendar in restart file: !! !! WRITE(READ) mode: !! kt : number of time step since the begining of the experiment at the !! end of the current(previous) run !! adatrj(0) : number of elapsed days since the begining of the experiment at the !! end of the current(previous) run (REAL -> keep fractions of day) !! ndastp : date at the end of the current(previous) run (coded as yyyymmdd integer) !! !! According to namelist parameter nrstdt, !! nn_rsttr = 0 no control on the date (nit000 is arbitrary). !! nn_rsttr = 1 we verify that nit000 is equal to the last !! time step of previous run + 1. !! In both those options, the exact duration of the experiment !! since the beginning (cumulated duration of all previous restart runs) !! is not stored in the restart and is assumed to be (nit000-1)*rdt. !! This is valid is the time step has remained constant. !! !! nn_rsttr = 2 the duration of the experiment in days (adatrj) !! has been stored in the restart file. !!---------------------------------------------------------------------- INTEGER , INTENT(in) :: kt ! ocean time-step CHARACTER(len=*), INTENT(in) :: cdrw ! "READ"/"WRITE" flag ! REAL(wp) :: zkt, zrdttrc1 REAL(wp) :: zndastp ! Time domain : restart ! --------------------- IF( TRIM(cdrw) == 'READ' ) THEN CALL iom_get ( numrtr, 'kt', zkt ) ! last time-step of previous run IF(lwp) THEN WRITE(numout,*) ' *** Info read in restart : ' WRITE(numout,*) ' previous time-step : ', NINT( zkt ) WRITE(numout,*) ' *** restart option' SELECT CASE ( nn_rsttr ) CASE ( 0 ) ; WRITE(numout,*) ' nn_rsttr = 0 : no control of nit000' CASE ( 1 ) ; WRITE(numout,*) ' nn_rsttr = 1 : no control the date at nit000 (use ndate0 read in the namelist)' CASE ( 2 ) ; WRITE(numout,*) ' nn_rsttr = 2 : calendar parameters read in restart' END SELECT WRITE(numout,*) ENDIF ! Control of date IF( nit000 - NINT( zkt ) /= 1 .AND. nn_rsttr /= 0 ) & & CALL ctl_stop( ' ===>>>> : problem with nit000 for the restart', & & ' verify the restart file or rerun with nn_rsttr = 0 (namelist)' ) IF( lk_offline ) THEN ! set the date in offline mode ! Check dynamics and tracer time-step consistency and force Euler restart if changed IF( iom_varid( numrtr, 'rdttrc1', ldstop = .FALSE. ) > 0 ) THEN CALL iom_get( numrtr, 'rdttrc1', zrdttrc1 ) IF( zrdttrc1 /= rdttrc(1) ) neuler = 0 ENDIF ! ! define ndastp and adatrj IF ( nn_rsttr == 2 ) THEN CALL iom_get( numrtr, 'ndastp', zndastp ) ndastp = NINT( zndastp ) CALL iom_get( numrtr, 'adatrj', adatrj ) ELSE ndastp = ndate0 - 1 ! ndate0 read in the namelist in dom_nam adatrj = ( REAL( nit000-1, wp ) * rdttra(1) ) / rday ! note this is wrong if time step has changed during run ENDIF ! IF(lwp) THEN WRITE(numout,*) ' *** Info used values : ' WRITE(numout,*) ' date ndastp : ', ndastp WRITE(numout,*) ' number of elapsed days since the begining of run : ', adatrj WRITE(numout,*) ENDIF ! CALL day_init ! compute calendar ! ENDIF ! ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN ! IF( kt == nitrst ) THEN IF(lwp) WRITE(numout,*) IF(lwp) WRITE(numout,*) 'trc_wri : write the TOP restart file (NetCDF) at it= ', kt, ' date= ', ndastp IF(lwp) WRITE(numout,*) '~~~~~~~' ENDIF CALL iom_rstput( kt, nitrst, numrtw, 'kt' , REAL( kt , wp) ) ! time-step CALL iom_rstput( kt, nitrst, numrtw, 'ndastp' , REAL( ndastp, wp) ) ! date CALL iom_rstput( kt, nitrst, numrtw, 'adatrj' , adatrj ) ! number of elapsed days since ! ! the begining of the run [s] ENDIF END SUBROUTINE trc_rst_cal SUBROUTINE trc_rst_stat !!---------------------------------------------------------------------- !! *** trc_rst_stat *** !! !! ** purpose : Compute tracers statistics !!---------------------------------------------------------------------- INTEGER :: jn REAL(wp) :: ztraf, zmin, zmax, zmean, zdrift !!---------------------------------------------------------------------- IF( lwp ) THEN WRITE(numout,*) WRITE(numout,*) ' ----TRACER STAT---- ' WRITE(numout,*) ENDIF DO jn = 1, jptra IF( .NOT.lk_degrad ) THEN ; ztraf = glob_sum( trn(:,:,:,jn) * cvol(:,:,:) ) ELSE ; ztraf = glob_sum( trn(:,:,:,jn) * cvol(:,:,:) * facvol(:,:,:) ) ENDIF zmin = MINVAL( trn(:,:,:,jn), mask= ((tmask*SPREAD(tmask_i,DIM=3,NCOPIES=jpk).NE.0.)) ) zmax = MAXVAL( trn(:,:,:,jn), mask= ((tmask*SPREAD(tmask_i,DIM=3,NCOPIES=jpk).NE.0.)) ) IF( lk_mpp ) THEN CALL mpp_min( zmin ) ! min over the global domain CALL mpp_max( zmax ) ! max over the global domain END IF zmean = ztraf / areatot zdrift = ( ( ztraf - trai(jn) ) / ( trai(jn) + 1.e-12 ) ) * 100._wp IF(lwp) WRITE(numout,*) ' tracer nb : ', jn,' ', TRIM( ctrcnm(jn) ) , & & ' mean = ', zmean, ' min = ', zmin, ' max = ', zmax, ' drift = ', zdrift, ' %' END DO WRITE(numout,*) END SUBROUTINE trc_rst_stat #else !!---------------------------------------------------------------------- !! Dummy module : No passive tracer !!---------------------------------------------------------------------- CONTAINS SUBROUTINE trc_rst_read ! Empty routines END SUBROUTINE trc_rst_read SUBROUTINE trc_rst_wri( kt ) INTEGER, INTENT ( in ) :: kt WRITE(*,*) 'trc_rst_wri: You should not have seen this print! error?', kt END SUBROUTINE trc_rst_wri #endif !!---------------------------------------------------------------------- !! NEMO/TOP 3.3 , NEMO Consortium (2010) !! $Id$ !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) !!====================================================================== END MODULE trcrst