 Timestamp:
 20200608T15:32:56+02:00 (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/UKMO/AddC3SDiags/NEMOGCM/NEMO/OPA_SRC/DIA/diahth.F90
r11101 r13063 37 37 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hd28 !: depth of 28 C isotherm [m] 38 38 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: htc3 !: heat content of first 300 m [W] 39 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tav3 !: Average temp of first 300 m [degC] 40 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: sav3 !: Average salinity of first 300 m [PSU] 39 41 40 42 !! * Substitutions … … 53 55 ! 54 56 ALLOCATE(hth(jpi,jpj), hd20(jpi,jpj), hd28(jpi,jpj), htc3(jpi,jpj), STAT=dia_hth_alloc) 57 ! 58 ALLOCATE(tav3(jpi,jpj), sav3(jpi,jpj), STAT=dia_hth_alloc) 55 59 ! 56 60 IF( lk_mpp ) CALL mpp_sum ( dia_hth_alloc ) … … 92 96 REAL(wp) :: zztmp, zzdep ! temporary scalars inside do loop 93 97 REAL(wp) :: zu, zv, zw, zut, zvt ! temporary workspace 98 REAL(wp) :: depth_remain ! temporary float 94 99 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zabs2 ! MLD: abs( tn  tn(10m) ) = ztem2 95 100 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: ztm2 ! Top of thermocline: tn = tn(10m)  ztem2 … … 300 305 CALL iom_put( "28d", hd28 ) ! depth of the 28 isotherm 301 306 307 !  ! 308 ! Average temp & salinity of first 300 m ! 309 ! Copied from the hc300 section below ! 310 !  ! 311 ! find ilevel with (ilevel+1) the deepest Wlevel above 300m (we assume we can use e3t_0 to do this search...) 312 ilevel = 0 313 zthick_0 = 0._wp 314 DO jk = 1, jpkm1 315 zthick_0 = zthick_0 + e3t_1d(jk) 316 IF( zthick_0 < 300. ) ilevel = jk 317 END DO 318 ! surface boundary condition 319 IF( lk_vvl ) THEN ; zthick(:,:) = 0._wp ; tav3(:,:) = 0._wp ; sav3(:,:) = 0._wp 320 ELSE 321 ! For the surface layer calculate the integrated temp and salinity 322 ! The SSH "layer" is considered to be at the same temp/salinity as the 323 ! first ocean level. 324 zthick(:,:) = sshn(:,:) * tmask(:,:,1) 325 tav3(:,:) = tsn(:,:,1,jp_tem) * tmask(:,:,1) * zthick(:,:) 326 sav3(:,:) = tsn(:,:,1,jp_sal) * tmask(:,:,1) * zthick(:,:) 327 ENDIF 328 329 ! integration down to ilevel 330 ! For each level above 300m (1 to ilevel) calculate the integrated 331 ! temp and salinty (TS * layer thickness), and sum the thickness 332 DO jk = 1, ilevel 333 zthick(:,:) = zthick(:,:) + fse3t(:,:,jk) * tmask(:,:,jk) 334 tav3 (:,:) = tav3 (:,:) + tsn(:,:,jk,jp_tem) * tmask(:,:,jk) * fse3t(:,:,jk) 335 sav3 (:,:) = sav3 (:,:) + tsn(:,:,jk,jp_sal) * tmask(:,:,jk) * fse3t(:,:,jk) 336 END DO 337 338 DO jj = 1, jpj 339 DO ji = 1, jpi 340 ! What distance is there remaining to get to 300m 341 depth_remain = 300.0  zthick(ji,jj) 342 343 ! add the value of reamining thickness * TS of the level beyond 344 ! 300m. 345 tav3(ji,jj) = tav3(ji,jj) + tsn(ji,jj,ilevel+1,jp_tem) * MIN( fse3t(ji,jj,ilevel+1), depth_remain ) & 346 * tmask(ji,jj,ilevel+1) 347 sav3(ji,jj) = sav3(ji,jj) + tsn(ji,jj,ilevel+1,jp_sal) * MIN( fse3t(ji,jj,ilevel+1), depth_remain ) & 348 * tmask(ji,jj,ilevel+1) 349 ! Add the remaining depth, this should take us to 300m (in most 350 ! places) 351 zthick(ji,jj) = zthick(ji,jj) + MIN( fse3t(ji,jj,ilevel+1), depth_remain ) * tmask(ji,jj,ilevel+1) 352 ! Divide by the depth to get the integrated tempsalinity 353 tav3(ji,jj) = tav3(ji,jj) / zthick(ji,jj) 354 sav3(ji,jj) = sav3(ji,jj) / zthick(ji,jj) 355 356 END DO 357 END DO 358 359 CALL iom_put( "tav300", tav3 ) ! first 300m average temperature 360 CALL iom_put( "sav300", sav3 ) ! first 300m average salinity 361 302 362 !  ! 303 363 ! Heat content of first 300 m !
Note: See TracChangeset
for help on using the changeset viewer.