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 3862 for branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/DIA – NEMO

Ignore:
Timestamp:
2013-04-08T13:07:43+02:00 (11 years ago)
Author:
acc
Message:

Branch 2013/dev_r3858_NOC_ZTC #863. Change names of gdept_0, gdepw_0, e3t_0, e3w_0 to gdept_1d, gdepw_1d, e3t_1d, e3w_1d respectively in preparation for the introduction of time invariant 3d arrays using the _0 suffix

Location:
branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/DIA
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90

    r3294 r3862  
    212212               ik = mbkt(ji,jj) 
    213213               IF( ik > 1 ) THEN 
    214                   zztmp = ( gdept_0(ik) - fsdept_0(ji,jj,ik) ) / ( gdept_0(ik) - gdept_0(ik-1) ) 
     214                  zztmp = ( gdept_1d(ik) - fsdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
    215215                  sn0(ji,jj,ik) = ( 1._wp - zztmp ) * sn0(ji,jj,ik) + zztmp * sn0(ji,jj,ik-1) 
    216216               ENDIF 

  • branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/DIA/diadimg.F90

    r3294 r3862  
    112112 
    113113    CASE ( 'T') 
    114        z4dep(:)=gdept_0(:) 
     114       z4dep(:)=gdept_1d(:) 
    115115 
    116116    CASE ( 'W' ) 
    117        z4dep(:)=gdepw_0(:) 
     117       z4dep(:)=gdepw_1d(:) 
    118118 
    119119    CASE ( '2' ) 

  • branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/DIA/diahth.F90

    r3764 r3862  
    304304      ! ----------------------------- ! 
    305305 
    306       ! find ilevel with (ilevel+1) the deepest W-level above 300m (we assume we can use e3t_0 to do this search...) 
     306      ! find ilevel with (ilevel+1) the deepest W-level above 300m (we assume we can use e3t_1d to do this search...) 
    307307      ilevel   = 0 
    308308      zthick_0 = 0._wp 
    309309      DO jk = 1, jpkm1                       
    310          zthick_0 = zthick_0 + e3t_0(jk) 
     310         zthick_0 = zthick_0 + e3t_1d(jk) 
    311311         IF( zthick_0 < 300. )   ilevel = jk 
    312312      END DO 

  • branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90

    r3764 r3862  
    662662            CALL histbeg(clhstnam, 1, zfoo, jpj, zphi,   & 
    663663               1, 1, 1, jpj, niter, zjulian, zdt*nn_fptr, nhoridz, numptr, domain_id=nidom_ptr) 
    664             ! Vertical grids : gdept_0, gdepw_0 
     664            ! Vertical grids : gdept_1d, gdepw_1d 
    665665            CALL histvert( numptr, "deptht", "Vertical T levels",   & 
    666                &                   "m", jpk, gdept_0, ndepidzt, "down" ) 
     666               &                   "m", jpk, gdept_1d, ndepidzt, "down" ) 
    667667            CALL histvert( numptr, "depthw", "Vertical W levels",   & 
    668                &                   "m", jpk, gdepw_0, ndepidzw, "down" ) 
     668               &                   "m", jpk, gdepw_1d, ndepidzw, "down" ) 
    669669            ! 
    670670            CALL wheneq ( jpj*jpk, MIN(sjk(:,:,1), 1._wp), 1, 1., ndex  , ndim  )      ! Lat-Depth 

  • branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r3704 r3862  
    323323            &          nit000-1, zjulian, zdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
    324324         CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept 
    325             &           "m", ipk, gdept_0, nz_T, "down" ) 
     325            &           "m", ipk, gdept_1d, nz_T, "down" ) 
    326326         !                                                            ! Index of ocean points 
    327327         CALL wheneq( jpi*jpj*ipk, tmask, 1, 1., ndex_T , ndim_T  )      ! volume 
    … …  
    359359            &          nit000-1, zjulian, zdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
    360360         CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept 
    361             &           "m", ipk, gdept_0, nz_U, "down" ) 
     361            &           "m", ipk, gdept_1d, nz_U, "down" ) 
    362362         !                                                            ! Index of ocean points 
    363363         CALL wheneq( jpi*jpj*ipk, umask, 1, 1., ndex_U , ndim_U  )      ! volume 
    … …  
    372372            &          nit000-1, zjulian, zdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
    373373         CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept 
    374             &          "m", ipk, gdept_0, nz_V, "down" ) 
     374            &          "m", ipk, gdept_1d, nz_V, "down" ) 
    375375         !                                                            ! Index of ocean points 
    376376         CALL wheneq( jpi*jpj*ipk, vmask, 1, 1., ndex_V , ndim_V  )      ! volume 
    … …  
    385385            &          nit000-1, zjulian, zdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 
    386386         CALL histvert( nid_W, "depthw", "Vertical W levels",      &  ! Vertical grid: gdepw 
    387             &          "m", ipk, gdepw_0, nz_W, "down" ) 
     387            &          "m", ipk, gdepw_1d, nz_W, "down" ) 
    388388 
    389389 
    … …  
    811811          1, jpi, 1, jpj, nit000-1, zjulian, zdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 
    812812      CALL histvert( id_i, "deptht", "Vertical T levels",   &    ! Vertical grid : gdept 
    813           "m", jpk, gdept_0, nz_i, "down") 
     813          "m", jpk, gdept_1d, nz_i, "down") 
    814814 
    815815      ! Declare all the output fields as NetCDF variables 
Note: See TracChangeset for help on using the changeset viewer.