- Timestamp:
- 2016-02-24T20:04:41+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90
r6204 r6353 80 80 INTEGER, INTENT(in) :: kt ! ocean time-step index 81 81 ! 82 INTEGER :: ji, jj, jk ! dummy loop indices83 INTEGER :: iikn, iiki, ikt , imkt! local integer84 REAL(wp) :: zN2_c ! local scalar82 INTEGER :: ji, jj, jk ! dummy loop indices 83 INTEGER :: iikn, iiki, ikt ! local integer 84 REAL(wp) :: zN2_c ! local scalar 85 85 INTEGER, POINTER, DIMENSION(:,:) :: imld ! 2D workspace 86 86 !!---------------------------------------------------------------------- … … 117 117 DO jj = 1, jpj 118 118 DO ji = 1, jpi 119 imkt = mikt(ji,jj) 120 IF( avt (ji,jj,jk) < avt_c ) imld(ji,jj) = MAX( imkt, jk ) ! Turbocline 119 IF( avt (ji,jj,jk) < avt_c * wmask(ji,jj,jk) ) imld(ji,jj) = jk ! Turbocline 121 120 END DO 122 121 END DO … … 127 126 iiki = imld(ji,jj) 128 127 iikn = nmln(ji,jj) 129 imkt = mikt(ji,jj) 130 hmld (ji,jj) = ( fsdepw(ji,jj,iiki ) - fsdepw(ji,jj,imkt ) ) * ssmask(ji,jj) ! Turbocline depth 131 hmlp (ji,jj) = ( fsdepw(ji,jj,iikn ) - fsdepw(ji,jj,imkt ) ) * ssmask(ji,jj) ! Mixed layer depth 132 hmlpt(ji,jj) = ( fsdept(ji,jj,iikn-1) - fsdepw(ji,jj,imkt ) ) * ssmask(ji,jj) ! depth of the last T-point inside the mixed layer 128 hmld (ji,jj) = fsdepw(ji,jj,iiki ) * ssmask(ji,jj) ! Turbocline depth 129 hmlp (ji,jj) = fsdepw(ji,jj,iikn ) * ssmask(ji,jj) ! Mixed layer depth 130 hmlpt(ji,jj) = fsdept(ji,jj,iikn-1) * ssmask(ji,jj) ! depth of the last T-point inside the mixed layer 133 131 END DO 134 132 END DO 135 IF( .NOT.lk_offline ) THEN ! no need to output in offline mode 136 CALL iom_put( "mldr10_1", hmlp ) ! mixed layer depth 137 CALL iom_put( "mldkz5" , hmld ) ! turbocline depth 133 ! no need to output in offline mode 134 IF( .NOT.lk_offline ) THEN 135 IF ( iom_use("mldr10_1") ) THEN 136 IF( ln_isfcav ) THEN 137 CALL iom_put( "mldr10_1", hmlp - risfdep) ! mixed layer thickness 138 ELSE 139 CALL iom_put( "mldr10_1", hmlp ) ! mixed layer depth 140 END IF 141 END IF 142 IF ( iom_use("mldkz5") ) THEN 143 IF( ln_isfcav ) THEN 144 CALL iom_put( "mldkz5" , hmld - risfdep ) ! turbocline thickness 145 ELSE 146 CALL iom_put( "mldkz5" , hmld ) ! turbocline depth 147 END IF 148 END IF 138 149 ENDIF 139 150
Note: See TracChangeset
for help on using the changeset viewer.