- Timestamp:
- 2010-05-04T10:39:48+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/DYN/dynzdf_exp.F90
r1537 r1858 54 54 !! ** Action : - Update (ua,va) with the vertical diffusive trend 55 55 !!--------------------------------------------------------------------- 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 64 62 !!---------------------------------------------------------------------- 65 63 … … 72 70 ! Local constant initialization 73 71 ! ----------------------------- 74 zrau0r = 1. / rau0 ! inverse of the reference density75 72 zlavmr = 1. / float( nn_zdfexp ) ! inverse of the number of sub time step 76 73 … … 81 78 ! Surface boundary condition 82 79 DO ji = 2, jpim1 83 zwy(ji,1) = utau(ji,jj) * zrau0r84 zww(ji,1) = vtau(ji,jj) * zrau0r80 zwy(ji,1) = utau(ji,jj) * r1_rau0 81 zww(ji,1) = vtau(ji,jj) * r1_rau0 85 82 END DO 86 83
Note: See TracChangeset
for help on using the changeset viewer.