Changeset 1317


Ignore:
Timestamp:
2009-02-17T17:31:42+01:00 (12 years ago)
Author:
smasson
Message:

nwrite = modulo referenced to nit000 in all ouputs, see ticket:339

Location:
trunk
Files:
15 edited

Legend:

Unmodified
Added
Removed
  • trunk/CONFIG/GYRE/EXP00/namelist

    r1312 r1317  
    3636   nleapy      =      30   !  Leap year calendar (1) or not (0) 
    3737   ninist      =       0   !  output the initial state (1) or not (0) 
    38    nstock      =    4320   !  frequency of creation of a restart file 
    39    nwrite      =      60   !  frequency of write in the output file  
     38   nstock      =    4320   !  frequency of creation of a restart file (modulo referenced to 1) 
     39   nwrite      =      60   !  frequency of write in the output file   (modulo referenced to nit000) 
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4141   ln_mskland  = .false.   !  mask land points (1.e+20) in NetCDF outputs (costly: + ~15%) 
  • trunk/CONFIG/GYRE_LOBSTER/EXP00/namelist

    r1312 r1317  
    3636   nleapy      =      30   !  Leap year calendar (1) or not (0) 
    3737   ninist      =       0   !  output the initial state (1) or not (0) 
    38    nstock      =    4320   !  frequency of creation of a restart file 
    39    nwrite      =     360   !  frequency of write in the output file  
     38   nstock      =    4320   !  frequency of creation of a restart file (modulo referenced to 1) 
     39   nwrite      =     360   !  frequency of write in the output file   (modulo referenced to nit000) 
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4141   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
  • trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist

    r1312 r1317  
    3636   nleapy      =       0   !  Leap year calendar (1) or not (0) 
    3737   ninist      =       0   !  output the initial state (1) or not (0) 
    38    nstock      =   10950   !  frequency of creation of a restart file 
    39    nwrite      =   10950   !  frequency of write in the output file  
     38   nstock      =   10950   !  frequency of creation of a restart file (modulo referenced to 1) 
     39   nwrite      =   10950   !  frequency of write in the output file   (modulo referenced to nit000) 
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4141   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
  • trunk/CONFIG/ORCA2_LIM/EXP00/namelist

    r1312 r1317  
    3636   nleapy      =       0   !  Leap year calendar (1) or not (0) 
    3737   ninist      =       0   !  output the initial state (1) or not (0) 
    38    nstock      =    5475   !  frequency of creation of a restart file 
    39    nwrite      =    5475   !  frequency of write in the output file  
     38   nstock      =    5475   !  frequency of creation of a restart file (modulo referenced to 1) 
     39   nwrite      =    5475   !  frequency of write in the output file   (modulo referenced to nit000) 
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4141   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
  • trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist

    r1312 r1317  
    3636   nleapy      =       0   !  Leap year calendar (1) or not (0) 
    3737   ninist      =       0   !  output the initial state (1) or not (0) 
    38    nstock      =    5475   !  frequency of creation of a restart file 
    39    nwrite      =    5475   !  frequency of write in the output file  
     38   nstock      =    5475   !  frequency of creation of a restart file (modulo referenced to 1) 
     39   nwrite      =    5475   !  frequency of write in the output file   (modulo referenced to nit000) 
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4141   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
  • trunk/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist

    r1312 r1317  
    2828   nleapy      =       0   !  Leap year calendar (1) or not (0) 
    2929   ninist      =       0   !  output the initial state (1) or not (0) 
    30    nstock      =     7300  !  frequency of creation of a restart file 
    31    nwrite      =     7300  !  frequency of write in the output file  
     30   nstock      =     7300  !  frequency of creation of a restart file (modulo referenced to 1) 
     31   nwrite      =     7300  !  frequency of write in the output file   (modulo referenced to nit000) 
    3232   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    3333/ 
  • trunk/NEMO/C1D_SRC/diawri_c1d.F90

    r1312 r1317  
    276276      ! donne le nombre d'elements, et ndex la liste des indices a sortir 
    277277 
    278       IF( lwp .AND. MOD( kt, nwrite ) == 0 ) THEN  
     278      IF( lwp .AND. MOD( it, nwrite ) == 0 ) THEN  
    279279         WRITE(numout,*) 'dia_wri : write model outputs in NetCDF files at ', kt, 'time-step' 
    280280         WRITE(numout,*) '~~~~~~ ' 
     
    345345      ! 3. Synchronise and close all files 
    346346      ! --------------------------------------- 
    347       IF( MOD( kt, nwrite ) == 0 .OR. kindic < 0 )   CALL histsync( nid_T ) 
     347      IF( MOD( it, nwrite ) == 0 .OR. kindic < 0 )   CALL histsync( nid_T ) 
    348348 
    349349      !  Create an output files (output.abort.nc) if S < 0 or u > 20 m/s 
  • trunk/NEMO/OPA_SRC/DIA/diagap.F90

    r1312 r1317  
    198198      ! ---------------------- 
    199199 
    200       IF( MOD( kt, ngap ) == 0 ) THEN 
     200      it = kt - nit000 + 1       ! define time axis 
     201      IF( MOD( it, ngap ) == 0 ) THEN 
    201202 
    202203         ! initialization 
     
    239240          ! ----------------------------====== 
    240241 
    241           IF( MOD( kt, nprg ) == 0 ) THEN 
     242          IF( MOD( it, nprg ) == 0 ) THEN 
    242243              IF(lwp) THEN 
    243244                  WRITE(numout,*) 'dia_gap: time step = ', kt, 'model - data' 
     
    259260 
    260261          IF (lwp ) THEN   ! Ok even in mpp after the call to mpp_sum 
    261              it = kt - nit000 + 1       ! define time axis 
    262262             CALL histwrite( numgap, "motempda", it, tdtag, jpk, ndex ) 
    263263             CALL histwrite( numgap, "mosalida", it, sdtag, jpk, ndex ) 
  • trunk/NEMO/OPA_SRC/DIA/diaptr.F90

    r1316 r1317  
    245245      !!---------------------------------------------------------------------- 
    246246 
    247       IF( kt == nit000 .OR. MOD( kt, nf_ptr ) == 0 )   THEN 
     247      IF( kt == nit000 .OR. MOD( kt - nit000 + 1, nf_ptr ) == 0 )   THEN 
    248248 
    249249         zsverdrup = 1.e-6 
     
    551551         ENDIF 
    552552          
    553  
    554553         CALL histend( numptr ) 
    555554 
    556555      ENDIF 
    557556 
    558       IF( MOD( kt, nf_ptr ) == 0 ) THEN 
     557      IF( MOD( it, nf_ptr ) == 0 ) THEN 
    559558 
    560559         IF(lwp) THEN 
     
    565564         ENDIF 
    566565 
    567          ! define time axis 
    568          it= kt - nit000 + 1 
    569566         ndex(1) = 0 
    570567         CALL histwrite( numptr, "zotemglo", it, tn_jk    , jpj*jpk, ndex ) 
  • trunk/NEMO/OPA_SRC/DIA/diawri.F90

    r1312 r1317  
    406406      ! donne le nombre d'elements, et ndex la liste des indices a sortir 
    407407 
    408       IF( lwp .AND. MOD( kt, nwrite ) == 0 ) THEN  
     408      IF( lwp .AND. MOD( it, nwrite ) == 0 ) THEN  
    409409         WRITE(numout,*) 'dia_wri : write model outputs in NetCDF files at ', kt, 'time-step' 
    410410         WRITE(numout,*) '~~~~~~ ' 
     
    513513      ! 3. Synchronise and close all files 
    514514      ! --------------------------------------- 
    515       IF( MOD( kt, nwrite ) == 0 .OR. kindic < 0 ) THEN 
     515      IF( MOD( it, nwrite ) == 0 .OR. kindic < 0 ) THEN 
    516516         CALL histsync( nid_T ) 
    517517         CALL histsync( nid_U ) 
  • trunk/NEMO/OPA_SRC/TRD/trdmld.F90

    r1312 r1317  
    389389      smltrd(:,:,:) = smltrd(:,:,:) * ucf   !  is no longer used, and is reset to 0. at next time step) 
    390390       
    391       MODULO_NTRD : IF( MOD( kt, ntrd ) == 0 ) THEN        ! nitend MUST be multiple of ntrd 
     391      it = kt - nit000 + 1 
     392 
     393      MODULO_NTRD : IF( MOD( it, ntrd ) == 0 ) THEN        ! nitend MUST be multiple of ntrd 
    392394         ! 
    393395         ztmltot (:,:) = 0.e0   ;   zsmltot (:,:) = 0.e0   ! reset arrays to zero 
     
    574576#if defined key_dimgout 
    575577 
    576       IF( MOD( kt, ntrd ) == 0 ) THEN 
     578      IF( MOD( it, ntrd ) == 0 ) THEN 
    577579         iyear =  ndastp/10000 
    578580         imon  = (ndastp-iyear*10000)/100 
     
    591593      ! ---------------------------------- 
    592594 
    593       IF( lwp .AND. MOD( kt , ntrd ) == 0 ) THEN 
     595      IF( lwp .AND. MOD( it , ntrd ) == 0 ) THEN 
    594596         WRITE(numout,*) ' ' 
    595597         WRITE(numout,*) 'trd_mld : write trends in the NetCDF file :' 
     
    600602      END IF 
    601603          
    602       it = kt - nit000 + 1 
    603  
    604604      !-- Write the trends for T/S instantaneous diagnostics  
    605605      IF( ln_trdmld_instant ) THEN            
     
    683683#endif 
    684684 
    685       IF( MOD( kt , ntrd ) == 0 ) THEN 
     685      IF( MOD( it, ntrd ) == 0 ) THEN 
    686686         ! 
    687687         ! III.5 Reset cumulative arrays to zero 
  • trunk/NEMO/OPA_SRC/TRD/trdvor.F90

    r1312 r1317  
    405405      !  ============================================= 
    406406 
    407       IF( MOD( kt - nit000+1, ntrd ) == 0 ) THEN 
     407      ! define time axis 
     408      it= kt - nit000 + 1 
     409 
     410      IF( MOD( it, ntrd ) == 0 ) THEN 
    408411 
    409412         ! III.1 compute total trend 
     
    452455      IF( kt >=  nit000+1 ) THEN 
    453456 
    454          ! define time axis 
    455          it= kt-nit000+1 
    456          IF( lwp .AND. MOD( kt, ntrd ) == 0 ) THEN 
     457         IF( lwp .AND. MOD( it, ntrd ) == 0 ) THEN 
    457458            WRITE(numout,*) '' 
    458459            WRITE(numout,*) 'trd_vor : write trends in the NetCDF file at kt = ', kt 
     
    481482      ENDIF 
    482483      ! 
    483       IF( MOD( kt - nit000+1, ntrd ) == 0 ) rotot(:,:)=0 
     484      IF( MOD( it, ntrd ) == 0 ) rotot(:,:)=0 
    484485      ! 
    485486      IF( kt == nitend )   CALL histclo( nidvor ) 
  • trunk/NEMO/TOP_SRC/SED/sedwri.F90

    r1312 r1317  
    251251         !! synchronise FILE 
    252252         ! 
    253          IF( MOD( kt, nwrised ) == 0 ) THEN 
     253         IF( MOD( it, nwrised ) == 0 ) THEN 
    254254              WRITE(numsed,*)  '**** sedwri : write NetCDF aditional arrays' 
    255255              CALL histsync( nised ) 
  • trunk/NEMO/TOP_SRC/TRP/trdmld_trc.F90

    r1312 r1317  
    575575      tmltrd_trc(:,:,:,:) = tmltrd_trc(:,:,:,:) * ucf_trc 
    576576 
    577       MODULO_NTRD : IF( MOD( kt, ntrd_trc ) == 0 ) THEN           ! nitend MUST be multiple of ntrd_trc 
     577      it = kt - nit000 + 1 
     578 
     579      MODULO_NTRD : IF( MOD( it, ntrd_trc ) == 0 ) THEN           ! nitend MUST be multiple of ntrd_trc 
    578580         ! 
    579581         ztmltot (:,:,:) = 0.e0                                   ! reset arrays to zero 
     
    817819      ! ---------------------------------- 
    818820 
    819       IF( lwp .AND. MOD( kt , ntrd_trc ) == 0 ) THEN 
     821      IF( lwp .AND. MOD( it , ntrd_trc ) == 0 ) THEN 
    820822         WRITE(numout,*) ' ' 
    821823         WRITE(numout,*) 'trd_mld_trc : write passive tracer trends in the NetCDF file :' 
     
    826828      ENDIF 
    827829          
    828       it = kt - nit000 + 1 
    829  
    830830      NETCDF_OUTPUT : IF( ln_trdmld_trc_instant ) THEN            ! <<< write the trends for passive tracer instant. diags 
    831831         ! 
     
    930930# endif /* key_dimgout */ 
    931931 
    932       IF( MOD( kt, ntrd_trc ) == 0 ) THEN 
     932      IF( MOD( it, ntrd_trc ) == 0 ) THEN 
    933933         ! 
    934934         ! Reset cumulative arrays to zero 
     
    11311131      ! ---------------------------------- 
    11321132 
    1133       IF( lwp .AND. MOD( kt , ntrd_trc ) == 0 ) THEN 
     1133      ! define time axis 
     1134      it = kt - nit000 + 1 
     1135 
     1136      IF( lwp .AND. MOD( it , ntrd_trc ) == 0 ) THEN 
    11341137         WRITE(numout,*) ' ' 
    11351138         WRITE(numout,*) 'trd_mld_bio : write ML bio trends in the NetCDF file :' 
     
    11411144 
    11421145 
    1143       ! define time axis 
    1144       it = kt - nit000 + 1 
    1145  
    1146  
    11471146      ! 2. Start writing data 
    11481147      ! --------------------- 
     
    11751174# endif /* key_dimgout */ 
    11761175 
    1177       IF( MOD( kt, ntrd_trc ) == 0 ) THEN 
     1176      IF( MOD( it, ntrd_trc ) == 0 ) THEN 
    11781177         ! 
    11791178         ! III.5 Reset cumulative arrays to zero 
  • trunk/NEMO/TOP_SRC/trcdia.F90

    r1316 r1317  
    201201      ! --------------------------------------- 
    202202 
    203       IF( lwp .AND. MOD( kt, nwritetrc ) == 0 ) THEN 
     203      IF( lwp .AND. MOD( it, nwritetrc ) == 0 ) THEN 
    204204         WRITE(numout,*) 'trcdit_wr : write NetCDF passive tracer concentrations at ', kt, 'time-step' 
    205205         WRITE(numout,*) '~~~~~~~~~ ' 
     
    214214 
    215215      ! synchronise file 
    216       IF( MOD( kt, nwritetrc ) == 0 .OR. kindic < 0 )   CALL histsync( nit5 ) 
     216      IF( MOD( it, nwritetrc ) == 0 .OR. kindic < 0 )   CALL histsync( nit5 ) 
    217217 
    218218 
     
    415415      ! trends for tracer concentrations 
    416416 
    417       IF( lwp .AND. MOD( kt, nwritetrd ) == 0 ) THEN 
     417      IF( lwp .AND. MOD( it, nwritetrd ) == 0 ) THEN 
    418418         WRITE(numout,*) 'trcdid_wr : write NetCDF dynamical trends at ', kt, 'time-step' 
    419419         WRITE(numout,*) '~~~~~~ ' 
     
    446446 
    447447      ! synchronise FILE 
    448       IF( MOD( kt, nwritetrd ) == 0 .OR. kindic < 0 ) THEN 
     448      IF( MOD( it, nwritetrd ) == 0 .OR. kindic < 0 ) THEN 
    449449         DO jn = 1, jptra 
    450450            IF (luttrd(jn))   CALL histsync( nit6(jn) ) 
     
    588588      ! --------------------- 
    589589 
    590       IF( lwp .AND. MOD( kt, nwritedia ) == 0 ) THEN 
     590      IF( lwp .AND. MOD( it, nwritedia ) == 0 ) THEN 
    591591         WRITE(numout,*) 'trcdii_wr : write NetCDF additional arrays at ', kt, 'time-step' 
    592592         WRITE(numout,*) '~~~~~~ ' 
     
    606606 
    607607      ! synchronise FILE 
    608       IF( MOD( kt, nwritedia ) == 0 .OR. kindic < 0 )   CALL histsync( nitd ) 
     608      IF( MOD( it, nwritedia ) == 0 .OR. kindic < 0 )   CALL histsync( nitd ) 
    609609 
    610610      ! Closing all files 
     
    724724 
    725725      ! biological trends 
    726       IF( lwp .AND. MOD( kt, nwritebio ) == 0 ) THEN 
     726      IF( lwp .AND. MOD( it, nwritebio ) == 0 ) THEN 
    727727         WRITE(numout,*) 'trcdit_wr : write NetCDF biological trends at ', kt, 'time-step' 
    728728         WRITE(numout,*) '~~~~~~ ' 
     
    735735 
    736736      ! synchronise FILE 
    737       IF( MOD( kt, nwritebio ) == 0 .OR. kindic < 0 )   CALL histsync( nitb ) 
     737      IF( MOD( it, nwritebio ) == 0 .OR. kindic < 0 )   CALL histsync( nitb ) 
    738738 
    739739      ! Closing all files 
Note: See TracChangeset for help on using the changeset viewer.