Changeset 11093
 Timestamp:
 20190611T14:58:46+02:00 (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

NEMO/releases/release4.0/src/OCE/DIA/diahth.F90
r11091 r11093 89 89 REAL(wp), DIMENSION(jpi,jpj) :: zhtc3 !: heat content of first 300 m [W] 90 90 91 <<<<<<< variant A 91 92 IF (iom_use("mlddzt").OR.iom_use("mldr0_3").OR.iom_use("mldr0_1")) THEN 92 93 !  ! … … 98 99 zmaxdzT(:,:) = 0._wp 99 100 IF( nla10 > 1 ) THEN 101 >>>>>>> variant B 102 IF(iom_use("mlddzt").OR.iom_use("mldr0_3").OR.iom_use("mldr0_1").OR.iom_use("mld_dt02") & 103 & .OR.iom_use("topthdep").OR.iom_use("mldr10_3").OR.iom_use("pycndep").OR.iom_use("tinv").OR.iom_use("depti")) THEN 104 ! initialization 105 ztinv (:,:) = 0._wp 106 zdepinv(:,:) = 0._wp 107 zmaxdzT(:,:) = 0._wp 108 DO jj = 1, jpj 109 DO ji = 1, jpi 110 zztmp = gdepw_n(ji,jj,mbkt(ji,jj)+1) 111 zhth (ji,jj) = zztmp 112 zabs2 (ji,jj) = zztmp 113 ztm2 (ji,jj) = zztmp 114 zrho10_3(ji,jj) = zztmp 115 zpycn (ji,jj) = zztmp 116 END DO 117 END DO 118 IF( nla10 > 1 ) THEN 119 ======= end 100 120 DO jj = 1, jpj 101 121 DO ji = 1, jpi … … 103 123 zrho0_3(ji,jj) = zztmp 104 124 zrho0_1(ji,jj) = zztmp 125 <<<<<<< variant A 105 126 zhth(ji,jj) = zztmp 106 END DO 107 END DO 127 >>>>>>> variant B 128 ======= end 129 END DO 130 END DO 131 <<<<<<< variant A 108 132 ELSE IF (iom_use("mlddzt")) THEN 109 133 DO jj = 1, jpj … … 115 139 ELSE 116 140 zhth(:,:) = 0._wp 141 >>>>>>> variant B 142 ======= end 117 143 ENDIF 118 144 <<<<<<< variant A 145 >>>>>>> variant B 146 ENDIF 147 ======= end 148 149 <<<<<<< variant A 150 >>>>>>> variant B 151 IF (iom_use("mlddzt").OR.iom_use("mldr0_3").OR.iom_use("mldr0_1")) THEN 152 !  ! 153 ! thermocline depth: strongest vertical gradient of temperature ! 154 ! turbocline depth (mixing layer depth): avt = zavt5 ! 155 ! MLD: rho = rho(1) + zrho3 ! 156 ! MLD: rho = rho(1) + zrho1 ! 157 !  ! 158 ======= end 119 159 DO jk = jpkm1, 2, 1 ! loop from bottom to 2 120 160 DO jj = 1, jpj … … 125 165 zzdep = zzdep * tmask(ji,jj,1) 126 166 167 <<<<<<< variant A 127 168 IF( zztmp > zmaxdzT(ji,jj) ) THEN 169 >>>>>>> variant B 170 IF( zztmp > zmaxdzT(ji,jj) ) THEN 171 ======= end 128 172 zmaxdzT(ji,jj) = zztmp ; zhth (ji,jj) = zzdep ! max and depth of dT/dz 129 173 ENDIF 130 174 175 <<<<<<< variant A 131 176 IF( nla10 > 1 ) THEN 177 >>>>>>> variant B 178 IF( nla10 > 1 ) THEN 179 ======= end 132 180 zztmp = rhop(ji,jj,jk)  rhop(ji,jj,1) ! delta rho(1) 133 181 IF( zztmp > zrho3 ) zrho0_3(ji,jj) = zzdep ! > 0.03 … … 140 188 141 189 IF (iom_use("mlddzt")) CALL iom_put( "mlddzt", zhth*tmask(:,:,1) ) ! depth of the thermocline 190 <<<<<<< variant A 142 191 IF( nla10 > 1 ) THEN 192 >>>>>>> variant B 193 IF( nla10 > 1 ) THEN 194 ======= end 143 195 IF (iom_use("mldr0_3")) CALL iom_put( "mldr0_3", zrho0_3*tmask(:,:,1) ) ! MLD delta rho(surf) = 0.03 144 196 IF (iom_use("mldr0_1")) CALL iom_put( "mldr0_1", zrho0_1*tmask(:,:,1) ) ! MLD delta rho(surf) = 0.01 … … 308 360 ENDIF 309 361 ! 362 IF( ln_timing ) CALL timing_stop('dia_hth') 363 ! 310 364 END SUBROUTINE dia_hth 311 365
Note: See TracChangeset
for help on using the changeset viewer.