New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7330 for branches/2016/dev_r7233_CMIP6_diags_trunk_version – NEMO

Ignore:
Timestamp:
2016-11-24T17:10:48+01:00 (7 years ago)
Author:
timgraham
Message:

Final bug fixes needed for compilation

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  
    4141         <field id="sstgrad2"     long_name="square of module of sst gradient"                                                        unit="degC2/m2" /> 
    4242         <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" /> 
    4344         <field id="sst_wl"       long_name="Delta SST of warm layer"                                                                 unit="degC"     /> 
    4445         <field id="sst_cs"       long_name="Delta SST of cool skin"                                                                  unit="degC"     /> 
    … …  
    426427         <field id="utbl"         long_name="zonal current in the Losh tbl"     unit="m/s" /> 
    427428 
    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" /> 
    430431         <field id="u_heattr"     long_name="ocean eulerian heat transport along i-axis"    standard_name="ocean_heat_x_transport"                          unit="W"                                /> 
    431432         <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  
    2222   USE dom_oce          ! ocean space and time domain 
    2323   USE phycst           ! physical constants 
    24    USE ldftra  
    2524   ! 
    2625   USE iom              ! IOM library 
    … …  
    372371         ENDIF 
    373372 
    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,:) 
    402391                  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 
    404401            ENDIF 
    405402         ENDIF 

  • branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r7240 r7330  
    306306         z2d(:,:) = 0.e0 
    307307         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) 
    309309            z2d(:,:) = z2d(:,:) + z3d(:,:,jk) 
    310310         END DO 

  • branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

    r7265 r7330  
    2525   USE ldfc1d_c2d      ! lateral diffusion: 1D & 2D cases  
    2626   USE diaar5, ONLY:   lk_diaar5 
     27   USE diaptr 
    2728   ! 
    2829   USE trc_oce, ONLY: lk_offline ! offline flag 
    … …  
    736737         ! 
    737738         zztmp = 0.5_wp * rau0 * rcp  
    738          IF( iom_use('ueiv_heattr') .OR. iom_use('ueiv_heattr3d') THEN 
     739         IF( iom_use('ueiv_heattr') .OR. iom_use('ueiv_heattr3d')) THEN 
    739740           zw2d(:,:) = 0._wp  
    740741           zw3d(:,:,:) = 0._wp  
    … …  
    742743              DO jj = 2, jpjm1 
    743744                 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                        &                              * ( tsn   (ji,jj,jk,jp_tem) + tsn   (ji+1,jj,jk,jp_tem) )  
     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) )  
    746747                    zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
    747748                 END DO 

  • branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

    r7250 r7330  
    8888      INTEGER ::   jk   ! dummy loop index 
    8989      REAL(wp), POINTER, DIMENSION(:,:,:) :: zun, zvn, zwn 
     90      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztrdt, ztrds   ! 3D workspace 
    9091      !!---------------------------------------------------------------------- 
    9192      ! 

  • branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_fct.F90

    r7241 r7330  
    2020   USE trdtra         ! tracers trends 
    2121   USE diaptr         ! poleward transport diagnostics 
     22   USE phycst, ONLY: rau0_rcp 
    2223   ! 
    2324   USE in_out_manager ! I/O manager 
     25   USE iom 
    2426   USE lib_mpp        ! MPP library 
    2527   USE lbclnk         ! ocean lateral boundary condition (or mpp link)  
    … …  
    8183      REAL(wp) ::   zfm_ui, zfm_vj, zfm_wk, zC2t_v, zC4t_v   !   -      - 
    8284      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zwi, zwx, zwy, zwz, ztu, ztv, zltu, zltv, ztw 
    83       REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztrdx, ztrdy, ztrd, zptry 
     85      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztrdx, ztrdy, ztrdz, zptry 
    8486      REAL(wp), POINTER, DIMENSION(:,:)   :: z2d 
    8587      !!---------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.