Changeset 5825


Ignore:
Timestamp:
2015-10-23T11:57:19+02:00 (6 years ago)
Author:
diovino
Message:
 
Location:
branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90

    r5282 r5825  
    567567      INTEGER            ::   inum                                    ! temporary logical unit 
    568568#endif 
    569       REAL(wp)           ::   zsto, zout, zdt, zjulian                ! temporary scalars 
     569      REAL(wp)           ::   zsto, zout, zjulian                ! temporary scalars 
    570570      !! 
    571571      REAL(wp), POINTER, DIMENSION(:)   ::   zphi, zfoo    ! 1D workspace 
     
    584584      IF( kt == nit000 ) THEN 
    585585         niter = ( nit000 - 1 ) / nn_fptr 
    586          zdt = rdt 
    587586         ! 
    588587         IF(lwp) THEN 
     
    640639            ! ---------------- 
    641640            ! Define frequency of output and means 
    642             zsto = nn_fptr * zdt 
     641            zsto = nn_fptr * rdt 
    643642            IF( ln_mskland )   THEN    ! put 1.e+20 on land (very expensive!!) 
    644643               clop      = "ave(only(x))" 
     
    649648            ENDIF 
    650649 
    651             zout = nn_fwri * zdt 
     650            zout = nn_fwri * rdt 
    652651            zfoo(1:jpj) = 0._wp 
    653652 
     
    670669            ! Horizontal grid : zphi() 
    671670            CALL histbeg(clhstnam, 1, zfoo, jpj, zphi,   & 
    672                1, 1, 1, jpj, niter, zjulian, zdt*nn_fptr, nhoridz, numptr, domain_id=nidom_ptr) 
     671               1, 1, 1, jpj, niter, zjulian, rdt*nn_fptr, nhoridz, numptr, domain_id=nidom_ptr) 
    673672            ! Vertical grids : gdept_1d, gdepw_1d 
    674673            CALL histvert( numptr, "deptht", "Vertical T levels",   & 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r5282 r5825  
    384384      INTEGER  ::   iimi, iima, ipk, it, itmod, ijmi, ijma   ! local integers 
    385385      INTEGER  ::   jn, ierror                               ! local integers 
    386       REAL(wp) ::   zsto, zout, zmax, zjulian, zdt           ! local scalars 
     386      REAL(wp) ::   zsto, zout, zmax, zjulian           ! local scalars 
    387387      !! 
    388388      REAL(wp), POINTER, DIMENSION(:,:)   :: zw2d       ! 2D workspace 
     
    409409 
    410410      ! Define frequency of output and means 
    411       zdt = rdt 
    412411      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    413412      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
    414413      ENDIF 
    415414#if defined key_diainstant 
    416       zsto = nwrite * zdt 
     415      zsto = nwrite * rdt 
    417416      clop = "inst("//TRIM(clop)//")" 
    418417#else 
    419       zsto=zdt 
     418      zsto=rdt 
    420419      clop = "ave("//TRIM(clop)//")" 
    421420#endif 
    422       zout = nwrite * zdt 
    423       zmax = ( nitend - nit000 + 1 ) * zdt 
     421      zout = nwrite * rdt 
     422      zmax = ( nitend - nit000 + 1 ) * rdt 
    424423 
    425424      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    463462         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    464463            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    465             &          nit000-1, zjulian, zdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
     464            &          nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
    466465         CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept 
    467466            &           "m", ipk, gdept_1d, nz_T, "down" ) 
     
    499498         CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu,           &  ! Horizontal grid: glamu and gphiu 
    500499            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    501             &          nit000-1, zjulian, zdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
     500            &          nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
    502501         CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept 
    503502            &           "m", ipk, gdept_1d, nz_U, "down" ) 
     
    512511         CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv,           &  ! Horizontal grid: glamv and gphiv 
    513512            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    514             &          nit000-1, zjulian, zdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
     513            &          nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
    515514         CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept 
    516515            &          "m", ipk, gdept_1d, nz_V, "down" ) 
     
    525524         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    526525            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    527             &          nit000-1, zjulian, zdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 
     526            &          nit000-1, zjulian, rdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 
    528527         CALL histvert( nid_W, "depthw", "Vertical W levels",      &  ! Vertical grid: gdepw 
    529528            &          "m", ipk, gdepw_1d, nz_W, "down" ) 
     
    929928      INTEGER  ::   id_i , nz_i, nh_i        
    930929      INTEGER, DIMENSION(1) ::   idex             ! local workspace 
    931       REAL(wp) ::   zsto, zout, zmax, zjulian, zdt 
     930      REAL(wp) ::   zsto, zout, zmax, zjulian 
    932931      !!---------------------------------------------------------------------- 
    933932      !  
     
    940939      clname = cdfile_name 
    941940      IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname) 
    942       zdt  = rdt 
    943941      zsto = rdt 
    944942      clop = "inst(x)"           ! no use of the mask value (require less cpu time) 
    945943      zout = rdt 
    946       zmax = ( nitend - nit000 + 1 ) * zdt 
     944      zmax = ( nitend - nit000 + 1 ) * rdt 
    947945 
    948946      IF(lwp) WRITE(numout,*) 
     
    959957      zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    960958      CALL histbeg( clname, jpi, glamt, jpj, gphit,   & 
    961           1, jpi, 1, jpj, nit000-1, zjulian, zdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 
     959          1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 
    962960      CALL histvert( id_i, "deptht", "Vertical T levels",   &    ! Vertical grid : gdept 
    963961          "m", jpk, gdept_1d, nz_i, "down") 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r5282 r5825  
    8181      REAL(wp) ::   zu, z0u, zzwx, zw         ! local scalars 
    8282      REAL(wp) ::   zv, z0v, zzwy, z0w        !   -      - 
    83       REAL(wp) ::   ztra, zbtr, zdt, zalpha   !   -      - 
     83      REAL(wp) ::   ztra, zbtr, zalpha   !   -      - 
    8484      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zslpx, zslpy   ! 3D workspace 
    8585      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zwx  , zwy     ! -      -  
     
    176176         !                                             !-- MUSCL horizontal advective fluxes 
    177177         DO jk = 1, jpkm1                                     ! interior values 
    178             zdt  = p2dt 
    179178            DO jj = 2, jpjm1 
    180179               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    182181                  z0u = SIGN( 0.5, pun(ji,jj,jk) ) 
    183182                  zalpha = 0.5 - z0u 
    184                   zu  = z0u - 0.5 * pun(ji,jj,jk) * zdt / ( e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk) ) 
     183                  zu  = z0u - 0.5 * pun(ji,jj,jk) * p2dt / ( e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk) ) 
    185184                  zzwx = ptb(ji+1,jj,jk,jn) + xind(ji,jj,jk) * zu * zslpx(ji+1,jj,jk) 
    186185                  zzwy = ptb(ji  ,jj,jk,jn) + xind(ji,jj,jk) * zu * zslpx(ji  ,jj,jk) 
     
    189188                  z0v = SIGN( 0.5, pvn(ji,jj,jk) ) 
    190189                  zalpha = 0.5 - z0v 
    191                   zv  = z0v - 0.5 * pvn(ji,jj,jk) * zdt / ( e1v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,jk) ) 
     190                  zv  = z0v - 0.5 * pvn(ji,jj,jk) * p2dt / ( e1v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,jk) ) 
    192191                  zzwx = ptb(ji,jj+1,jk,jn) + xind(ji,jj,jk) * zv * zslpy(ji,jj+1,jk) 
    193192                  zzwy = ptb(ji,jj  ,jk,jn) + xind(ji,jj,jk) * zv * zslpy(ji,jj  ,jk) 
     
    259258         ! 
    260259         DO jk = 1, jpkm1                                     ! interior values 
    261             zdt  = p2dt 
    262             DO jj = 2, jpjm1       
     260            DO jj = 2, jpjm1 
    263261               DO ji = fs_2, fs_jpim1   ! vector opt. 
    264262                  zbtr = 1. / ( e1t(ji,jj) * e2t(ji,jj) * fse3w(ji,jj,jk+1) ) 
    265263                  z0w = SIGN( 0.5, pwn(ji,jj,jk+1) ) 
    266264                  zalpha = 0.5 + z0w 
    267                   zw  = z0w - 0.5 * pwn(ji,jj,jk+1) * zdt * zbtr  
     265                  zw  = z0w - 0.5 * pwn(ji,jj,jk+1) * p2dt * zbtr 
    268266                  zzwx = ptb(ji,jj,jk+1,jn) + xind(ji,jj,jk) * zw * zslpx(ji,jj,jk+1) 
    269267                  zzwy = ptb(ji,jj,jk  ,jn) + xind(ji,jj,jk) * zw * zslpx(ji,jj,jk  ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r5282 r5825  
    7171      REAL(wp) ::   zu, z0u, zzwx, zw         ! local scalars 
    7272      REAL(wp) ::   zv, z0v, zzwy, z0w        !   -      - 
    73       REAL(wp) ::   ztra, zbtr, zdt, zalpha   !   -      - 
     73      REAL(wp) ::   ztra, zbtr, zalpha   !   -      - 
    7474      REAL(wp), POINTER, DIMENSION(:,:,:) :: zslpx, zslpy , zwx, zwy 
    7575      !!---------------------------------------------------------------------- 
     
    133133        !                                             !-- MUSCL horizontal advective fluxes 
    134134         DO jk = 1, jpkm1                                     ! interior values 
    135             zdt  = p2dt 
    136135            DO jj = 2, jpjm1 
    137136               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    139138                  z0u = SIGN( 0.5, pun(ji,jj,jk) ) 
    140139                  zalpha = 0.5 - z0u 
    141                   zu  = z0u - 0.5 * pun(ji,jj,jk) * zdt / ( e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk) ) 
     140                  zu  = z0u - 0.5 * pun(ji,jj,jk) * p2dt / ( e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk) ) 
    142141                  zzwx = ptb(ji+1,jj,jk,jn) + zu * zslpx(ji+1,jj,jk) 
    143142                  zzwy = ptb(ji  ,jj,jk,jn) + zu * zslpx(ji  ,jj,jk) 
     
    146145                  z0v = SIGN( 0.5, pvn(ji,jj,jk) ) 
    147146                  zalpha = 0.5 - z0v 
    148                   zv  = z0v - 0.5 * pvn(ji,jj,jk) * zdt / ( e1v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,jk) ) 
     147                  zv  = z0v - 0.5 * pvn(ji,jj,jk) * p2dt / ( e1v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,jk) ) 
    149148                  zzwx = ptb(ji,jj+1,jk,jn) + zv * zslpy(ji,jj+1,jk) 
    150149                  zzwy = ptb(ji,jj  ,jk,jn) + zv * zslpy(ji,jj  ,jk) 
     
    240239         ! 
    241240         DO jk = 1, jpkm1                                     ! interior values 
    242             zdt  = p2dt 
    243241            DO jj = 2, jpjm1 
    244242               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    246244                  z0w = SIGN( 0.5, pwn(ji,jj,jk+1) ) 
    247245                  zalpha = 0.5 + z0w 
    248                   zw  = z0w - 0.5 * pwn(ji,jj,jk+1) * zdt * zbtr 
     246                  zw  = z0w - 0.5 * pwn(ji,jj,jk+1) * p2dt * zbtr 
    249247                  zzwx = ptb(ji,jj,jk+1,jn) + zw * zslpx(ji,jj,jk+1) 
    250248                  zzwy = ptb(ji,jj,jk  ,jn) + zw * zslpx(ji,jj,jk  ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90

    r5282 r5825  
    131131      !! 
    132132      INTEGER  :: ji, jj, jk, jn   ! dummy loop indices 
    133       REAL(wp) :: ztra, zbtr, zdir, zdx, zdt, zmsk   ! local scalars 
     133      REAL(wp) :: ztra, zbtr, zdir, zdx, zmsk   ! local scalars 
    134134      REAL(wp), POINTER, DIMENSION(:,:,:) :: zwx, zfu, zfc, zfd 
    135135      !---------------------------------------------------------------------- 
     
    170170         ! 
    171171         DO jk = 1, jpkm1   
    172             zdt =  p2dt 
    173172            DO jj = 2, jpjm1 
    174173               DO ji = fs_2, fs_jpim1   ! vector opt.    
    175174                  zdir = 0.5 + SIGN( 0.5, pun(ji,jj,jk) )   ! if pun > 0 : zdir = 1 otherwise zdir = 0  
    176175                  zdx = ( zdir * e1t(ji,jj) + ( 1. - zdir ) * e1t(ji+1,jj) ) * e2u(ji,jj) * fse3u(ji,jj,jk) 
    177                   zwx(ji,jj,jk)  = ABS( pun(ji,jj,jk) ) * zdt / zdx    ! (0<zc_cfl<1 : Courant number on x-direction) 
     176                  zwx(ji,jj,jk)  = ABS( pun(ji,jj,jk) ) * p2dt / zdx    ! (0<zc_cfl<1 : Courant number on x-direction) 
    178177                  zfc(ji,jj,jk)  = zdir * ptb(ji  ,jj,jk,jn) + ( 1. - zdir ) * ptb(ji+1,jj,jk,jn)  ! FC in the x-direction for T 
    179178                  zfd(ji,jj,jk)  = zdir * ptb(ji+1,jj,jk,jn) + ( 1. - zdir ) * ptb(ji  ,jj,jk,jn)  ! FD in the x-direction for T 
     
    252251      !! 
    253252      INTEGER  :: ji, jj, jk, jn   ! dummy loop indices 
    254       REAL(wp) :: ztra, zbtr, zdir, zdx, zdt, zmsk   ! local scalars 
     253      REAL(wp) :: ztra, zbtr, zdir, zdx, zmsk   ! local scalars 
    255254      REAL(wp), POINTER, DIMENSION(:,:,:) :: zwy, zfu, zfc, zfd 
    256255      !---------------------------------------------------------------------- 
     
    293292         ! 
    294293         DO jk = 1, jpkm1   
    295             zdt =  p2dt 
    296294            DO jj = 2, jpjm1 
    297295               DO ji = fs_2, fs_jpim1   ! vector opt.    
    298296                  zdir = 0.5 + SIGN( 0.5, pvn(ji,jj,jk) )   ! if pun > 0 : zdir = 1 otherwise zdir = 0  
    299297                  zdx = ( zdir * e2t(ji,jj) + ( 1. - zdir ) * e2t(ji,jj+1) ) * e1v(ji,jj) * fse3v(ji,jj,jk) 
    300                   zwy(ji,jj,jk)  = ABS( pvn(ji,jj,jk) ) * zdt / zdx    ! (0<zc_cfl<1 : Courant number on x-direction) 
     298                  zwy(ji,jj,jk)  = ABS( pvn(ji,jj,jk) ) * p2dt / zdx    ! (0<zc_cfl<1 : Courant number on x-direction) 
    301299                  zfc(ji,jj,jk)  = zdir * ptb(ji,jj  ,jk,jn) + ( 1. - zdir ) * ptb(ji,jj+1,jk,jn)  ! FC in the x-direction for T 
    302300                  zfd(ji,jj,jk)  = zdir * ptb(ji,jj+1,jk,jn) + ( 1. - zdir ) * ptb(ji,jj  ,jk,jn)  ! FD in the x-direction for T 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r5282 r5825  
    8181      INTEGER  ::   ji, jj, jk, jn           ! dummy loop indices 
    8282      INTEGER  ::   ik   
    83       REAL(wp) ::   z2dtt, zbtr, ztra        ! local scalar 
     83      REAL(wp) ::   zbtr, ztra        ! local scalar 
    8484      REAL(wp) ::   zfp_ui, zfp_vj, zfp_wk   !   -      - 
    8585      REAL(wp) ::   zfm_ui, zfm_vj, zfm_wk   !   -      - 
     
    161161         ! total advective trend 
    162162         DO jk = 1, jpkm1 
    163             z2dtt = p2dt 
    164163            DO jj = 2, jpjm1 
    165164               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    171170                  ! update and guess with monotonic sheme 
    172171                  pta(ji,jj,jk,jn) =   pta(ji,jj,jk,jn)         + ztra   * tmask(ji,jj,jk) 
    173                   zwi(ji,jj,jk)    = ( ptb(ji,jj,jk,jn) + z2dtt * ztra ) * tmask(ji,jj,jk) 
     172                  zwi(ji,jj,jk)    = ( ptb(ji,jj,jk,jn) + p2dt * ztra ) * tmask(ji,jj,jk) 
    174173               END DO 
    175174            END DO 
     
    299298      INTEGER  ::   jtaken          ! toggle for collecting appropriate fluxes from sub timesteps 
    300299      REAL(wp) ::   z_rzts          ! Fractional length of Euler forward sub-timestep for vertical advection 
    301       REAL(wp) ::   z2dtt, zbtr, ztra        ! local scalar 
     300      REAL(wp) ::   zbtr, ztra        ! local scalar 
    302301      REAL(wp) ::   zfp_ui, zfp_vj, zfp_wk   !   -      - 
    303302      REAL(wp) ::   zfm_ui, zfm_vj, zfm_wk   !   -      - 
     
    375374         ! total advective trend 
    376375         DO jk = 1, jpkm1 
    377             z2dtt = p2dt 
    378376            DO jj = 2, jpjm1 
    379377               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    385383                  ! update and guess with monotonic sheme 
    386384                  pta(ji,jj,jk,jn) =   pta(ji,jj,jk,jn)         + ztra 
    387                   zwi(ji,jj,jk)    = ( ptb(ji,jj,jk,jn) + z2dtt * ztra ) * tmask(ji,jj,jk) 
     385                  zwi(ji,jj,jk)    = ( ptb(ji,jj,jk,jn) + p2dt * ztra ) * tmask(ji,jj,jk) 
    388386               END DO 
    389387            END DO 
     
    559557      INTEGER  ::   ji, jj, jk   ! dummy loop indices 
    560558      INTEGER  ::   ikm1         ! local integer 
    561       REAL(wp) ::   zpos, zneg, zbt, za, zb, zc, zbig, zrtrn, z2dtt   ! local scalars 
     559      REAL(wp) ::   zpos, zneg, zbt, za, zb, zc, zbig, zrtrn   ! local scalars 
    562560      REAL(wp) ::   zau, zbu, zcu, zav, zbv, zcv, zup, zdo            !   -      - 
    563561      REAL(wp), POINTER, DIMENSION(:,:,:) :: zbetup, zbetdo, zbup, zbdo 
     
    584582            DO jk = mikt(ji,jj), jpkm1 
    585583               ikm1 = MAX(jk-1,mikt(ji,jj)) 
    586                z2dtt = p2dt 
    587584                
    588585               ! search maximum in neighbourhood 
     
    609606 
    610607               ! up & down beta terms 
    611                zbt = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) / z2dtt 
     608               zbt = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) / p2dt 
    612609               zbetup(ji,jj,jk) = ( zup            - paft(ji,jj,jk) ) / ( zpos + zrtrn ) * zbt 
    613610               zbetdo(ji,jj,jk) = ( paft(ji,jj,jk) - zdo            ) / ( zneg + zrtrn ) * zbt 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_ubs.F90

    r5282 r5825  
    8787      ! 
    8888      INTEGER  ::   ji, jj, jk, jn   ! dummy loop indices 
    89       REAL(wp) ::   ztra, zbtr, zcoef, z2dtt                       ! local scalars 
     89      REAL(wp) ::   ztra, zbtr, zcoef                      ! local scalars 
    9090      REAL(wp) ::   zfp_ui, zfm_ui, zcenut, ztak, zfp_wk, zfm_wk   !   -      - 
    9191      REAL(wp) ::   zfp_vj, zfm_vj, zcenvt, zeeu, zeev, z_hdivn    !   -      - 
     
    207207         ! update and guess with monotonic sheme 
    208208         DO jk = 1, jpkm1 
    209             z2dtt = p2dt 
    210209            DO jj = 2, jpjm1 
    211210               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    213212                  ztak = - ( ztw(ji,jj,jk) - ztw(ji,jj,jk+1) ) * zbtr 
    214213                  pta(ji,jj,jk,jn) =   pta(ji,jj,jk,jn) +  ztak  
    215                   zti(ji,jj,jk)    = ( ptb(ji,jj,jk,jn) + z2dtt * ( ztak + zltu(ji,jj,jk) ) ) * tmask(ji,jj,jk) 
     214                  zti(ji,jj,jk)    = ( ptb(ji,jj,jk,jn) + p2dt * ( ztak + zltu(ji,jj,jk) ) ) * tmask(ji,jj,jk) 
    216215               END DO 
    217216            END DO 
     
    288287      INTEGER  ::   ji, jj, jk   ! dummy loop indices 
    289288      INTEGER  ::   ikm1         ! local integer 
    290       REAL(wp) ::   zpos, zneg, zbt, za, zb, zc, zbig, zrtrn, z2dtt   ! local scalars 
     289      REAL(wp) ::   zpos, zneg, zbt, za, zb, zc, zbig, zrtrn   ! local scalars 
    291290      REAL(wp), POINTER, DIMENSION(:,:,:) :: zbetup, zbetdo 
    292291      !!---------------------------------------------------------------------- 
     
    340339 
    341340      DO jk = 1, jpkm1 
    342          z2dtt = p2dt 
    343341         DO jj = 2, jpjm1 
    344342            DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    347345               zneg = MAX( 0., pcc(ji  ,jj  ,jk  ) ) - MIN( 0., pcc(ji  ,jj  ,jk+1) ) 
    348346               ! up & down beta terms 
    349                zbt = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) / z2dtt 
     347               zbt = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) / p2dt 
    350348               zbetup(ji,jj,jk) = ( zbetup(ji,jj,jk) - paft(ji,jj,jk) ) / (zpos+zrtrn) * zbt 
    351349               zbetdo(ji,jj,jk) = ( paft(ji,jj,jk) - zbetdo(ji,jj,jk) ) / (zneg+zrtrn) * zbt 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/SAS_SRC/diawri.F90

    r5282 r5825  
    136136      INTEGER  ::   ierr                                     ! error code return from allocation 
    137137      INTEGER  ::   iimi, iima, ipk, it, itmod, ijmi, ijma   ! local integers 
    138       REAL(wp) ::   zsto, zout, zmax, zjulian, zdt           ! local scalars 
     138      REAL(wp) ::   zsto, zout, zmax, zjulian                ! local scalars 
    139139      !!---------------------------------------------------------------------- 
    140140      !  
     
    155155 
    156156      ! Define frequency of output and means 
    157       zdt = rdt 
    158157      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    159158      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
    160159      ENDIF 
    161160#if defined key_diainstant 
    162       zsto = nwrite * zdt 
     161      zsto = nwrite * rdt 
    163162      clop = "inst("//TRIM(clop)//")" 
    164163#else 
    165       zsto=zdt 
     164      zsto=rdt 
    166165      clop = "ave("//TRIM(clop)//")" 
    167166#endif 
    168       zout = nwrite * zdt 
    169       zmax = ( nitend - nit000 + 1 ) * zdt 
     167      zout = nwrite * rdt 
     168      zmax = ( nitend - nit000 + 1 ) * rdt 
    170169 
    171170      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    209208         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    210209            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    211             &          nit000-1, zjulian, zdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
     210            &          nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
    212211         CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept 
    213212            &           "m", ipk, gdept_1d, nz_T, "down" ) 
     
    221220         CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu,           &  ! Horizontal grid: glamu and gphiu 
    222221            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    223             &          nit000-1, zjulian, zdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
     222            &          nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
    224223         CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept 
    225224            &           "m", ipk, gdept_1d, nz_U, "down" ) 
     
    233232         CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv,           &  ! Horizontal grid: glamv and gphiv 
    234233            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    235             &          nit000-1, zjulian, zdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
     234            &          nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
    236235         CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept 
    237236            &          "m", ipk, gdept_1d, nz_V, "down" ) 
     
    353352      INTEGER  ::   id_i , nz_i, nh_i        
    354353      INTEGER, DIMENSION(1) ::   idex             ! local workspace 
    355       REAL(wp) ::   zsto, zout, zmax, zjulian, zdt 
     354      REAL(wp) ::   zsto, zout, zmax, zjulian 
    356355      !!---------------------------------------------------------------------- 
    357356      !  
     
    364363      clname = cdfile_name 
    365364      IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname) 
    366       zdt  = rdt 
    367365      zsto = rdt 
    368366      clop = "inst(x)"           ! no use of the mask value (require less cpu time) 
    369367      zout = rdt 
    370       zmax = ( nitend - nit000 + 1 ) * zdt 
     368      zmax = ( nitend - nit000 + 1 ) * rdt 
    371369 
    372370      IF(lwp) WRITE(numout,*) 
     
    383381      zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    384382      CALL histbeg( clname, jpi, glamt, jpj, gphit,   & 
    385           1, jpi, 1, jpj, nit000-1, zjulian, zdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 
     383          1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_i, id_i, domain_id=nidom, snc4chunks=snc4set ) ! Horizontal grid : glamt and gphit 
    386384      CALL histvert( id_i, "deptht", "Vertical T levels",   &    ! Vertical grid : gdept 
    387385          "m", jpk, gdept_1d, nz_i, "down") 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/trcdia.F90

    r4292 r5825  
    108108      CHARACTER (len=20) :: cltra, cltrau 
    109109      CHARACTER (len=80) :: cltral 
    110       REAL(wp) :: zsto, zout, zdt 
     110      REAL(wp) :: zsto, zout 
    111111      INTEGER  :: iimi, iima, ijmi, ijma, ipk, it, itmod, iiter 
    112112      !!---------------------------------------------------------------------- 
     
    120120 
    121121      ! Define frequency of output and means 
    122       zdt = rdt 
    123122      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    124123      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     
    128127      clop = "inst("//TRIM(clop)//")" 
    129128# else 
    130       zsto = zdt 
     129      zsto = rdt 
    131130      clop = "ave("//TRIM(clop)//")" 
    132131# endif 
    133       zout = nn_writetrc * zdt 
     132      zout = nn_writetrc * rdt 
    134133 
    135134      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    184183         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,     & 
    185184            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,         &  
    186             &          iiter, zjulian, zdt, nhorit5, nit5 , domain_id=nidom, snc4chunks=snc4set) 
     185            &          iiter, zjulian, rdt, nhorit5, nit5 , domain_id=nidom, snc4chunks=snc4set) 
    187186 
    188187         ! Vertical grid for tracer : gdept 
     
    252251      INTEGER  ::   jl 
    253252      INTEGER  ::   iimi, iima, ijmi, ijma, ipk, it, itmod, iiter 
    254       REAL(wp) ::   zsto, zout, zdt 
     253      REAL(wp) ::   zsto, zout 
    255254      !!---------------------------------------------------------------------- 
    256255 
     
    263262      ! 
    264263      ! Define frequency of output and means 
    265       zdt = rdt 
    266264      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    267265      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
    268266      ENDIF 
    269267#  if defined key_diainstant 
    270       zsto = nn_writedia * zdt 
     268      zsto = nn_writedia * rdt 
    271269      clop = "inst("//TRIM(clop)//")" 
    272270#  else 
    273       zsto = zdt 
     271      zsto = rdt 
    274272      clop = "ave("//TRIM(clop)//")" 
    275273#  endif 
    276       zout = nn_writedia * zdt 
     274      zout = nn_writedia * rdt 
    277275 
    278276      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    304302         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,             & 
    305303            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,         & 
    306             &          iiter, zjulian, zdt, nhoritd, nitd , domain_id=nidom, snc4chunks=snc4set ) 
     304            &          iiter, zjulian, rdt, nhoritd, nitd , domain_id=nidom, snc4chunks=snc4set ) 
    307305 
    308306         ! Vertical grid for 2d and 3d arrays 
     
    389387      INTEGER  ::   ji, jj, jk, jl 
    390388      INTEGER  ::   iimi, iima, ijmi, ijma, ipk, it, itmod, iiter 
    391       REAL(wp) ::   zsto, zout, zdt 
     389      REAL(wp) ::   zsto, zout 
    392390      !!---------------------------------------------------------------------- 
    393391 
     
    400398 
    401399      ! Define frequency of output and means 
    402       zdt = rdt 
    403400      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    404401      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
    405402      ENDIF 
    406403#        if defined key_diainstant 
    407       zsto = nn_writebio * zdt 
     404      zsto = nn_writebio * rdt 
    408405      clop = "inst("//TRIM(clop)//")" 
    409406#        else 
    410       zsto = zdt 
     407      zsto = rdt 
    411408      clop = "ave("//TRIM(clop)//")" 
    412409#        endif 
    413       zout = nn_writebio * zdt 
     410      zout = nn_writebio * rdt 
    414411 
    415412      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    437434         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,      & 
    438435            &    iimi, iima-iimi+1, ijmi, ijma-ijmi+1,          & 
    439             &    iiter, zjulian, zdt, nhoritb, nitb , domain_id=nidom, snc4chunks=snc4set ) 
     436            &    iiter, zjulian, rdt, nhoritb, nitb , domain_id=nidom, snc4chunks=snc4set ) 
    440437         ! Vertical grid for biological trends 
    441438         CALL histvert(nitb, 'deptht', 'Vertical T levels', 'm', ipk, gdept_1d, ndepitb) 
Note: See TracChangeset for help on using the changeset viewer.