Changeset 7330
- Timestamp:
- 2016-11-24T17:10:48+01:00 (8 years ago)
- Location:
- branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/CONFIG/SHARED/field_def.xml
r7262 r7330 41 41 <field id="sstgrad2" long_name="square of module of sst gradient" unit="degC2/m2" /> 42 42 <field id="sbt" long_name="sea bottom temperature" unit="degC" /> 43 <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" /> 43 44 <field id="sst_wl" long_name="Delta SST of warm layer" unit="degC" /> 44 45 <field id="sst_cs" long_name="Delta SST of cool skin" unit="degC" /> … … 426 427 <field id="utbl" long_name="zonal current in the Losh tbl" unit="m/s" /> 427 428 428 < !-- variables available with key_diaar5 -->429 <field id="u_masstr " long_name="ocean eulerian mass transport along i-axis" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" />429 <field id="u_masstr" long_name="Ocean Mass X Transport" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" /> 430 <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" /> 430 431 <field id="u_heattr" long_name="ocean eulerian heat transport along i-axis" standard_name="ocean_heat_x_transport" unit="W" /> 431 432 <field id="u_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_x_transport" unit="1e-3*kg/s" /> -
branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90
r7236 r7330 22 22 USE dom_oce ! ocean space and time domain 23 23 USE phycst ! physical constants 24 USE ldftra25 24 ! 26 25 USE iom ! IOM library … … 372 371 ENDIF 373 372 374 IF(ln_ldfeiv) THEN 375 IF( iom_use("sophteiv") .OR. iom_use("sopsteiv") ) THEN 376 z2d(1,:) = htr_eiv(:,1) * rc_pwatt ! (conversion in PW) 377 DO ji = 1, jpi 378 z2d(ji,:) = z2d(1,:) 379 ENDDO 380 cl1 = 'sophteiv' 381 CALL iom_put( TRIM(cl1), z2d ) 382 z2d(1,:) = str_eiv(:,1) * rc_ggram ! (conversion in Gg) 383 DO ji = 1, jpi 384 z2d(ji,:) = z2d(1,:) 385 ENDDO 386 cl1 = 'sopsteiv' 387 CALL iom_put( TRIM(cl1), z2d ) 388 IF( ln_subbas ) THEN 389 DO jn=2,nptr 390 z2d(1,:) = htr_eiv(:,jn) * rc_pwatt ! (conversion in PW) 391 DO ji = 1, jpi 392 z2d(ji,:) = z2d(1,:) 393 ENDDO 394 cl1 = TRIM('sophteiv_'//clsubb(jn)) 395 CALL iom_put( cl1, z2d ) 396 z2d(1,:) = str_eiv(:,jn) * rc_ggram ! (conversion in Gg) 397 DO ji = 1, jpi 398 z2d(ji,:) = z2d(1,:) 399 ENDDO 400 cl1 = TRIM('sopsteiv_'//clsubb(jn)) 401 CALL iom_put( cl1, z2d ) 373 IF( iom_use("sophteiv") .OR. iom_use("sopsteiv") ) THEN 374 z2d(1,:) = htr_eiv(:,1) * rc_pwatt ! (conversion in PW) 375 DO ji = 1, jpi 376 z2d(ji,:) = z2d(1,:) 377 ENDDO 378 cl1 = 'sophteiv' 379 CALL iom_put( TRIM(cl1), z2d ) 380 z2d(1,:) = str_eiv(:,1) * rc_ggram ! (conversion in Gg) 381 DO ji = 1, jpi 382 z2d(ji,:) = z2d(1,:) 383 ENDDO 384 cl1 = 'sopsteiv' 385 CALL iom_put( TRIM(cl1), z2d ) 386 IF( ln_subbas ) THEN 387 DO jn=2,nptr 388 z2d(1,:) = htr_eiv(:,jn) * rc_pwatt ! (conversion in PW) 389 DO ji = 1, jpi 390 z2d(ji,:) = z2d(1,:) 402 391 ENDDO 403 ENDIF 392 cl1 = TRIM('sophteiv_'//clsubb(jn)) 393 CALL iom_put( cl1, z2d ) 394 z2d(1,:) = str_eiv(:,jn) * rc_ggram ! (conversion in Gg) 395 DO ji = 1, jpi 396 z2d(ji,:) = z2d(1,:) 397 ENDDO 398 cl1 = TRIM('sopsteiv_'//clsubb(jn)) 399 CALL iom_put( cl1, z2d ) 400 ENDDO 404 401 ENDIF 405 402 ENDIF -
branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90
r7240 r7330 306 306 z2d(:,:) = 0.e0 307 307 DO jk = 1, jpkm1 308 z3d(:,:,jk) = rau0 * un(:,:,jk) * e2u(:,:) * fse3u(:,:,jk) * umask(:,:,jk)308 z3d(:,:,jk) = rau0 * un(:,:,jk) * e2u(:,:) * e3u_n(:,:,jk) * umask(:,:,jk) 309 309 z2d(:,:) = z2d(:,:) + z3d(:,:,jk) 310 310 END DO -
branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90
r7265 r7330 25 25 USE ldfc1d_c2d ! lateral diffusion: 1D & 2D cases 26 26 USE diaar5, ONLY: lk_diaar5 27 USE diaptr 27 28 ! 28 29 USE trc_oce, ONLY: lk_offline ! offline flag … … 736 737 ! 737 738 zztmp = 0.5_wp * rau0 * rcp 738 IF( iom_use('ueiv_heattr') .OR. iom_use('ueiv_heattr3d') THEN739 IF( iom_use('ueiv_heattr') .OR. iom_use('ueiv_heattr3d')) THEN 739 740 zw2d(:,:) = 0._wp 740 741 zw3d(:,:,:) = 0._wp … … 742 743 DO jj = 2, jpjm1 743 744 DO ji = fs_2, fs_jpim1 ! vector opt. 744 zw3d(ji,jj,jk) = zw3d(ji,jj,jk) + *( psi_uw(ji,jj,jk+1) - psi_uw(ji,jj,jk) ) &745 & 745 zw3d(ji,jj,jk) = zw3d(ji,jj,jk) + ( psi_uw(ji,jj,jk+1) - psi_uw(ji,jj,jk) ) & 746 & * ( tsn (ji,jj,jk,jp_tem) + tsn (ji+1,jj,jk,jp_tem) ) 746 747 zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 747 748 END DO -
branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90
r7250 r7330 88 88 INTEGER :: jk ! dummy loop index 89 89 REAL(wp), POINTER, DIMENSION(:,:,:) :: zun, zvn, zwn 90 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdt, ztrds ! 3D workspace 90 91 !!---------------------------------------------------------------------- 91 92 ! -
branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_fct.F90
r7241 r7330 20 20 USE trdtra ! tracers trends 21 21 USE diaptr ! poleward transport diagnostics 22 USE phycst, ONLY: rau0_rcp 22 23 ! 23 24 USE in_out_manager ! I/O manager 25 USE iom 24 26 USE lib_mpp ! MPP library 25 27 USE lbclnk ! ocean lateral boundary condition (or mpp link) … … 81 83 REAL(wp) :: zfm_ui, zfm_vj, zfm_wk, zC2t_v, zC4t_v ! - - 82 84 REAL(wp), POINTER, DIMENSION(:,:,:) :: zwi, zwx, zwy, zwz, ztu, ztv, zltu, zltv, ztw 83 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdx, ztrdy, ztrd , zptry85 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdx, ztrdy, ztrdz, zptry 84 86 REAL(wp), POINTER, DIMENSION(:,:) :: z2d 85 87 !!----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.