Ignore:
Timestamp:
2020-02-28T16:55:11+01:00 (14 months ago)
Author:
davestorkey
Message:

Preparation for new timestepping scheme #2390.
Main changes:

  1. Initial euler timestep now handled in stp and not in TRA/DYN routines.
  2. Renaming of all timestep parameters. In summary, the namelist parameter is now rn_Dt and the current timestep is rDt (and rDt_ice, rDt_trc etc).
  3. Renaming of a few miscellaneous parameters, eg. atfp → rn_atfp (namelist parameter used everywhere) and rau0 → rho0.

This version gives bit-comparable results to the previous version of the trunk.

Location:
NEMO/trunk/src/OCE/LDF
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/LDF/ldfdyn.F90

    r12377 r12489  
    407407            zcmsmag   = (rn_csmc/rpi)**2                                            ! (C_smag/pi)^2 
    408408            zstabf_lo = rn_minfac * rn_minfac / ( 2._wp * 12._wp * 12._wp * zcmsmag ) ! lower limit stability factor scaling 
    409             zstabf_up = rn_maxfac / ( 4._wp * zcmsmag * 2._wp * rdt )               ! upper limit stability factor scaling 
     409            zstabf_up = rn_maxfac / ( 4._wp * zcmsmag * 2._wp * rn_Dt )               ! upper limit stability factor scaling 
    410410            IF( ln_dynldf_blp ) zstabf_lo = ( 16._wp / 9._wp ) * zstabf_lo          ! provide |U|L^3/12 lower limit instead  
    411411            !                                                                       ! of |U|L^3/16 in blp case 
  • NEMO/trunk/src/OCE/LDF/ldftra.F90

    r12377 r12489  
    820820      ! 
    821821      IF( iom_use('weiv_masstr') ) THEN   ! vertical mass transport & its square value 
    822          zw2d(:,:) = rau0 * e1e2t(:,:) 
     822         zw2d(:,:) = rho0 * e1e2t(:,:) 
    823823         DO jk = 1, jpk 
    824824            zw3d(:,:,jk) = zw3d(:,:,jk) * zw2d(:,:) 
     
    830830         zw3d(:,:,:) = 0.e0 
    831831         DO jk = 1, jpkm1 
    832             zw3d(:,:,jk) = rau0 * ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) )  
     832            zw3d(:,:,jk) = rho0 * ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) )  
    833833         END DO 
    834834         CALL iom_put( "ueiv_masstr", zw3d )                  ! mass transport in i-direction 
    835835      ENDIF 
    836836      ! 
    837       zztmp = 0.5_wp * rau0 * rcp  
     837      zztmp = 0.5_wp * rho0 * rcp  
    838838      IF( iom_use('ueiv_heattr') .OR. iom_use('ueiv_heattr3d') ) THEN 
    839839        zw2d(:,:)   = 0._wp  
     
    853853         zw3d(:,:,:) = 0.e0 
    854854         DO jk = 1, jpkm1 
    855             zw3d(:,:,jk) = rau0 * ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) )  
     855            zw3d(:,:,jk) = rho0 * ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) )  
    856856         END DO 
    857857         CALL iom_put( "veiv_masstr", zw3d )                  ! mass transport in i-direction 
Note: See TracChangeset for help on using the changeset viewer.