Changeset 12252
- Timestamp:
- 2019-12-14T14:57:23+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019
- Files:
-
- 2 deleted
- 55 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/1_namelist_cfg
r12205 r12252 382 382 !! namflo float parameters (default: OFF) 383 383 !! nam_diadct transports through some sections (default: OFF) 384 !! nam_diatmb Top Middle Bottom Output (default: OFF)385 384 !! nam_dia25h 25h Mean Output (default: OFF) 386 385 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/2_namelist_cfg
r12205 r12252 331 331 !! namflo float parameters (default: OFF) 332 332 !! nam_diadct transports through some sections (default: OFF) 333 !! nam_diatmb Top Middle Bottom Output (default: OFF)334 333 !! nam_dia25h 25h Mean Output (default: OFF) 335 334 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/3_namelist_cfg
r12205 r12252 313 313 !! namflo float parameters (default: OFF) 314 314 !! nam_diadct transports through some sections (default: OFF) 315 !! nam_diatmb Top Middle Bottom Output (default: OFF)316 315 !! nam_dia25h 25h Mean Output (default: OFF) 317 316 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/namelist_cfg
r12205 r12252 382 382 !! namflo float parameters (default: OFF) 383 383 !! nam_diadct transports through some sections (default: OFF) 384 !! nam_diatmb Top Middle Bottom Output (default: OFF)385 384 !! nam_dia25h 25h Mean Output (default: OFF) 386 385 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AMM12/EXPREF/file_def_nemo-oce.xml
r12205 r12252 13 13 <file_group id="1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files --> 14 14 15 <!--old <file_group id="1h" output_freq="1h" output_level="10" enabled=".TRUE."/> old --> <!-- 1h files -->16 <!-- TMB files -->15 <!--old <file_group id="1h" output_freq="1h" output_level="10" enabled=".TRUE."/> old --> <!-- 1h files --> 16 <!-- TB files --> 17 17 <file_group id="1h" output_freq="1h" output_level="10" enabled=".TRUE." > 18 18 19 <file id="file9" name_suffix="_shelftmb_grid_T" description="ocean T grid variables" enabled=".TRUE." > 20 <field field_ref="top_temp" name="votemper_top" operation="instant" enabled=".TRUE." /> 21 <field field_ref="mid_temp" name="votemper_mid" operation="instant" enabled=".TRUE." /> 22 <field field_ref="bot_temp" name="votemper_bot" operation="instant" enabled=".TRUE." /> 23 <field field_ref="sshnmasked" name="sossheig" operation="instant" enabled=".TRUE." /> 24 <field field_ref="top_sal" name="vosaline_top" operation="instant" enabled=".TRUE." /> 25 <field field_ref="mid_sal" name="vosaline_mid" operation="instant" enabled=".TRUE." /> 26 <field field_ref="bot_sal" name="vosaline_bot" operation="instant" enabled=".TRUE." /> 27 </file> 28 29 <file id="file10" name_suffix="_shelftmb_grid_U" description="TMB ocean U grid variables" enabled=".TRUE." > 30 <field field_ref="top_u" name="vozocrtx_top" operation="instant" /> 31 <field field_ref="mid_u" name="vozocrtx_mid" operation="instant" /> 32 <field field_ref="bot_u" name="vozocrtx_bot" operation="instant" /> 33 <field field_ref="baro_u" name="vobtcrtx" operation="instant" /> 34 </file> 35 36 <file id="file11" name_suffix="_shelftmb_grid_V" description="TMB ocean V grid variables" enabled=".TRUE." > 37 <field field_ref="top_v" name="vomecrty_top" operation="instant" /> 38 <field field_ref="mid_v" name="vomecrty_mid" operation="instant" /> 39 <field field_ref="bot_v" name="vomecrty_bot" operation="instant" /> 40 <field field_ref="baro_v" name="vobtcrty" operation="instant" /> 41 </file> 42 19 <file id="file9" name_suffix="_shelftb_grid_T" description="TB ocean T grid variables" enabled=".TRUE." > 20 <field field_ref="sst" operation="instant" enabled=".TRUE." /> 21 <field field_ref="sbt" operation="instant" enabled=".TRUE." /> 22 <field field_ref="sss" operation="instant" enabled=".TRUE." /> 23 <field field_ref="sbs" operation="instant" enabled=".TRUE." /> 24 </file> 25 26 <file id="file10" name_suffix="_shelftb_grid_U" description="TB ocean U grid variables" enabled=".TRUE." > 27 <field field_ref="ssu" operation="instant" enabled=".TRUE." /> 28 <field field_ref="sbu" operation="instant" enabled=".TRUE." /> 29 </file> 30 31 <file id="file11" name_suffix="_shelftb_grid_V" description="TB ocean V grid variables" enabled=".TRUE." > 32 <field field_ref="ssv" operation="instant" enabled=".TRUE." /> 33 <field field_ref="sbv" operation="instant" enabled=".TRUE." /> 34 </file> 35 43 36 </file_group> 44 37 -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AMM12/EXPREF/namelist_cfg
r12205 r12252 351 351 !! namflo float parameters (default: OFF) 352 352 !! nam_diadct transports through some sections (default: OFF) 353 !! nam_diatmb Top Middle Bottom Output (default: OFF)354 353 !! nam_dia25h 25h Mean Output (default: OFF) 355 354 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") 356 355 !!====================================================================== 357 356 ! 358 !-----------------------------------------------------------------------359 &nam_diatmb ! Top Middle Bottom Output (default F)360 !-----------------------------------------------------------------------361 ln_diatmb = .true. ! Choose Top Middle and Bottom output or not362 /363 357 !----------------------------------------------------------------------- 364 358 &nam_dia25h ! 25h Mean Output (default F) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/C1D_PAPA/EXPREF/namelist_cfg
r12205 r12252 421 421 !! namflo float parameters (default: OFF) 422 422 !! nam_diadct transports through some sections (default: OFF) 423 !! nam_diatmb Top Middle Bottom Output (default: OFF)424 423 !! nam_dia25h 25h Mean Output (default: OFF) 425 424 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") … … 448 447 !----------------------------------------------------------------------- 449 448 &nam_diadct ! transports through some sections (default: OFF) 450 !-----------------------------------------------------------------------451 /452 !-----------------------------------------------------------------------453 &nam_diatmb ! Top Middle Bottom Output (default: OFF)454 449 !----------------------------------------------------------------------- 455 450 / -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/GYRE_BFM/EXPREF/namelist_cfg
r12205 r12252 229 229 !! namflo float parameters (default: OFF) 230 230 !! nam_diadct transports through some sections (default: OFF) 231 !! nam_diatmb Top Middle Bottom Output (default: OFF)232 231 !! nam_dia25h 25h Mean Output (default: OFF) 233 232 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/GYRE_PISCES/EXPREF/namelist_cfg
r12205 r12252 223 223 !! namflo float parameters (default: OFF) 224 224 !! nam_diadct transports through some sections (default: OFF) 225 !! nam_diatmb Top Middle Bottom Output (default: OFF)226 225 !! nam_dia25h 25h Mean Output (default: OFF) 227 226 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg
r12205 r12252 392 392 !! namflo float parameters (default: OFF) 393 393 !! nam_diadct transports through some sections (default: OFF) 394 !! nam_diatmb Top Middle Bottom Output (default: OFF)395 394 !! nam_dia25h 25h Mean Output (default: OFF) 396 395 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_cfg
r12205 r12252 378 378 !! namflo float parameters (default: OFF) 379 379 !! nam_diadct transports through some sections (default: OFF) 380 !! nam_diatmb Top Middle Bottom Output (default: OFF)381 380 !! nam_dia25h 25h Mean Output (default: OFF) 382 381 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") … … 405 404 !----------------------------------------------------------------------- 406 405 &nam_diadct ! transports through some sections (default: OFF) 407 !-----------------------------------------------------------------------408 /409 !-----------------------------------------------------------------------410 &nam_diatmb ! Top Middle Bottom Output (default: OFF)411 406 !----------------------------------------------------------------------- 412 407 / -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_OFF_TRC/EXPREF/namelist_cfg
r12205 r12252 376 376 !! namflo float parameters (default: OFF) 377 377 !! namdct transports through some sections ("key_diadct") 378 !! nam_diatmb Top Middle Bottom Output (default: OFF)379 378 !! nam_dia25h 25h Mean Output (default: OFF) 380 379 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") … … 403 402 !----------------------------------------------------------------------- 404 403 &namdct ! transports through some sections ("key_diadct") 405 !-----------------------------------------------------------------------406 /407 !-----------------------------------------------------------------------408 &nam_diatmb ! Top Middle Bottom Output (default: OFF)409 404 !----------------------------------------------------------------------- 410 405 / -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_SAS_ICE/EXPREF/namelist_cfg
r12205 r12252 197 197 !! namflo float parameters (default: OFF) 198 198 !! nam_diadct transports through some sections (default: OFF) 199 !! nam_diatmb Top Middle Bottom Output (default: OFF)200 199 !! nam_dia25h 25h Mean Output (default: OFF) 201 200 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/field_def_nemo-oce.xml
r12229 r12252 1061 1061 </field_group> 1062 1062 1063 <!-- TMB diagnostic output -->1064 <field_group id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant">1065 <field id="top_temp" name="votemper_top" unit="degC" />1066 <field id="mid_temp" name="votemper_mid" unit="degC" />1067 <field id="bot_temp" name="votemper_bot" unit="degC" />1068 <field id="top_sal" name="vosaline_top" unit="psu" />1069 <field id="mid_sal" name="vosaline_mid" unit="psu" />1070 <field id="bot_sal" name="vosaline_bot" unit="psu" />1071 <field id="sshnmasked" name="sossheig" unit="m" />1072 </field_group>1073 1074 1063 <field_group id="1h_grid_U_tmb" grid_ref="grid_U_2D" operation="instant"> 1075 1064 <field id="top_u" name="vozocrtx_top" unit="m/s" /> -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/grid_def_nemo.xml
r12205 r12252 24 24 <!-- --> 25 25 <grid id="grid_T_3DS" > 26 <domain id="grid_T" />27 <axis id="profsed" />26 <domain domain_ref="grid_T" /> 27 <axis axis_ref="profsed" /> 28 28 </grid> 29 29 <!-- --> -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/namelist_ice_ref
r12182 r12252 60 60 rn_icebfr = 15. ! maximum bottom stress per unit volume [N/m3] 61 61 rn_lfrelax = 1.e-5 ! relaxation time scale to reach static friction [s-1] 62 rn_tensile = 0. 2! isotropic tensile strength [0-0.5??]62 rn_tensile = 0.05 ! isotropic tensile strength [0-0.5??] 63 63 / 64 64 !------------------------------------------------------------------------------ … … 91 91 !------------------------------------------------------------------------------ 92 92 ln_rhg_EVP = .true. ! EVP rheology 93 ln_aEVP = . false.! adaptive rheology (Kimmritz et al. 2016 & 2017)93 ln_aEVP = .true. ! adaptive rheology (Kimmritz et al. 2016 & 2017) 94 94 rn_creepl = 2.0e-9 ! creep limit [1/s] 95 95 rn_ecc = 2.0 ! eccentricity of the elliptical yield curve -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/namelist_ref
r12236 r12252 1199 1199 !! namflo float parameters (default: OFF) 1200 1200 !! nam_diadct transports through some sections (default: OFF) 1201 !! nam_diatmb Top Middle Bottom Output (default: OFF)1202 1201 !! nam_dia25h 25h Mean Output (default: OFF) 1203 1202 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") … … 1260 1259 ! ! -1 : debug all section 1261 1260 ! ! 0 < n : debug section number n 1262 /1263 !-----------------------------------------------------------------------1264 &nam_diatmb ! Top Middle Bottom Output (default: OFF)1265 !-----------------------------------------------------------------------1266 ln_diatmb = .true. ! Choose Top Middle and Bottom output or not1267 1261 / 1268 1262 !----------------------------------------------------------------------- -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SPITZ12/EXPREF/namelist_cfg
r12205 r12252 349 349 !! namflo float parameters (default: OFF) 350 350 !! nam_diadct transports through some sections (default: OFF) 351 !! nam_diatmb Top Middle Bottom Output (default: OFF)352 351 !! nam_dia25h 25h Mean Output (default: OFF) 353 352 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/doc/latex/NEMO/subfiles/chap_DIA.tex
r12205 r12252 1976 1976 1977 1977 %% ================================================================================================= 1978 \subsection{Top middle and bed hourly output}1979 1980 \begin{listing}1981 \nlst{nam_diatmb}1982 \caption{\forcode{&nam_diatmb}}1983 \label{lst:nam_diatmb}1984 \end{listing}1985 1986 A module is available to output the surface (top), mid water and bed diagnostics of a set of standard variables.1987 This can be a useful diagnostic when hourly or sub-hourly output is required in high resolution tidal outputs.1988 The tidal signal is retained but the overall data usage is cut to just three vertical levels.1989 Also the bottom level is calculated for each cell.1990 This diagnostic is actived with the logical $ln\_diatmb$.1991 1992 %% =================================================================================================1993 1978 \subsection{Courant numbers} 1994 1979 -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/ICE/icedyn_adv.F90
r11960 r12252 88 88 CASE( np_advPRA ) ! PRATHER scheme ! 89 89 ! !-----------------------! 90 CALL ice_dyn_adv_pra( kt, u_ice, v_ice, &90 CALL ice_dyn_adv_pra( kt, u_ice, v_ice, h_i, h_s, h_ip, & 91 91 & ato_i, v_i, v_s, sv_i, oa_i, a_i, a_ip, v_ip, e_s, e_i ) 92 92 END SELECT -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/ICE/icedyn_adv_pra.F90
r11812 r12252 54 54 CONTAINS 55 55 56 SUBROUTINE ice_dyn_adv_pra( kt, pu_ice, pv_ice, &56 SUBROUTINE ice_dyn_adv_pra( kt, pu_ice, pv_ice, ph_i, ph_s, ph_ip, & 57 57 & pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, pe_s, pe_i ) 58 58 !!---------------------------------------------------------------------- … … 70 70 REAL(wp), DIMENSION(:,:) , INTENT(in ) :: pu_ice ! ice i-velocity 71 71 REAL(wp), DIMENSION(:,:) , INTENT(in ) :: pv_ice ! ice j-velocity 72 REAL(wp), DIMENSION(:,:,:) , INTENT(in ) :: ph_i ! ice thickness 73 REAL(wp), DIMENSION(:,:,:) , INTENT(in ) :: ph_s ! snw thickness 74 REAL(wp), DIMENSION(:,:,:) , INTENT(in ) :: ph_ip ! ice pond thickness 72 75 REAL(wp), DIMENSION(:,:) , INTENT(inout) :: pato_i ! open water area 73 76 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: pv_i ! ice volume … … 87 90 REAL(wp), DIMENSION(jpi,jpj) :: zati1, zati2 88 91 REAL(wp), DIMENSION(jpi,jpj) :: zudy, zvdx 92 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zhi_max, zhs_max, zhip_max 89 93 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zarea 90 94 REAL(wp), DIMENSION(jpi,jpj,jpl) :: z0ice, z0snw, z0ai, z0smi, z0oi … … 95 99 ! 96 100 IF( kt == nit000 .AND. lwp ) WRITE(numout,*) '-- ice_dyn_adv_pra: Prather advection scheme' 101 ! 102 ! --- Record max of the surrounding 9-pts ice thick. (for call Hbig) --- ! 103 DO jl = 1, jpl 104 DO jj = 2, jpjm1 105 DO ji = fs_2, fs_jpim1 106 zhip_max(ji,jj,jl) = MAX( epsi20, ph_ip(ji,jj,jl), ph_ip(ji+1,jj ,jl), ph_ip(ji ,jj+1,jl), & 107 & ph_ip(ji-1,jj ,jl), ph_ip(ji ,jj-1,jl), & 108 & ph_ip(ji+1,jj+1,jl), ph_ip(ji-1,jj-1,jl), & 109 & ph_ip(ji+1,jj-1,jl), ph_ip(ji-1,jj+1,jl) ) 110 zhi_max (ji,jj,jl) = MAX( epsi20, ph_i (ji,jj,jl), ph_i (ji+1,jj ,jl), ph_i (ji ,jj+1,jl), & 111 & ph_i (ji-1,jj ,jl), ph_i (ji ,jj-1,jl), & 112 & ph_i (ji+1,jj+1,jl), ph_i (ji-1,jj-1,jl), & 113 & ph_i (ji+1,jj-1,jl), ph_i (ji-1,jj+1,jl) ) 114 zhs_max (ji,jj,jl) = MAX( epsi20, ph_s (ji,jj,jl), ph_s (ji+1,jj ,jl), ph_s (ji ,jj+1,jl), & 115 & ph_s (ji-1,jj ,jl), ph_s (ji ,jj-1,jl), & 116 & ph_s (ji+1,jj+1,jl), ph_s (ji-1,jj-1,jl), & 117 & ph_s (ji+1,jj-1,jl), ph_s (ji-1,jj+1,jl) ) 118 END DO 119 END DO 120 END DO 121 CALL lbc_lnk_multi( 'icedyn_adv_pra', zhi_max, 'T', 1., zhs_max, 'T', 1., zhip_max, 'T', 1. ) 97 122 ! 98 123 ! --- If ice drift is too fast, use subtime steps for advection (CFL test for stability) --- ! … … 239 264 ! (because advected fields are not perfectly bounded and tiny negative values can occur, e.g. -1.e-20) 240 265 CALL ice_var_zapneg( zdt, pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, pe_s, pe_i ) 266 ! 267 ! --- Make sure ice thickness is not too big --- ! 268 ! (because ice thickness can be too large where ice concentration is very small) 269 CALL Hbig( zdt, zhi_max, zhs_max, zhip_max, pv_i, pv_s, pa_i, pa_ip, pv_ip, pe_s ) 241 270 ! 242 271 ! --- Ensure snow load is not too big --- ! … … 588 617 ! 589 618 END SUBROUTINE adv_y 619 620 621 SUBROUTINE Hbig( pdt, phi_max, phs_max, phip_max, pv_i, pv_s, pa_i, pa_ip, pv_ip, pe_s ) 622 !!------------------------------------------------------------------- 623 !! *** ROUTINE Hbig *** 624 !! 625 !! ** Purpose : Thickness correction in case advection scheme creates 626 !! abnormally tick ice or snow 627 !! 628 !! ** Method : 1- check whether ice thickness is larger than the surrounding 9-points 629 !! (before advection) and reduce it by adapting ice concentration 630 !! 2- check whether snow thickness is larger than the surrounding 9-points 631 !! (before advection) and reduce it by sending the excess in the ocean 632 !! 633 !! ** input : Max thickness of the surrounding 9-points 634 !!------------------------------------------------------------------- 635 REAL(wp) , INTENT(in ) :: pdt ! tracer time-step 636 REAL(wp), DIMENSION(:,:,:) , INTENT(in ) :: phi_max, phs_max, phip_max ! max ice thick from surrounding 9-pts 637 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: pv_i, pv_s, pa_i, pa_ip, pv_ip 638 REAL(wp), DIMENSION(:,:,:,:), INTENT(inout) :: pe_s 639 ! 640 INTEGER :: ji, jj, jl ! dummy loop indices 641 REAL(wp) :: z1_dt, zhip, zhi, zhs, zfra 642 !!------------------------------------------------------------------- 643 ! 644 z1_dt = 1._wp / pdt 645 ! 646 DO jl = 1, jpl 647 648 DO jj = 1, jpj 649 DO ji = 1, jpi 650 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 651 ! 652 ! ! -- check h_ip -- ! 653 ! if h_ip is larger than the surrounding 9 pts => reduce h_ip and increase a_ip 654 IF( ln_pnd_H12 .AND. pv_ip(ji,jj,jl) > 0._wp ) THEN 655 zhip = pv_ip(ji,jj,jl) / MAX( epsi20, pa_ip(ji,jj,jl) ) 656 IF( zhip > phip_max(ji,jj,jl) .AND. pa_ip(ji,jj,jl) < 0.15 ) THEN 657 pa_ip(ji,jj,jl) = pv_ip(ji,jj,jl) / phip_max(ji,jj,jl) 658 ENDIF 659 ENDIF 660 ! 661 ! ! -- check h_i -- ! 662 ! if h_i is larger than the surrounding 9 pts => reduce h_i and increase a_i 663 zhi = pv_i(ji,jj,jl) / pa_i(ji,jj,jl) 664 IF( zhi > phi_max(ji,jj,jl) .AND. pa_i(ji,jj,jl) < 0.15 ) THEN 665 pa_i(ji,jj,jl) = pv_i(ji,jj,jl) / MIN( phi_max(ji,jj,jl), hi_max(jpl) ) !-- bound h_i to hi_max (99 m) 666 ENDIF 667 ! 668 ! ! -- check h_s -- ! 669 ! if h_s is larger than the surrounding 9 pts => put the snow excess in the ocean 670 zhs = pv_s(ji,jj,jl) / pa_i(ji,jj,jl) 671 IF( pv_s(ji,jj,jl) > 0._wp .AND. zhs > phs_max(ji,jj,jl) .AND. pa_i(ji,jj,jl) < 0.15 ) THEN 672 zfra = phs_max(ji,jj,jl) / MAX( zhs, epsi20 ) 673 ! 674 wfx_res(ji,jj) = wfx_res(ji,jj) + ( pv_s(ji,jj,jl) - pa_i(ji,jj,jl) * phs_max(ji,jj,jl) ) * rhos * z1_dt 675 hfx_res(ji,jj) = hfx_res(ji,jj) - SUM( pe_s(ji,jj,1:nlay_s,jl) ) * ( 1._wp - zfra ) * z1_dt ! W.m-2 <0 676 ! 677 pe_s(ji,jj,1:nlay_s,jl) = pe_s(ji,jj,1:nlay_s,jl) * zfra 678 pv_s(ji,jj,jl) = pa_i(ji,jj,jl) * phs_max(ji,jj,jl) 679 ENDIF 680 ! 681 ENDIF 682 END DO 683 END DO 684 END DO 685 ! 686 END SUBROUTINE Hbig 590 687 591 688 -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/ICE/icedyn_adv_umx.F90
r11627 r12252 352 352 CALL ice_var_zapneg( zdt, pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, pe_s, pe_i ) 353 353 ! 354 ! Make sure ice thickness is not too big 355 ! (because ice thickness can be too large where ice concentration is very small) 356 CALL Hbig( zdt, zhi_max, zhs_max, zhip_max, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, pe_s, pe_i ) 357 354 ! --- Make sure ice thickness is not too big --- ! 355 ! (because ice thickness can be too large where ice concentration is very small) 356 CALL Hbig( zdt, zhi_max, zhs_max, zhip_max, pv_i, pv_s, pa_i, pa_ip, pv_ip, pe_s ) 357 ! 358 ! --- Ensure snow load is not too big --- ! 359 CALL Hsnow( zdt, pv_i, pv_s, pa_i, pa_ip, pe_s ) 360 ! 358 361 END DO 359 362 ! … … 1514 1517 1515 1518 1516 SUBROUTINE Hbig( pdt, phi_max, phs_max, phip_max, pv_i, pv_s, p sv_i, poa_i, pa_i, pa_ip, pv_ip, pe_s, pe_i)1519 SUBROUTINE Hbig( pdt, phi_max, phs_max, phip_max, pv_i, pv_s, pa_i, pa_ip, pv_ip, pe_s ) 1517 1520 !!------------------------------------------------------------------- 1518 1521 !! *** ROUTINE Hbig *** … … 1525 1528 !! 2- check whether snow thickness is larger than the surrounding 9-points 1526 1529 !! (before advection) and reduce it by sending the excess in the ocean 1527 !! 3- check whether snow load deplets the snow-ice interface below sea level$1528 !! and reduce it by sending the excess in the ocean1529 !! 4- correct pond concentration to avoid a_ip > a_i1530 1530 !! 1531 1531 !! ** input : Max thickness of the surrounding 9-points … … 1533 1533 REAL(wp) , INTENT(in ) :: pdt ! tracer time-step 1534 1534 REAL(wp), DIMENSION(:,:,:) , INTENT(in ) :: phi_max, phs_max, phip_max ! max ice thick from surrounding 9-pts 1535 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: pv_i, pv_s, p sv_i, poa_i, pa_i, pa_ip, pv_ip1535 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: pv_i, pv_s, pa_i, pa_ip, pv_ip 1536 1536 REAL(wp), DIMENSION(:,:,:,:), INTENT(inout) :: pe_s 1537 REAL(wp), DIMENSION(:,:,:,:), INTENT(inout) :: pe_i 1538 ! 1539 INTEGER :: ji, jj, jk, jl ! dummy loop indices 1540 REAL(wp) :: z1_dt, zhip, zhi, zhs, zvs_excess, zfra 1541 REAL(wp), DIMENSION(jpi,jpj) :: zswitch 1537 ! 1538 INTEGER :: ji, jj, jl ! dummy loop indices 1539 REAL(wp) :: z1_dt, zhip, zhi, zhs, zfra 1542 1540 !!------------------------------------------------------------------- 1543 1541 ! … … 1578 1576 pv_s(ji,jj,jl) = pa_i(ji,jj,jl) * phs_max(ji,jj,jl) 1579 1577 ENDIF 1578 ! 1579 ENDIF 1580 END DO 1581 END DO 1582 END DO 1583 ! 1584 END SUBROUTINE Hbig 1585 1586 1587 SUBROUTINE Hsnow( pdt, pv_i, pv_s, pa_i, pa_ip, pe_s ) 1588 !!------------------------------------------------------------------- 1589 !! *** ROUTINE Hsnow *** 1590 !! 1591 !! ** Purpose : 1- Check snow load after advection 1592 !! 2- Correct pond concentration to avoid a_ip > a_i 1593 !! 1594 !! ** Method : If snow load makes snow-ice interface to deplet below the ocean surface 1595 !! then put the snow excess in the ocean 1596 !! 1597 !! ** Notes : This correction is crucial because of the call to routine icecor afterwards 1598 !! which imposes a mini of ice thick. (rn_himin). This imposed mini can artificially 1599 !! make the snow very thick (if concentration decreases drastically) 1600 !! This behavior has been seen in Ultimate-Macho and supposedly it can also be true for Prather 1601 !!------------------------------------------------------------------- 1602 REAL(wp) , INTENT(in ) :: pdt ! tracer time-step 1603 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: pv_i, pv_s, pa_i, pa_ip 1604 REAL(wp), DIMENSION(:,:,:,:), INTENT(inout) :: pe_s 1605 ! 1606 INTEGER :: ji, jj, jl ! dummy loop indices 1607 REAL(wp) :: z1_dt, zvs_excess, zfra 1608 !!------------------------------------------------------------------- 1609 ! 1610 z1_dt = 1._wp / pdt 1611 ! 1612 ! -- check snow load -- ! 1613 DO jl = 1, jpl 1614 DO jj = 1, jpj 1615 DO ji = 1, jpi 1616 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 1580 1617 ! 1581 ! ! -- check snow load -- !1582 ! if snow load makes snow-ice interface to deplet below the ocean surface => put the snow excess in the ocean1583 ! this correction is crucial because of the call to routine icecor afterwards which imposes a mini of ice thick. (rn_himin)1584 ! this imposed mini can artificially make the snow very thick (if concentration decreases drastically)1585 1618 zvs_excess = MAX( 0._wp, pv_s(ji,jj,jl) - pv_i(ji,jj,jl) * (rau0-rhoi) * r1_rhos ) 1586 IF( zvs_excess > 0._wp ) THEN 1619 ! 1620 IF( zvs_excess > 0._wp ) THEN ! snow-ice interface deplets below the ocean surface 1621 ! put snow excess in the ocean 1587 1622 zfra = ( pv_s(ji,jj,jl) - zvs_excess ) / MAX( pv_s(ji,jj,jl), epsi20 ) 1588 1623 wfx_res(ji,jj) = wfx_res(ji,jj) + zvs_excess * rhos * z1_dt 1589 1624 hfx_res(ji,jj) = hfx_res(ji,jj) - SUM( pe_s(ji,jj,1:nlay_s,jl) ) * ( 1._wp - zfra ) * z1_dt ! W.m-2 <0 1590 ! 1625 ! correct snow volume and heat content 1591 1626 pe_s(ji,jj,1:nlay_s,jl) = pe_s(ji,jj,1:nlay_s,jl) * zfra 1592 1627 pv_s(ji,jj,jl) = pv_s(ji,jj,jl) - zvs_excess 1593 1628 ENDIF 1594 1629 ! 1595 1630 ENDIF 1596 1631 END DO 1597 1632 END DO 1598 END DO 1599 ! !-- correct pond concentration to avoid a_ip > a_i 1633 END DO 1634 ! 1635 !-- correct pond concentration to avoid a_ip > a_i -- ! 1600 1636 WHERE( pa_ip(:,:,:) > pa_i(:,:,:) ) pa_ip(:,:,:) = pa_i(:,:,:) 1601 1637 ! 1602 !1603 END SUBROUTINE Hbig 1604 1638 END SUBROUTINE Hsnow 1639 1640 1605 1641 #else 1606 1642 !!---------------------------------------------------------------------- -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdyini.F90
r12205 r12252 395 395 jpbdta = MAXVAL(nblendta(1:jpbgrd,1:nb_bdy)) 396 396 ALLOCATE( nbidta(jpbdta, jpbgrd, nb_bdy), nbjdta(jpbdta, jpbgrd, nb_bdy), nbrdta(jpbdta, jpbgrd, nb_bdy) ) 397 397 nbrdta(:,:,:) = 0 ! initialize nbrdta as it may not be completely defined for each bdy 398 398 399 ! Calculate global boundary index arrays or read in from file 399 400 !------------------------------------------------------------ -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdyvol.F90
r12150 r12252 143 143 ! Check the cumulated transport through unstructured OBC once barotropic velocities corrected 144 144 ! ------------------------------------------------------ 145 IF( MOD( kt, nn_write) == 0 .AND. ( kc == 1 ) ) THEN145 IF( MOD( kt, MAX(nn_write,1) ) == 0 .AND. ( kc == 1 ) ) THEN 146 146 ! 147 147 ! compute residual transport across boundary -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DIA/diawri.F90
r12182 r12252 49 49 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 50 50 USE in_out_manager ! I/O manager 51 USE diatmb ! Top,middle,bottom output52 51 USE dia25h ! 25h Mean output 53 52 USE iom ! … … 399 398 CALL iom_put( "bn2", rn2 ) ! Brunt-Vaisala buoyancy frequency (N^2) 400 399 ! 401 402 IF (ln_diatmb) CALL dia_tmb( Kmm ) ! tmb values 403 400 404 401 IF (ln_dia25h) CALL dia_25h( kt, Kmm ) ! 25h averaging 405 402 -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DYN/divhor.F90
r12150 r12252 73 73 IF(lwp) WRITE(numout,*) 'div_hor : horizontal velocity divergence ' 74 74 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 75 hdiv(:,:,:) = 0._wp ! initialize hdiv for the halos at the first time step 75 76 ENDIF 76 77 ! -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DYN/dynspg_ts.F90
r12250 r12252 46 46 USE tide_mod ! 47 47 USE sbcwave ! surface wave 48 USE diatmb ! Top,middle,bottom output49 48 #if defined key_agrif 50 49 USE agrif_oce_interp ! agrif … … 61 60 USE iom ! IOM library 62 61 USE restart ! only for lrst_oce 63 USE diatmb ! Top,middle,bottom output64 62 65 63 USE iom ! to remove … … 154 152 REAL(wp) :: r1_2dt_b, z1_hu, z1_hv ! local scalars 155 153 REAL(wp) :: za0, za1, za2, za3 ! - - 156 REAL(wp) :: z mdi, zztmp, zldg ! - -154 REAL(wp) :: zztmp, zldg ! - - 157 155 REAL(wp) :: zhu_bck, zhv_bck, zhdiv ! - - 158 156 REAL(wp) :: zun_save, zvn_save ! - - … … 178 176 ! !* Allocate temporary arrays 179 177 IF( ln_wd_dl ) ALLOCATE( ztwdmask(jpi,jpj), zuwdmask(jpi,jpj), zvwdmask(jpi,jpj), zuwdav2(jpi,jpj), zvwdav2(jpi,jpj)) 180 !181 zmdi=1.e+20 ! missing data indicator for masking182 178 ! 183 179 zwdramp = r_rn_wdmin1 ! simplest ramp … … 835 831 IF( ln_wd_dl ) DEALLOCATE( ztwdmask, zuwdmask, zvwdmask, zuwdav2, zvwdav2 ) 836 832 ! 837 IF( ln_diatmb ) THEN 838 CALL iom_put( "baro_u" , puu_b(:,:,Kmm)*ssumask(:,:)+zmdi*(1.-ssumask(:,:) ) ) ! Barotropic U Velocity 839 CALL iom_put( "baro_v" , pvv_b(:,:,Kmm)*ssvmask(:,:)+zmdi*(1.-ssvmask(:,:) ) ) ! Barotropic V Velocity 840 ENDIF 833 CALL iom_put( "baro_u" , puu_b(:,:,Kmm) ) ! Barotropic U Velocity 834 CALL iom_put( "baro_v" , pvv_b(:,:,Kmm) ) ! Barotropic V Velocity 841 835 ! 842 836 END SUBROUTINE dyn_spg_ts -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/IOM/iom.F90
r12205 r12252 837 837 CHARACTER(LEN=100) :: clinfo ! info character 838 838 !--------------------------------------------------------------------- 839 ! 840 IF( iom_open_init == 0 ) RETURN ! avoid to use iom_file(jf)%nfid that us not yet initialized 839 841 ! 840 842 clinfo = ' iom_close ~~~ ' -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/LDF/ldfdyn.F90
r11960 r12252 416 416 ! 417 417 zcmsmag = (rn_csmc/rpi)**2 ! (C_smag/pi)^2 418 zstabf_lo = rn_minfac * rn_minfac / ( 2._wp * 4._wp * zcmsmag )! lower limit stability factor scaling418 zstabf_lo = rn_minfac * rn_minfac / ( 2._wp * 12._wp * 12._wp * zcmsmag ) ! lower limit stability factor scaling 419 419 zstabf_up = rn_maxfac / ( 4._wp * zcmsmag * 2._wp * rdt ) ! upper limit stability factor scaling 420 420 IF( ln_dynldf_blp ) zstabf_lo = ( 16._wp / 9._wp ) * zstabf_lo ! provide |U|L^3/12 lower limit instead -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/LDF/ldfslp.F90
r12236 r12252 209 209 zbu = MIN( zbu, - z1_slpmax * ABS( zau ) , -7.e+3_wp/e3u(ji,jj,jk,Kmm)* ABS( zau ) ) 210 210 zbv = MIN( zbv, - z1_slpmax * ABS( zav ) , -7.e+3_wp/e3v(ji,jj,jk,Kmm)* ABS( zav ) ) 211 ! ! Fred Dupont: add a correction for bottom partial steps: 212 ! ! max slope = 1/2 * e3 / e1 213 IF (ln_zps .AND. jk==mbku(ji,jj)) & 214 zbu = MIN( zbu, - z1_slpmax * ABS( zau ) , - 2._wp * e1u(ji,jj) / e3u(ji,jj,jk,Kmm)* ABS( zau ) ) 215 IF (ln_zps .AND. jk==mbkv(ji,jj)) & 216 zbv = MIN( zbv, - z1_slpmax * ABS( zav ) , - 2._wp * e2v(ji,jj) / e3v(ji,jj,jk,Kmm)* ABS( zav ) ) 211 217 ! ! uslp and vslp output in zwz and zww, resp. 212 218 zfi = MAX( omlmask(ji,jj,jk), omlmask(ji+1,jj,jk) ) … … 405 411 REAL(wp) :: zbeta0, ze3_e1, ze3_e2 406 412 REAL(wp), DIMENSION(jpi,jpj) :: z1_mlbw 407 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zalbet408 413 REAL(wp), DIMENSION(jpi,jpj,jpk,0:1) :: zdxrho , zdyrho, zdzrho ! Horizontal and vertical density gradients 409 414 REAL(wp), DIMENSION(jpi,jpj,0:1,0:1) :: zti_mlb, ztj_mlb ! for Griffies operator only … … 461 466 zdks = 0._wp 462 467 ENDIF 463 zdzrho_raw = ( - rab_b(ji,jj,jk +kp,jp_tem) * zdkt &464 & + rab_b(ji,jj,jk +kp,jp_sal) * zdks &468 zdzrho_raw = ( - rab_b(ji,jj,jk ,jp_tem) * zdkt & 469 & + rab_b(ji,jj,jk ,jp_sal) * zdks & 465 470 & ) / e3w(ji,jj,jk+kp,Kmm) 466 471 zdzrho(ji,jj,jk,kp) = - MIN( - repsln , zdzrho_raw ) ! force zdzrho >= repsln -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/SBC/cpl_oasis3.F90
r12236 r12252 306 306 ! End of definition phase 307 307 !------------------------------------------------------------------ 308 308 ! 309 #if defined key_agrif 310 IF( agrif_fixed() == Agrif_Nb_Fine_Grids() ) THEN 311 #endif 309 312 CALL oasis_enddef(nerror) 310 313 IF( nerror /= OASIS_Ok ) CALL oasis_abort ( ncomp_id, 'cpl_define', 'Failure in oasis_enddef') 314 #if defined key_agrif 315 ENDIF 316 #endif 311 317 ! 312 318 IF( ltmp_wapatch ) THEN … … 357 363 WRITE(numout,*) 'oasis_put: kstep ', kstep 358 364 WRITE(numout,*) 'oasis_put: info ', kinfo 359 WRITE(numout,*) ' - Minimum value is ', MINVAL(pdata( :,:,jc))360 WRITE(numout,*) ' - Maximum value is ', MAXVAL(pdata( :,:,jc))361 WRITE(numout,*) ' - Sum value is ', SUM(pdata( :,:,jc))365 WRITE(numout,*) ' - Minimum value is ', MINVAL(pdata(nldi:nlei,nldj:nlej,jc)) 366 WRITE(numout,*) ' - Maximum value is ', MAXVAL(pdata(nldi:nlei,nldj:nlej,jc)) 367 WRITE(numout,*) ' - Sum value is ', SUM(pdata(nldi:nlei,nldj:nlej,jc)) 362 368 WRITE(numout,*) '****************' 363 369 ENDIF -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/SBC/sbc_oce.F90
r12199 r12252 107 107 !! Ocean Surface Boundary Condition fields 108 108 !!---------------------------------------------------------------------- 109 INTEGER , PUBLIC :: ncpl_qsr_freq !: qsr coupling frequency per days from atmosphere109 INTEGER , PUBLIC :: ncpl_qsr_freq = 0 !: qsr coupling frequency per days from atmosphere (used by top) 110 110 ! 111 111 !! !! now ! before !! -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/SBC/sbccpl.F90
r12193 r12252 571 571 IF( TRIM( sn_rcv_emp%clcat ) == 'yes' ) srcv(jpr_ievp)%nct = nn_cats_cpl 572 572 573 #if defined key_si3 574 IF( ln_cndflx .AND. .NOT.ln_cndemulate ) THEN 575 IF( .NOT.srcv(jpr_ts_ice)%laction ) & 576 & CALL ctl_stop( 'sbc_cpl_init: srcv(jpr_ts_ice)%laction should be set to true when ln_cndflx=T' ) 577 ENDIF 578 #endif 573 579 ! ! ------------------------- ! 574 580 ! ! Wave breaking ! … … 860 866 ELSE IF( sn_snd_crt%clvgrd /= 'T' ) THEN 861 867 CALL ctl_stop( 'sn_snd_crt%clvgrd must be equal to T' ) 862 ssnd(jps_ocx1:jps_ivz1)%clgrid = 'T' ! all oce and ice components on the same unique grid863 868 ENDIF 864 869 ssnd(jps_ocx1:jps_ivz1)%laction = .TRUE. ! default: all are send … … 1038 1043 ENDIF 1039 1044 xcplmask(:,:,0) = 1. - SUM( xcplmask(:,:,1:nn_cplmodel), dim = 3 ) 1040 !1041 ncpl_qsr_freq = cpl_freq( 'O_QsrOce' ) + cpl_freq( 'O_QsrMix' ) + cpl_freq( 'I_QsrOce' ) + cpl_freq( 'I_QsrMix' )1042 IF( ln_dm2dc .AND. ln_cpl .AND. ncpl_qsr_freq /= 86400 ) &1043 & CALL ctl_stop( 'sbc_cpl_init: diurnal cycle reconstruction (ln_dm2dc) needs daily couping for solar radiation' )1044 IF( ln_dm2dc .AND. ln_cpl ) ncpl_qsr_freq = 86400 / ncpl_qsr_freq1045 1045 ! 1046 1046 END SUBROUTINE sbc_cpl_init … … 1109 1109 REAL(wp), DIMENSION(jpi,jpj) :: ztx, zty, zmsk, zemp, zqns, zqsr 1110 1110 !!---------------------------------------------------------------------- 1111 ! 1112 IF( kt == nit000 ) THEN 1113 ! cannot be done in the init phase when we use agrif as cpl_freq requires that oasis_enddef is done 1114 ncpl_qsr_freq = cpl_freq( 'O_QsrOce' ) + cpl_freq( 'O_QsrMix' ) + cpl_freq( 'I_QsrOce' ) + cpl_freq( 'I_QsrMix' ) 1115 IF( ln_dm2dc .AND. ncpl_qsr_freq /= 86400 ) & 1116 & CALL ctl_stop( 'sbc_cpl_rcv: diurnal cycle reconstruction (ln_dm2dc) needs daily couping for solar radiation' ) 1117 ncpl_qsr_freq = 86400 / ncpl_qsr_freq ! used by top 1118 ENDIF 1111 1119 ! 1112 1120 IF( ln_mixcpl ) zmsk(:,:) = 1. - xcplmask(:,:,0) … … 1242 1250 IF( srcv(jpr_co2)%laction ) atm_co2(:,:) = frcv(jpr_co2)%z3(:,:,1) 1243 1251 ! 1244 ! ! ================== !1245 ! ! ice skin temp. !1246 ! ! ================== !1247 #if defined key_si31248 ! needed by Met Office1249 IF( srcv(jpr_ts_ice)%laction ) THEN1250 WHERE ( frcv(jpr_ts_ice)%z3(:,:,:) > 0.0 ) ; tsfc_ice(:,:,:) = 0.01251 ELSEWHERE( frcv(jpr_ts_ice)%z3(:,:,:) < -60. ) ; tsfc_ice(:,:,:) = -60.1252 ELSEWHERE ; tsfc_ice(:,:,:) = frcv(jpr_ts_ice)%z3(:,:,:)1253 END WHERE1254 ENDIF1255 #endif1256 1252 ! ! ========================= ! 1257 1253 ! ! Mean Sea Level Pressure ! (taum) … … 1633 1629 !! sprecip solid precipitation over the ocean 1634 1630 !!---------------------------------------------------------------------- 1635 REAL(wp), INTENT(in) , DIMENSION(:,:) :: picefr ! ice fraction [0 to 1]1636 ! !! ! optional arguments, used only in 'mixed oce-ice' case1637 REAL(wp), INTENT(in) , DIMENSION(:,:,:), OPTIONAL :: palbi ! all skies ice albedo1638 REAL(wp), INTENT(in) , DIMENSION(:,: ), OPTIONAL :: psst ! sea surface temperature [Celsius]1639 REAL(wp), INTENT(in ), DIMENSION(:,:,:), OPTIONAL :: pist ! ice surface temperature [Kelvin]1640 REAL(wp), INTENT(in) , DIMENSION(:,:,:), OPTIONAL :: phs ! snow depth [m]1641 REAL(wp), INTENT(in) , DIMENSION(:,:,:), OPTIONAL :: phi ! ice thickness [m]1631 REAL(wp), INTENT(in) , DIMENSION(:,:) :: picefr ! ice fraction [0 to 1] 1632 ! !! ! optional arguments, used only in 'mixed oce-ice' case or for Met-Office coupling 1633 REAL(wp), INTENT(in) , DIMENSION(:,:,:), OPTIONAL :: palbi ! all skies ice albedo 1634 REAL(wp), INTENT(in) , DIMENSION(:,: ), OPTIONAL :: psst ! sea surface temperature [Celsius] 1635 REAL(wp), INTENT(inout), DIMENSION(:,:,:), OPTIONAL :: pist ! ice surface temperature [Kelvin] => inout for Met-Office 1636 REAL(wp), INTENT(in) , DIMENSION(:,:,:), OPTIONAL :: phs ! snow depth [m] 1637 REAL(wp), INTENT(in) , DIMENSION(:,:,:), OPTIONAL :: phi ! ice thickness [m] 1642 1638 ! 1643 1639 INTEGER :: ji, jj, jl ! dummy loop index … … 1646 1642 REAL(wp), DIMENSION(jpi,jpj) :: zemp_tot, zemp_ice, zemp_oce, ztprecip, zsprecip , zevap_oce, zdevap_ice 1647 1643 REAL(wp), DIMENSION(jpi,jpj) :: zqns_tot, zqns_oce, zqsr_tot, zqsr_oce, zqprec_ice, zqemp_oce, zqemp_ice 1648 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zqns_ice, zqsr_ice, zdqns_ice, zqevap_ice, zevap_ice !!gm , zfrqsr_tr_i1644 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zqns_ice, zqsr_ice, zdqns_ice, zqevap_ice, zevap_ice, zqtr_ice_top, ztsu 1649 1645 !!---------------------------------------------------------------------- 1650 1646 ! … … 1814 1810 ! ** NEED TO SORT OUT HOW THIS SHOULD WORK IN THE MULTI-CATEGORY CASE - CURRENTLY NOT ALLOWED WHEN INTERFACE INITIALISED ** 1815 1811 zqns_tot(:,: ) = frcv(jpr_qnsmix)%z3(:,:,1) 1816 zqns_ice(:,:,1) = frcv(jpr_qnsmix)%z3(:,:,1) & 1817 & + frcv(jpr_dqnsdt)%z3(:,:,1) * ( pist(:,:,1) - ( (rt0 + psst(:,: ) ) * ziceld(:,:) & 1818 & + pist(:,:,1) * picefr(:,:) ) ) 1812 IF ( TRIM(sn_rcv_qsr%clcat) == 'yes' ) THEN 1813 DO jl = 1, jpl 1814 zqns_ice(:,:,jl) = frcv(jpr_qnsmix)%z3(:,:,jl) & 1815 & + frcv(jpr_dqnsdt)%z3(:,:,jl) * ( pist(:,:,jl) - ( ( rt0 + psst(:,:) ) * ziceld(:,:) & 1816 & + pist(:,:,jl) * picefr(:,:) ) ) 1817 END DO 1818 ELSE 1819 DO jl = 1, jpl 1820 zqns_ice(:,:,jl) = frcv(jpr_qnsmix)%z3(:,:, 1) & 1821 & + frcv(jpr_dqnsdt)%z3(:,:, 1) * ( pist(:,:,jl) - ( ( rt0 + psst(:,:) ) * ziceld(:,:) & 1822 & + pist(:,:,jl) * picefr(:,:) ) ) 1823 END DO 1824 ENDIF 1819 1825 END SELECT 1820 1826 ! … … 1930 1936 END DO 1931 1937 ENDIF 1932 zqsr_tot(:,: ) = frcv(jpr_qsrmix)%z3(:,:,1)1933 zqsr_ice(:,:,1) = frcv(jpr_qsrice)%z3(:,:,1)1934 1938 CASE( 'oce and ice' ) 1935 1939 zqsr_tot(:,: ) = ziceld(:,:) * frcv(jpr_qsroce)%z3(:,:,1) … … 1951 1955 ! Create solar heat flux over ice using incoming solar heat flux and albedos 1952 1956 ! ( see OASIS3 user guide, 5th edition, p39 ) 1953 zqsr_ice(:,:,1) = frcv(jpr_qsrmix)%z3(:,:,1) * ( 1.- palbi(:,:,1) ) & 1954 & / ( 1.- ( alb_oce_mix(:,: ) * ziceld(:,:) & 1955 & + palbi (:,:,1) * picefr(:,:) ) ) 1957 IF ( TRIM(sn_rcv_qsr%clcat) == 'yes' ) THEN 1958 DO jl = 1, jpl 1959 zqsr_ice(:,:,jl) = frcv(jpr_qsrmix)%z3(:,:,jl) * ( 1.- palbi(:,:,jl) ) & 1960 & / ( 1.- ( alb_oce_mix(:,: ) * ziceld(:,:) & 1961 & + palbi (:,:,jl) * picefr(:,:) ) ) 1962 END DO 1963 ELSE 1964 DO jl = 1, jpl 1965 zqsr_ice(:,:,jl) = frcv(jpr_qsrmix)%z3(:,:, 1) * ( 1.- palbi(:,:,jl) ) & 1966 & / ( 1.- ( alb_oce_mix(:,: ) * ziceld(:,:) & 1967 & + palbi (:,:,jl) * picefr(:,:) ) ) 1968 END DO 1969 ENDIF 1956 1970 CASE( 'none' ) ! Not available as for now: needs additional coding 1957 1971 ! ! since fields received, here zqsr_tot, are not defined with none option … … 2013 2027 ! ! ========================= ! 2014 2028 CASE ('coupled') 2015 qml_ice(:,:,:) = frcv(jpr_topm)%z3(:,:,:) 2016 qcn_ice(:,:,:) = frcv(jpr_botm)%z3(:,:,:) 2029 IF( ln_mixcpl ) THEN 2030 DO jl=1,jpl 2031 qml_ice(:,:,jl) = qml_ice(:,:,jl) * xcplmask(:,:,0) + frcv(jpr_topm)%z3(:,:,jl) * zmsk(:,:) 2032 qcn_ice(:,:,jl) = qcn_ice(:,:,jl) * xcplmask(:,:,0) + frcv(jpr_botm)%z3(:,:,jl) * zmsk(:,:) 2033 ENDDO 2034 ELSE 2035 qml_ice(:,:,:) = frcv(jpr_topm)%z3(:,:,:) 2036 qcn_ice(:,:,:) = frcv(jpr_botm)%z3(:,:,:) 2037 ENDIF 2017 2038 END SELECT 2018 !2019 2039 ! ! ========================= ! 2020 2040 ! ! Transmitted Qsr ! [W/m2] … … 2023 2043 ! 2024 2044 ! ! ===> used prescribed cloud fraction representative for polar oceans in summer (0.81) 2025 ztri = 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ! surface transmission parameter(Grenfell Maykut 77)2045 ztri = 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ! surface transmission when hi>10cm (Grenfell Maykut 77) 2026 2046 ! 2027 qtr_ice_top(:,:,:) = ztri * qsr_ice(:,:,:) 2028 WHERE( phs(:,:,:) >= 0.0_wp ) qtr_ice_top(:,:,:) = 0._wp ! snow fully opaque 2029 WHERE( phi(:,:,:) <= 0.1_wp ) qtr_ice_top(:,:,:) = qsr_ice(:,:,:) ! thin ice transmits all solar radiation 2047 WHERE ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 2048 zqtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ( ztri + ( 1._wp - ztri ) * ( 1._wp - phi(:,:,:) * 10._wp ) ) 2049 ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp ) ! constant (ztri) when hi>10cm 2050 zqtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ztri 2051 ELSEWHERE ! zero when hs>0 2052 zqtr_ice_top(:,:,:) = 0._wp 2053 END WHERE 2030 2054 ! 2031 2055 ELSEIF( ln_cndflx .AND. .NOT.ln_cndemulate ) THEN !== conduction flux as surface forcing ==! … … 2033 2057 ! ! ===> here we must receive the qtr_ice_top array from the coupler 2034 2058 ! for now just assume zero (fully opaque ice) 2035 qtr_ice_top(:,:,:) = 0._wp 2059 zqtr_ice_top(:,:,:) = 0._wp 2060 ! 2061 ENDIF 2062 ! 2063 IF( ln_mixcpl ) THEN 2064 DO jl=1,jpl 2065 qtr_ice_top(:,:,jl) = qtr_ice_top(:,:,jl) * xcplmask(:,:,0) + zqtr_ice_top(:,:,jl) * zmsk(:,:) 2066 ENDDO 2067 ELSE 2068 qtr_ice_top(:,:,:) = zqtr_ice_top(:,:,:) 2069 ENDIF 2070 ! ! ================== ! 2071 ! ! ice skin temp. ! 2072 ! ! ================== ! 2073 ! needed by Met Office 2074 IF( srcv(jpr_ts_ice)%laction ) THEN 2075 WHERE ( frcv(jpr_ts_ice)%z3(:,:,:) > 0.0 ) ; ztsu(:,:,:) = 0. + rt0 2076 ELSEWHERE( frcv(jpr_ts_ice)%z3(:,:,:) < -60. ) ; ztsu(:,:,:) = -60. + rt0 2077 ELSEWHERE ; ztsu(:,:,:) = frcv(jpr_ts_ice)%z3(:,:,:) + rt0 2078 END WHERE 2079 ! 2080 IF( ln_mixcpl ) THEN 2081 DO jl=1,jpl 2082 pist(:,:,jl) = pist(:,:,jl) * xcplmask(:,:,0) + ztsu(:,:,jl) * zmsk(:,:) 2083 ENDDO 2084 ELSE 2085 pist(:,:,:) = ztsu(:,:,:) 2086 ENDIF 2036 2087 ! 2037 2088 ENDIF … … 2197 2248 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_thick%clcat' ) 2198 2249 END SELECT 2199 IF( ssnd(jps_fice)%laction )CALL cpl_snd( jps_fice, isec, ztmp3, info )2250 CALL cpl_snd( jps_fice, isec, ztmp3, info ) 2200 2251 ENDIF 2201 2252 … … 2257 2308 ! ! Ice melt ponds ! 2258 2309 ! ! ------------------------- ! 2259 ! needed by Met Office 2310 ! needed by Met Office: 1) fraction of ponded ice 2) local/actual pond depth 2260 2311 IF( ssnd(jps_a_p)%laction .OR. ssnd(jps_ht_p)%laction ) THEN 2261 2312 SELECT CASE( sn_snd_mpnd%cldes) … … 2263 2314 SELECT CASE( sn_snd_mpnd%clcat ) 2264 2315 CASE( 'yes' ) 2265 ztmp3(:,:,1:jpl) = a_ip (:,:,1:jpl)2266 ztmp4(:,:,1:jpl) = v_ip(:,:,1:jpl)2316 ztmp3(:,:,1:jpl) = a_ip_frac(:,:,1:jpl) 2317 ztmp4(:,:,1:jpl) = h_ip(:,:,1:jpl) 2267 2318 CASE( 'no' ) 2268 2319 ztmp3(:,:,:) = 0.0 2269 2320 ztmp4(:,:,:) = 0.0 2270 2321 DO jl=1,jpl 2271 ztmp3(:,:,1) = ztmp3(:,:,1) + a_ip (:,:,jpl)2272 ztmp4(:,:,1) = ztmp4(:,:,1) + v_ip(:,:,jpl)2322 ztmp3(:,:,1) = ztmp3(:,:,1) + a_ip_frac(:,:,jpl) 2323 ztmp4(:,:,1) = ztmp4(:,:,1) + h_ip(:,:,jpl) 2273 2324 ENDDO 2274 2325 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_mpnd%clcat' ) -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/SBC/sbcrnf.F90
r12193 r12252 366 366 IF( h_rnf(ji,jj) > 0._wp ) THEN 367 367 jk = 2 368 DO WHILE ( jk /=mbkt(ji,jj) .AND. gdept_0(ji,jj,jk) < h_rnf(ji,jj) ) ; jk = jk + 1368 DO WHILE ( jk < mbkt(ji,jj) .AND. gdept_0(ji,jj,jk) < h_rnf(ji,jj) ) ; jk = jk + 1 369 369 END DO 370 370 nk_rnf(ji,jj) = jk … … 423 423 IF( zrnfcl(ji,jj,1) > 0._wp ) THEN 424 424 jk = 2 425 DO WHILE ( jk /=mbkt(ji,jj) .AND. gdept_0(ji,jj,jk) < h_rnf(ji,jj) ) ; jk = jk + 1425 DO WHILE ( jk < mbkt(ji,jj) .AND. gdept_0(ji,jj,jk) < h_rnf(ji,jj) ) ; jk = jk + 1 426 426 END DO 427 427 nk_rnf(ji,jj) = jk -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/nemogcm.F90
r12236 r12252 73 73 USE step_diu ! diurnal bulk SST timestepping (called from here if run offline) 74 74 USE crsini ! initialise grid coarsening utility 75 USE diatmb ! Top,middle,bottom output76 75 USE dia25h ! 25h mean output 77 76 USE diadetide ! Weights computation for daily detiding of model diagnostics … … 497 496 CALL trd_init( Nnn ) ! Mixed-layer/Vorticity/Integral constraints trends 498 497 CALL dia_obs_init( Nnn ) ! Initialize observational data 499 CALL dia_tmb_init ! TMB outputs500 498 CALL dia_25h_init( Nbb ) ! 25h mean outputs 501 499 CALL dia_detide_init ! Weights computation for daily detiding of model diagnostics -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p4zmeso.F90
r12236 r12252 72 72 REAL(wp) :: zgrazn, zgrazpoc, zgraznf, zgrazf 73 73 REAL(wp) :: zgrazfffp, zgrazfffg, zgrazffep, zgrazffeg 74 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing2, zfezoo2, zz2ligprod 74 75 CHARACTER (len=25) :: charout 75 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing, zfezoo276 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zw3d, zz2ligprod77 76 !!--------------------------------------------------------------------- 78 77 ! 79 78 IF( ln_timing ) CALL timing_start('p4z_meso') 80 !81 zgrazing(:,:,:) = 0._wp82 zfezoo2 (:,:,:) = 0._wp83 !84 IF (ln_ligand) THEN85 ALLOCATE( zz2ligprod(jpi,jpj,jpk) )86 zz2ligprod(:,:,:) = 0._wp87 ENDIF88 79 ! 89 80 DO jk = 1, jpkm1 … … 163 154 164 155 ! Total grazing ( grazing by microzoo is already computed in p4zmicro ) 165 zgrazing (ji,jj,jk) = zgraztotc156 zgrazing2(ji,jj,jk) = zgraztotc 166 157 167 158 ! Mesozooplankton efficiency … … 234 225 ! 235 226 IF( lk_iomput .AND. knt == nrdttrc ) THEN 236 ALLOCATE( zw3d(jpi,jpj,jpk) ) 237 IF( iom_use( "GRAZ2" ) ) THEN 238 zw3d(:,:,:) = zgrazing(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ! Total grazing of phyto by zooplankton 239 CALL iom_put( "GRAZ2", zw3d ) 227 zgrazing2(:,:,jpk) = 0._wp 228 zfezoo2 (:,:,jpk) = 0._wp 229 CALL iom_put( "GRAZ2" , zgrazing2(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ) ! Total grazing of phyto by zooplankton 230 CALL iom_put( "PCAL" , prodcal (:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ) ! Calcite production 231 CALL iom_put( "FEZOO2", zfezoo2 (:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 232 IF( ln_ligand ) THEN 233 zz2ligprod(:,:,jpk) = 0._wp 234 CALL iom_put( "LPRODZ2", zz2ligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 240 235 ENDIF 241 IF( iom_use( "PCAL" ) ) THEN242 zw3d(:,:,:) = prodcal(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ! Calcite production243 CALL iom_put( "PCAL", zw3d )244 ENDIF245 IF( iom_use( "FEZOO2" ) ) THEN246 zw3d(:,:,:) = zfezoo2(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) !247 CALL iom_put( "FEZOO2", zw3d )248 ENDIF249 IF( iom_use( "LPRODZ2" ) .AND. ln_ligand ) THEN250 zw3d(:,:,:) = zz2ligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:)251 CALL iom_put( "LPRODZ2" , zw3d )252 ENDIF253 DEALLOCATE( zw3d )254 236 ENDIF 255 !256 IF (ln_ligand) DEALLOCATE( zz2ligprod )257 237 ! 258 238 IF(sn_cfctl%l_prttrc) THEN ! print mean trends (used for debugging) -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p4zmicro.F90
r12236 r12252 70 70 REAL(wp) :: zgrazp, zgrazm, zgrazsd 71 71 REAL(wp) :: zgrazmf, zgrazsf, zgrazpf 72 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing, zfezoo 73 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zw3d, zzligprod 72 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing, zfezoo, zzligprod 74 73 CHARACTER (len=25) :: charout 75 74 !!--------------------------------------------------------------------- 76 75 ! 77 76 IF( ln_timing ) CALL timing_start('p4z_micro') 78 !79 IF (ln_ligand) THEN80 ALLOCATE( zzligprod(jpi,jpj,jpk) )81 zzligprod(:,:,:) = 0._wp82 ENDIF83 77 ! 84 78 DO jk = 1, jpkm1 … … 187 181 END DO 188 182 ! 189 IF( lk_iomput ) THEN 190 IF( knt == nrdttrc ) THEN 191 ALLOCATE( zw3d(jpi,jpj,jpk) ) 192 IF( iom_use( "GRAZ1" ) ) THEN 193 zw3d(:,:,:) = zgrazing(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ! Total grazing of phyto by zooplankton 194 CALL iom_put( "GRAZ1", zw3d ) 195 ENDIF 196 IF( iom_use( "FEZOO" ) ) THEN 197 zw3d(:,:,:) = zfezoo(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ! 198 CALL iom_put( "FEZOO", zw3d ) 199 ENDIF 200 IF( iom_use( "LPRODZ" ) .AND. ln_ligand ) THEN 201 zw3d(:,:,:) = zzligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) 202 CALL iom_put( "LPRODZ" , zw3d ) 203 ENDIF 204 DEALLOCATE( zw3d ) 183 IF( lk_iomput .AND. knt == nrdttrc ) THEN 184 zgrazing(:,:,jpk) = 0._wp 185 zfezoo (:,:,jpk) = 0._wp 186 CALL iom_put( "GRAZ1", zgrazing(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ) ! Total grazing of phyto by zooplankton 187 CALL iom_put( "FEZOO", zfezoo (:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 188 IF( ln_ligand ) THEN 189 zzligprod(:,:,jpk) = 0._wp 190 CALL iom_put( "LPRODZ", zzligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 205 191 ENDIF 206 192 ENDIF 207 !208 IF (ln_ligand) DEALLOCATE( zzligprod )209 193 ! 210 194 IF(sn_cfctl%l_prttrc) THEN ! print mean trends (used for debugging) -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/trcini.F90
r12236 r12252 58 58 IF(lwp) WRITE(numout,*) '~~~~~~~~' 59 59 ! 60 CALL trc_ini_ctl ! control61 60 CALL trc_nam ! read passive tracers namelists 62 61 CALL top_alloc() ! allocate TOP arrays … … 84 83 ! 85 84 END SUBROUTINE trc_init 86 87 88 SUBROUTINE trc_ini_ctl89 !!----------------------------------------------------------------------90 !! *** ROUTINE trc_ini_ctl ***91 !! ** Purpose : Control + ocean volume92 !!----------------------------------------------------------------------93 INTEGER :: jk ! dummy loop indices94 !95 ! Define logical parameter ton control dirunal cycle in TOP96 l_trcdm2dc = ln_dm2dc .OR. ( ln_cpl .AND. ncpl_qsr_freq /= 1 )97 l_trcdm2dc = l_trcdm2dc .AND. .NOT. l_offline98 IF( l_trcdm2dc .AND. lwp ) CALL ctl_warn( 'Coupling with passive tracers and used of diurnal cycle.', &99 & 'Computation of a daily mean shortwave for some biogeochemical models ' )100 !101 END SUBROUTINE trc_ini_ctl102 85 103 86 -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/trcstp.F90
r12236 r12252 72 72 ll_trcstat = ( sn_cfctl%l_trcstat ) .AND. & 73 73 & ( ( MOD( kt, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) ) 74 75 IF( kt == nittrc000 ) CALL trc_stp_ctl ! control 74 76 IF( kt == nittrc000 .AND. lk_trdmxl_trc ) CALL trd_mxl_trc_init ! trends: Mixed-layer 75 77 ! … … 131 133 ! 132 134 END SUBROUTINE trc_stp 135 136 137 SUBROUTINE trc_stp_ctl 138 !!---------------------------------------------------------------------- 139 !! *** ROUTINE trc_stp_ctl *** 140 !! ** Purpose : Control + ocean volume 141 !!---------------------------------------------------------------------- 142 ! 143 ! Define logical parameter ton control dirunal cycle in TOP 144 l_trcdm2dc = ln_dm2dc .OR. ( ln_cpl .AND. ncpl_qsr_freq /= 1 ) 145 l_trcdm2dc = l_trcdm2dc .AND. .NOT. l_offline 146 IF( l_trcdm2dc .AND. lwp ) CALL ctl_warn( 'Coupling with passive tracers and used of diurnal cycle.', & 147 & 'Computation of a daily mean shortwave for some biogeochemical models ' ) 148 ! 149 END SUBROUTINE trc_stp_ctl 150 133 151 134 152 SUBROUTINE trc_mean_qsr( kt ) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/CANAL/EXPREF/field_def_nemo-oce.xml
r12205 r12252 1 1 <?xml version="1.0"?> 2 2 <!-- $id$ --> 3 4 <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined --> 5 6 <!-- 7 ===================================================================================================== 8 = Configurable diagnostics = 9 ===================================================================================================== 10 --> 11 12 <field_group id="diamlr_fields"> 13 14 <!-- 15 ===================================================================================================== 16 Configuration of multiple-linear-regression analysis (diamlr) 17 ===================================================================================================== 18 19 This field group configures diamlr for tidal harmonic analysis of field 20 ssh: in addition to a regressor for fitting the mean value (diamlr_r101), 21 it includes the regressors for the analysis of the tidal constituents 22 that are available in the tidal-forcing implementation (see 23 ./src/OCE/SBC/tide.h90). 24 25 --> 26 27 <!-- Time --> 28 <field id="diamlr_time" grid_ref="diamlr_grid_T_2D" prec="8" /> 29 30 <!-- Regressors for tidal harmonic analysis --> 31 <field id="diamlr_r001" field_ref="diamlr_time" expr="sin( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M2" /> 32 <field id="diamlr_r002" field_ref="diamlr_time" expr="cos( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M2" /> 33 <field id="diamlr_r003" field_ref="diamlr_time" expr="sin( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:N2" /> 34 <field id="diamlr_r004" field_ref="diamlr_time" expr="cos( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:N2" /> 35 <field id="diamlr_r005" field_ref="diamlr_time" expr="sin( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:2N2" /> 36 <field id="diamlr_r006" field_ref="diamlr_time" expr="cos( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:2N2" /> 37 <field id="diamlr_r007" field_ref="diamlr_time" expr="sin( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S2" /> 38 <field id="diamlr_r008" field_ref="diamlr_time" expr="cos( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S2" /> 39 <field id="diamlr_r009" field_ref="diamlr_time" expr="sin( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K2" /> 40 <field id="diamlr_r010" field_ref="diamlr_time" expr="cos( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K2" /> 41 <field id="diamlr_r011" field_ref="diamlr_time" expr="sin( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K1" /> 42 <field id="diamlr_r012" field_ref="diamlr_time" expr="cos( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K1" /> 43 <field id="diamlr_r013" field_ref="diamlr_time" expr="sin( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:O1" /> 44 <field id="diamlr_r014" field_ref="diamlr_time" expr="cos( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:O1" /> 45 <field id="diamlr_r015" field_ref="diamlr_time" expr="sin( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Q1" /> 46 <field id="diamlr_r016" field_ref="diamlr_time" expr="cos( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Q1" /> 47 <field id="diamlr_r017" field_ref="diamlr_time" expr="sin( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:P1" /> 48 <field id="diamlr_r018" field_ref="diamlr_time" expr="cos( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:P1" /> 49 <field id="diamlr_r019" field_ref="diamlr_time" expr="sin( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M4" /> 50 <field id="diamlr_r020" field_ref="diamlr_time" expr="cos( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M4" /> 51 <field id="diamlr_r021" field_ref="diamlr_time" expr="sin( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mf" /> 52 <field id="diamlr_r022" field_ref="diamlr_time" expr="cos( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mf" /> 53 <field id="diamlr_r023" field_ref="diamlr_time" expr="sin( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mm" /> 54 <field id="diamlr_r024" field_ref="diamlr_time" expr="cos( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mm" /> 55 <field id="diamlr_r025" field_ref="diamlr_time" expr="sin( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Msqm" /> 56 <field id="diamlr_r026" field_ref="diamlr_time" expr="cos( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Msqm" /> 57 <field id="diamlr_r027" field_ref="diamlr_time" expr="sin( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mtm" /> 58 <field id="diamlr_r028" field_ref="diamlr_time" expr="cos( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mtm" /> 59 <field id="diamlr_r029" field_ref="diamlr_time" expr="sin( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S1" /> 60 <field id="diamlr_r030" field_ref="diamlr_time" expr="cos( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S1" /> 61 <field id="diamlr_r031" field_ref="diamlr_time" expr="sin( __TDE_MU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:MU2" /> 62 <field id="diamlr_r032" field_ref="diamlr_time" expr="cos( __TDE_MU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:MU2" /> 63 <field id="diamlr_r033" field_ref="diamlr_time" expr="sin( __TDE_NU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:NU2" /> 64 <field id="diamlr_r034" field_ref="diamlr_time" expr="cos( __TDE_NU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:NU2" /> 65 <field id="diamlr_r035" field_ref="diamlr_time" expr="sin( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:L2" /> 66 <field id="diamlr_r036" field_ref="diamlr_time" expr="cos( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:L2" /> 67 <field id="diamlr_r037" field_ref="diamlr_time" expr="sin( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:T2" /> 68 <field id="diamlr_r038" field_ref="diamlr_time" expr="cos( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:T2" /> 69 <field id="diamlr_r101" field_ref="diamlr_time" expr="diamlr_time^0.0" enabled=".TRUE." comment="mean" /> 70 71 <!-- Fields selected for regression analysis --> 72 <field id="diamlr_f001" field_ref="ssh" enabled=".TRUE." /> 73 74 </field_group> 3 75 4 76 <!-- … … 8 80 ============================================================================================================ 9 81 --> 10 <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined --> 82 83 <field_group id="diadetide_fields"> 84 85 <!-- 86 ===================================================================================================== 87 Weight fields for the computation of daily detided model diagnostics (diadetide) 88 ===================================================================================================== 89 90 --> 91 92 <field id="diadetide_weight" grid_ref="diadetide_grid_T_2D" enabled=".TRUE." /> 93 <field id="diadetide_weight_grid_T_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_T_2D" enabled=".TRUE." > this </field> 94 <field id="diadetide_weight_grid_U_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_U_2D" enabled=".TRUE." > this </field> 95 <field id="diadetide_weight_grid_V_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_V_2D" enabled=".TRUE." > this </field> 96 <field id="diadetide_weight_grid_T_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_T_3D" enabled=".TRUE." > this </field> 97 <field id="diadetide_weight_grid_U_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_U_3D" enabled=".TRUE." > this </field> 98 <field id="diadetide_weight_grid_V_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_V_3D" enabled=".TRUE." > this </field> 99 <field id="diadetide_weight_grid_W_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_W_3D" enabled=".TRUE." > this </field> 100 101 </field_group> 11 102 12 103 <!-- … … 19 110 20 111 <field_group id="grid_T" grid_ref="grid_T_2D" > 21 <field id="salgrad" long_name="module of salinity gradient" unit="psu/m" grid_ref="grid_T_3D"/> 22 <field id="salgrad2" long_name="square of module of salinity gradient" unit="psu2/m2" grid_ref="grid_T_3D"/> 23 <field id="ke" long_name="kinetic energy" unit="m2/s2" grid_ref="grid_T_3D"/> 24 <field id="ke_zint" long_name="vertical integration of kinetic energy" unit="m3/s2" /> 25 <field id="relvor" long_name="relative vorticity" unit="s-1" grid_ref="grid_T_3D"/> 26 <field id="absvor" long_name="absolute vorticity" unit="s-1" grid_ref="grid_T_3D"/> 27 <field id="potvor" long_name="potential vorticity" unit="s-1" grid_ref="grid_T_3D"/> 28 <field id="e3t" long_name="T-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_T_3D"/> 29 <field id="e3t_surf" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_SFC"/> 30 <field id="e3t_0" long_name="Initial T-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_T_3D"/> 31 32 <field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/> 33 <field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field > 34 <field id="soce" long_name="salinity" standard_name="sea_water_practical_salinity" unit="1e-3" grid_ref="grid_T_3D"/> 35 <field id="soce_e3t" long_name="salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce * e3t </field > 36 37 <!-- t-eddy viscosity coefficients (ldfdyn) --> 38 <field id="ahmt_2d" long_name=" surface t-eddy viscosity coefficient" unit="m2/s or m4/s" /> 39 <field id="ahmt_3d" long_name=" 3D t-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 40 41 <field id="sst" long_name="sea surface temperature" standard_name="sea_surface_temperature" unit="degC" /> 42 <field id="sst2" long_name="square of sea surface temperature" standard_name="square_of_sea_surface_temperature" unit="degC2" > sst * sst </field > 43 <field id="sstmax" long_name="max of sea surface temperature" field_ref="sst" operation="maximum" /> 44 <field id="sstmin" long_name="min of sea surface temperature" field_ref="sst" operation="minimum" /> 45 <field id="sstgrad" long_name="module of sst gradient" unit="degC/m" /> 46 <field id="sstgrad2" long_name="square of module of sst gradient" unit="degC2/m2" /> 47 <field id="sbt" long_name="sea bottom temperature" unit="degC" /> 48 <field id="tosmint" long_name="vertical integral of temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C" /> 49 <field id="sst_wl" long_name="Delta SST of warm layer" unit="degC" /> 50 <field id="sst_cs" long_name="Delta SST of cool skin" unit="degC" /> 51 <field id="temp_3m" long_name="temperature at 3m" unit="degC" /> 52 53 <field id="sss" long_name="sea surface salinity" standard_name="sea_surface_salinity" unit="1e-3" /> 54 <field id="sss2" long_name="square of sea surface salinity" unit="1e-6" > sss * sss </field > 55 <field id="sssmax" long_name="max of sea surface salinity" field_ref="sss" operation="maximum" /> 56 <field id="sssmin" long_name="min of sea surface salinity" field_ref="sss" operation="minimum" /> 57 <field id="sbs" long_name="sea bottom salinity" unit="0.001" /> 58 <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" /> 59 60 <field id="taubot" long_name="bottom stress module" unit="N/m2" /> 61 62 <field id="ssh" long_name="sea surface height" standard_name="sea_surface_height_above_geoid" unit="m" /> 63 <field id="ssh2" long_name="square of sea surface height" standard_name="square_of_sea_surface_height_above_geoid" unit="m2" > ssh * ssh </field > 64 <field id="wetdep" long_name="wet depth" standard_name="wet_depth" unit="m" /> 65 <field id="sshmax" long_name="max of sea surface height" field_ref="ssh" operation="maximum" /> 66 67 <field id="mldkz5" long_name="Turbocline depth (Kz = 5e-4)" standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity" unit="m" /> 68 <field id="mldr10_1" long_name="Mixed Layer Depth (dsigma = 0.01 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 69 <field id="mldr10_1max" long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="maximum" /> 70 <field id="mldr10_1min" long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="minimum" /> 71 <field id="heatc" long_name="Heat content vertically integrated" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" /> 72 <field id="saltc" long_name="Salt content vertically integrated" unit="1e-3*kg/m2" /> 73 <field id="salt2c" long_name="Salt content vertically integrated" unit="1e-3*kg/m2" /> 74 75 <!-- EOS --> 76 <field id="alpha" long_name="thermal expansion" unit="degC-1" grid_ref="grid_T_3D" /> 77 <field id="beta" long_name="haline contraction" unit="1e3" grid_ref="grid_T_3D" /> 78 <field id="bn2" long_name="squared Brunt-Vaisala frequency" unit="s-1" grid_ref="grid_T_3D" /> 79 <field id="rhop" long_name="potential density (sigma0)" standard_name="sea_water_sigma_theta" unit="kg/m3" grid_ref="grid_T_3D" /> 80 81 <!-- Energy - horizontal divergence --> 82 <field id="eken" long_name="kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_3D" /> 83 <field id="hdiv" long_name="horizontal divergence" unit="s-1" grid_ref="grid_T_3D" /> 84 85 <!-- variables available with MLE --> 86 <field id="Lf_NHpf" long_name="MLE: Lf = N H / f" unit="m" /> 87 88 <!-- next variables available with key_diahth --> 89 <field id="mlddzt" long_name="Thermocline Depth (depth of max dT/dz)" standard_name="depth_at_maximum_upward_derivative_of_sea_water_potential_temperature" unit="m" /> 90 <field id="mldr10_3" long_name="Mixed Layer Depth (dsigma = 0.03 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 91 <field id="mldr0_1" long_name="Mixed Layer Depth (dsigma = 0.01 wrt sfc)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 92 <field id="mldr0_3" long_name="Mixed Layer Depth (dsigma = 0.03 wrt sfc)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 93 <field id="mld_dt02" long_name="Mixed Layer Depth (|dT| = 0.2 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_temperature" unit="m" /> 94 <field id="topthdep" long_name="Top of Thermocline Depth (dT = -0.2 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_temperature" unit="m" /> 95 <field id="pycndep" long_name="Pycnocline Depth (dsigma[dT=-0.2] wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 96 <field id="BLT" long_name="Barrier Layer Thickness" unit="m" > topthdep - pycndep </field> 97 <field id="tinv" long_name="Max of vertical invertion of temperature" unit="degC" /> 98 <field id="depti" long_name="Depth of max. vert. inv. of temperature" unit="m" /> 99 <field id="20d" long_name="Depth of 20C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_20C" /> 100 <field id="28d" long_name="Depth of 28C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_28C" /> 101 <field id="hc300" long_name="Heat content 0-300m" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" /> 102 103 <!-- variables available with diaar5 --> 104 <field id="botpres" long_name="Sea Water Pressure at Sea Floor" standard_name="sea_water_pressure_at_sea_floor" unit="dbar" /> 105 <field id="sshdyn" long_name="dynamic sea surface height" standard_name="dynamic_sea_surface_height_above_geoid" unit="m" /> 106 <field id="sshdyn2" long_name="square of dynamic sea surface height" standard_name="dynamic_sea_surface_height_above_geoid_squared" unit="m2" > sshdyn * sshdyn </field> 107 <field id="tnpeo" long_name="Tendency of ocean potential energy content" unit="W/m2" /> 108 109 <!-- variables available ln_linssh=.FALSE. --> 110 <field id="tpt_dep" long_name="T-point depth" standard_name="depth_below_geoid" unit="m" grid_ref="grid_T_3D" /> 111 <field id="e3tdef" long_name="T-cell thickness deformation" unit="%" grid_ref="grid_T_3D" /> 112 113 <field id="salgrad" long_name="module of salinity gradient" unit="psu/m" grid_ref="grid_T_3D"/> 114 <field id="salgrad2" long_name="square of module of salinity gradient" unit="psu2/m2" grid_ref="grid_T_3D"/> 115 <field id="ke" long_name="kinetic energy" unit="m2/s2" grid_ref="grid_T_3D"/> 116 <field id="ke_zint" long_name="vertical integration of kinetic energy" unit="m3/s2" /> 117 <field id="relvor" long_name="relative vorticity" unit="s-1" grid_ref="grid_T_3D"/> 118 <field id="absvor" long_name="absolute vorticity" unit="s-1" grid_ref="grid_T_3D"/> 119 <field id="potvor" long_name="potential vorticity" unit="s-1" grid_ref="grid_T_3D"/> 120 <field id="sstgrad2" long_name="square of module of sst gradient" unit="degC2/m2" /> 121 122 <field id="e3t" long_name="T-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_T_3D" /> 123 <field id="e3ts" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_SFC"/> 124 <field id="e3t_0" long_name="Initial T-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_T_3D" /> 125 <field id="e3tb" long_name="bottom T-cell thickness" standard_name="bottom_cell_thickness" unit="m" grid_ref="grid_T_2D"/> 126 <field id="e3t_300" field_ref="e3t" grid_ref="grid_T_zoom_300" detect_missing_value="true" /> 127 <field id="e3t_vsum300" field_ref="e3t_300" grid_ref="grid_T_vsum" detect_missing_value="true" /> 128 <field id="masscello" long_name="Sea Water Mass per unit area" standard_name="sea_water_mass_per_unit_area" unit="kg/m2" grid_ref="grid_T_3D"/> 129 <field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D"/> 130 <field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/> 131 <field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field > 132 <field id="soce" long_name="salinity" standard_name="sea_water_practical_salinity" unit="1e-3" grid_ref="grid_T_3D"/> 133 <field id="soce_e3t" long_name="salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce * e3t </field > 134 135 <field id="toce_e3t_300" field_ref="toce_e3t" unit="degree_C" grid_ref="grid_T_zoom_300" detect_missing_value="true" /> 136 <field id="toce_e3t_vsum300" field_ref="toce_e3t_300" unit="degress_C*m" grid_ref="grid_T_vsum" detect_missing_value="true" /> 137 <field id="toce_vmean300" field_ref="toce_e3t_vsum300" unit="degree_C" grid_ref="grid_T_vsum" detect_missing_value="true" > toce_e3t_vsum300/e3t_vsum300 </field> 138 139 <!-- AGRIF sponge --> 140 <field id="agrif_spt" long_name=" AGRIF t-sponge coefficient" unit=" " /> 141 142 <!-- t-eddy viscosity coefficients (ldfdyn) --> 143 <field id="ahmt_2d" long_name=" surface t-eddy viscosity coefficient" unit="m2/s or m4/s" /> 144 <field id="ahmt_3d" long_name=" 3D t-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 145 146 <field id="sst" long_name="Bulk sea surface temperature" standard_name="bulk_sea_surface_temperature" unit="degC" /> 147 <field id="t_skin" long_name="Skin temperature aka SSST" standard_name="skin_temperature" unit="degC" /> 148 <field id="sst2" long_name="square of sea surface temperature" standard_name="square_of_sea_surface_temperature" unit="degC2" > sst * sst </field > 149 <field id="sstmax" long_name="max of sea surface temperature" field_ref="sst" operation="maximum" /> 150 <field id="sstmin" long_name="min of sea surface temperature" field_ref="sst" operation="minimum" /> 151 <field id="sstgrad" long_name="module of sst gradient" unit="degC/m" /> 152 <field id="sstgrad2" long_name="square of module of sst gradient" unit="degC2/m2" /> 153 <field id="sbt" long_name="sea bottom temperature" unit="degC" /> 154 <field id="tosmint" long_name="vertical integral of temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C" /> 155 <field id="sst_wl" long_name="Delta SST of warm layer" unit="degC" /> 156 <field id="sst_cs" long_name="Delta SST of cool skin" unit="degC" /> 157 <field id="temp_3m" long_name="temperature at 3m" unit="degC" /> 158 159 <field id="sss" long_name="sea surface salinity" standard_name="sea_surface_salinity" unit="1e-3" /> 160 <field id="sss2" long_name="square of sea surface salinity" unit="1e-6" > sss * sss </field > 161 <field id="sssmax" long_name="max of sea surface salinity" field_ref="sss" operation="maximum" /> 162 <field id="sssmin" long_name="min of sea surface salinity" field_ref="sss" operation="minimum" /> 163 <field id="sbs" long_name="sea bottom salinity" unit="0.001" /> 164 <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" /> 165 166 <field id="taubot" long_name="bottom stress module" unit="N/m2" /> 167 168 <!-- Case EOS = TEOS-10 : output potential temperature --> 169 <field id="toce_pot" long_name="Sea Water Potential Temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/> 170 <field id="sst_pot" long_name="potential sea surface temperature" standard_name="sea_surface_temperature" unit="degC" /> 171 <field id="tosmint_pot" long_name="vertical integral of potential temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C" /> 172 173 174 <field id="ssh" long_name="sea surface height" standard_name="sea_surface_height_above_geoid" unit="m" /> 175 <field id="ssh2" long_name="square of sea surface height" standard_name="square_of_sea_surface_height_above_geoid" unit="m2" > ssh * ssh </field > 176 <field id="wetdep" long_name="wet depth" standard_name="wet_depth" unit="m" /> 177 <field id="sshmax" long_name="max of sea surface height" field_ref="ssh" operation="maximum" /> 178 179 <field id="mldkz5" long_name="Turbocline depth (Kz = 5e-4)" standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity" unit="m" /> 180 <field id="mldr10_1" long_name="Mixed Layer Depth (dsigma = 0.01 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 181 <field id="mldr10_1max" long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="maximum" /> 182 <field id="mldr10_1min" long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="minimum" /> 183 <field id="heatc" long_name="Heat content vertically integrated" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" /> 184 <field id="saltc" long_name="Salt content vertically integrated" unit="1e-3*kg/m2" /> 185 186 <!-- EOS --> 187 <field id="alpha" long_name="thermal expansion" unit="degC-1" grid_ref="grid_T_3D" /> 188 <field id="beta" long_name="haline contraction" unit="1e3" grid_ref="grid_T_3D" /> 189 <field id="bn2" long_name="squared Brunt-Vaisala frequency" unit="s-1" grid_ref="grid_T_3D" /> 190 <field id="rhop" long_name="potential density (sigma0)" standard_name="sea_water_sigma_theta" unit="kg/m3" grid_ref="grid_T_3D" /> 191 192 <!-- Energy - horizontal divergence --> 193 <field id="eken" long_name="kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_3D" /> 194 <field id="hdiv" long_name="horizontal divergence" unit="s-1" grid_ref="grid_T_3D" /> 195 196 <!-- variables available with MLE --> 197 <field id="Lf_NHpf" long_name="MLE: Lf = N H / f" unit="m" /> 198 199 <!-- next variables available with ln_zad_Aimp=.true. --> 200 <field id="Courant" long_name="Courant number" unit="#" grid_ref="grid_T_3D" /> 201 <field id="wimp" long_name="Implicit vertical velocity" unit="m/s" grid_ref="grid_T_3D" /> 202 <field id="wexp" long_name="Explicit vertical velocity" unit="m/s" grid_ref="grid_T_3D" /> 203 <field id="wi_cff" long_name="Fraction of implicit vertical velocity" unit="#" grid_ref="grid_T_3D" /> 204 205 <!-- next variables available with key_diahth --> 206 <field id="mlddzt" long_name="Thermocline Depth (depth of max dT/dz)" standard_name="depth_at_maximum_upward_derivative_of_sea_water_potential_temperature" unit="m" /> 207 <field id="mldr10_3" long_name="Mixed Layer Depth (dsigma = 0.03 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 208 <field id="mldr0_1" long_name="Mixed Layer Depth (dsigma = 0.01 wrt sfc)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 209 <field id="mldr0_3" long_name="Mixed Layer Depth (dsigma = 0.03 wrt sfc)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 210 <field id="mld_dt02" long_name="Mixed Layer Depth (|dT| = 0.2 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_temperature" unit="m" /> 211 <field id="topthdep" long_name="Top of Thermocline Depth (dT = -0.2 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_temperature" unit="m" /> 212 <field id="pycndep" long_name="Pycnocline Depth (dsigma[dT=-0.2] wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" /> 213 <field id="BLT" long_name="Barrier Layer Thickness" unit="m" > topthdep - pycndep </field> 214 <field id="tinv" long_name="Max of vertical invertion of temperature" unit="degC" /> 215 <field id="depti" long_name="Depth of max. vert. inv. of temperature" unit="m" /> 216 <field id="20d" long_name="Depth of 20C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_20C" /> 217 <field id="26d" long_name="Depth of 26C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_26C" /> 218 <field id="28d" long_name="Depth of 28C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_28C" /> 219 <field id="hc300" long_name="Heat content 0-300m" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" /> 220 <field id="hc700" long_name="Heat content 0-700m" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" /> 221 <field id="hc2000" long_name="Heat content 0-2000m" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" /> 222 223 <!-- variables available with diaar5 --> 224 <field id="botpres" long_name="Sea Water Pressure at Sea Floor" standard_name="sea_water_pressure_at_sea_floor" unit="dbar" /> 225 <field id="sshdyn" long_name="dynamic sea surface height" standard_name="dynamic_sea_surface_height_above_geoid" unit="m" /> 226 <field id="sshdyn2" long_name="square of dynamic sea surface height" standard_name="dynamic_sea_surface_height_above_geoid_squared" unit="m2" > sshdyn * sshdyn </field> 227 <field id="tnpeo" long_name="Tendency of ocean potential energy content" unit="W/m2" /> 228 229 <!-- variables available ln_linssh=.FALSE. --> 230 <field id="tpt_dep" long_name="T-point depth" standard_name="depth_below_geoid" unit="m" grid_ref="grid_T_3D" /> 231 <field id="e3tdef" long_name="T-cell thickness deformation" unit="%" grid_ref="grid_T_3D" /> 232 233 <!-- variables available with ln_diacfl=.true. --> 234 <field id="cfl_cu" long_name="u-courant number" unit="#" /> 235 <field id="cfl_cv" long_name="v-courant number" unit="#" /> 236 <field id="cfl_cw" long_name="w-courant number" unit="#" /> 237 238 </field_group> <!-- grid_T --> 239 240 <!-- Tides --> 241 242 <field_group id="Tides_T" grid_ref="grid_T_2D" operation="once" > 243 <!-- Tidal potential --> 244 <field id="tide_pot" long_name="Total tidal potential" unit="m" /> 245 <field id="tide_pot_M2" long_name="M2 tidal potential" unit="m" /> 246 <field id="tide_pot_N2" long_name="N2 tidal potential" unit="m" /> 247 <field id="tide_pot_2N2" long_name="2N2 tidal potential" unit="m" /> 248 <field id="tide_pot_S2" long_name="S2 tidal potential" unit="m" /> 249 <field id="tide_pot_K2" long_name="K2 tidal potential" unit="m" /> 250 <field id="tide_pot_K1" long_name="K1 tidal potential" unit="m" /> 251 <field id="tide_pot_O1" long_name="O1 tidal potential" unit="m" /> 252 <field id="tide_pot_Q1" long_name="Q1 tidal potential" unit="m" /> 253 <field id="tide_pot_P1" long_name="P1 tidal potential" unit="m" /> 254 <field id="tide_pot_M4" long_name="M4 tidal potential" unit="m" /> 255 <field id="tide_pot_Mf" long_name="Mf tidal potential" unit="m" /> 256 <field id="tide_pot_Mm" long_name="Mm tidal potential" unit="m" /> 257 <field id="tide_pot_Msqm" long_name="Msqm tidal potential" unit="m" /> 258 <field id="tide_pot_Mtm" long_name="Mtm tidal potential" unit="m" /> 259 <field id="tide_pot_S1" long_name="S1 tidal potential" unit="m" /> 260 <field id="tide_pot_MU2" long_name="MU2 tidal potential" unit="m" /> 261 <field id="tide_pot_NU2" long_name="NU2 tidal potential" unit="m" /> 262 <field id="tide_pot_L2" long_name="L2 tidal potential" unit="m" /> 263 <field id="tide_pot_T2" long_name="T2 tidal potential" unit="m" /> 112 264 </field_group> 113 265 … … 115 267 116 268 <field_group id="OSMOSIS_T" grid_ref="grid_T_2D"> 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 269 <field id="zwth0" long_name="surface non-local temperature flux" unit="deg m/s" /> 270 <field id="zws0" long_name="surface non-local salinity flux" unit="psu m/s" /> 271 <field id="hbl" long_name="boundary layer depth" unit="m" /> 272 <field id="hbli" long_name="initial boundary layer depth" unit="m" /> 273 <field id="dstokes" long_name="stokes drift depth scale" unit="m" /> 274 <field id="zustke" long_name="magnitude of stokes drift at T-points" unit="m/s" /> 275 <field id="zwstrc" long_name="convective velocity scale" unit="m/s" /> 276 <field id="zwstrl" long_name="langmuir velocity scale" unit="m/s" /> 277 <field id="zustar" long_name="friction velocity" unit="m/s" /> 278 <field id="zhbl" long_name="boundary layer depth" unit="m" /> 279 <field id="zhml" long_name="mixed layer depth" unit="m" /> 280 <field id="wind_wave_abs_power" long_name="\rho |U_s| x u*^2" unit="mW" /> 281 <field id="wind_wave_power" long_name="U_s \dot tau" unit="mW" /> 282 <field id="wind_power" long_name="\rho u*^3" unit="mW" /> 283 284 <!-- extra OSMOSIS diagnostics --> 285 <field id="zwthav" long_name="av turb flux of T in ml" unit="deg m/s" /> 286 <field id="zt_ml" long_name="av T in ml" unit="deg" /> 287 <field id="zwth_ent" long_name="entrainment turb flux of T" unit="deg m/s" /> 288 <field id="zhol" long_name="Hoenekker number" unit="#" /> 289 <field id="zdh" long_name="Pycnocline depth - grid" unit=" m" /> 138 290 </field_group> 139 291 140 292 <field_group id="OSMOSIS_W" grid_ref="grid_W_3D" operation="instant" > 141 142 143 293 <field id="ghamt" long_name="non-local temperature flux" unit="deg m/s" /> 294 <field id="ghams" long_name="non-local salinity flux" unit="psu m/s" /> 295 <field id="zdtdz_pyc" long_name="Pycnocline temperature gradient" unit=" deg/m" /> 144 296 </field_group> 145 297 146 298 <field_group id="OSMOSIS_U" grid_ref="grid_U_2D" > 147 148 299 <field id="ghamu" long_name="non-local u-momentum flux" grid_ref="grid_U_3D" unit="m^2/s^2" /> 300 <field id="us_x" long_name="i component of Stokes drift" unit="m/s" /> 149 301 </field_group> 150 302 151 303 <field_group id="OSMOSIS_V" grid_ref="grid_V_2D" > 152 153 304 <field id="ghamv" long_name="non-local v-momentum flux" grid_ref="grid_V_3D" unit="m^2/s^2" /> 305 <field id="us_y" long_name="j component of Stokes drift" unit="m/s" /> 154 306 </field_group> 155 307 156 308 <!-- SBC --> 157 158 <field_group id="SBC" grid_ref="grid_T_2D" > <!-- time step automaticaly defined based on nn_fsbc --> 159 <field id="empmr" long_name="Net Upward Water Flux" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" /> 160 <field id="empbmr" long_name="Net Upward Water Flux at pre. tstep" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" /> 161 <field id="emp_oce" long_name="Evap minus Precip over ocean" standard_name="evap_minus_precip_over_sea_water" unit="kg/m2/s" /> 162 <field id="emp_ice" long_name="Evap minus Precip over ice" standard_name="evap_minus_precip_over_sea_ice" unit="kg/m2/s" /> 163 <field id="saltflx" long_name="Downward salt flux" unit="1e-3/m2/s" /> 164 <field id="fmmflx" long_name="Water flux due to freezing/melting" unit="kg/m2/s" /> 165 <field id="snowpre" long_name="Snow precipitation" standard_name="snowfall_flux" unit="kg/m2/s" /> 166 <field id="runoffs" long_name="River Runoffs" standard_name="water_flux_into_sea_water_from_rivers" unit="kg/m2/s" /> 167 <field id="precip" long_name="Total precipitation" standard_name="precipitation_flux" unit="kg/m2/s" /> 168 169 <field id="qt" long_name="Net Downward Heat Flux" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" /> 170 <field id="qns" long_name="non solar Downward Heat Flux" unit="W/m2" /> 171 <field id="qsr" long_name="Shortwave Radiation" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" /> 172 <field id="qsr3d" long_name="Shortwave Radiation 3D distribution" standard_name="downwelling_shortwave_flux_in_sea_water" unit="W/m2" grid_ref="grid_T_3D" /> 173 <field id="qrp" long_name="Surface Heat Flux: Damping" standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation" unit="W/m2" /> 174 <field id="erp" long_name="Surface Water Flux: Damping" standard_name="water_flux_out_of_sea_water_due_to_newtonian_relaxation" unit="kg/m2/s" /> 175 <field id="taum" long_name="wind stress module" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" /> 176 <field id="wspd" long_name="wind speed module" standard_name="wind_speed" unit="m/s" /> 177 178 <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 179 <field id="ssh_ib" long_name="Inverse barometer sea surface height" standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency" unit="m" /> 180 181 <!-- * variable related to ice shelf forcing * --> 182 <field id="fwfisf" long_name="Ice shelf melting" unit="kg/m2/s" /> 183 <field id="fwfisf3d" long_name="Ice shelf melting" unit="kg/m2/s" grid_ref="grid_T_3D" /> 184 <field id="qlatisf" long_name="Ice shelf latent heat flux" unit="W/m2" /> 185 <field id="qlatisf3d" long_name="Ice shelf latent heat flux" unit="W/m2" grid_ref="grid_T_3D" /> 186 <field id="qhcisf" long_name="Ice shelf heat content flux" unit="W/m2" /> 187 <field id="qhcisf3d" long_name="Ice shelf heat content flux" unit="W/m2" grid_ref="grid_T_3D" /> 188 <field id="isfgammat" long_name="transfert coefficient for isf (temperature) " unit="m/s" /> 189 <field id="isfgammas" long_name="transfert coefficient for isf (salinity) " unit="m/s" /> 190 <field id="stbl" long_name="salinity in the Losh tbl " unit="PSU" /> 191 <field id="ttbl" long_name="temperature in the Losh tbl " unit="C" /> 192 <field id="utbl" long_name="zonal current in the Losh tbl at T point " unit="m/s" /> 193 <field id="vtbl" long_name="merid current in the Losh tbl at T point " unit="m/s" /> 194 <field id="thermald" long_name="thermal driving of ice shelf melting " unit="C" /> 195 <field id="tfrz" long_name="top freezing point (used to compute melt) " unit="C" /> 196 <field id="tinsitu" long_name="top insitu temperature (used to cmpt melt) " unit="C" /> 197 <field id="ustar" long_name="ustar at T point used in ice shelf melting " unit="m/s" /> 198 199 <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 200 <field id="qlw_oce" long_name="Longwave Downward Heat Flux over open ocean" standard_name="surface_net_downward_longwave_flux" unit="W/m2" /> 201 <field id="qsb_oce" long_name="Sensible Downward Heat Flux over open ocean" standard_name="surface_downward_sensible_heat_flux" unit="W/m2" /> 202 <field id="qla_oce" long_name="Latent Downward Heat Flux over open ocean" standard_name="surface_downward_latent_heat_flux" unit="W/m2" /> 203 <field id="qt_oce" long_name="total flux at ocean surface" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" /> 204 <field id="qsr_oce" long_name="solar heat flux at ocean surface" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" /> 205 <field id="qns_oce" long_name="non-solar heat flux at ocean surface (including E-P)" unit="W/m2" /> 206 <field id="qemp_oce" long_name="Downward Heat Flux from E-P over open ocean" unit="W/m2" /> 207 <field id="taum_oce" long_name="wind stress module over open ocean" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" /> 208 209 <!-- available key_oasis3 --> 210 <field id="snow_ao_cea" long_name="Snow over ice-free ocean (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" /> 211 <field id="snow_ai_cea" long_name="Snow over sea-ice (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" /> 212 <field id="subl_ai_cea" long_name="Sublimation over sea-ice (cell average)" standard_name="surface_snow_and_ice_sublimation_flux" unit="kg/m2/s" /> 213 <field id="icealb_cea" long_name="Ice albedo (cell average)" standard_name="sea_ice_albedo" unit="1" /> 214 <field id="calving_cea" long_name="Calving" standard_name="water_flux_into_sea_water_from_icebergs" unit="kg/m2/s" /> 215 <field id="iceberg_cea" long_name="Iceberg" standard_name="water_flux_into_sea_water_from_icebergs" unit="kg/m2/s" /> 216 <field id="iceshelf_cea" long_name="Iceshelf" standard_name="water_flux_into_sea_water_from_iceshelf" unit="kg/m2/s" /> 217 218 219 <!-- available if key_oasis3 + conservative method --> 220 <field id="rain" long_name="Liquid precipitation" standard_name="rainfall_flux" unit="kg/m2/s" /> 221 <field id="evap_ao_cea" long_name="Evaporation over ice-free ocean (cell average)" standard_name="water_evaporation_flux" unit="kg/m2/s" /> 222 <field id="isnwmlt_cea" long_name="Snow over Ice melting (cell average)" standard_name="surface_snow_melt_flux" unit="kg/m2/s" /> 223 <field id="fsal_virt_cea" long_name="Virtual salt flux due to ice formation (cell average)" standard_name="virtual_salt_flux_into_sea_water_due_to_sea_ice_thermodynamics" unit="kg/m2/s" /> 224 <field id="fsal_real_cea" long_name="Real salt flux due to ice formation (cell average)" standard_name="downward_sea_ice_basal_salt_flux" unit="kg/m2/s" /> 225 <field id="hflx_rain_cea" long_name="heat flux due to rainfall" standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water" unit="W/m2" /> 226 <field id="hflx_evap_cea" long_name="heat flux due to evaporation" standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water" unit="W/m2" /> 227 <field id="hflx_snow_cea" long_name="heat flux due to snow falling" standard_name="heat_flux_onto_ocean_and_ice_due_to_snow_thermodynamics" unit="W/m2" /> 228 <field id="hflx_snow_ai_cea" long_name="heat flux due to snow falling over ice" standard_name="heat_flux_onto_ice_due_to_snow_thermodynamics" unit="W/m2" /> 229 <field id="hflx_snow_ao_cea" long_name="heat flux due to snow falling over ice-free ocean" standard_name="heat_flux_onto_sea_water_due_to_snow_thermodynamics" unit="W/m2" /> 230 <field id="hflx_ice_cea" long_name="heat flux due to ice thermodynamics" standard_name="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics" unit="W/m2" /> 231 <field id="hflx_rnf_cea" long_name="heat flux due to runoffs" standard_name="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water" unit="W/m2" /> 232 <field id="hflx_cal_cea" long_name="heat flux due to calving" standard_name="heat_flux_into_sea_water_due_to_calving" unit="W/m2" /> 233 <field id="hflx_icb_cea" long_name="heat flux due to iceberg" standard_name="heat_flux_into_sea_water_due_to_icebergs" unit="W/m2" /> 234 <field id="hflx_isf_cea" long_name="heat flux due to iceshelf" standard_name="heat_flux_into_sea_water_due_to_iceshelf" unit="W/m2" /> 235 <field id="bicemel_cea" long_name="Rate of Melt at Sea Ice Base (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_basal_melting" unit="kg/m2/s" /> 236 <field id="licepro_cea" long_name="Lateral Sea Ice Growth Rate (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_lateral_growth_of_ice_floes" unit="kg/m2/s" /> 237 <field id="snowmel_cea" long_name="Snow Melt Rate (cell average)" standard_name="surface_snow_melt_flux" unit="kg/m2/s" /> 238 <field id="sntoice_cea" long_name="Snow-Ice Formation Rate (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_snow_conversion" unit="kg/m2/s" /> 239 <field id="ticemel_cea" long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_surface_melting" unit="kg/m2/s" /> 240 241 <!-- ice field (nn_ice=1) --> 242 <field id="ice_cover" long_name="Ice fraction" standard_name="sea_ice_area_fraction" unit="1" /> 243 244 <!-- dilution --> 245 <field id="emp_x_sst" long_name="Concentration/Dilution term on SST" unit="kg*degC/m2/s" /> 246 <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" /> 247 <field id="rnf_x_sst" long_name="Runoff term on SST" unit="kg*degC/m2/s" /> 248 <field id="rnf_x_sss" long_name="Runoff term on SSS" unit="kg*1e-3/m2/s" /> 249 250 <!-- sbcssm variables --> 251 <field id="sst_m" unit="degC" /> 252 <field id="sss_m" unit="psu" /> 253 <field id="ssu_m" unit="m/s" /> 254 <field id="ssv_m" unit="m/s" /> 255 <field id="ssh_m" unit="m" /> 256 <field id="e3t_m" unit="m" /> 257 <field id="frq_m" unit="-" /> 258 309 <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc --> 310 311 <field_group id="SBC_2D" grid_ref="grid_T_2D" > 312 313 <field id="empmr" long_name="Net Upward Water Flux" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" /> 314 <field id="empbmr" long_name="Net Upward Water Flux at pre. tstep" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" /> 315 <field id="emp_oce" long_name="Evap minus Precip over ocean" standard_name="evap_minus_precip_over_sea_water" unit="kg/m2/s" /> 316 <field id="emp_ice" long_name="Evap minus Precip over ice" standard_name="evap_minus_precip_over_sea_ice" unit="kg/m2/s" /> 317 <field id="saltflx" long_name="Downward salt flux" unit="1e-3/m2/s" /> 318 <field id="fmmflx" long_name="Water flux due to freezing/melting" unit="kg/m2/s" /> 319 <field id="snowpre" long_name="Snow precipitation" standard_name="snowfall_flux" unit="kg/m2/s" /> 320 <field id="runoffs" long_name="River Runoffs" standard_name="water_flux_into_sea_water_from_rivers" unit="kg/m2/s" /> 321 <field id="precip" long_name="Total precipitation" standard_name="precipitation_flux" unit="kg/m2/s" /> 322 <field id="wclosea" long_name="closed sea empmr correction" standard_name="closea_empmr" unit="kg/m2/s" /> 323 324 <field id="qt" long_name="Net Downward Heat Flux" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" /> 325 <field id="qns" long_name="non solar Downward Heat Flux" unit="W/m2" /> 326 <field id="qsr" long_name="Shortwave Radiation" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" /> 327 <field id="qsr3d" long_name="Shortwave Radiation 3D distribution" standard_name="downwelling_shortwave_flux_in_sea_water" unit="W/m2" grid_ref="grid_T_3D" /> 328 <field id="qrp" long_name="Surface Heat Flux: Damping" standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation" unit="W/m2" /> 329 <field id="qclosea" long_name="closed sea heat content flux" standard_name="closea_heat_content_downward_flux" unit="W/m2" /> 330 <field id="erp" long_name="Surface Water Flux: Damping" standard_name="water_flux_out_of_sea_water_due_to_newtonian_relaxation" unit="kg/m2/s" /> 331 <field id="taum" long_name="wind stress module" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" /> 332 <field id="wspd" long_name="wind speed module" standard_name="wind_speed" unit="m/s" /> 333 334 <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 335 <field id="ssh_ib" long_name="Inverse barometer sea surface height" standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency" unit="m" /> 336 337 <!-- * variable related to ice shelf forcing * --> 338 <field id="isftfrz_cav" long_name="freezing point temperature at ocean/isf interface" unit="degC" /> 339 <field id="isftfrz_par" long_name="freezing point temperature in the parametrization boundary layer" unit="degC" /> 340 <field id="fwfisf_cav" long_name="Ice shelf melt rate" unit="kg/m2/s" /> 341 <field id="fwfisf_par" long_name="Ice shelf melt rate" unit="kg/m2/s" /> 342 <field id="qoceisf_cav" long_name="Ice shelf ocean heat flux" unit="W/m2" /> 343 <field id="qoceisf_par" long_name="Ice shelf ocean heat flux" unit="W/m2" /> 344 <field id="qlatisf_cav" long_name="Ice shelf latent heat flux" unit="W/m2" /> 345 <field id="qlatisf_par" long_name="Ice shelf latent heat flux" unit="W/m2" /> 346 <field id="qhcisf_cav" long_name="Ice shelf heat content flux of injected water" unit="W/m2" /> 347 <field id="qhcisf_par" long_name="Ice shelf heat content flux of injected water" unit="W/m2" /> 348 <field id="fwfisf3d_cav" long_name="Ice shelf melt rate" unit="kg/m2/s" grid_ref="grid_T_3D" /> 349 <field id="fwfisf3d_par" long_name="Ice shelf melt rate" unit="kg/m2/s" grid_ref="grid_T_3D" /> 350 <field id="qoceisf3d_cav" long_name="Ice shelf ocean heat flux" unit="W/m2" grid_ref="grid_T_3D" /> 351 <field id="qoceisf3d_par" long_name="Ice shelf ocean heat flux" unit="W/m2" grid_ref="grid_T_3D" /> 352 <field id="qlatisf3d_cav" long_name="Ice shelf latent heat flux" unit="W/m2" grid_ref="grid_T_3D" /> 353 <field id="qlatisf3d_par" long_name="Ice shelf latent heat flux" unit="W/m2" grid_ref="grid_T_3D" /> 354 <field id="qhcisf3d_cav" long_name="Ice shelf heat content flux of injected water" unit="W/m2" grid_ref="grid_T_3D" /> 355 <field id="qhcisf3d_par" long_name="Ice shelf heat content flux of injected water" unit="W/m2" grid_ref="grid_T_3D" /> 356 <field id="ttbl_cav" long_name="temperature in Losch tbl" unit="degC" /> 357 <field id="ttbl_par" long_name="temperature in the parametrisation boundary layer" unit="degC" /> 358 <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting" unit="degC" /> 359 <field id="isfthermald_par" long_name="thermal driving of ice shelf melting" unit="degC" /> 360 <field id="isfgammat" long_name="Ice shelf heat-transfert velocity" unit="m/s" /> 361 <field id="isfgammas" long_name="Ice shelf salt-transfert velocity" unit="m/s" /> 362 <field id="stbl" long_name="salinity in the Losh tbl" unit="1e-3" /> 363 <field id="utbl" long_name="zonal current in the Losh tbl at T point" unit="m/s" /> 364 <field id="vtbl" long_name="merid current in the Losh tbl at T point" unit="m/s" /> 365 <field id="isfustar" long_name="ustar at T point used in ice shelf melting" unit="m/s" /> 366 <field id="qconisf" long_name="Conductive heat flux through the ice shelf" unit="W/m2" /> 367 368 <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 369 <field id="rho_air" long_name="Air density at 10m above sea surface" standard_name="rho_air_10m" unit="kg/m3" /> 370 <field id="dt_skin" long_name="SSST-SST temperature difference" standard_name="SSST-SST" unit="K" /> 371 <field id="qlw_oce" long_name="Longwave Downward Heat Flux over open ocean" standard_name="surface_net_downward_longwave_flux" unit="W/m2" /> 372 <field id="qsb_oce" long_name="Sensible Downward Heat Flux over open ocean" standard_name="surface_downward_sensible_heat_flux" unit="W/m2" /> 373 <field id="qla_oce" long_name="Latent Downward Heat Flux over open ocean" standard_name="surface_downward_latent_heat_flux" unit="W/m2" /> 374 <field id="evap_oce" long_name="Evaporation over open ocean" standard_name="evaporation" unit="kg/m2/s" /> 375 <field id="qt_oce" long_name="total flux at ocean surface" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" /> 376 <field id="qsr_oce" long_name="solar heat flux at ocean surface" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" /> 377 <field id="qns_oce" long_name="non-solar heat flux at ocean surface (including E-P)" unit="W/m2" /> 378 <field id="qemp_oce" long_name="Downward Heat Flux from E-P over open ocean" unit="W/m2" /> 379 <field id="taum_oce" long_name="wind stress module over open ocean" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" /> 380 381 <!-- available key_oasis3 --> 382 <field id="snow_ao_cea" long_name="Snow over ice-free ocean (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" /> 383 <field id="snow_ai_cea" long_name="Snow over sea-ice (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" /> 384 <field id="subl_ai_cea" long_name="Sublimation over sea-ice (cell average)" standard_name="surface_snow_and_ice_sublimation_flux" unit="kg/m2/s" /> 385 <field id="icealb_cea" long_name="Ice albedo (cell average)" standard_name="sea_ice_albedo" unit="1" /> 386 <field id="calving_cea" long_name="Calving" standard_name="water_flux_into_sea_water_from_icebergs" unit="kg/m2/s" /> 387 <field id="iceberg_cea" long_name="Iceberg" standard_name="water_flux_into_sea_water_from_icebergs" unit="kg/m2/s" /> 388 <field id="iceshelf_cea" long_name="Iceshelf" standard_name="water_flux_into_sea_water_from_iceshelf" unit="kg/m2/s" /> 389 390 391 <!-- available if key_oasis3 + conservative method --> 392 <field id="rain" long_name="Liquid precipitation" standard_name="rainfall_flux" unit="kg/m2/s" /> 393 <field id="rain_ao_cea" long_name="Liquid precipitation over ice-free ocean (cell average)" standard_name="rainfall_flux" unit="kg/m2/s" /> 394 <field id="evap_ao_cea" long_name="Evaporation over ice-free ocean (cell average)" standard_name="water_evaporation_flux" unit="kg/m2/s" /> 395 <field id="isnwmlt_cea" long_name="Snow over Ice melting (cell average)" standard_name="surface_snow_melt_flux" unit="kg/m2/s" /> 396 <field id="fsal_virt_cea" long_name="Virtual salt flux due to ice formation (cell average)" standard_name="virtual_salt_flux_into_sea_water_due_to_sea_ice_thermodynamics" unit="kg/m2/s" /> 397 <field id="fsal_real_cea" long_name="Real salt flux due to ice formation (cell average)" standard_name="downward_sea_ice_basal_salt_flux" unit="kg/m2/s" /> 398 <field id="hflx_rain_cea" long_name="heat flux due to rainfall" standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water" unit="W/m2" /> 399 <field id="hflx_evap_cea" long_name="heat flux due to evaporation" standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water" unit="W/m2" /> 400 <field id="hflx_snow_cea" long_name="heat flux due to snow falling" standard_name="heat_flux_onto_ocean_and_ice_due_to_snow_thermodynamics" unit="W/m2" /> 401 <field id="hflx_snow_ai_cea" long_name="heat flux due to snow falling over ice" standard_name="heat_flux_onto_ice_due_to_snow_thermodynamics" unit="W/m2" /> 402 <field id="hflx_snow_ao_cea" long_name="heat flux due to snow falling over ice-free ocean" standard_name="heat_flux_onto_sea_water_due_to_snow_thermodynamics" unit="W/m2" /> 403 <field id="hflx_ice_cea" long_name="heat flux due to ice thermodynamics" standard_name="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics" unit="W/m2" /> 404 <field id="hflx_rnf_cea" long_name="heat flux due to runoffs" standard_name="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water" unit="W/m2" /> 405 <field id="hflx_cal_cea" long_name="heat flux due to calving" standard_name="heat_flux_into_sea_water_due_to_calving" unit="W/m2" /> 406 <field id="hflx_icb_cea" long_name="heat flux due to iceberg" standard_name="heat_flux_into_sea_water_due_to_icebergs" unit="W/m2" /> 407 <field id="hflx_isf_cea" long_name="heat flux due to iceshelf" standard_name="heat_flux_into_sea_water_due_to_iceshelf" unit="W/m2" /> 408 <field id="bicemel_cea" long_name="Rate of Melt at Sea Ice Base (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_basal_melting" unit="kg/m2/s" /> 409 <field id="licepro_cea" long_name="Lateral Sea Ice Growth Rate (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_lateral_growth_of_ice_floes" unit="kg/m2/s" /> 410 <field id="snowmel_cea" long_name="Snow Melt Rate (cell average)" standard_name="surface_snow_melt_flux" unit="kg/m2/s" /> 411 <field id="sntoice_cea" long_name="Snow-Ice Formation Rate (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_snow_conversion" unit="kg/m2/s" /> 412 <field id="ticemel_cea" long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_surface_melting" unit="kg/m2/s" /> 413 414 <!-- ice field (nn_ice=1) --> 415 <field id="ice_cover" long_name="Ice fraction" standard_name="sea_ice_area_fraction" unit="1" /> 416 417 <!-- dilution --> 418 <field id="emp_x_sst" long_name="Concentration/Dilution term on SST" unit="kg*degC/m2/s" /> 419 <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" /> 420 <field id="rnf_x_sst" long_name="Runoff term on SST" unit="kg*degC/m2/s" /> 421 <field id="rnf_x_sss" long_name="Runoff term on SSS" unit="kg*1e-3/m2/s" /> 422 423 <!-- sbcssm variables --> 424 <field id="sst_m" unit="degC" /> 425 <field id="sss_m" unit="psu" /> 426 <field id="ssu_m" unit="m/s" /> 427 <field id="ssv_m" unit="m/s" /> 428 <field id="ssh_m" unit="m" /> 429 <field id="e3t_m" unit="m" /> 430 <field id="frq_m" unit="-" /> 431 432 </field_group> 433 434 435 </field_group> <!-- SBC --> 436 437 <!-- ABL --> 438 <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> 439 440 <!-- variables available with ABL on atmospheric T grid--> 441 <field_group id="grid_ABL3D" grid_ref="grid_TA_3D" > 442 <field id="u_abl" long_name="ABL i-horizontal velocity" standard_name="abl_x_velocity" unit="m/s" /> 443 <field id="v_abl" long_name="ABL j-horizontal velocity" standard_name="abl_y_velocity" unit="m/s" /> 444 <field id="t_abl" long_name="ABL potential temperature" standard_name="abl_theta" unit="K" /> 445 <field id="q_abl" long_name="ABL specific humidity" standard_name="abl_qspe" unit="kg/kg" /> 446 <!-- debug (to be removed) --> 447 <field id="u_dta" long_name="DTA i-horizontal velocity" standard_name="dta_x_velocity" unit="m/s" /> 448 <field id="v_dta" long_name="DTA j-horizontal velocity" standard_name="dta_y_velocity" unit="m/s" /> 449 <field id="t_dta" long_name="DTA potential temperature" standard_name="dta_theta" unit="K" /> 450 <field id="q_dta" long_name="DTA specific humidity" standard_name="dta_qspe" unit="kg/kg" /> 451 <field id="coeft" long_name="ABL nudging coefficient" standard_name="coeft" unit="" /> 452 <field id="tke_abl" long_name="ABL turbulent kinetic energy" standard_name="abl_tke" unit="m2/s2" /> 453 <field id="avm_abl" long_name="ABL turbulent viscosity" standard_name="abl_avm" unit="m2/s" /> 454 <field id="avt_abl" long_name="ABL turbulent diffusivity" standard_name="abl_avt" unit="m2/s" /> 455 <field id="mxl_abl" long_name="ABL mixing length" standard_name="abl_mxl" unit="m" /> 456 </field_group> 457 458 <field_group id="grid_ABL2D" grid_ref="grid_TA_2D" > 459 <field id="pblh" long_name="ABL height" standard_name="abl_height" unit="m" /> 460 <field id="uz1_abl" long_name="ABL i-horizontal velocity" standard_name="abl_x_velocity" unit="m/s" /> 461 <field id="vz1_abl" long_name="ABL j-horizontal velocity" standard_name="abl_y_velocity" unit="m/s" /> 462 <field id="uvz1_abl" long_name="ABL wind speed module" standard_name="abl_wind_speed" unit="m/s" > sqrt( uz1_abl^2 + vz1_abl^2 ) </field> 463 <field id="tz1_abl" long_name="ABL potential temperature" standard_name="abl_theta" unit="K" /> 464 <field id="qz1_abl" long_name="ABL specific humidity" standard_name="abl_qspe" unit="kg/kg" /> 465 <field id="uz1_dta" long_name="DTA i-horizontal velocity" standard_name="dta_x_velocity" unit="m/s" /> 466 <field id="vz1_dta" long_name="DTA j-horizontal velocity" standard_name="dta_y_velocity" unit="m/s" /> 467 <field id="uvz1_dta" long_name="DTA wind speed module" standard_name="dta_wind_speed" unit="m/s" > sqrt( uz1_dta^2 + vz1_dta^2 ) </field> 468 <field id="tz1_dta" long_name="DTA potential temperature" standard_name="dta_theta" unit="K" /> 469 <field id="qz1_dta" long_name="DTA specific humidity" standard_name="dta_qspe" unit="kg/kg" /> 470 <!-- debug (to be removed) --> 471 <field id="uz1_geo" long_name="GEO i-horizontal velocity" standard_name="geo_x_velocity" unit="m/s" /> 472 <field id="vz1_geo" long_name="GEO j-horizontal velocity" standard_name="geo_y_velocity" unit="m/s" /> 473 <field id="uvz1_geo" long_name="GEO wind speed module" standard_name="geo_wind_speed" unit="m/s" > sqrt( uz1_geo^2 + vz1_geo^2 ) </field> 474 </field_group> 475 476 </field_group> <!-- ABL --> 477 478 479 <!-- U grid --> 480 481 <field_group id="grid_U" grid_ref="grid_U_2D"> 482 <field id="e2u" long_name="U-cell width in meridional direction" standard_name="cell_width" unit="m" /> 483 <field id="e3u" long_name="U-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_U_3D" /> 484 <field id="e3u_0" long_name="Initial U-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_U_3D"/> 485 <field id="utau" long_name="Wind Stress along i-axis" standard_name="surface_downward_x_stress" unit="N/m2" /> 486 <field id="uoce" long_name="ocean current along i-axis" standard_name="sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 487 <field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field> 488 <field id="uoce_e3u_vsum" long_name="ocean current along i-axis * e3u summed on the vertical" field_ref="uoce_e3u" unit="m3/s" grid_ref="grid_U_vsum"/> 489 <field id="uocetr_vsum" long_name="ocean transport along i-axis summed on the vertical" field_ref="e2u" unit="m3/s"> this * uoce_e3u_vsum </field> 490 491 <field id="uocetr_vsum_op" long_name="ocean current along i-axis * e3u * e2u summed on the vertical" read_access="true" freq_op="1mo" field_ref="e2u" unit="m3/s"> @uocetr_vsum </field> 492 <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo" unit="m3/s" /> 493 <field id="msftbarot" long_name="ocean_barotropic_mass_streamfunction" unit="kg s-1" > uocetr_vsum_cumul * $rau0 </field> 494 495 496 <field id="ssu" long_name="ocean surface current along i-axis" unit="m/s" /> 497 <field id="sbu" long_name="ocean bottom current along i-axis" unit="m/s" /> 498 <field id="ubar" long_name="ocean barotropic current along i-axis" unit="m/s" /> 499 <field id="uocetr_eff" long_name="Effective ocean transport along i-axis" standard_name="ocean_volume_x_transport" unit="m3/s" grid_ref="grid_U_3D" /> 500 <field id="uocet" long_name="ocean transport along i-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_U_3D" /> 501 <field id="uoces" long_name="ocean transport along i-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_U_3D" /> 502 <field id="ssuww" long_name="ocean surface wind work along i-axis" standard_name="surface_x_wind_work" unit="N/m*s" > utau * ssu </field> 503 <!-- AGRIF sponge --> 504 <field id="agrif_spu" long_name=" AGRIF u-sponge coefficient" unit=" " /> 505 <!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 506 <field id="ahtu_2d" long_name=" surface u-eddy diffusivity coefficient" unit="m2/s or m4/s" /> 507 <field id="ahtu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s or m4/s" grid_ref="grid_U_3D"/> 508 <!-- u-eiv diffusivity coefficients (available if ln_ldfeiv=F) --> 509 <field id="aeiu_2d" long_name=" surface u-EIV coefficient" unit="m2/s" /> 510 <field id="aeiu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s" grid_ref="grid_U_3D"/> 511 512 <!-- variables available with MLE (ln_mle=T) --> 513 <field id="psiu_mle" long_name="MLE streamfunction along i-axis" unit="m3/s" grid_ref="grid_U_3D" /> 514 515 <!-- uoce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 516 <field id="uoce_eiv" long_name="EIV ocean current along i-axis" standard_name="bolus_sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 517 <field id="ueiv_masstr" long_name="EIV Ocean Mass X Transport" standard_name="bolus_ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" /> 518 <field id="ueiv_heattr" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" /> 519 <field id="ueiv_salttr" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" /> 520 <field id="ueiv_heattr3d" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" grid_ref="grid_U_3D" /> 521 <field id="ueiv_salttr3d" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_U_3D" /> 522 523 <!-- uoce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 524 <field id="uoce_bbl" long_name="BBL ocean current along i-axis" unit="m/s" /> 525 <!-- ahu_bbl : available with ln_trabbl=T and nn_bbl_ldf=1 --> 526 <field id="ahu_bbl" long_name="BBL diffusive flux along i-axis" unit="m3/s" /> 527 528 <!-- variables available with WAVE (ln_wave=T) --> 529 <field id="ustokes" long_name="Stokes Drift Velocity i-axis" standard_name="StokesDrift_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 530 <field id="ustokes_e3u" long_name="Stokes Drift Velocity i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > ustokes * e3u </field> 531 532 <!-- variable for ice shelves --> 533 <field id="utbl" long_name="zonal current in the Losh tbl" unit="m/s" /> 534 535 <!-- variables available with diaar5 --> 536 <field id="u_masstr" long_name="Ocean Mass X Transport" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" /> 537 <field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport" unit="kg/s" /> 538 <field id="u_heattr" long_name="ocean eulerian heat transport along i-axis" standard_name="ocean_heat_x_transport" unit="W" /> 539 <field id="u_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_x_transport" unit="1e-3*kg/s" /> 540 <field id="uadv_heattr" long_name="ocean advective heat transport along i-axis" standard_name="advectice_ocean_heat_x_transport" unit="W" /> 541 <field id="uadv_salttr" long_name="ocean advective salt transport along i-axis" standard_name="advectice_ocean_salt_x_transport" unit="1e-3*kg/s" /> 542 <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_diffusion" unit="W" /> 543 <field id="udiff_salttr" long_name="ocean diffusion salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_diffusion" unit="1e-3*kg/s" /> 259 544 </field_group> 260 261 <!-- U grid --> 262 263 <field_group id="grid_U" grid_ref="grid_U_2D"> 264 <field id="e3u" long_name="U-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_U_3D" /> 265 <field id="e3u_0" long_name="Initial U-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_U_3D"/> 266 <field id="utau" long_name="Wind Stress along i-axis" standard_name="surface_downward_x_stress" unit="N/m2" /> 267 <field id="uoce" long_name="ocean current along i-axis" standard_name="sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 268 <field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field> 269 <field id="ssu" long_name="ocean surface current along i-axis" unit="m/s" /> 270 <field id="sbu" long_name="ocean bottom current along i-axis" unit="m/s" /> 271 <field id="ubar" long_name="ocean barotropic current along i-axis" unit="m/s" /> 272 <field id="uocetr_eff" long_name="Effective ocean transport along i-axis" standard_name="ocean_volume_x_transport" unit="m3/s" grid_ref="grid_U_3D" /> 273 <field id="uocet" long_name="ocean transport along i-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_U_3D" /> 274 <field id="uoces" long_name="ocean transport along i-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_U_3D" /> 275 276 <!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 277 <field id="ahtu_2d" long_name=" surface u-eddy diffusivity coefficient" unit="m2/s or m4/s" /> 278 <field id="ahtu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s or m4/s" grid_ref="grid_U_3D"/> 279 <!-- u-eiv diffusivity coefficients (available if ln_ldfeiv=F) --> 280 <field id="aeiu_2d" long_name=" surface u-EIV coefficient" unit="m2/s" /> 281 <field id="aeiu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s" grid_ref="grid_U_3D"/> 282 283 <!-- variables available with MLE (ln_mle=T) --> 284 <field id="psiu_mle" long_name="MLE streamfunction along i-axis" unit="m3/s" grid_ref="grid_U_3D" /> 285 286 <!-- uoce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 287 <field id="uoce_eiv" long_name="EIV ocean current along i-axis" standard_name="bolus_sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 288 289 <!-- uoce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 290 <field id="uoce_bbl" long_name="BBL ocean current along i-axis" unit="m/s" /> 291 <!-- ahu_bbl : available with ln_trabbl=T and nn_bbl_ldf=1 --> 292 <field id="ahu_bbl" long_name="BBL diffusive flux along i-axis" unit="m3/s" /> 293 294 <!-- variable for ice shelves --> 295 <field id="utbl" long_name="zonal current in the Losh tbl" unit="m/s" /> 296 297 <field id="u_masstr" long_name="Ocean Mass X Transport" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" /> 298 <field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport" unit="kg/s" /> 299 <field id="u_heattr" long_name="ocean eulerian heat transport along i-axis" standard_name="ocean_heat_x_transport" unit="W" /> 300 <field id="u_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_x_transport" unit="1e-3*kg/s" /> 301 <field id="uadv_heattr" long_name="ocean advective heat transport along i-axis" standard_name="advectice_ocean_heat_x_transport" unit="W" /> 302 <field id="uadv_salttr" long_name="ocean advective salt transport along i-axis" standard_name="advectice_ocean_salt_x_transport" unit="1e-3*kg/s" /> 303 <field id="ueiv_heattr" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" /> 304 <field id="ueiv_salttr" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" /> 305 <field id="ueiv_heattr3d" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" grid_ref="grid_U_3D" /> 306 <field id="ueiv_salttr3d" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_U_3D" /> 307 <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_diffusion" unit="W" /> 308 <field id="udiff_salttr" long_name="ocean diffusion salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_diffusion" unit="1e-3*kg/s" /> 545 546 <!-- V grid --> 547 548 <field_group id="grid_V" grid_ref="grid_V_2D"> 549 <field id="e1v" long_name="V-cell width in longitudinal direction" standard_name="cell_width" unit="m" /> 550 <field id="e3v" long_name="V-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_V_3D" /> 551 <field id="e3v_0" long_name="Initial V-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_V_3D" /> 552 <field id="vtau" long_name="Wind Stress along j-axis" standard_name="surface_downward_y_stress" unit="N/m2" /> 553 <field id="voce" long_name="ocean current along j-axis" standard_name="sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 554 <field id="voce_e3v" long_name="ocean current along j-axis (thickness weighted)" unit="m/s" grid_ref="grid_V_3D" > voce * e3v </field> 555 <field id="ssv" long_name="ocean surface current along j-axis" unit="m/s" /> 556 <field id="sbv" long_name="ocean bottom current along j-axis" unit="m/s" /> 557 <field id="vbar" long_name="ocean barotropic current along j-axis" unit="m/s" /> 558 <field id="vocetr_eff" long_name="Effective ocean transport along j-axis" standard_name="ocean_volume_y_transport" unit="m3/s" grid_ref="grid_V_3D" /> 559 <field id="vocet" long_name="ocean transport along j-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_V_3D" /> 560 <field id="voces" long_name="ocean transport along j-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_V_3D" /> 561 <field id="ssvww" long_name="ocean surface wind work along j-axis" standard_name="surface_y_wind_work" unit="N/m*s" > vtau * ssv </field> 562 <!-- AGRIF sponge --> 563 <field id="agrif_spv" long_name=" AGRIF v-sponge coefficient" unit=" " /> 564 <!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 565 <field id="ahtv_2d" long_name=" surface v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" /> 566 <field id="ahtv_3d" long_name=" 3D v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" grid_ref="grid_V_3D"/> 567 <!-- v-eiv diffusivity coefficients (available if ln_ldfeiv=F) --> 568 <field id="aeiv_2d" long_name=" surface v-EIV coefficient" unit="m2/s" /> 569 <field id="aeiv_3d" long_name=" 3D v-EIV coefficient" unit="m2/s" grid_ref="grid_V_3D" /> 570 571 <!-- variables available with MLE (ln_mle=T) --> 572 <field id="psiv_mle" long_name="MLE streamfunction along j-axis" unit="m3/s" grid_ref="grid_V_3D" /> 573 574 <!-- voce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 575 <field id="voce_eiv" long_name="EIV ocean current along j-axis" standard_name="bolus_sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 576 <field id="veiv_masstr" long_name="EIV Ocean Mass Y Transport" standard_name="bolus_ocean_mass_y_transport" unit="kg/s" grid_ref="grid_V_3D" /> 577 <field id="veiv_heattr" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" /> 578 <field id="veiv_salttr" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" /> 579 <field id="veiv_heattr3d" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" grid_ref="grid_V_3D" /> 580 <field id="veiv_salttr3d" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_V_3D" /> 581 582 583 <!-- voce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 584 <field id="voce_bbl" long_name="BBL ocean current along j-axis" unit="m/s" /> 585 <!-- ahu_bbl : available with ln_trabbl=T and nn_bbl_ldf=1 --> 586 <field id="ahv_bbl" long_name="BBL diffusive flux along j-axis" unit="m3/s" /> 587 588 <!-- variables available with WAVE (ln_wave=T) --> 589 <field id="vstokes" long_name="Stokes Drift Velocity j-axis" standard_name="StokesDrift_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 590 <field id="vstokes_e3v" long_name="Stokes Drift Velocity j-axis (thickness weighted)" unit="m/s" grid_ref="grid_V_3D" > vstokes * e3v </field> 591 592 <!-- variable for ice shelves --> 593 <field id="vtbl" long_name="meridional current in the Losh tbl" unit="m/s" /> 594 595 <!-- variables available with diaar5 --> 596 <field id="v_masstr" long_name="ocean eulerian mass transport along j-axis" standard_name="ocean_mass_y_transport" unit="kg/s" grid_ref="grid_V_3D" /> 597 <field id="v_heattr" long_name="ocean eulerian heat transport along j-axis" standard_name="ocean_heat_y_transport" unit="W" /> 598 <field id="v_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_y_transport" unit="1e-3*kg/s" /> 599 <field id="vadv_heattr" long_name="ocean advective heat transport along j-axis" standard_name="advectice_ocean_heat_y_transport" unit="W" /> 600 <field id="vadv_salttr" long_name="ocean advective salt transport along j-axis" standard_name="advectice_ocean_salt_y_transport" unit="1e-3*kg/s" /> 601 <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_diffusion" unit="W" /> 602 <field id="vdiff_salttr" long_name="ocean diffusion salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_diffusion" unit="1e-3*kg/s" /> 309 603 </field_group> 310 604 311 <!-- V grid --> 312 313 <field_group id="grid_V" grid_ref="grid_V_2D"> 314 <field id="e3v" long_name="V-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_V_3D" /> 315 <field id="e3v_0" long_name="Initial V-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_V_3D"/> 316 <field id="vtau" long_name="Wind Stress along j-axis" standard_name="surface_downward_y_stress" unit="N/m2" /> 317 <field id="voce" long_name="ocean current along j-axis" standard_name="sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 318 <field id="voce_e3v" long_name="ocean current along j-axis (thickness weighted)" unit="m/s" grid_ref="grid_V_3D" > voce * e3v </field> 319 <field id="ssv" long_name="ocean surface current along j-axis" unit="m/s" /> 320 <field id="sbv" long_name="ocean bottom current along j-axis" unit="m/s" /> 321 <field id="vbar" long_name="ocean barotropic current along j-axis" unit="m/s" /> 322 <field id="vocetr_eff" long_name="Effective ocean transport along j-axis" standard_name="ocean_volume_y_transport" unit="m3/s" grid_ref="grid_V_3D" /> 323 <field id="vocet" long_name="ocean transport along j-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_V_3D" /> 324 <field id="voces" long_name="ocean transport along j-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_V_3D" /> 325 326 <!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 327 <field id="ahtv_2d" long_name=" surface v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" /> 328 <field id="ahtv_3d" long_name=" 3D v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" grid_ref="grid_V_3D"/> 329 <!-- v-eiv diffusivity coefficients (available if ln_ldfeiv=F) --> 330 <field id="aeiv_2d" long_name=" surface v-EIV coefficient" unit="m2/s" /> 331 <field id="aeiv_3d" long_name=" 3D v-EIV coefficient" unit="m2/s" grid_ref="grid_V_3D" /> 332 333 <!-- variables available with MLE (ln_mle=T) --> 334 <field id="psiv_mle" long_name="MLE streamfunction along j-axis" unit="m3/s" grid_ref="grid_V_3D" /> 335 336 <!-- voce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 337 <field id="voce_eiv" long_name="EIV ocean current along j-axis" standard_name="bolus_sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 338 339 <!-- voce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 340 <field id="voce_bbl" long_name="BBL ocean current along j-axis" unit="m/s" /> 341 <!-- ahu_bbl : available with ln_trabbl=T and nn_bbl_ldf=1 --> 342 <field id="ahv_bbl" long_name="BBL diffusive flux along j-axis" unit="m3/s" /> 343 344 <!-- variable for ice shelves --> 345 <field id="vtbl" long_name="meridional current in the Losh tbl" unit="m/s" /> 346 347 <!-- variables available with diaar5 --> 348 <field id="v_masstr" long_name="ocean eulerian mass transport along j-axis" standard_name="ocean_mass_y_transport" unit="kg/s" grid_ref="grid_V_3D" /> 349 <field id="v_heattr" long_name="ocean eulerian heat transport along j-axis" standard_name="ocean_heat_y_transport" unit="W" /> 350 <field id="v_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_y_transport" unit="1e-3*kg/s" /> 351 <field id="vadv_heattr" long_name="ocean advective heat transport along j-axis" standard_name="advectice_ocean_heat_y_transport" unit="W" /> 352 <field id="vadv_salttr" long_name="ocean advective salt transport along j-axis" standard_name="advectice_ocean_salt_y_transport" unit="1e-3*kg/s" /> 353 <field id="veiv_heattr" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" /> 354 <field id="veiv_salttr" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" /> 355 <field id="veiv_heattr3d" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" grid_ref="grid_V_3D" /> 356 <field id="veiv_salttr3d" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_V_3D" /> 357 <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_diffusion" unit="W" /> 358 <field id="vdiff_salttr" long_name="ocean diffusion salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_diffusion" unit="1e-3*kg/s" /> 605 <!-- W grid --> 606 607 <field_group id="grid_W" grid_ref="grid_W_3D"> 608 <field id="e3w" long_name="W-cell thickness" standard_name="cell_thickness" unit="m" /> 609 <field id="woce" long_name="ocean vertical velocity" standard_name="upward_sea_water_velocity" unit="m/s" /> 610 <field id="woce_e3w" long_name="ocean vertical velocity * e3w" unit="m2/s" > woce * e3w </field> 611 <field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s" /> 612 613 <!-- woce_eiv: available with EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 614 <field id="woce_eiv" long_name="EIV ocean vertical velocity" standard_name="bolus_upward_sea_water_velocity" unit="m/s" /> 615 <field id="weiv_masstr" long_name="EIV Upward Ocean Mass Transport" standard_name="bolus_upward_ocean_mass_transport" unit="kg/s" /> 616 <field id="weiv_heattr3d" long_name="ocean bolus heat transport" standard_name="ocean_heat_z_transport_due_to_bolus_advection" unit="W" /> 617 <field id="weiv_salttr3d" long_name="ocean bolus salt transport" standard_name="ocean_salt_z_transport_due_to_bolus_advection" unit="kg" /> 618 619 <field id="avt" long_name="vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 620 <field id="avt_e3w" long_name="vertical heat diffusivity * e3w" unit="m3/s" > avt * e3w </field> 621 <field id="logavt" long_name="logarithm of vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 622 <field id="avm" long_name="vertical eddy viscosity" standard_name="ocean_vertical_momentum_diffusivity" unit="m2/s" /> 623 <field id="avm_e3w" long_name="vertical eddy viscosity * e3w" unit="m3/s" > avm * e3w </field> 624 625 <!-- avs: /= avt with ln_zdfddm=T --> 626 <field id="avs" long_name="salt vertical eddy diffusivity" standard_name="ocean_vertical_salt_diffusivity" unit="m2/s" /> 627 <field id="avs_e3w" long_name="vertical salt diffusivity * e3w" unit="m3/s" > avs * e3w </field> 628 <field id="logavs" long_name="logarithm of salt vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 629 630 <!-- avt_evd and avm_evd: available with ln_zdfevd --> 631 <field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" /> 632 <field id="avt_evd_e3w" long_name="convective enhancement to vertical diffusivity * e3w " unit="m3/s" > avt_evd * e3w </field> 633 <field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" /> 634 635 <!-- avt_tide: available with ln_zdfiwm=T --> 636 <field id="av_ratio" long_name="S over T diffusivity ratio" standard_name="salinity_over_temperature_diffusivity_ratio" unit="1" /> 637 <field id="av_wave" long_name="internal wave-induced vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves" unit="m2/s" /> 638 <field id="bflx_iwm" long_name="internal wave-induced buoyancy flux" standard_name="buoyancy_flux_due_to_internal_waves" unit="W/kg" /> 639 <field id="pcmap_iwm" long_name="power consumed by wave-driven mixing" standard_name="vertically_integrated_power_consumption_by_wave_driven_mixing" unit="W/m2" grid_ref="grid_W_2D" /> 640 <field id="emix_iwm" long_name="power density available for mixing" standard_name="power_available_for_mixing_from_breaking_internal_waves" unit="W/kg" /> 641 642 <!-- variables available with WAVE (ln_wave=T) --> 643 <field id="wstokes" long_name="Stokes Drift vertical velocity" standard_name="upward_StokesDrift_velocity" unit="m/s" /> 644 645 <!-- variables available with diaar5 --> 646 <field id="w_masstr" long_name="vertical mass transport" standard_name="upward_ocean_mass_transport" unit="kg/s" /> 647 <field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" /> 648 359 649 </field_group> 360 650 361 <!-- W grid -->362 363 <field_group id="grid_W" grid_ref="grid_W_3D">364 <field id="e3w" long_name="W-cell thickness" standard_name="cell_thickness" unit="m" />365 <field id="woce" long_name="ocean vertical velocity" standard_name="upward_sea_water_velocity" unit="m/s" />366 <field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s" />367 368 <!-- woce_eiv: available with EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) -->369 <field id="woce_eiv" long_name="EIV ocean vertical velocity" standard_name="bolus_upward_sea_water_velocity" unit="m/s" />370 371 <field id="avt" long_name="vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />372 <field id="logavt" long_name="logarithm of vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />373 <field id="avm" long_name="vertical eddy viscosity" standard_name="ocean_vertical_momentum_diffusivity" unit="m2/s" />374 375 <!-- avs: /= avt with ln_zdfddm=T -->376 <field id="avs" long_name="salt vertical eddy diffusivity" standard_name="ocean_vertical_salt_diffusivity" unit="m2/s" />377 <field id="logavs" long_name="logarithm of salt vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />378 379 <!-- avt_evd and avm_evd: available with ln_zdfevd -->380 <field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" />381 <field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" />382 383 <!-- avt_tide: available with ln_zdfiwm=T -->384 <field id="av_ratio" long_name="S over T diffusivity ratio" standard_name="salinity_over_temperature_diffusivity_ratio" unit="1" />385 <field id="av_wave" long_name="internal wave-induced vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves" unit="m2/s" />386 <field id="bflx_iwm" long_name="internal wave-induced buoyancy flux" standard_name="buoyancy_flux_due_to_internal_waves" unit="W/kg" />387 <field id="pcmap_iwm" long_name="power consumed by wave-driven mixing" standard_name="vertically_integrated_power_consumption_by_wave_driven_mixing" unit="W/m2" grid_ref="grid_W_2D" />388 <field id="emix_iwm" long_name="power density available for mixing" standard_name="power_available_for_mixing_from_breaking_internal_waves" unit="W/kg" />389 390 <!-- variables available with diaar5 -->391 <field id="w_masstr" long_name="vertical mass transport" standard_name="upward_ocean_mass_transport" unit="kg/s" />392 <field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" />393 394 </field_group>395 396 651 <!-- F grid --> 652 <!-- AGRIF sponge --> 653 <field id="agrif_spf" long_name=" AGRIF f-sponge coefficient" unit=" " /> 397 654 <!-- f-eddy viscosity coefficients (ldfdyn) --> 398 655 <field id="ahmf_2d" long_name=" surface f-eddy viscosity coefficient" unit="m2/s or m4/s" /> 399 <field id="ahmf_3d" long_name=" 3D f-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 400 401 <field_group id="scalar" grid_ref="grid_T_2D" > 402 <!-- Need to have a salinity reference climatological file : sali_ref_clim_monthly --> 403 <field id="voltot" long_name="global total volume" standard_name="sea_water_volume" unit="m3" /> 404 <field id="sshtot" long_name="global mean ssh" standard_name="global_average_sea_level_change" unit="m" /> 405 <field id="sshsteric" long_name="global mean ssh steric" standard_name="global_average_steric_sea_level_change" unit="m" /> 406 <field id="sshthster" long_name="global mean ssh thermosteric" standard_name="global_average_thermosteric_sea_level_change" unit="m" /> 407 <field id="masstot" long_name="global total mass" standard_name="sea_water_mass" unit="kg" /> 408 <field id="temptot" long_name="global mean temperature" standard_name="sea_water_potential_temperature" unit="degC" /> 409 <field id="saltot" long_name="global mean salinity" standard_name="sea_water_salinity" unit="1e-3" /> 410 <field id="fram_trans" long_name="Sea Ice Mass Transport Through Fram Strait" standard_name="sea_ice_transport_across_line" unit="kg/s" /> 411 412 <!-- available with ln_diahsb --> 413 <field id="bgtemper" long_name="drift in global mean temperature wrt timestep 1" standard_name="change_over_time_in_sea_water_potential_temperature" unit="degC" /> 414 <field id="bgsaline" long_name="drift in global mean salinity wrt timestep 1" standard_name="change_over_time_in_sea_water_practical_salinity" unit="1e-3" /> 415 <field id="bgheatco" long_name="drift in global mean heat content wrt timestep 1" unit="1.e20J" /> 416 <field id="bgheatfx" long_name="drift in global mean heat flux wrt timestep 1" unit="W/m2" /> 417 <field id="bgsaltco" long_name="drift in global mean salt content wrt timestep 1" unit="1e-3*km3" /> 418 <field id="bgvolssh" long_name="drift in global mean ssh volume wrt timestep 1" unit="km3" /> 419 <field id="bgvole3t" long_name="drift in global mean volume variation (e3t) wrt timestep 1" unit="km3" /> 420 <field id="bgfrcvol" long_name="global mean volume from forcing" unit="km3" /> 421 <field id="bgfrctem" long_name="global mean heat content from forcing" unit="1.e20J" /> 422 <field id="bgfrchfx" long_name="global mean heat flux from forcing" unit="W/m2" /> 423 <field id="bgfrcsal" long_name="global mean salt content from forcing" unit="1e-3*km3" /> 424 <field id="bgmistem" long_name="global mean temperature error due to free surface (linssh true)" unit="degC" /> 425 <field id="bgmissal" long_name="global mean salinity error due to free surface (linssh true)" unit="1e-3" /> 656 <field id="ahmf_3d" long_name=" 3D f-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 657 658 <field_group id="scalar" grid_ref="grid_scalar" > 659 <!-- Need to have a salinity reference climatological file : sali_ref_clim_monthly --> 660 <field id="voltot" long_name="global total volume" standard_name="sea_water_volume" unit="m3" /> 661 <field id="sshtot" long_name="global mean ssh" standard_name="global_average_sea_level_change" unit="m" /> 662 <field id="sshsteric" long_name="global mean ssh steric" standard_name="global_average_steric_sea_level_change" unit="m" /> 663 <field id="sshthster" long_name="global mean ssh thermosteric" standard_name="global_average_thermosteric_sea_level_change" unit="m" /> 664 <field id="masstot" long_name="global total mass" standard_name="sea_water_mass" unit="kg" /> 665 <field id="temptot" long_name="global mean temperature" standard_name="sea_water_potential_temperature" unit="degC" /> 666 <field id="saltot" long_name="global mean salinity" standard_name="sea_water_salinity" unit="1e-3" /> 667 <field id="ssttot" long_name="global mean sea surface temperature" standard_name="sea_water_conservative_temperature" unit="degC" /> 668 <!-- EOS = TEOS-10 --> 669 <field id="temptot_pot" long_name="global mean potential temperature" standard_name="sea_water_potential_temperature" unit="degC" /> 670 671 <!-- available with ln_diahsb --> 672 <field id="bgtemper" long_name="drift in global mean temperature wrt timestep 1" standard_name="change_over_time_in_sea_water_potential_temperature" unit="degC" /> 673 <field id="bgsaline" long_name="drift in global mean salinity wrt timestep 1" standard_name="change_over_time_in_sea_water_practical_salinity" unit="1e-3" /> 674 <field id="bgheatco" long_name="drift in global mean heat content wrt timestep 1" unit="1.e20J" /> 675 <field id="bgheatfx" long_name="drift in global mean heat flux wrt timestep 1" unit="W/m2" /> 676 <field id="bgsaltco" long_name="drift in global mean salt content wrt timestep 1" unit="1e-3*km3" /> 677 <field id="bgvolssh" long_name="drift in global mean ssh volume wrt timestep 1" unit="km3" /> 678 <field id="bgvole3t" long_name="drift in global mean volume variation (e3t) wrt timestep 1" unit="km3" /> 679 <field id="bgfrcvol" long_name="global mean volume from forcing" unit="km3" /> 680 <field id="bgfrctem" long_name="global mean heat content from forcing" unit="1.e20J" /> 681 <field id="bgfrchfx" long_name="global mean heat flux from forcing" unit="W/m2" /> 682 <field id="bgfrcsal" long_name="global mean salt content from forcing" unit="1e-3*km3" /> 683 <field id="bgmistem" long_name="global mean temperature error due to free surface (linssh true)" unit="degC" /> 684 <field id="bgmissal" long_name="global mean salinity error due to free surface (linssh true)" unit="1e-3" /> 426 685 </field_group> 427 428 <!-- variables available with key_float --> 686 687 688 <!-- transects --> 689 <field_group id="oce_straits"> 690 <field id="uoce_e3u_ave" long_name="Monthly average of u*e3u" field_ref="uoce_e3u" freq_op="1mo" freq_offset="_reset_" > @uoce_e3u </field> 691 <field id="uoce_e3u_ave_vsum" long_name="Vertical sum of u*e3u" field_ref="uoce_e3u_ave" grid_ref="grid_U_vsum" /> 692 <field id="uocetr_vsum_section" long_name="Total 2D transport in i-direction" field_ref="uoce_e3u_ave_vsum" grid_ref="grid_U_scalar" detect_missing_value="true"> this * e2u </field> 693 <field id="uocetr_strait" long_name="Total transport across lines in i-direction" field_ref="uocetr_vsum_section" grid_ref="grid_U_4strait" /> 694 <field id="u_masstr_strait" long_name="Sea water transport across line in i-direction" field_ref="uocetr_strait" grid_ref="grid_U_4strait_hsum" unit="kg/s"> this * maskMFO_u * $rau0 </field> 695 696 <field id="voce_e3v_ave" long_name="Monthly average of v*e3v" field_ref="voce_e3v" freq_op="1mo" freq_offset="_reset_" > @voce_e3v </field> 697 <field id="voce_e3v_ave_vsum" long_name="Vertical sum of v*e3v" field_ref="voce_e3v_ave" grid_ref="grid_V_vsum" /> 698 <field id="vocetr_vsum_section" long_name="Total 2D transport of in j-direction" field_ref="voce_e3v_ave_vsum" grid_ref="grid_V_scalar" detect_missing_value="true"> this * e1v </field> 699 <field id="vocetr_strait" long_name="Total transport across lines in j-direction" field_ref="vocetr_vsum_section" grid_ref="grid_V_4strait" /> 700 <field id="v_masstr_strait" long_name="Sea water transport across line in j-direction" field_ref="vocetr_strait" grid_ref="grid_V_4strait_hsum" unit="kg/s"> this * maskMFO_v * $rau0 </field> 701 702 <field id="masstr_strait" long_name="Sea water transport across line" grid_ref="grid_4strait" > u_masstr_strait + v_masstr_strait </field> 703 </field_group> 704 705 706 <!-- variables available with ln_floats --> 429 707 430 708 <field_group id="floatvar" grid_ref="grid_T_nfloat" operation="instant" > 431 432 433 434 435 436 437 709 <field id="traj_lon" long_name="floats longitude" unit="degrees_east" /> 710 <field id="traj_lat" long_name="floats latitude" unit="degrees_north" /> 711 <field id="traj_dep" long_name="floats depth" unit="m" /> 712 <field id="traj_temp" long_name="floats temperature" standard_name="sea_water_potential_temperature" unit="degC" /> 713 <field id="traj_salt" long_name="floats salinity" standard_name="sea_water_practical_salinity" unit="1e-3" /> 714 <field id="traj_dens" long_name="floats in-situ density" standard_name="sea_water_density" unit="kg/m3" /> 715 <field id="traj_group" long_name="floats group" unit="1" /> 438 716 </field_group> 439 717 … … 441 719 442 720 <field_group id="icbvar" domain_ref="grid_T" > 443 <field id="berg_melt" long_name="icb melt rate of icebergs" unit="kg/m2/s" /> 444 <field id="berg_buoy_melt" long_name="icb buoyancy component of iceberg melt rate" unit="kg/m2/s" /> 445 <field id="berg_eros_melt" long_name="icb erosion component of iceberg melt rate" unit="kg/m2/s" /> 446 <field id="berg_conv_melt" long_name="icb convective component of iceberg melt rate" unit="kg/m2/s" /> 447 <field id="berg_virtual_area" long_name="icb virtual coverage by icebergs" unit="m2" /> 448 <field id="bits_src" long_name="icb mass source of bergy bits" unit="kg/m2/s" /> 449 <field id="bits_melt" long_name="icb melt rate of bergy bits" unit="kg/m2/s" /> 450 <field id="bits_mass" long_name="icb bergy bit density field" unit="kg/m2" /> 451 <field id="berg_mass" long_name="icb iceberg density field" unit="kg/m2" /> 452 <field id="calving" long_name="icb calving mass input" unit="kg/s" /> 453 <field id="berg_floating_melt" long_name="icb melt rate of icebergs + bits" unit="kg/m2/s" /> 454 <field id="berg_real_calving" long_name="icb calving into iceberg class" unit="kg/s" axis_ref="icbcla" /> 455 <field id="berg_stored_ice" long_name="icb accumulated ice mass by class" unit="kg" axis_ref="icbcla" /> 721 <field id="berg_melt" long_name="icb melt rate of icebergs" unit="kg/m2/s" /> 722 <field id="berg_melt_hcflx" long_name="icb heat flux to ocean due to melting heat content" unit="J/m2/s" /> 723 <field id="berg_melt_qlat" long_name="icb heat flux to ocean due to melting latent heat" unit="J/m2/s" /> 724 <field id="berg_buoy_melt" long_name="icb buoyancy component of iceberg melt rate" unit="kg/m2/s" /> 725 <field id="berg_eros_melt" long_name="icb erosion component of iceberg melt rate" unit="kg/m2/s" /> 726 <field id="berg_conv_melt" long_name="icb convective component of iceberg melt rate" unit="kg/m2/s" /> 727 <field id="berg_virtual_area" long_name="icb virtual coverage by icebergs" unit="m2" /> 728 <field id="bits_src" long_name="icb mass source of bergy bits" unit="kg/m2/s" /> 729 <field id="bits_melt" long_name="icb melt rate of bergy bits" unit="kg/m2/s" /> 730 <field id="bits_mass" long_name="icb bergy bit density field" unit="kg/m2" /> 731 <field id="berg_mass" long_name="icb iceberg density field" unit="kg/m2" /> 732 <field id="calving" long_name="icb calving mass input" unit="kg/s" /> 733 <field id="berg_floating_melt" long_name="icb melt rate of icebergs + bits" unit="kg/m2/s" /> 734 <field id="berg_real_calving" long_name="icb calving into iceberg class" unit="kg/s" axis_ref="icbcla" /> 735 <field id="berg_stored_ice" long_name="icb accumulated ice mass by class" unit="kg" axis_ref="icbcla" /> 456 736 </field_group> 457 737 458 738 <!-- Poleward transport : ptr --> 459 739 <field_group id="diaptr" > 460 <field id="zomsfglo" long_name="Meridional Stream-Function: Global" unit="Sv" grid_ref="gznl_W_3D" /> 461 <field id="zomsfatl" long_name="Meridional Stream-Function: Atlantic" unit="Sv" grid_ref="gznl_W_3D" /> 462 <field id="zomsfpac" long_name="Meridional Stream-Function: Pacific" unit="Sv" grid_ref="gznl_W_3D" /> 463 <field id="zomsfind" long_name="Meridional Stream-Function: Indian" unit="Sv" grid_ref="gznl_W_3D" /> 464 <field id="zomsfipc" long_name="Meridional Stream-Function: Pacific+Indian" unit="Sv" grid_ref="gznl_W_3D" /> 465 <field id="zotemglo" long_name="Zonal Mean Temperature : Global" unit="degree_C" grid_ref="gznl_T_3D" /> 466 <field id="zotematl" long_name="Zonal Mean Temperature : Atlantic" unit="degree_C" grid_ref="gznl_T_3D" /> 467 <field id="zotempac" long_name="Zonal Mean Temperature : Pacific" unit="degree_C" grid_ref="gznl_T_3D" /> 468 <field id="zotemind" long_name="Zonal Mean Temperature : Indian" unit="degree_C" grid_ref="gznl_T_3D" /> 469 <field id="zotemipc" long_name="Zonal Mean Temperature : Pacific+Indian" unit="degree_C" grid_ref="gznl_T_3D" /> 470 <field id="zosalglo" long_name="Zonal Mean Salinity : Global" unit="0.001" grid_ref="gznl_T_3D" /> 471 <field id="zosalatl" long_name="Zonal Mean Salinity : Atlantic" unit="0.001" grid_ref="gznl_T_3D" /> 472 <field id="zosalpac" long_name="Zonal Mean Salinity : Pacific" unit="0.001" grid_ref="gznl_T_3D" /> 473 <field id="zosalind" long_name="Zonal Mean Salinity : Indian" unit="0.001" grid_ref="gznl_T_3D" /> 474 <field id="zosalipc" long_name="Zonal Mean Salinity : Pacific+Indian" unit="0.001" grid_ref="gznl_T_3D" /> 475 <field id="zosrfglo" long_name="Zonal Mean Surface" unit="m2" grid_ref="gznl_T_3D" /> 476 <field id="zosrfatl" long_name="Zonal Mean Surface : Atlantic" unit="m2" grid_ref="gznl_T_3D" /> 477 <field id="zosrfpac" long_name="Zonal Mean Surface : Pacific" unit="m2" grid_ref="gznl_T_3D" /> 478 <field id="zosrfind" long_name="Zonal Mean Surface : Indian" unit="m2" grid_ref="gznl_T_3D" /> 479 <field id="zosrfipc" long_name="Zonal Mean Surface : Pacific+Indian" unit="m2" grid_ref="gznl_T_3D" /> 480 <field id="sophtadv" long_name="Advective Heat Transport" unit="PW" grid_ref="gznl_T_2D" /> 481 <field id="sophtadv_atl" long_name="Advective Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 482 <field id="sophtadv_pac" long_name="Advective Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 483 <field id="sophtadv_ind" long_name="Advective Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 484 <field id="sophtadv_ipc" long_name="Advective Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 485 <field id="sophtldf" long_name="Diffusive Heat Transport" unit="PW" grid_ref="gznl_T_2D" /> 486 <field id="sophtldf_atl" long_name="Diffusive Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 487 <field id="sophtldf_pac" long_name="Diffusive Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 488 <field id="sophtldf_ind" long_name="Diffusive Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 489 <field id="sophtldf_ipc" long_name="Diffusive Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 490 <field id="sophtove" long_name="Overturning Heat Transport" unit="PW" grid_ref="gznl_T_2D" /> 491 <field id="sophtove_atl" long_name="Overturning Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 492 <field id="sophtove_pac" long_name="Overturning Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 493 <field id="sophtove_ind" long_name="Overturning Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 494 <field id="sophtove_ipc" long_name="Overturning Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 495 <field id="sophtbtr" long_name="Barotropic Heat Transport" unit="PW" grid_ref="gznl_T_2D" /> 496 <field id="sophtbtr_atl" long_name="Barotropic Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 497 <field id="sophtbtr_pac" long_name="Barotropic Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 498 <field id="sophtbtr_ind" long_name="Barotropic Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 499 <field id="sophtbtr_ipc" long_name="Barotropic Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 500 <field id="sophteiv" long_name="Heat Transport from mesoscale eddy advection" unit="PW" grid_ref="gznl_T_2D" /> 501 <field id="sophteiv_atl" long_name="Heat Transport from mesoscale eddy advection: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 502 <field id="sophteiv_pac" long_name="Heat Transport from mesoscale eddy advection: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 503 <field id="sophteiv_ind" long_name="Heat Transport from mesoscale eddy advection: Indian" unit="PW" grid_ref="gznl_T_2D" /> 504 <field id="sophteiv_ipc" long_name="Heat Transport from mesoscale eddy advection: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 505 <field id="sopstadv" long_name="Advective Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" /> 506 <field id="sopstadv_atl" long_name="Advective Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 507 <field id="sopstadv_pac" long_name="Advective Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 508 <field id="sopstadv_ind" long_name="Advective Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 509 <field id="sopstadv_ipc" long_name="Advective Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 510 <field id="sopstove" long_name="Overturning Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" /> 511 <field id="sopstove_atl" long_name="Overturning Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 512 <field id="sopstove_pac" long_name="Overturning Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 513 <field id="sopstove_ind" long_name="Overturning Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 514 <field id="sopstove_ipc" long_name="Overturning Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 515 <field id="sopstbtr" long_name="Barotropic Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" /> 516 <field id="sopstbtr_atl" long_name="Barotropic Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 517 <field id="sopstbtr_pac" long_name="Barotropic Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 518 <field id="sopstbtr_ind" long_name="Barotropic Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 519 <field id="sopstbtr_ipc" long_name="Barotropic Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 520 <field id="sopstldf" long_name="Diffusive Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" /> 521 <field id="sopstldf_atl" long_name="Diffusive Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 522 <field id="sopstldf_pac" long_name="Diffusive Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 523 <field id="sopstldf_ind" long_name="Diffusive Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 524 <field id="sopstldf_ipc" long_name="Diffusive Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 525 <field id="sopsteiv" long_name="Salt Transport from mesoscale eddy advection" unit="Giga g/s" grid_ref="gznl_T_2D" /> 526 <field id="sopsteiv_atl" long_name="Salt Transport from mesoscale eddy advection: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 527 <field id="sopsteiv_pac" long_name="Salt Transport from mesoscale eddy advection: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 528 <field id="sopsteiv_ind" long_name="Salt Transport from mesoscale eddy advection: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 529 <field id="sopsteiv_ipc" long_name="Salt Transport from mesoscale eddy advection: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 530 </field_group> 740 <field id="zomsf" long_name="Overturning Stream-Function : All basins" unit="Sv" grid_ref="grid_znl_W_3D" /> 741 <field id="zotem" long_name="Zonal Mean Temperature : All basins" unit="degree_C" grid_ref="grid_znl_T_3D" /> 742 <field id="zosal" long_name="Zonal Mean Salinity : All basins" unit="0.001" grid_ref="grid_znl_T_3D" /> 743 <field id="zosrf" long_name="Zonal Mean Surface : All basins" unit="m2" grid_ref="grid_znl_T_3D" /> 744 <field id="sophtove" long_name="Overturning Heat Transport: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 745 <field id="sopstove" long_name="Overturning Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 746 <field id="sophtbtr" long_name="Barotropic Heat Transport: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 747 <field id="sopstbtr" long_name="Barotropic Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 748 <field id="sophtadv" long_name="Advective Heat Transport: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 749 <field id="sopstadv" long_name="Advective Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 750 <field id="sophtldf" long_name="Diffusive Heat Transport: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 751 <field id="sopstldf" long_name="Diffusive Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 752 <field id="sophtvtr" long_name="Heat Transport : All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 753 <field id="sopstvtr" long_name="Salt Transport : All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 754 <field id="sophteiv" long_name="Heat Transport from mesoscale eddy advection: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 755 <field id="sopsteiv" long_name="Salt Transport from mesoscale eddy advection : All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 756 <field id="sopstadv" long_name="Advective Salt Transport" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 757 <field id="sophtgyre" long_name="Overturning heat transport due to gyre" field_ref="sophtove" unit="W" grid_ref="grid_znl_T_2D" > sophtvtr - sophtove </field> 758 <field id="sopstgyre" long_name="Overturning salt transport due to gyre" field_ref="sopstove" unit="kg/s" grid_ref="grid_znl_T_2D" > sophtvtr - sopstove </field> 759 </field_group> 760 761 <field_group id="constant_fields" grid_ref="grid_T_2D" operation="once" > 762 <field id="bathy" long_name="Sea floor depth below geoid" standard_name="sea_floor_depth_below_geoid" unit="m"/> 763 <field id="areacello" long_name="Horizontal area of ocean grid cells" standard_name="cell_area" unit="m2" /> 764 <field id="hfgeou" long_name="Upward geothermal heat flux at sea floor" standard_name="upward_geothermal_heat_flux_at_sea_floor" unit="W/m2"/> 765 <field id="basins" long_name="ocean tracer region masks" standard_name="ocean_tracer_region_masks" unit="none" grid_ref="grid_basin" /> 766 </field_group> 767 531 768 532 769 <!-- … … 539 776 <!-- Asselin trends calculated on odd time steps--> 540 777 <field_group id="trendT_odd" grid_ref="grid_T_3D"> 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 <field id="ttrd_sad" long_name="temperature-trend: surface adv. (linssh true)"unit="degC/s" grid_ref="grid_T_2D" />564 <field id="strd_sad" long_name="salinity -trend: surface adv. (linssh true)"unit="1e-3/s" grid_ref="grid_T_2D" />565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 <field id="ttrd_qns_e3t" unit="degC/s * m" > ttrd_qns * e3t_surf</field>622 <field id="strd_cdt_e3t" unit="degC/s * m" > strd_cdt * e3t_surf</field>623 624 625 626 627 <field id="ttrd_totad_li" long_name="layer integrated heat-trend : total advection"unit="W/m^2" > ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field>628 <field id="strd_totad_li" long_name="layer integrated salt -trend : total advection"unit="kg/(m^2 s)" > strd_totad_e3t * 1026.0 * 0.001 </field>629 <field id="ttrd_evd_li" long_name="layer integrated heat-trend : EVD convection"unit="W/m^2" > ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field>630 <field id="strd_evd_li" long_name="layer integrated salt -trend : EVD convection"unit="kg/(m^2 s)" > strd_evd_e3t * 1026.0 * 0.001 </field>631 <field id="ttrd_iso_li" long_name="layer integrated heat-trend : isopycnal diffusion"unit="W/m^2" > ttrd_iso_e3t * 1026.0 * 3991.86795711963 </field>632 <field id="strd_iso_li" long_name="layer integrated salt -trend : isopycnal diffusion"unit="kg/(m^2 s)" > strd_iso_e3t * 1026.0 * 0.001 </field>633 <field id="ttrd_zdfp_li" long_name="layer integrated heat-trend : pure vert. diffusion"unit="W/m^2" > ttrd_zdfp_e3t * 1026.0 * 3991.86795711963 </field>634 <field id="strd_zdfp_li" long_name="layer integrated salt -trend : pure vert. diffusion"unit="kg/(m^2 s)" > strd_zdfp_e3t * 1026.0 * 0.001 </field>635 <field id="ttrd_qns_li" long_name="layer integrated heat-trend : non-solar flux + runoff"unit="W/m^2" grid_ref="grid_T_2D"> ttrd_qns_e3t * 1026.0 * 3991.86795711963 </field>636 <field id="ttrd_qsr_li" long_name="layer integrated heat-trend : solar flux"unit="W/m^2" grid_ref="grid_T_3D"> ttrd_qsr_e3t * 1026.0 * 3991.86795711963 </field>637 <field id="ttrd_bbl_li" long_name="layer integrated heat-trend: bottom boundary layer "unit="W/m^2" > ttrd_bbl_e3t * 1026.0 * 3991.86795711963 </field>638 <field id="strd_bbl_li" long_name="layer integrated salt -trend: bottom boundary layer "unit="kg/(m^2 s)" > strd_bbl_e3t * 1026.0 * 0.001 </field>639 <field id="ttrd_evd_li" long_name="layer integrated heat -trend: evd convection "unit="W/m^2" >ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field>640 <field id="strd_evd_li" long_name="layer integrated salt -trend: evd convection "unit="kg/(m^2 s)" > strd_evd_e3t * 1026.0 * 0.001 </field>641 642 643 644 645 646 647 648 649 650 651 652 653 <field id="strd_tot_li" long_name="layer integrated salt-trend: total model trend :" unit="kg/(m^2 s)" > strd_tot_e3t * 1026.0 * 0.001 </field>654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 778 <field id="ttrd_atf" long_name="temperature-trend: asselin time filter" unit="degree_C/s" /> 779 <field id="strd_atf" long_name="salinity -trend: asselin time filter" unit="0.001/s" /> 780 <!-- Thickness weighted versions: --> 781 <field id="ttrd_atf_e3t" unit="degC/s * m" > ttrd_atf * e3t </field> 782 <field id="strd_atf_e3t" unit="1e-3/s * m" > strd_atf * e3t </field> 783 <!-- OMIP layer-integrated trends --> 784 <field id="ttrd_atf_li" long_name="layer integrated heat-trend: asselin time filter " unit="W/m^2" > ttrd_atf_e3t * 1026.0 * 3991.86795711963 </field> 785 <field id="strd_atf_li" long_name="layer integrated salt -trend: asselin time filter " unit="kg/(m^2 s)" > strd_atf_e3t * 1026.0 * 0.001 </field> 786 </field_group> 787 788 <!-- Other trends calculated on even time steps--> 789 <field_group id="trendT_even" grid_ref="grid_T_3D"> 790 <field id="ttrd_xad" long_name="temperature-trend: i-advection" unit="degC/s" /> 791 <field id="strd_xad" long_name="salinity -trend: i-advection" unit="1e-3/s" /> 792 <field id="ttrd_yad" long_name="temperature-trend: j-advection" unit="degC/s" /> 793 <field id="strd_yad" long_name="salinity -trend: j-advection" unit="1e-3/s" /> 794 <field id="ttrd_zad" long_name="temperature-trend: k-advection" unit="degC/s" /> 795 <field id="strd_zad" long_name="salinity -trend: k-advection" unit="1e-3/s" /> 796 <field id="ttrd_ad" long_name="temperature-trend: advection" standard_name="tendency_of_sea_water_temperature_due_to_advection" unit="degC/s" > sqrt( ttrd_xad^2 + ttrd_yad^2 + ttrd_zad^2 ) </field> 797 <field id="strd_ad" long_name="salinity -trend: advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="1e-3/s" > sqrt( strd_xad^2 + strd_yad^2 + strd_zad^2 ) </field> 798 <field id="ttrd_totad" long_name="temperature-trend: total advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="degC/s" /> 799 <field id="strd_totad" long_name="salinity -trend: total advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="1e-3/s" /> 800 <field id="ttrd_sad" long_name="temperature-trend: surface adv. (linssh true)" unit="degC/s" grid_ref="grid_T_2D" /> 801 <field id="strd_sad" long_name="salinity -trend: surface adv. (linssh true)" unit="1e-3/s" grid_ref="grid_T_2D" /> 802 <field id="ttrd_ldf" long_name="temperature-trend: lateral diffusion" standard_name="tendency_of_sea_water_temperature_due_to_horizontal_mixing" unit="degC/s" /> 803 <field id="strd_ldf" long_name="salinity -trend: lateral diffusion" standard_name="tendency_of_sea_water_salinity_due_to_horizontal_mixing" unit="1e-3/s" /> 804 <field id="ttrd_zdf" long_name="temperature-trend: vertical diffusion" standard_name="tendency_of_sea_water_temperature_due_to_vertical_mixing" unit="degC/s" /> 805 <field id="strd_zdf" long_name="salinity -trend: vertical diffusion" standard_name="tendency_of_sea_water_salinity_due_to_vertical_mixing" unit="1e-3/s" /> 806 <field id="ttrd_evd" long_name="temperature-trend: EVD convection" unit="degC/s" /> 807 <field id="strd_evd" long_name="salinity -trend: EVD convection" unit="1e-3/s" /> 808 809 <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 810 <field id="ttrd_iso" long_name="temperature-trend: isopycnal diffusion" unit="degC/s" > ttrd_ldf + ttrd_zdf - ttrd_zdfp </field> 811 <field id="strd_iso" long_name="salinity -trend: isopycnal diffusion" unit="1e-3/s" > strd_ldf + strd_zdf - strd_zdfp </field> 812 <field id="ttrd_zdfp" long_name="temperature-trend: pure vert. diffusion" unit="degC/s" /> 813 <field id="strd_zdfp" long_name="salinity -trend: pure vert. diffusion" unit="1e-3/s" /> 814 815 <!-- --> 816 <field id="ttrd_dmp" long_name="temperature-trend: interior restoring" unit="degC/s" /> 817 <field id="strd_dmp" long_name="salinity -trend: interior restoring" unit="1e-3/s" /> 818 <field id="ttrd_bbl" long_name="temperature-trend: bottom boundary layer" unit="degC/s" /> 819 <field id="strd_bbl" long_name="salinity -trend: bottom boundary layer" unit="1e-3/s" /> 820 <field id="ttrd_npc" long_name="temperature-trend: non-penetrative conv." unit="degC/s" /> 821 <field id="strd_npc" long_name="salinity -trend: non-penetrative conv." unit="1e-3/s" /> 822 <field id="ttrd_qns" long_name="temperature-trend: non-solar flux + runoff" unit="degC/s" grid_ref="grid_T_2D" /> 823 <field id="strd_cdt" long_name="salinity -trend: C/D term + runoff" unit="degC/s" grid_ref="grid_T_2D" /> 824 <field id="ttrd_qsr" long_name="temperature-trend: solar penetr. heating" unit="degC/s" /> 825 <field id="ttrd_bbc" long_name="temperature-trend: geothermal heating" unit="degC/s" /> 826 827 <!-- Thickness weighted versions: --> 828 <field id="ttrd_xad_e3t" unit="degC/s * m" > ttrd_xad * e3t </field> 829 <field id="strd_xad_e3t" unit="1e-3/s * m" > strd_xad * e3t </field> 830 <field id="ttrd_yad_e3t" unit="degC/s * m" > ttrd_yad * e3t </field> 831 <field id="strd_yad_e3t" unit="1e-3/s * m" > strd_yad * e3t </field> 832 <field id="ttrd_zad_e3t" unit="degC/s * m" > ttrd_zad * e3t </field> 833 <field id="strd_zad_e3t" unit="1e-3/s * m" > strd_zad * e3t </field> 834 <field id="ttrd_ad_e3t" unit="degC/s * m" > ttrd_ad * e3t </field> 835 <field id="strd_ad_e3t" unit="1e-3/s * m" > strd_ad * e3t </field> 836 <field id="ttrd_totad_e3t" unit="degC/s * m" > ttrd_totad * e3t </field> 837 <field id="strd_totad_e3t" unit="1e-3/s * m" > strd_totad * e3t </field> 838 <field id="ttrd_ldf_e3t" unit="degC/s * m" > ttrd_ldf * e3t </field> 839 <field id="strd_ldf_e3t" unit="1e-3/s * m" > strd_ldf * e3t </field> 840 <field id="ttrd_zdf_e3t" unit="degC/s * m" > ttrd_zdf * e3t </field> 841 <field id="strd_zdf_e3t" unit="1e-3/s * m" > strd_zdf * e3t </field> 842 <field id="ttrd_evd_e3t" unit="degC/s * m" > ttrd_evd * e3t </field> 843 <field id="strd_evd_e3t" unit="1e-3/s * m" > strd_evd * e3t </field> 844 845 <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 846 <field id="ttrd_iso_e3t" unit="degC/s * m" > ttrd_iso * e3t </field> 847 <field id="strd_iso_e3t" unit="1e-3/s * m" > strd_iso * e3t </field> 848 <field id="ttrd_zdfp_e3t" unit="degC/s * m" > ttrd_zdfp * e3t </field> 849 <field id="strd_zdfp_e3t" unit="1e-3/s * m" > strd_zdfp * e3t </field> 850 851 <!-- --> 852 <field id="ttrd_dmp_e3t" unit="degC/s * m" > ttrd_dmp * e3t </field> 853 <field id="strd_dmp_e3t" unit="1e-3/s * m" > strd_dmp * e3t </field> 854 <field id="ttrd_bbl_e3t" unit="degC/s * m" > ttrd_bbl * e3t </field> 855 <field id="strd_bbl_e3t" unit="1e-3/s * m" > strd_bbl * e3t </field> 856 <field id="ttrd_npc_e3t" unit="degC/s * m" > ttrd_npc * e3t </field> 857 <field id="strd_npc_e3t" unit="1e-3/s * m" > strd_npc * e3t </field> 858 <field id="ttrd_qns_e3t" unit="degC/s * m" > ttrd_qns * e3ts </field> 859 <field id="strd_cdt_e3t" unit="degC/s * m" > strd_cdt * e3ts </field> 860 <field id="ttrd_qsr_e3t" unit="degC/s * m" > ttrd_qsr * e3t </field> 861 <field id="ttrd_bbc_e3t" unit="degC/s * m" > ttrd_bbc * e3t </field> 862 863 <!-- OMIP layer-integrated trends --> 864 <field id="ttrd_totad_li" long_name="layer integrated heat-trend: total advection" unit="W/m^2" > ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field> 865 <field id="strd_totad_li" long_name="layer integrated salt-trend: total advection" unit="kg/(m^2 s)" > strd_totad_e3t * 1026.0 * 0.001 </field> 866 <field id="ttrd_evd_li" long_name="layer integrated heat-trend: EVD convection" unit="W/m^2" > ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field> 867 <field id="strd_evd_li" long_name="layer integrated salt-trend: EVD convection" unit="kg/(m^2 s)" > strd_evd_e3t * 1026.0 * 0.001 </field> 868 <field id="ttrd_iso_li" long_name="layer integrated heat-trend: isopycnal diffusion" unit="W/m^2" > ttrd_iso_e3t * 1026.0 * 3991.86795711963 </field> 869 <field id="strd_iso_li" long_name="layer integrated salt-trend: isopycnal diffusion" unit="kg/(m^2 s)" > strd_iso_e3t * 1026.0 * 0.001 </field> 870 <field id="ttrd_zdfp_li" long_name="layer integrated heat-trend: pure vert. diffusion" unit="W/m^2" > ttrd_zdfp_e3t * 1026.0 * 3991.86795711963 </field> 871 <field id="strd_zdfp_li" long_name="layer integrated salt-trend: pure vert. diffusion" unit="kg/(m^2 s)" > strd_zdfp_e3t * 1026.0 * 0.001 </field> 872 <field id="ttrd_qns_li" long_name="layer integrated heat-trend: non-solar flux + runoff" unit="W/m^2" grid_ref="grid_T_2D"> ttrd_qns_e3t * 1026.0 * 3991.86795711963 </field> 873 <field id="ttrd_qsr_li" long_name="layer integrated heat-trend: solar flux" unit="W/m^2" grid_ref="grid_T_3D"> ttrd_qsr_e3t * 1026.0 * 3991.86795711963 </field> 874 <field id="ttrd_bbl_li" long_name="layer integrated heat-trend: bottom boundary layer " unit="W/m^2" > ttrd_bbl_e3t * 1026.0 * 3991.86795711963 </field> 875 <field id="strd_bbl_li" long_name="layer integrated salt-trend: bottom boundary layer " unit="kg/(m^2 s)" > strd_bbl_e3t * 1026.0 * 0.001 </field> 876 <field id="ttrd_evd_li" long_name="layer integrated heat-trend: evd convection " unit="W/m^2" >ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field> 877 <field id="strd_evd_li" long_name="layer integrated salt-trend: evd convection " unit="kg/(m^2 s)" > strd_evd_e3t * 1026.0 * 0.001 </field> 878 879 </field_group> 880 881 <!-- Total trends calculated every time step--> 882 <field_group id="trendT" grid_ref="grid_T_3D"> 883 <field id="ttrd_tot" long_name="temperature-trend: total model trend" unit="degC/s" /> 884 <field id="strd_tot" long_name="salinity -trend: total model trend" unit="1e-3/s" /> 885 <!-- Thickness weighted versions: --> 886 <field id="ttrd_tot_e3t" unit="degC/s * m" > ttrd_tot * e3t </field> 887 <field id="strd_tot_e3t" unit="1e-3/s * m" > strd_tot * e3t </field> 888 <!-- OMIP layer-integrated total trends --> 889 <field id="ttrd_tot_li" long_name="layer integrated heat-trend: total model trend :" unit="W/m^2" > ttrd_tot_e3t * 1026.0 * 3991.86795711963 </field> 890 <field id="strd_tot_li" long_name="layer integrated salt-trend: total model trend :" unit="kg/(m^2 s)" > strd_tot_e3t * 1026.0 * 0.001 </field> 891 892 <!-- **** these trends have not been apportioned to all/even/odd ts yet **** --> 893 <!-- variables available with ln_KE_trd --> 894 <field id="ketrd_hpg" long_name="ke-trend: hydrostatic pressure gradient" unit="W/s^3" /> 895 <field id="ketrd_spg" long_name="ke-trend: surface pressure gradient" unit="W/s^3" /> 896 <field id="ketrd_spgexp" long_name="ke-trend: surface pressure gradient (explicit)" unit="W/s^3" /> 897 <field id="ketrd_spgflt" long_name="ke-trend: surface pressure gradient (filter)" unit="W/s^3" /> 898 <field id="ssh_flt" long_name="filtered contribution to ssh (dynspg_flt)" unit="m" grid_ref="grid_T_2D" /> 899 <field id="w0" long_name="surface vertical velocity" unit="m/s" grid_ref="grid_T_2D" /> 900 <field id="pw0_exp" long_name="surface pressure flux due to ssh" unit="W/s^2" grid_ref="grid_T_2D" /> 901 <field id="pw0_flt" long_name="surface pressure flux due to filtered ssh" unit="W/s^2" grid_ref="grid_T_2D" /> 902 <field id="ketrd_keg" long_name="ke-trend: KE gradient or hor. adv." unit="W/s^3" /> 903 <field id="ketrd_rvo" long_name="ke-trend: relative vorticity or metric term" unit="W/s^3" /> 904 <field id="ketrd_pvo" long_name="ke-trend: planetary vorticity" unit="W/s^3" /> 905 <field id="ketrd_zad" long_name="ke-trend: vertical advection" unit="W/s^3" /> 906 <field id="ketrd_udx" long_name="ke-trend: U.dx[U]" unit="W/s^3" /> 907 <field id="ketrd_ldf" long_name="ke-trend: lateral diffusion" unit="W/s^3" /> 908 <field id="ketrd_zdf" long_name="ke-trend: vertical diffusion" unit="W/s^3" /> 909 <field id="ketrd_tau" long_name="ke-trend: wind stress " unit="W/s^3" grid_ref="grid_T_2D" /> 910 <field id="ketrd_bfr" long_name="ke-trend: bottom friction (explicit)" unit="W/s^3" /> 911 <field id="ketrd_bfri" long_name="ke-trend: bottom friction (implicit)" unit="W/s^3" /> 912 <field id="ketrd_atf" long_name="ke-trend: asselin time filter trend" unit="W/s^3" /> 913 <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)" unit="W/s^3" /> 914 <field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" /> 915 916 <!-- variables available when explicit lateral mixing is used (ln_dynldf_OFF=F) --> 917 <field id="dispkexyfo" long_name="KE-trend: lateral mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 918 <field id="dispkevfo" long_name="KE-trend: vertical mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 919 <!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T --> 920 <field id="eketrd_eiv" long_name="EKE-trend due to parameterized eddy advection" standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection" unit="W/m^2" grid_ref="grid_T_2D" /> 921 922 <!-- variables available with ln_PE_trd --> 923 <field id="petrd_xad" long_name="pe-trend: i-advection" unit="W/m^3" /> 924 <field id="petrd_yad" long_name="pe-trend: j-advection" unit="W/m^3" /> 925 <field id="petrd_zad" long_name="pe-trend: k-advection" unit="W/m^3" /> 926 <field id="petrd_sad" long_name="pe-trend: surface adv. (linssh true)" unit="W/m^3" grid_ref="grid_T_2D" /> 927 <field id="petrd_ldf" long_name="pe-trend: lateral diffusion" unit="W/m^3" /> 928 <field id="petrd_zdf" long_name="pe-trend: vertical diffusion" unit="W/m^3" /> 929 <field id="petrd_zdfp" long_name="pe-trend: pure vert. diffusion" unit="W/m^3" /> 930 <field id="petrd_dmp" long_name="pe-trend: interior restoring" unit="W/m^3" /> 931 <field id="petrd_bbl" long_name="pe-trend: bottom boundary layer" unit="W/m^3" /> 932 <field id="petrd_npc" long_name="pe-trend: non-penetrative conv." unit="W/m^3" /> 933 <field id="petrd_nsr" long_name="pe-trend: surface forcing + runoff" unit="W/m^3" /> 934 <field id="petrd_qsr" long_name="pe-trend: solar penetr. heating" unit="W/m^3" /> 935 <field id="petrd_bbc" long_name="pe-trend: geothermal heating" unit="W/m^3" /> 936 <field id="petrd_atf" long_name="pe-trend: asselin time filter" unit="W/m^3" /> 937 <field id="PEanom" long_name="potential energy anomaly" unit="1" /> 938 <field id="alphaPE" long_name="partial deriv. of PEanom wrt T" unit="degC-1" /> 939 <field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3" /> 940 </field_group> 941 942 <field_group id="trendU" grid_ref="grid_U_3D"> 943 <!-- variables available with ln_dyn_trd --> 944 <field id="utrd_hpg" long_name="i-trend: hydrostatic pressure gradient" unit="m/s^2" /> 945 <field id="utrd_spg" long_name="i-trend: surface pressure gradient" unit="m/s^2" /> 946 <field id="utrd_spgexp" long_name="i-trend: surface pressure gradient (explicit)" unit="m/s^2" /> 947 <field id="utrd_spgflt" long_name="i-trend: surface pressure gradient (filtered)" unit="m/s^2" /> 948 <field id="utrd_keg" long_name="i-trend: KE gradient or hor. adv." unit="m/s^2" /> 949 <field id="utrd_rvo" long_name="i-trend: relative vorticity or metric term" unit="m/s^2" /> 950 <field id="utrd_pvo" long_name="i-trend: planetary vorticity" unit="m/s^2" /> 951 <field id="utrd_zad" long_name="i-trend: vertical advection" unit="m/s^2" /> 952 <field id="utrd_udx" long_name="i-trend: U.dx[U]" unit="m/s^2" /> 953 <field id="utrd_ldf" long_name="i-trend: lateral diffusion" unit="m/s^2" /> 954 <field id="utrd_zdf" long_name="i-trend: vertical diffusion" unit="m/s^2" /> 955 <field id="utrd_tau" long_name="i-trend: wind stress " unit="m/s^2" grid_ref="grid_U_2D" /> 956 <field id="utrd_bfr" long_name="i-trend: bottom friction (explicit)" unit="m/s^2" /> 957 <field id="utrd_bfri" long_name="i-trend: bottom friction (implicit)" unit="m/s^2" /> 958 <field id="utrd_tot" long_name="i-trend: total momentum trend before atf" unit="m/s^2" /> 959 <field id="utrd_atf" long_name="i-trend: asselin time filter trend" unit="m/s^2" /> 960 </field_group> 961 962 <field_group id="trendV" grid_ref="grid_V_3D"> 963 <!-- variables available with ln_dyn_trd --> 964 <field id="vtrd_hpg" long_name="j-trend: hydrostatic pressure gradient" unit="m/s^2" /> 965 <field id="vtrd_spg" long_name="j-trend: surface pressure gradient" unit="m/s^2" /> 966 <field id="vtrd_spgexp" long_name="j-trend: surface pressure gradient (explicit)" unit="m/s^2" /> 967 <field id="vtrd_spgflt" long_name="j-trend: surface pressure gradient (filtered)" unit="m/s^2" /> 968 <field id="vtrd_keg" long_name="j-trend: KE gradient or hor. adv." unit="m/s^2" /> 969 <field id="vtrd_rvo" long_name="j-trend: relative vorticity or metric term" unit="m/s^2" /> 970 <field id="vtrd_pvo" long_name="j-trend: planetary vorticity" unit="m/s^2" /> 971 <field id="vtrd_zad" long_name="j-trend: vertical advection" unit="m/s^2" /> 972 <field id="vtrd_vdy" long_name="i-trend: V.dx[V]" unit="m/s^2" /> 973 <field id="vtrd_ldf" long_name="j-trend: lateral diffusion" unit="m/s^2" /> 974 <field id="vtrd_zdf" long_name="j-trend: vertical diffusion" unit="m/s^2" /> 975 <field id="vtrd_tau" long_name="j-trend: wind stress " unit="m/s^2" grid_ref="grid_V_2D" /> 976 <field id="vtrd_bfr" long_name="j-trend: bottom friction (explicit)" unit="m/s^2" /> 977 <field id="vtrd_bfri" long_name="j-trend: bottom friction (implicit)" unit="m/s^2" /> 978 <field id="vtrd_tot" long_name="j-trend: total momentum trend before atf" unit="m/s^2" /> 979 <field id="vtrd_atf" long_name="j-trend: asselin time filter trend" unit="m/s^2" /> 980 </field_group> 744 981 745 982 … … 750 987 --> 751 988 752 753 754 755 756 757 758 759 760 761 762 989 <field_group id="TRD" > 990 <field field_ref="ttrd_totad_li" name="opottempadvect" /> 991 <field field_ref="ttrd_iso_li" name="opottemppmdiff" /> 992 <field field_ref="ttrd_zdfp_li" name="opottempdiff" /> 993 <field field_ref="ttrd_evd_li" name="opottempevd" /> 994 <field field_ref="strd_evd_li" name="osaltevd" /> 995 <field field_ref="ttrd_qns_li" name="opottempqns" /> 996 <field field_ref="ttrd_qsr_li" name="rsdoabsorb" operation="accumulate" /> 997 <field field_ref="strd_totad_li" name="osaltadvect" /> 998 <field field_ref="strd_iso_li" name="osaltpmdiff" /> 999 <field field_ref="strd_zdfp_li" name="osaltdiff" /> 763 1000 </field_group> 764 1001 765 1002 <field_group id="mooring" > 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 1003 <field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature" /> 1004 <field field_ref="soce" name="so" long_name="sea_water_salinity" /> 1005 <field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" /> 1006 <field field_ref="voce" name="vo" long_name="sea_water_y_velocity" /> 1007 <field field_ref="woce" name="wo" long_name="sea_water_z_velocity" /> 1008 <field field_ref="avt" name="difvho" long_name="ocean_vertical_heat_diffusivity" /> 1009 <field field_ref="avm" name="difvmo" long_name="ocean_vertical_momentum_diffusivity" /> 1010 1011 <field field_ref="sst" name="tos" long_name="sea_surface_temperature" /> 1012 <field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature" /> 1013 <field field_ref="sstgrad" name="tosgrad" long_name="module_of_sea_surface_temperature_gradient" /> 1014 <field field_ref="sss" name="sos" long_name="sea_surface_salinity" /> 1015 <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" /> 1016 <field field_ref="empmr" name="wfo" long_name="water_flux_into_sea_water" /> 1017 <field field_ref="qsr" name="rsntds" long_name="surface_net_downward_shortwave_flux" /> 1018 <field field_ref="qt" name="tohfls" long_name="surface_net_downward_total_heat_flux" /> 1019 <field field_ref="taum" /> 1020 <field field_ref="20d" /> 1021 <field field_ref="mldkz5" /> 1022 <field field_ref="mldr10_1" /> 1023 <field field_ref="mldr10_3" /> 1024 <field field_ref="mldr0_1" /> 1025 <field field_ref="mldr0_3" /> 1026 <field field_ref="mld_dt02" /> 1027 <field field_ref="topthdep" /> 1028 <field field_ref="pycndep" /> 1029 <field field_ref="tinv" /> 1030 <field field_ref="depti" /> 1031 <field field_ref="BLT" name="blt" long_name="barrier_layer_thickness" /> 1032 <field field_ref="utau" name="tauuo" long_name="surface_downward_x_stress" /> 1033 <field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" /> 797 1034 </field_group> 798 1035 799 1036 <field_group id="groupT" > 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 1037 <field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature" /> 1038 <field field_ref="soce" name="so" long_name="sea_water_salinity" /> 1039 <field field_ref="sst" name="tos" long_name="sea_surface_temperature" /> 1040 <field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature" /> 1041 <field field_ref="sss" name="sos" long_name="sea_surface_salinity" /> 1042 <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" /> 1043 <field field_ref="empmr" name="wfo" long_name="water_flux_into_sea_water" /> 1044 <field field_ref="qsr" name="rsntds" long_name="surface_net_downward_shortwave_flux" /> 1045 <field field_ref="qt" name="tohfls" long_name="surface_net_downward_total_heat_flux" /> 1046 <field field_ref="taum" /> 1047 <field field_ref="20d" /> 1048 <field field_ref="mldkz5" /> 1049 <field field_ref="mldr10_1" /> 1050 <field field_ref="mldr10_3" /> 1051 <field field_ref="mld_dt02" /> 1052 <field field_ref="topthdep" /> 1053 <field field_ref="pycndep" /> 1054 <field field_ref="tinv" /> 1055 <field field_ref="depti" /> 1056 <field field_ref="BLT" name="blt" long_name="Barrier Layer Thickness" /> 820 1057 </field_group> 821 1058 822 1059 <field_group id="groupU" > 823 <field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" /> 824 <field field_ref="ssu" name="uos" long_name="sea_surface_x_velocity" /> 825 <field field_ref="utau" name="tauuo" long_name="surface_downward_x_stress" /> 1060 <field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" /> 1061 <field field_ref="utau" name="tauuo" long_name="surface_downward_x_stress" /> 826 1062 </field_group> 827 1063 828 1064 <field_group id="groupV" > 829 <field field_ref="voce" name="vo" long_name="sea_water_y_velocity" /> 830 <field field_ref="ssv" name="vos" long_name="sea_surface_y_velocity" /> 831 <field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" /> 1065 <field field_ref="voce" name="vo" long_name="sea_water_y_velocity" /> 1066 <field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" /> 832 1067 </field_group> 833 1068 834 1069 <field_group id="groupW" > 835 <field field_ref="woce" name="wo" long_name="ocean vertical velocity" /> 836 </field_group> 837 838 <!-- TMB diagnostic output --> 839 <field_group id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant"> 840 <field id="top_temp" name="votemper_top" unit="degC" /> 841 <field id="mid_temp" name="votemper_mid" unit="degC" /> 842 <field id="bot_temp" name="votemper_bot" unit="degC" /> 843 <field id="top_sal" name="vosaline_top" unit="psu" /> 844 <field id="mid_sal" name="vosaline_mid" unit="psu" /> 845 <field id="bot_sal" name="vosaline_bot" unit="psu" /> 846 <field id="sshnmasked" name="sossheig" unit="m" /> 1070 <field field_ref="woce" name="wo" long_name="ocean vertical velocity" /> 847 1071 </field_group> 848 1072 849 1073 <field_group id="1h_grid_U_tmb" grid_ref="grid_U_2D" operation="instant"> 850 851 852 853 1074 <field id="top_u" name="vozocrtx_top" unit="m/s" /> 1075 <field id="mid_u" name="vozocrtx_mid" unit="m/s" /> 1076 <field id="bot_u" name="vozocrtx_bot" unit="m/s" /> 1077 <field id="baro_u" name="vobtcrtx" unit="m/s" /> 854 1078 </field_group> 855 1079 856 1080 <field_group id="1h_grid_V_tmb" grid_ref="grid_V_2D" operation="instant"> 857 858 859 860 1081 <field id="top_v" name="vomecrty_top" unit="m/s" /> 1082 <field id="mid_v" name="vomecrty_mid" unit="m/s" /> 1083 <field id="bot_v" name="vomecrty_bot" unit="m/s" /> 1084 <field id="baro_v" name="vobtcrty" unit="m/s" /> 861 1085 </field_group> 862 1086 863 1087 <!-- 25h diagnostic output --> 864 1088 <field_group id="25h_grid_T" grid_ref="grid_T_3D" operation="instant"> 865 866 867 868 1089 <field id="temper25h" name="potential temperature 25h mean" unit="degC" /> 1090 <field id="tempis25h" name="insitu temperature 25h mean" unit="degC" /> 1091 <field id="salin25h" name="salinity 25h mean" unit="psu" /> 1092 <field id="ssh25h" name="sea surface height 25h mean" grid_ref="grid_T_2D" unit="m" /> 869 1093 </field_group> 870 1094 871 1095 <field_group id="25h_grid_U" grid_ref="grid_U_3D" operation="instant" > 872 1096 <field id="vozocrtx25h" name="i current 25h mean" unit="m/s" /> 873 1097 </field_group> 874 1098 875 1099 <field_group id="25h_grid_V" grid_ref="grid_V_3D" operation="instant"> 876 1100 <field id="vomecrty25h" name="j current 25h mean" unit="m/s" /> 877 1101 </field_group> 878 1102 879 1103 <field_group id="25h_grid_W" grid_ref="grid_W_3D" operation="instant"> 880 <field id="vomecrtz25h" name="k current 25h mean" unit="m/s" />881 <field id="avt25h"name="vertical diffusivity25h mean" unit="m2/s" />882 883 884 885 </field_group> 886 887 1104 <field id="vovecrtz25h" name="k current 25h mean" unit="m/s" /> 1105 <field id="avt25h" name="vertical diffusivity25h mean" unit="m2/s" /> 1106 <field id="avm25h" name="vertical viscosity 25h mean" unit="m2/s" /> 1107 <field id="tke25h" name="turbulent kinetic energy 25h mean" /> 1108 <field id="mxln25h" name="mixing length 25h mean" unit="m" /> 1109 </field_group> 1110 1111 <!-- 888 1112 ============================================================================================================ 889 890 891 1113 --> 1114 <!-- output variables for my configuration (example) --> 1115 892 1116 <field_group id="myvarOCE" > 893 894 895 896 897 898 899 900 901 902 903 904 905 1117 <!-- grid T --> 1118 <field field_ref="e3t" name="e3t" long_name="vertical scale factor" /> 1119 <field field_ref="sst" name="tos" long_name="sea_surface_temperature" /> 1120 <field field_ref="sss" name="sos" long_name="sea_surface_salinity" /> 1121 <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" /> 1122 1123 <!-- grid U --> 1124 <field field_ref="e3u" name="e3u" long_name="vertical scale factor" /> 1125 <field field_ref="ssu" name="uos" long_name="sea_surface_x_velocity" /> 1126 1127 <!-- grid V --> 1128 <field field_ref="e3v" name="e3v" long_name="vertical scale factor" /> 1129 <field field_ref="ssv" name="vos" long_name="sea_surface_y_velocity" /> 906 1130 </field_group> 907 1131 -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/CANAL/EXPREF/namelist_cfg
r12205 r12252 277 277 !! namflo float parameters (default: OFF) 278 278 !! nam_diadct transports through some sections (default: OFF) 279 !! nam_diatmb Top Middle Bottom Output (default: OFF)280 279 !! nam_dia25h 25h Mean Output (default: OFF) 281 280 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/CANAL/MY_SRC/diawri.F90
r11949 r12252 46 46 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 47 47 USE in_out_manager ! I/O manager 48 USE diatmb ! Top,middle,bottom output49 48 USE dia25h ! 25h Mean output 50 49 USE iom ! … … 493 492 CALL iom_put( "bn2", rn2 ) ! Brunt-Vaisala buoyancy frequency (N^2) 494 493 ! 495 496 IF (ln_diatmb) CALL dia_tmb( Kmm ) ! tmb values497 494 498 495 IF (ln_dia25h) CALL dia_25h( kt, Kmm ) ! 25h averaging -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg
r12205 r12252 203 203 !! namflo float parameters (default: OFF) 204 204 !! nam_diadct transports through some sections (default: OFF) 205 !! nam_diatmb Top Middle Bottom Output (default: OFF)206 205 !! nam_dia25h 25h Mean Output (default: OFF) 207 206 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_120pts
r12205 r12252 203 203 !! namflo float parameters (default: OFF) 204 204 !! nam_diadct transports through some sections (default: OFF) 205 !! nam_diatmb Top Middle Bottom Output (default: OFF)206 205 !! nam_dia25h 25h Mean Output (default: OFF) 207 206 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_240pts
r12205 r12252 203 203 !! namflo float parameters (default: OFF) 204 204 !! nam_diadct transports through some sections (default: OFF) 205 !! nam_diatmb Top Middle Bottom Output (default: OFF)206 205 !! nam_dia25h 25h Mean Output (default: OFF) 207 206 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_60pts
r12205 r12252 203 203 !! namflo float parameters (default: OFF) 204 204 !! nam_diadct transports through some sections (default: OFF) 205 !! nam_diatmb Top Middle Bottom Output (default: OFF)206 205 !! nam_dia25h 25h Mean Output (default: OFF) 207 206 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV2D/EXPREF/namelist_cfg
r12205 r12252 203 203 !! namflo float parameters (default: OFF) 204 204 !! nam_diadct transports through some sections (default: OFF) 205 !! nam_diatmb Top Middle Bottom Output (default: OFF)206 205 !! nam_dia25h 25h Mean Output (default: OFF) 207 206 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_AGRIF/EXPREF/1_namelist_cfg
r12205 r12252 203 203 !! namflo float parameters (default: OFF) 204 204 !! nam_diadct transports through some sections (default: OFF) 205 !! nam_diatmb Top Middle Bottom Output (default: OFF)206 205 !! nam_dia25h 25h Mean Output (default: OFF) 207 206 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_AGRIF/EXPREF/namelist_cfg
r12205 r12252 203 203 !! namflo float parameters (default: OFF) 204 204 !! nam_diadct transports through some sections (default: OFF) 205 !! nam_diatmb Top Middle Bottom Output (default: OFF)206 205 !! nam_dia25h 25h Mean Output (default: OFF) 207 206 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ISOMIP/EXPREF/namelist_cfg
r12205 r12252 429 429 !! namflo float parameters (default: OFF) 430 430 !! nam_diadct transports through some sections (default: OFF) 431 !! nam_diatmb Top Middle Bottom Output (default: OFF)432 431 !! nam_dia25h 25h Mean Output (default: OFF) 433 432 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_ubs_cfg
r12205 r12252 277 277 !! namflo float parameters (default: OFF) 278 278 !! nam_diadct transports through some sections (default: OFF) 279 !! nam_diatmb Top Middle Bottom Output (default: OFF)280 279 !! nam_dia25h 25h Mean Output (default: OFF) 281 280 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_flux_ubs_cfg
r12205 r12252 295 295 !! namflo float parameters (default: OFF) 296 296 !! nam_diadct transports through some sections (default: OFF) 297 !! nam_diatmb Top Middle Bottom Output (default: OFF)298 297 !! nam_dia25h 25h Mean Output (default: OFF) 299 298 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/VORTEX/EXPREF/1_namelist_cfg
r12205 r12252 274 274 !! namflo float parameters (default: OFF) 275 275 !! nam_diadct transports through some sections (default: OFF) 276 !! nam_diatmb Top Middle Bottom Output (default: OFF)277 276 !! nam_dia25h 25h Mean Output (default: OFF) 278 277 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/VORTEX/EXPREF/namelist_cfg
r12205 r12252 265 265 !! namflo float parameters (default: OFF) 266 266 !! nam_diadct transports through some sections (default: OFF) 267 !! nam_diatmb Top Middle Bottom Output (default: OFF)268 267 !! nam_dia25h 25h Mean Output (default: OFF) 269 268 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/WAD/EXPREF/namelist_cfg
r12205 r12252 423 423 !! namflo float parameters (default: OFF) 424 424 !! nam_diadct transports through some sections (default: OFF) 425 !! nam_diatmb Top Middle Bottom Output (default: OFF)426 425 !! nam_dia25h 25h Mean Output (default: OFF) 427 426 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")
Note: See TracChangeset
for help on using the changeset viewer.