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 1858 for branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/DYN/dynzdf_exp.F90 – NEMO

Ignore:
Timestamp:
2010-05-04T10:39:48+02:00 (14 years ago)
Author:
gm
Message:

ticket:#665 : step 1 - heat content of freezing-melting ice

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/DYN/dynzdf_exp.F90

    r1537 r1858  
    5454      !! ** Action : - Update (ua,va) with the vertical diffusive trend 
    5555      !!--------------------------------------------------------------------- 
    56       !! * Arguments 
    57       INTEGER , INTENT( in ) ::   kt                           ! ocean time-step index 
    58       REAL(wp), INTENT( in ) ::   p2dt                         ! time-step  
    59  
    60       !! * Local declarations 
    61       INTEGER ::   ji, jj, jk, jl                              ! dummy loop indices 
    62       REAL(wp) ::   zrau0r, zlavmr, zua, zva                   ! temporary scalars 
    63       REAL(wp), DIMENSION(jpi,jpk) ::   zwx, zwy, zwz, zww     ! temporary workspace arrays 
     56      INTEGER , INTENT( in ) ::   kt     ! ocean time-step index 
     57      REAL(wp), INTENT( in ) ::   p2dt   ! time-step  
     58      !!  
     59      INTEGER  ::   ji, jj, jk, jl     ! dummy loop indices 
     60      REAL(wp) ::   zlavmr, zua, zva   ! temporary scalars 
     61      REAL(wp), DIMENSION(jpi,jpk) ::   zwx, zwy, zwz, zww   ! temporary workspace arrays 
    6462      !!---------------------------------------------------------------------- 
    6563 
     
    7270      ! Local constant initialization 
    7371      ! ----------------------------- 
    74       zrau0r = 1. / rau0                                   ! inverse of the reference density 
    7572      zlavmr = 1. / float( nn_zdfexp )                      ! inverse of the number of sub time step 
    7673 
     
    8178         ! Surface boundary condition 
    8279         DO ji = 2, jpim1 
    83             zwy(ji,1) = utau(ji,jj) * zrau0r 
    84             zww(ji,1) = vtau(ji,jj) * zrau0r 
     80            zwy(ji,1) = utau(ji,jj) * r1_rau0 
     81            zww(ji,1) = vtau(ji,jj) * r1_rau0 
    8582         END DO   
    8683 
Note: See TracChangeset for help on using the changeset viewer.