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 6225 for branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/SAS_SRC/diawri.F90 – NEMO

Ignore:
Timestamp:
2016-01-08T10:35:19+01:00 (8 years ago)
Author:
jamesharle
Message:

Update MPP_BDY_UPDATE branch to be consistent with head of trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/SAS_SRC/diawri.F90

    • Property svn:keywords set to Id
    r4292 r6225  
    2626   USE dom_oce         ! ocean space and time domain 
    2727   USE zdf_oce         ! ocean vertical physics 
    28    USE ldftra_oce      ! ocean active tracers: lateral physics 
    29    USE ldfdyn_oce      ! ocean dynamics: lateral physics 
    30    USE traldf_iso_grif, ONLY : psix_eiv, psiy_eiv 
    31    USE sol_oce         ! solver variables 
    3228   USE sbc_oce         ! Surface boundary condition: ocean fields 
    3329   USE sbc_ice         ! Surface boundary condition: ice fields 
     
    4036   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    4137   USE in_out_manager  ! I/O manager 
    42    USE diadimg         ! dimg direct access file format output 
    4338   USE diaar5, ONLY :   lk_diaar5 
    4439   USE iom 
     
    6560 
    6661   !! * Substitutions 
    67 #  include "zdfddm_substitute.h90" 
    68 #  include "domzgr_substitute.h90" 
    6962#  include "vectopt_loop_substitute.h90" 
    7063   !!---------------------------------------------------------------------- 
    7164   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    72    !! $Id $ 
     65   !! $Id$ 
    7366   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    7467   !!---------------------------------------------------------------------- 
     
    8578  END FUNCTION dia_wri_alloc 
    8679 
    87 #if defined key_dimgout 
    88    !!---------------------------------------------------------------------- 
    89    !!   'key_dimgout'                                      DIMG output file 
    90    !!---------------------------------------------------------------------- 
    91 #   include "diawri_dimg.h90" 
    92  
    93 #else 
    9480   !!---------------------------------------------------------------------- 
    9581   !!   Default option                                   NetCDF output file 
    9682   !!---------------------------------------------------------------------- 
    97 # if defined key_iomput 
     83#if defined key_iomput 
    9884   !!---------------------------------------------------------------------- 
    9985   !!   'key_iomput'                                        use IOM library 
     
    11096      !! ** Method  :  use iom_put 
    11197      !!---------------------------------------------------------------------- 
    112       !! 
    113       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
     98      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
    11499      !!---------------------------------------------------------------------- 
    115100      !  
     
    144129      INTEGER  ::   ierr                                     ! error code return from allocation 
    145130      INTEGER  ::   iimi, iima, ipk, it, itmod, ijmi, ijma   ! local integers 
    146       REAL(wp) ::   zsto, zout, zmax, zjulian, zdt           ! local scalars 
     131      REAL(wp) ::   zsto, zout, zmax, zjulian                ! local scalars 
    147132      !!---------------------------------------------------------------------- 
    148133      !  
     
    163148 
    164149      ! Define frequency of output and means 
    165       zdt = rdt 
    166       IF( nacc == 1 ) zdt = rdtmin 
    167150      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    168151      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
    169152      ENDIF 
    170153#if defined key_diainstant 
    171       zsto = nwrite * zdt 
     154      zsto = nwrite * rdt 
    172155      clop = "inst("//TRIM(clop)//")" 
    173156#else 
    174       zsto=zdt 
     157      zsto=rdt 
    175158      clop = "ave("//TRIM(clop)//")" 
    176159#endif 
    177       zout = nwrite * zdt 
    178       zmax = ( nitend - nit000 + 1 ) * zdt 
     160      zout = nwrite * rdt 
     161      zmax = ( nitend - nit000 + 1 ) * rdt 
    179162 
    180163      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    218201         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    219202            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    220             &          nit000-1, zjulian, zdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
     203            &          nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
    221204         CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept 
    222205            &           "m", ipk, gdept_1d, nz_T, "down" ) 
     
    230213         CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu,           &  ! Horizontal grid: glamu and gphiu 
    231214            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    232             &          nit000-1, zjulian, zdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
     215            &          nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
    233216         CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept 
    234217            &           "m", ipk, gdept_1d, nz_U, "down" ) 
     
    242225         CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv,           &  ! Horizontal grid: glamv and gphiv 
    243226            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    244             &          nit000-1, zjulian, zdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
     227            &          nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
    245228         CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept 
    246229            &          "m", ipk, gdept_1d, nz_V, "down" ) 
     
    338321      ! 
    339322   END SUBROUTINE dia_wri 
    340 # endif 
    341  
    342323#endif 
    343324 
     
    362343      INTEGER  ::   id_i , nz_i, nh_i        
    363344      INTEGER, DIMENSION(1) ::   idex             ! local workspace 
    364       REAL(wp) ::   zsto, zout, zmax, zjulian, zdt 
     345      REAL(wp) ::   zsto, zout, zmax, zjulian 
    365346      !!---------------------------------------------------------------------- 
    366347      !  
     
    373354      clname = cdfile_name 
    374355      IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname) 
    375       zdt  = rdt 
    376356      zsto = rdt 
    377357      clop = "inst(x)"           ! no use of the mask value (require less cpu time) 
    378358      zout = rdt 
    379       zmax = ( nitend - nit000 + 1 ) * zdt 
     359      zmax = ( nitend - nit000 + 1 ) * rdt 
    380360 
    381361      IF(lwp) WRITE(numout,*) 
     
    392372      zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    393373      CALL histbeg( clname, jpi, glamt, jpj, gphit,   & 
    394           1, jpi, 1, jpj, nit000-1, zjulian, zdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 
     374          1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 
    395375      CALL histvert( id_i, "deptht", "Vertical T levels",   &    ! Vertical grid : gdept 
    396376          "m", jpk, gdept_1d, nz_i, "down") 
     
    435415      ! ----------------- 
    436416      CALL histclo( id_i ) 
    437 #if ! defined key_iomput && ! defined key_dimgout 
     417#if ! defined key_iomput 
    438418      IF( ninist /= 1  ) THEN 
    439419         CALL histclo( nid_T ) 
Note: See TracChangeset for help on using the changeset viewer.