Changeset 7063
- Timestamp:
- 2016-10-21T15:11:17+02:00 (7 years ago)
- Location:
- branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/CONFIG/SHARED/field_def.xml
r7040 r7063 677 677 <field id="sophtbtr_ind" long_name="Barotropic Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 678 678 <field id="sophtbtr_ipc" long_name="Barotropic Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 679 <field id="sophteiv" long_name="Heat Transport from mesoscale eddy oveection" unit="PW" grid_ref="gznl_T_2D" />679 <field id="sophteiv" long_name="Heat Transport from mesoscale eddy advection" unit="PW" grid_ref="gznl_T_2D" /> 680 680 <field id="sophteiv_atl" long_name="Heat Transport from mesoscale eddy advection: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 681 681 <field id="sophteiv_pac" long_name="Heat Transport from mesoscale eddy advection: Pacific" unit="PW" grid_ref="gznl_T_2D" /> … … 782 782 <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)" unit="W/s^3" /> 783 783 <field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" /> 784 <field id="ketrd_eiv" long_name="Kinetic energy trend from mesoscale eddy advection" unit="W/s^2" 784 <field id="ketrd_eiv" long_name="Kinetic energy trend from mesoscale eddy advection" unit="W/s^2" grid_ref="grid_T_2D" /> 785 785 786 786 <!-- variables available with ln_PE_trd --> … … 802 802 <field id="alphaPE" long_name="partial deriv. of PEanom wrt T" unit="degree_C-1" /> 803 803 <field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3" /> 804 <field id="petrd_zdfp_e3t" unit="W/m^2" grid_ref="grid_T_3D" > petrd_zdfp * e3t </field> 804 805 </field_group> 805 806 … … 1014 1015 <field id="INTPBSI" long_name="Vertically integrated of biogenic Si production" unit="mol/m2/s" /> 1015 1016 <field id="INTPCAL" long_name="Vertically integrated of calcite production" unit="mol/m2/s" /> 1017 <field id="pe_trd" long_name="Tendency of ocean potential energy content" unit="W/m2" /> 1016 1018 1017 1019 <!-- PISCES light : variables available with key_pisces_reduced --> -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90
r6437 r7063 24 24 USE phycst ! physical constant 25 25 USE in_out_manager ! I/O manager 26 USE zdfddm 27 USE zdf_oce 26 28 27 29 IMPLICIT NONE … … 42 44 !! * Substitutions 43 45 # include "domzgr_substitute.h90" 46 # include "zdfddm_substitute.h90" 44 47 !!---------------------------------------------------------------------- 45 48 !! NEMO/OPA 3.3 , NEMO Consortium (2010) … … 75 78 INTEGER :: ji, jj, jk ! dummy loop arguments 76 79 REAL(wp) :: zvolssh, zvol, zssh_steric, zztmp, zarho, ztemp, zsal, zmass 80 REAL(wp) :: zaw, zbw, zrw 77 81 ! 78 82 REAL(wp), POINTER, DIMENSION(:,:) :: zarea_ssh , zbotpres ! 2D workspace 83 REAL(wp), POINTER, DIMENSION(:,:) :: pe ! 2D workspace 79 84 REAL(wp), POINTER, DIMENSION(:,:,:) :: zrhd , zrhop ! 3D workspace 80 85 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztsn ! 4D workspace … … 82 87 IF( nn_timing == 1 ) CALL timing_start('dia_ar5') 83 88 84 CALL wrk_alloc( jpi , jpj , zarea_ssh , zbotpres )89 CALL wrk_alloc( jpi , jpj , zarea_ssh , zbotpres, pe ) 85 90 CALL wrk_alloc( jpi , jpj , jpk , zrhd , zrhop ) 86 91 CALL wrk_alloc( jpi , jpj , jpk , jpts , ztsn ) -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90
r6427 r7063 166 166 ! 167 167 ENDIF 168 IF( l_trddyn ) THEN ! Put here so code doesn't crash when doing KE trend but needs to be done properly 169 CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) 170 ENDIF 168 171 ! 169 172 ELSE ! fixed volume (add the surface pressure gradient + unweighted time stepping) -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90
r7040 r7063 29 29 USE timing ! Timing 30 30 USE diaptr ! Heat/Salt transport diagnostics 31 USE trddyn 32 USE trd_oce 31 33 32 34 IMPLICIT NONE … … 241 243 ENDIF 242 244 243 IF( ln_ trddyn .AND. ln_KE_trd ) CALL trd_dyn(u_eiv, v_eiv, jptrd_eivke, kt )245 IF( ln_KE_trd ) CALL trd_dyn(u_eiv, v_eiv, jpdyn_eivke, kt ) 244 246 # endif 245 247 -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdini.F90
r6427 r7063 91 91 !!gm end 92 92 ! 93 IF( lk_vvl .AND. ( l_trdtra .OR. l_trddyn ) ) CALL ctl_stop( 'trend diagnostics with variable volume not validated' )93 ! IF( lk_vvl .AND. ( l_trdtra .OR. l_trddyn ) ) CALL ctl_stop( 'trend diagnostics with variable volume not validated' ) 94 94 95 95 !!gm : Potential BUG : 3D output only for vector invariant form! add a ctl_stop or code the flux form case -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90
r7040 r7063 105 105 ENDIF 106 106 ! 107 zke(:,:,jpk) = 0._wp 108 zke(1,:, : ) = 0._wp 107 109 zke(:,1, : ) = 0._wp 108 110 DO jk = 1, jpkm1 … … 200 202 ! S = isoneutral slope. 201 203 ! Most terms are on W grid so work on this grid 204 CALL wrk_alloc( jpi, jpj, zke2d ) 202 205 zke2d(:,:) = 0._wp 203 206 DO jk = 1,jpk … … 207 210 & * ( wslpi(ji, jj, jk) * wslpi(ji,jj,jk) & 208 211 & + wslpj(ji, jj, jk) * wslpj(ji,jj,jk) ) & 209 & * bn2(ji,jj,jk) * fse3w(ji, jj, jk)212 & * rn2(ji,jj,jk) * fse3w(ji, jj, jk) 210 213 ENDDO 211 214 ENDDO 212 215 ENDDO 213 iom_put("ketrd_eiv", zke2d) 216 CALL iom_put("ketrd_eiv", zke2d) 217 CALL wrk_dealloc( jpi, jpj, zke2d ) 214 218 ! 215 219 END SELECT -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdpen.F90
r6427 r7063 150 150 rab_pe(:,:,:,:) = 0._wp 151 151 ! 152 IF ( lk_vvl ) CALL ctl_stop('trd_pen_init : PE trends not coded for variable volume')152 ! IF ( lk_vvl ) CALL ctl_stop('trd_pen_init : PE trends not coded for variable volume') 153 153 ! 154 154 nkstp = nit000 - 1
Note: See TracChangeset
for help on using the changeset viewer.