- Timestamp:
- 2011-10-22T15:46:41+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/OPA_SRC/DIA/diahth.F90
r2715 r2977 160 160 DO ji = 1, jpi 161 161 IF( tmask(ji,jj,nla10) == 1. ) THEN 162 zu = 1779.50 + 11.250*tn(ji,jj,nla10) - 3.80*sn(ji,jj,nla10) - 0.0745*tn(ji,jj,nla10)*tn(ji,jj,nla10) & 163 & - 0.0100*tn(ji,jj,nla10)*sn(ji,jj,nla10) 164 zv = 5891.00 + 38.000*tn(ji,jj,nla10) + 3.00*sn(ji,jj,nla10) - 0.3750*tn(ji,jj,nla10)*tn(ji,jj,nla10) 165 zut = 11.25 - 0.149*tn(ji,jj,nla10) - 0.01*sn(ji,jj,nla10) 166 zvt = 38.00 - 0.750*tn(ji,jj,nla10) 162 zu = 1779.50 + 11.250 * tsn(ji,jj,nla10,jp_tem) - 3.80 * tsn(ji,jj,nla10,jp_sal) & 163 & - 0.0745 * tsn(ji,jj,nla10,jp_tem) * tsn(ji,jj,nla10,jp_tem) & 164 & - 0.0100 * tsn(ji,jj,nla10,jp_tem) * tsn(ji,jj,nla10,jp_sal) 165 zv = 5891.00 + 38.000 * tsn(ji,jj,nla10,jp_tem) + 3.00 * tsn(ji,jj,nla10,jp_sal) & 166 & - 0.3750 * tsn(ji,jj,nla10,jp_tem) * tsn(ji,jj,nla10,jp_tem) 167 zut = 11.25 - 0.149 * tsn(ji,jj,nla10,jp_tem) - 0.01 * tsn(ji,jj,nla10,jp_sal) 168 zvt = 38.00 - 0.750 * tsn(ji,jj,nla10,jp_tem) 167 169 zw = (zu + 0.698*zv) * (zu + 0.698*zv) 168 170 zdelr(ji,jj) = ztem2 * (1000.*(zut*zv - zvt*zu)/zw) … … 184 186 ! 185 187 zzdep = fsdepw(ji,jj,jk) 186 zztmp = ( t n(ji,jj,jk-1) - tn(ji,jj,jk) ) / zzdep * tmask(ji,jj,jk) ! vertical gradient of temperature (dT/dz)188 zztmp = ( tsn(ji,jj,jk-1,jp_tem) - tsn(ji,jj,jk,jp_tem) ) / zzdep * tmask(ji,jj,jk) ! vertical gradient of temperature (dT/dz) 187 189 zzdep = zzdep * tmask(ji,jj,1) 188 190 … … 221 223 zzdep = fsdepw(ji,jj,jk) * tmask(ji,jj,1) 222 224 ! 223 zztmp = t n(ji,jj,nla10) - tn(ji,jj,jk)! - delta T(10m)225 zztmp = tsn(ji,jj,nla10,jp_tem) - tsn(ji,jj,jk,jp_tem) ! - delta T(10m) 224 226 IF( ABS(zztmp) > ztem2 ) zabs2 (ji,jj) = zzdep ! abs > 0.2 225 227 IF( zztmp > ztem2 ) ztm2 (ji,jj) = zzdep ! > 0.2 … … 254 256 DO jj = 1, jpj 255 257 DO ji = 1, jpi 256 zztmp = t n(ji,jj,jk)258 zztmp = tsn(ji,jj,jk,jp_tem) 257 259 IF( zztmp >= 20. ) ik20(ji,jj) = jk 258 260 IF( zztmp >= 28. ) ik28(ji,jj) = jk … … 273 275 zztmp = fsdept(ji,jj,iid ) & ! linear interpolation 274 276 & + ( fsdept(ji,jj,iid+1) - fsdept(ji,jj,iid) ) & 275 & * ( 20.*tmask(ji,jj,iid+1) - tn(ji,jj,iid) ) &276 & / ( tn(ji,jj,iid+1) - tn(ji,jj,iid) + (1.-tmask(ji,jj,1)) )277 & * ( 20.*tmask(ji,jj,iid+1) - tsn(ji,jj,iid,jp_tem) ) & 278 & / ( tsn(ji,jj,iid+1,jp_tem) - tsn(ji,jj,iid,jp_tem) + (1.-tmask(ji,jj,1)) ) 277 279 hd20(ji,jj) = MIN( zztmp , zzdep) * tmask(ji,jj,1) ! bound by the ocean depth 278 280 ELSE … … 284 286 zztmp = fsdept(ji,jj,iid ) & ! linear interpolation 285 287 & + ( fsdept(ji,jj,iid+1) - fsdept(ji,jj,iid) ) & 286 & * ( 28.*tmask(ji,jj,iid+1) - tn(ji,jj,iid) ) &287 & / ( tn(ji,jj,iid+1) - tn(ji,jj,iid) + (1.-tmask(ji,jj,1)) )288 & * ( 28.*tmask(ji,jj,iid+1) - tsn(ji,jj,iid,jp_tem) ) & 289 & / ( tsn(ji,jj,iid+1,jp_tem) - tsn(ji,jj,iid,jp_tem) + (1.-tmask(ji,jj,1)) ) 288 290 hd28(ji,jj) = MIN( zztmp , zzdep ) * tmask(ji,jj,1) ! bound by the ocean depth 289 291 ELSE … … 309 311 ! surface boundary condition 310 312 IF( lk_vvl ) THEN ; zthick(:,:) = 0._wp ; htc3(:,:) = 0._wp 311 ELSE ; zthick(:,:) = sshn(:,:) ; htc3(:,:) = t n(:,:,jk) * sshn(:,:) * tmask(:,:,jk)313 ELSE ; zthick(:,:) = sshn(:,:) ; htc3(:,:) = tsn(:,:,jk,jp_tem) * sshn(:,:) * tmask(:,:,jk) 312 314 ENDIF 313 315 ! integration down to ilevel 314 316 DO jk = 1, ilevel 315 317 zthick(:,:) = zthick(:,:) + fse3t(:,:,jk) 316 htc3 (:,:) = htc3 (:,:) + fse3t(:,:,jk) * t n(:,:,jk) * tmask(:,:,jk)318 htc3 (:,:) = htc3 (:,:) + fse3t(:,:,jk) * tsn(:,:,jk,jp_tem) * tmask(:,:,jk) 317 319 END DO 318 320 ! deepest layer … … 320 322 DO jj = 1, jpj 321 323 DO ji = 1, jpi 322 htc3(ji,jj) = htc3(ji,jj) + t n(ji,jj,ilevel+1) * MIN( fse3t(ji,jj,ilevel+1), zthick(ji,jj) ) * tmask(ji,jj,ilevel+1)324 htc3(ji,jj) = htc3(ji,jj) + tsn(ji,jj,ilevel+1,jp_tem) * MIN( fse3t(ji,jj,ilevel+1), zthick(ji,jj) ) * tmask(ji,jj,ilevel+1) 323 325 END DO 324 326 END DO
Note: See TracChangeset
for help on using the changeset viewer.