Changeset 6352 for trunk/NEMOGCM/NEMO/OPA_SRC/ZDF
- Timestamp:
- 2016-02-24T19:52:54+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90
r6140 r6352 31 31 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hmld !: mixing layer depth (turbocline) [m] 32 32 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hmlp !: mixed layer depth (rho=rho0+zdcrit) [m] 33 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hmlpt !: depth of the last T-point inside the mixed layer 33 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hmlpt !: depth of the last T-point inside the mixed layer [m] 34 34 35 35 REAL(wp), PUBLIC :: rho_c = 0.01_wp !: density criterion for mixed layer depth … … 77 77 INTEGER, INTENT(in) :: kt ! ocean time-step index 78 78 ! 79 INTEGER :: ji, jj, jk ! dummy loop indices80 INTEGER :: iikn, iiki, ikt , imkt! local integer81 REAL(wp) :: zN2_c ! local scalar79 INTEGER :: ji, jj, jk ! dummy loop indices 80 INTEGER :: iikn, iiki, ikt ! local integer 81 REAL(wp) :: zN2_c ! local scalar 82 82 INTEGER, POINTER, DIMENSION(:,:) :: imld ! 2D workspace 83 83 !!---------------------------------------------------------------------- … … 123 123 iiki = imld(ji,jj) 124 124 iikn = nmln(ji,jj) 125 imkt = mikt(ji,jj)126 125 hmld (ji,jj) = gdepw_n(ji,jj,iiki ) * ssmask(ji,jj) ! Turbocline depth 127 126 hmlp (ji,jj) = gdepw_n(ji,jj,iikn ) * ssmask(ji,jj) ! Mixed layer depth … … 129 128 END DO 130 129 END DO 131 IF( .NOT.lk_offline ) THEN ! no need to output in offline mode 130 ! no need to output in offline mode 131 IF( .NOT.lk_offline ) THEN 132 132 IF ( iom_use("mldr10_1") ) THEN 133 IF( .NOT. ln_isfcav ) CALL iom_put( "mldr10_1", hmlp ) ! mixed layer depth 134 IF( ln_isfcav ) CALL iom_put( "mldr10_1", hmlp - risfdep) ! mixed layer thickness 133 IF( ln_isfcav ) THEN 134 CALL iom_put( "mldr10_1", hmlp - risfdep) ! mixed layer thickness 135 ELSE 136 CALL iom_put( "mldr10_1", hmlp ) ! mixed layer depth 137 END IF 135 138 END IF 136 139 IF ( iom_use("mldkz5") ) THEN 137 IF( .NOT. ln_isfcav ) CALL iom_put( "mldkz5" , hmld ) ! turbocline depth 138 IF( ln_isfcav ) CALL iom_put( "mldkz5" , hmld - risfdep ) ! turbocline thickness 140 IF( ln_isfcav ) THEN 141 CALL iom_put( "mldkz5" , hmld - risfdep ) ! turbocline thickness 142 ELSE 143 CALL iom_put( "mldkz5" , hmld ) ! turbocline depth 144 END IF 139 145 END IF 140 146 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.