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.
Changeset 12969 for NEMO/branches/2020/dev_12905_xios_restart/src/TOP – NEMO

Ignore:
Timestamp:
2020-05-26T12:05:09+02:00 (4 years ago)
Author:
andmirek
Message:

ticket #2462: read restart with XIOS independently for each component

Location:
NEMO/branches/2020/dev_12905_xios_restart/src/TOP
Files:
11 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/C14/trcini_c14.F90

    r12961 r12969  
    6868        IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 
    6969        !  
    70         IF(lrxios) CALL iom_swap(crtxios_context) 
    71         CALL iom_get( numrtr, 'co2sbc', co2sbc, ldxios = lrxios )  
    72         CALL iom_get( numrtr, jpdom_autoglo, 'c14sbc', c14sbc, ldxios = lrxios )  
    73         CALL iom_get( numrtr, jpdom_autoglo, 'exch_co2', exch_co2, ldxios = lrxios )  
    74         CALL iom_get( numrtr, jpdom_autoglo, 'exch_c14', exch_c14, ldxios = lrxios )  
    75         CALL iom_get( numrtr, jpdom_autoglo, 'qtr_c14', qtr_c14, ldxios = lrxios ) 
    76         IF(lrxios) CALL iom_swap(cxios_context) 
     70        IF(lrtxios) CALL iom_swap(crtxios_context) 
     71        CALL iom_get( numrtr, 'co2sbc', co2sbc, ldxios = lrtxios )  
     72        CALL iom_get( numrtr, jpdom_autoglo, 'c14sbc', c14sbc, ldxios = lrtxios )  
     73        CALL iom_get( numrtr, jpdom_autoglo, 'exch_co2', exch_co2, ldxios = lrtxios )  
     74        CALL iom_get( numrtr, jpdom_autoglo, 'exch_c14', exch_c14, ldxios = lrtxios )  
     75        CALL iom_get( numrtr, jpdom_autoglo, 'qtr_c14', qtr_c14, ldxios = lrtxios ) 
     76        IF(lrtxios) CALL iom_swap(cxios_context) 
    7777        ! 
    7878      END IF 
    … …  
    8787      ELSE 
    8888        ! 
    89         IF(lrxios) CALL iom_swap(crtxios_context) 
    90         CALL iom_get( numrtr, jpdom_autoglo, 'qint_c14', qint_c14, ldxios = lrxios )  
    91         IF(lrxios) CALL iom_swap(cxios_context) 
     89        IF(lrtxios) CALL iom_swap(crtxios_context) 
     90        CALL iom_get( numrtr, jpdom_autoglo, 'qint_c14', qint_c14, ldxios = lrtxios )  
     91        IF(lrtxios) CALL iom_swap(cxios_context) 
    9292        ! 
    9393      ENDIF 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/CFC/trcsms_cfc.F90

    r12961 r12969  
    297297         ! 
    298298         jl = 0 
    299          IF(lrxios) CALL iom_swap(crtxios_context) 
     299         IF(lrtxios) CALL iom_swap(crtxios_context) 
    300300         DO jn = jp_cfc0, jp_cfc1 
    301301            jl = jl + 1 
    302             CALL iom_get( numrtr, jpdom_autoglo, 'qint_'//ctrcnm(jn), qint_cfc(:,:,jl), ldxios = lrxios )  
     302            CALL iom_get( numrtr, jpdom_autoglo, 'qint_'//ctrcnm(jn), qint_cfc(:,:,jl), ldxios = lrtxios )  
    303303         END DO 
    304          IF(lrxios) CALL iom_swap(cxios_context) 
     304         IF(lrtxios) CALL iom_swap(cxios_context) 
    305305      ENDIF 
    306306      IF(lwp) WRITE(numout,*) 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/PISCES/P2Z/p2zexp.F90

    r12961 r12969  
    215215      ! 
    216216      IF( ln_rsttr ) THEN 
    217          IF(lrxios) CALL iom_swap(crtxios_context) 
    218          CALL iom_get( numrtr, jpdom_autoglo, 'SEDB'//ctrcnm(jpdet), sedpocb(:,:), ldxios = lrxios ) 
    219          CALL iom_get( numrtr, jpdom_autoglo, 'SEDN'//ctrcnm(jpdet), sedpocn(:,:), ldxios = lrxios ) 
    220          IF(lrxios) CALL iom_swap(cxios_context) 
     217         IF(lrtxios) CALL iom_swap(crtxios_context) 
     218         CALL iom_get( numrtr, jpdom_autoglo, 'SEDB'//ctrcnm(jpdet), sedpocb(:,:), ldxios = lrtxios ) 
     219         CALL iom_get( numrtr, jpdom_autoglo, 'SEDN'//ctrcnm(jpdet), sedpocn(:,:), ldxios = lrtxios ) 
     220         IF(lrtxios) CALL iom_swap(cxios_context) 
    221221      ELSE 
    222222         sedpocb(:,:) = 0._wp 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/PISCES/P4Z/p4zsms.F90

    r12961 r12969  
    339339         IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 
    340340         !  
    341          IF(lrxios) CALL iom_swap(crtxios_context) 
     341         IF(lrtxios) CALL iom_swap(crtxios_context) 
    342342         IF( iom_varid( numrtr, 'PH', ldstop = .FALSE. ) > 0 ) THEN 
    343             CALL iom_get( numrtr, jpdom_autoglo, 'PH' , hi(:,:,:), ldxios = lrxios  ) 
     343            CALL iom_get( numrtr, jpdom_autoglo, 'PH' , hi(:,:,:), ldxios = lrtxios  ) 
    344344         ELSE 
    345345            CALL p4z_che( Kbb, Kmm )                  ! initialize the chemical constants 
    346346            CALL ahini_for_at( hi, Kbb ) 
    347347         ENDIF 
    348          CALL iom_get( numrtr, jpdom_autoglo, 'Silicalim', xksi(:,:), ldxios = lrxios ) 
     348         CALL iom_get( numrtr, jpdom_autoglo, 'Silicalim', xksi(:,:), ldxios = lrtxios ) 
    349349         IF( iom_varid( numrtr, 'Silicamax', ldstop = .FALSE. ) > 0 ) THEN 
    350             CALL iom_get( numrtr, jpdom_autoglo, 'Silicamax' , xksimax(:,:), ldxios = lrxios  ) 
     350            CALL iom_get( numrtr, jpdom_autoglo, 'Silicamax' , xksimax(:,:), ldxios = lrtxios  ) 
    351351         ELSE 
    352352            xksimax(:,:) = xksi(:,:) 
    … …  
    354354         ! 
    355355         IF( iom_varid( numrtr, 'tcflxcum', ldstop = .FALSE. ) > 0 ) THEN  ! cumulative total flux of carbon 
    356             CALL iom_get( numrtr, 'tcflxcum' , t_oce_co2_flx_cum, ldxios = lrxios  ) 
     356            CALL iom_get( numrtr, 'tcflxcum' , t_oce_co2_flx_cum, ldxios = lrtxios  ) 
    357357         ELSE 
    358358            t_oce_co2_flx_cum = 0._wp 
    … …  
    361361         IF( ln_p5z ) THEN 
    362362            IF( iom_varid( numrtr, 'sized', ldstop = .FALSE. ) > 0 ) THEN 
    363                CALL iom_get( numrtr, jpdom_autoglo, 'sizep' , sizep(:,:,:), ldxios = lrxios  ) 
    364                CALL iom_get( numrtr, jpdom_autoglo, 'sizen' , sizen(:,:,:), ldxios = lrxios  ) 
    365                CALL iom_get( numrtr, jpdom_autoglo, 'sized' , sized(:,:,:), ldxios = lrxios  ) 
     363               CALL iom_get( numrtr, jpdom_autoglo, 'sizep' , sizep(:,:,:), ldxios = lrtxios  ) 
     364               CALL iom_get( numrtr, jpdom_autoglo, 'sizen' , sizen(:,:,:), ldxios = lrtxios  ) 
     365               CALL iom_get( numrtr, jpdom_autoglo, 'sized' , sized(:,:,:), ldxios = lrtxios  ) 
    366366            ELSE 
    367367               sizep(:,:,:) = 1. 
    … …  
    370370            ENDIF 
    371371        ENDIF 
    372         IF(lrxios) CALL iom_swap(cxios_context) 
     372        IF(lrtxios) CALL iom_swap(cxios_context) 
    373373        ! 
    374374      ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/PISCES/SED/sed.F90

    r10425 r12969  
    5353   INTEGER             , PUBLIC   ::  nn_rstsed      !: control of the time step ( 0 or 1 ) for pass. tr. 
    5454   INTEGER , PUBLIC               ::  nn_dtsed = 1   !: frequency of step on passive tracers 
     55   LOGICAL,  PUBLIC               ::  lrsxios        !: use xios to read restart 
    5556   CHARACTER(len = 80) , PUBLIC   ::  cn_sedrst_in   !: suffix of pass. tracer restart name (input) 
    5657   CHARACTER(len = 256), PUBLIC   ::  cn_sedrst_indir  !: restart input directory 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/PISCES/SED/sedrst.F90

    r12961 r12969  
    139139      zdta2 = 0. 
    140140 
    141       IF(lrxios) CALL iom_swap(crsxios_context) 
     141      IF(lrsxios) CALL iom_swap(crsxios_context) 
    142142      DO jn = 1, jptrased 
    143143         cltra = TRIM(sedtrcd(jn)) 
    144144         IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN 
    145             CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta(:,:,:,jn), ldxios = lrxios ) 
     145            CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta(:,:,:,jn), ldxios = lrsxios ) 
    146146         ELSE 
    147147            zdta(:,:,:,jn) = 0.0 
    … …  
    162162         cltra = TRIM(seddia3d(jn)) 
    163163         IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN 
    164             CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta1(:,:,:,jn), ldxios = lrxios ) 
     164            CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta1(:,:,:,jn), ldxios = lrsxios ) 
    165165         ELSE 
    166166            zdta1(:,:,:,jn) = 0.0 
    … …  
    189189      cltra = "dbioturb" 
    190190      IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN 
    191          CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:), ldxios = lrxios ) 
     191         CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:), ldxios = lrsxios ) 
    192192      ELSE 
    193193         zdta2(:,:,:) = 0.0 
    … …  
    199199      cltra = "irrig" 
    200200      IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN 
    201          CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:), ldxios = lrxios ) 
     201         CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:), ldxios = lrsxios ) 
    202202      ELSE 
    203203         zdta2(:,:,:) = 0.0 
    … …  
    209209      cltra = "sedligand" 
    210210      IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN 
    211          CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:), ldxios = lrxios ) 
     211         CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:), ldxios = lrsxios ) 
    212212      ELSE 
    213213         zdta2(:,:,:) = 0.0 
    … …  
    216216      CALL pack_arr( jpoce, sedligand(1:jpoce,1:jpksed), & 
    217217         &             zdta2(1:jpi,1:jpj,1:jpksed), iarroce(1:jpoce) ) 
    218       IF(lrxios) CALL iom_swap(cxios_context) 
     218      IF(lrsxios) CALL iom_swap(cxios_context) 
    219219      IF( ln_timing )  CALL timing_stop('sed_rst_read') 
    220220      
    … …  
    382382            CALL iom_open( TRIM(cn_sedrst_indir)//'/'//cn_sedrst_in, numrsr ) 
    383383 
    384       IF( lrxios .AND. .NOT. lxios_sini) THEN 
    385           CALL ctl_stop('OCE and SED restart must be in a single file when XIOS is used to read restart') 
    386       ENDIF 
    387       IF( lrxios) THEN 
    388           crsxios_context = 'sed_rst' 
    389           IF(lwp) WRITE(numout,*) 'Enable restart reading by XIOS for SED' 
    390           IF( TRIM(Agrif_CFixed()) == '0' ) THEN 
    391              clpname = cn_sedrst_in 
    392           ELSE 
    393              clpname = TRIM(Agrif_CFixed())//"_"//cn_sedrst_in    
    394           ENDIF 
    395           CALL iom_init( crsxios_context, fname = TRIM(cn_sedrst_indir)//'/'//TRIM(clpname), & 
    396                                           idfp = iom_file(numrsr)%nfid, ld_tmppatch = .TRUE. ) 
    397       ENDIF 
    398             IF(lrxios) CALL iom_swap(crsxios_context) 
    399             CALL iom_get ( numrsr, 'kt', zkt, ldxios = lrxios )   ! last time-step of previous run 
    400             IF(lrxios) CALL iom_swap(cxios_context) 
     384            lrsxios = lrxios.AND.lxios_sini 
     385            IF( lrsxios) THEN 
     386                crsxios_context = 'sed_rst' 
     387                IF(lwp) WRITE(numout,*) 'Enable restart reading by XIOS for SED' 
     388                IF( TRIM(Agrif_CFixed()) == '0' ) THEN 
     389                   clpname = cn_sedrst_in 
     390                ELSE 
     391                   clpname = TRIM(Agrif_CFixed())//"_"//cn_sedrst_in    
     392                ENDIF 
     393                CALL iom_init( crsxios_context, fname = TRIM(cn_sedrst_indir)//'/'//TRIM(clpname), & 
     394                                                idfp = iom_file(numrsr)%nfid, ld_tmppatch = .TRUE. ) 
     395            ENDIF 
     396            IF(lrsxios) CALL iom_swap(crsxios_context) 
     397            CALL iom_get ( numrsr, 'kt', zkt, ldxios = lrsxios )   ! last time-step of previous run 
     398            IF(lrsxios) CALL iom_swap(cxios_context) 
    401399            IF(lwp) THEN 
    402400               WRITE(numsed,*) ' *** Info read in restart : ' 
    … …  
    419417            !                                          ! set the date in offline mode 
    420418            IF( ln_rst_sed .AND. nn_rstsed == 2 ) THEN 
    421                IF(lrxios) CALL iom_swap(crsxios_context) 
    422                CALL iom_get( numrsr, 'ndastp', zndastp, ldxios = lrxios ) 
     419               IF(lrsxios) CALL iom_swap(crsxios_context) 
     420               CALL iom_get( numrsr, 'ndastp', zndastp, ldxios = lrsxios ) 
    423421               ndastp = NINT( zndastp ) 
    424                CALL iom_get( numrsr, 'adatrj', adatrj, ldxios = lrxios  ) 
    425                IF(lrxios) CALL iom_swap(crxios_context) 
     422               CALL iom_get( numrsr, 'adatrj', adatrj, ldxios = lrsxios  ) 
     423               IF(lrsxios) CALL iom_swap(crxios_context) 
    426424             ELSE 
    427425               ndastp = ndate0 - 1     ! ndate0 read in the namelist in dom_nam 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/PISCES/SED/sedstp.F90

    r12961 r12969  
    8686      IF( kt == nitsed000 ) THEN 
    8787          CALL iom_close( numrsr )       ! close input tracer restart file 
    88           IF(lrxios) CALL iom_context_finalize(      crsxios_context  ) 
     88          IF(lrsxios) CALL iom_context_finalize(      crsxios_context  ) 
    8989!          IF(lwm) CALL FLUSH( numont )   ! flush namelist output 
    9090      ENDIF 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/TRP/trcsbc.F90

    r12961 r12969  
    8686            IF(lwp) WRITE(numout,*) '          nittrc000-1 surface tracer content forcing fields read in the restart file' 
    8787            zfact = 0.5_wp 
    88             IF(lrxios) CALL iom_swap(crtxios_context) 
     88            IF(lrtxios) CALL iom_swap(crtxios_context) 
    8989            DO jn = 1, jptra 
    90                CALL iom_get( numrtr, jpdom_autoglo, 'sbc_'//TRIM(ctrcnm(jn))//'_b', sbc_trc_b(:,:,jn), ldxios = lrxios )   ! before tracer content sbc 
     90               CALL iom_get( numrtr, jpdom_autoglo, 'sbc_'//TRIM(ctrcnm(jn))//'_b', sbc_trc_b(:,:,jn), ldxios = lrtxios )   ! before tracer content sbc 
    9191            END DO 
    92             IF(lrxios) CALL iom_swap(cxios_context) 
     92            IF(lrtxios) CALL iom_swap(cxios_context) 
    9393         ELSE                                         ! No restart or restart not found: Euler forward time stepping 
    9494           zfact = 1._wp 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/trc.F90

    r12489 r12969  
    6969   INTEGER             , PUBLIC ::   nittrc000          !: first time step of passive tracers model 
    7070   LOGICAL             , PUBLIC ::   l_trcdm2dc         !: Diurnal cycle for TOP 
     71   LOGICAL             , PUBLIC ::   lrtxios            !: read TOP restart with XIOS 
    7172 
    7273   !! Information for the ice module for tracers 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/trcrst.F90

    r12961 r12969  
    130130 
    131131      ! READ prognostic variables and computes diagnostic variable 
    132       IF(lrxios) CALL iom_swap(crtxios_context) 
     132      IF(lrtxios) CALL iom_swap(crtxios_context) 
    133133      DO jn = 1, jptra 
    134          CALL iom_get( numrtr, jpdom_autoglo, 'TRN'//ctrcnm(jn), tr(:,:,:,jn,Kmm), ldxios = lrxios ) 
     134         CALL iom_get( numrtr, jpdom_autoglo, 'TRN'//ctrcnm(jn), tr(:,:,:,jn,Kmm), ldxios = lrtxios ) 
    135135      END DO 
    136136 
    137137      DO jn = 1, jptra 
    138          CALL iom_get( numrtr, jpdom_autoglo, 'TRB'//ctrcnm(jn), tr(:,:,:,jn,Kbb), ldxios = lrxios ) 
     138         CALL iom_get( numrtr, jpdom_autoglo, 'TRB'//ctrcnm(jn), tr(:,:,:,jn,Kbb), ldxios = lrtxios ) 
    139139      END DO 
    140140      ! 
    141141      CALL iom_delay_rst( 'READ', 'TOP', numrtr )   ! read only TOP delayed global communication variables 
    142       IF(lrxios) CALL iom_swap(cxios_context)  
     142      IF(lrtxios) CALL iom_swap(cxios_context)  
    143143   END SUBROUTINE trc_rst_read 
    144144 
    … …  
    234234            lxios_sini = .FALSE. 
    235235            CALL iom_open( TRIM(cn_trcrst_indir)//'/'//cn_trcrst_in, numrtr ) 
    236             IF( lrxios .AND. .NOT. lxios_sini) THEN 
    237                CALL ctl_stop('OCE and TOP restart must be in a single file when XIOS is used to read restart') 
    238             ENDIF 
    239             IF( lrxios) THEN 
     236            lrtxios = lrxios.AND.lxios_sini 
     237            IF( lrtxios) THEN 
    240238                crtxios_context = 'top_rst' 
    241239                IF(lwp) WRITE(numout,*) 'Enable restart reading by XIOS for TOP' 
    … …  
    249247            ENDIF 
    250248 
    251             IF(lrxios) CALL iom_swap(crtxios_context) 
    252             CALL iom_get ( numrtr, 'kt', zkt, ldxios = lrxios )   ! last time-step of previous run 
    253             IF(lrxios) CALL iom_swap(cxios_context) 
     249            IF(lrtxios) CALL iom_swap(crtxios_context) 
     250            CALL iom_get ( numrtr, 'kt', zkt, ldxios = lrtxios )   ! last time-step of previous run 
     251            IF(lrtxios) CALL iom_swap(cxios_context) 
    254252 
    255253            IF(lwp) THEN 
    … …  
    273271            !                                          ! set the date in offline mode 
    274272            IF( ln_rsttr .AND. nn_rsttr == 2 ) THEN 
    275                IF(lrxios) CALL iom_swap(crtxios_context) 
    276                CALL iom_get( numrtr, 'ndastp', zndastp, ldxios = lrxios ) 
     273               IF(lrtxios) CALL iom_swap(crtxios_context) 
     274               CALL iom_get( numrtr, 'ndastp', zndastp, ldxios = lrtxios ) 
    277275               ndastp = NINT( zndastp ) 
    278                CALL iom_get( numrtr, 'adatrj', adatrj, ldxios = lrxios  ) 
    279                CALL iom_get( numrtr, 'ntime' , ktime, ldxios = lrxios   ) 
    280                IF(lrxios) CALL iom_swap(cxios_context) 
     276               CALL iom_get( numrtr, 'adatrj', adatrj, ldxios = lrtxios  ) 
     277               CALL iom_get( numrtr, 'ntime' , ktime, ldxios = lrtxios   ) 
     278               IF(lrtxios) CALL iom_swap(cxios_context) 
    281279               nn_time0=INT(ktime) 
    282280               ! calculate start time in hours and minutes 

  • NEMO/branches/2020/dev_12905_xios_restart/src/TOP/trcstp.F90

    r12961 r12969  
    109109      IF( kt == nittrc000 ) THEN 
    110110         CALL iom_close( numrtr )                         ! close input tracer restart file 
    111          IF(lrxios) CALL iom_context_finalize(      crtxios_context          ) 
     111         IF(lrtxios) CALL iom_context_finalize(      crtxios_context          ) 
    112112         IF(lwm) CALL FLUSH( numont )                     ! flush namelist output 
    113113      ENDIF 
    … …  
    196196           &                              .AND. iom_varid( numrtr, 'ktdcy'    , ldstop = .FALSE. ) > 0  & 
    197197           &                              .AND. iom_varid( numrtr, 'nrdcy'    , ldstop = .FALSE. ) > 0  ) THEN 
    198             IF(lrxios) CALL iom_swap(crtxios_context) 
    199             CALL iom_get( numrtr, 'ktdcy', zkt, ldxios = lrxios )   
     198            IF(lrtxios) CALL iom_swap(crtxios_context) 
     199            CALL iom_get( numrtr, 'ktdcy', zkt, ldxios = lrtxios )   
    200200            rsecfst = INT( zkt ) * rn_Dt 
    201201            IF(lwp) WRITE(numout,*) 'trc_qsr_mean:   qsr_mean read in the restart file at time-step rsecfst =', rsecfst, ' s ' 
    202             CALL iom_get( numrtr, jpdom_autoglo, 'qsr_mean', qsr_mean, ldxios = lrxios )   !  A mean of qsr 
    203             CALL iom_get( numrtr, 'nrdcy', zrec, ldxios = lrxios )   !  Number of record per days 
     202            CALL iom_get( numrtr, jpdom_autoglo, 'qsr_mean', qsr_mean, ldxios = lrtxios )   !  A mean of qsr 
     203            CALL iom_get( numrtr, 'nrdcy', zrec, ldxios = lrtxios )   !  Number of record per days 
    204204            IF( INT( zrec ) == nb_rec_per_day ) THEN 
    205205               DO jn = 1, nb_rec_per_day  
    206206                  IF( jn <= 9 )  THEN 
    207207                    WRITE(cl1,'(i1)') jn 
    208                     CALL iom_get( numrtr, jpdom_autoglo, 'qsr_arr_'//cl1, qsr_arr(:,:,jn), ldxios = lrxios )   !  A mean of qsr 
     208                    CALL iom_get( numrtr, jpdom_autoglo, 'qsr_arr_'//cl1, qsr_arr(:,:,jn), ldxios = lrtxios )   !  A mean of qsr 
    209209                  ELSE 
    210210                    WRITE(cl2,'(i2.2)') jn 
    211                     CALL iom_get( numrtr, jpdom_autoglo, 'qsr_arr_'//cl2, qsr_arr(:,:,jn), ldxios = lrxios )   !  A mean of qsr 
     211                    CALL iom_get( numrtr, jpdom_autoglo, 'qsr_arr_'//cl2, qsr_arr(:,:,jn), ldxios = lrtxios )   !  A mean of qsr 
    212212                  ENDIF 
    213213              END DO 
    … …  
    217217               ENDDO 
    218218            ENDIF 
    219             IF(lrxios) CALL iom_swap(cxios_context) 
     219            IF(lrtxios) CALL iom_swap(cxios_context) 
    220220         ELSE                                         !* no restart: set from nit000 values 
    221221            IF(lwp) WRITE(numout,*) 'trc_qsr_mean:   qsr_mean set to nit000 values' 
Note: See TracChangeset for help on using the changeset viewer.