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 5034 for branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z – NEMO

Ignore:
Timestamp:
2015-01-15T14:48:42+01:00 (9 years ago)
Author:
andrewryan
Message:

merge with trunk

Location:
branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90

    r4624 r5034  
    2121   USE lbclnk          !  
    2222   USE prtctl_trc      ! Print control for debbuging 
    23    USE trdmod_oce 
    24    USE trdmod_trc 
     23   USE trd_oce 
     24   USE trdtrc 
    2525   USE iom 
    2626    
    … …  
    110110      IF( nn_timing == 1 )  CALL timing_start('p2z_bio') 
    111111      ! 
    112       IF( ln_diatrc ) THEN 
     112      IF( ln_diatrc .OR. lk_iomput ) THEN 
    113113         CALL wrk_alloc( jpi, jpj,     17, zw2d ) 
    114114         CALL wrk_alloc( jpi, jpj, jpk, 3, zw3d ) 
    … …  
    122122 
    123123      xksi(:,:) = 0.e0        ! zooplakton closure ( fbod) 
    124       IF( ln_diatrc ) THEN 
     124      IF( ln_diatrc .OR. lk_iomput ) THEN 
    125125         zw2d  (:,:,:) = 0.e0 
    126126         zw3d(:,:,:,:) = 0.e0 
    … …  
    239239                  !  trend number 17 in p2zexp 
    240240                ENDIF 
    241                 IF( ln_diatrc ) THEN 
     241                IF( ln_diatrc .OR. lk_iomput ) THEN 
    242242                  ! convert fluxes in per day 
    243243                  ze3t = fse3t(ji,jj,jk) * 86400. 
    … …  
    360360                  !  trend number 17 in p2zexp  
    361361                ENDIF 
    362                 IF( ln_diatrc ) THEN 
     362                IF( ln_diatrc .OR. lk_iomput ) THEN 
    363363                  ! convert fluxes in per day 
    364364                  ze3t = fse3t(ji,jj,jk) * 86400. 
    … …  
    390390      END DO 
    391391 
    392       IF( ln_diatrc ) THEN 
    393          ! 
     392      IF( ln_diatrc .OR. lk_iomput ) THEN 
    394393         DO jl = 1, 17  
    395394            CALL lbc_lnk( zw2d(:,:,jl),'T', 1. ) 
    … …  
    398397            CALL lbc_lnk( zw3d(:,:,:,jl),'T', 1. ) 
    399398         END DO 
    400          IF( lk_iomput ) THEN 
     399      ENDIF 
     400      IF( lk_iomput ) THEN 
    401401            ! Save diagnostics 
    402             CALL iom_put( "TNO3PHY", zw2d(:,:,1) ) 
    403             CALL iom_put( "TNH4PHY", zw2d(:,:,2) ) 
    404             CALL iom_put( "TPHYDOM", zw2d(:,:,3) ) 
    405             CALL iom_put( "TPHYNH4", zw2d(:,:,4) ) 
    406             CALL iom_put( "TPHYZOO", zw2d(:,:,5) ) 
    407             CALL iom_put( "TPHYDET", zw2d(:,:,6) ) 
    408             CALL iom_put( "TDETZOO", zw2d(:,:,7) ) 
    409             CALL iom_put( "TZOODET", zw2d(:,:,8) ) 
    410             CALL iom_put( "TZOOBOD", zw2d(:,:,9) ) 
    411             CALL iom_put( "TZOONH4", zw2d(:,:,10) ) 
    412             CALL iom_put( "TZOODOM", zw2d(:,:,11) ) 
    413             CALL iom_put( "TNH4NO3", zw2d(:,:,12) ) 
    414             CALL iom_put( "TDOMNH4", zw2d(:,:,13) ) 
    415             CALL iom_put( "TDETNH4", zw2d(:,:,14) ) 
    416             CALL iom_put( "TPHYTOT", zw2d(:,:,15) ) 
    417             CALL iom_put( "TZOOTOT", zw2d(:,:,16) ) 
    418             !  
    419             CALL iom_put( "FNO3PHY", zw3d(:,:,:,1) ) 
    420             CALL iom_put( "FNH4PHY", zw3d(:,:,:,2) ) 
    421             CALL iom_put( "FNH4NO3", zw3d(:,:,:,3) ) 
    422             ! 
    423          ELSE 
     402        CALL iom_put( "TNO3PHY", zw2d(:,:,1) ) 
     403        CALL iom_put( "TNH4PHY", zw2d(:,:,2) ) 
     404        CALL iom_put( "TPHYDOM", zw2d(:,:,3) ) 
     405        CALL iom_put( "TPHYNH4", zw2d(:,:,4) ) 
     406        CALL iom_put( "TPHYZOO", zw2d(:,:,5) ) 
     407        CALL iom_put( "TPHYDET", zw2d(:,:,6) ) 
     408        CALL iom_put( "TDETZOO", zw2d(:,:,7) ) 
     409        CALL iom_put( "TZOODET", zw2d(:,:,8) ) 
     410        CALL iom_put( "TZOOBOD", zw2d(:,:,9) ) 
     411        CALL iom_put( "TZOONH4", zw2d(:,:,10) ) 
     412        CALL iom_put( "TZOODOM", zw2d(:,:,11) ) 
     413        CALL iom_put( "TNH4NO3", zw2d(:,:,12) ) 
     414        CALL iom_put( "TDOMNH4", zw2d(:,:,13) ) 
     415        CALL iom_put( "TDETNH4", zw2d(:,:,14) ) 
     416        CALL iom_put( "TPHYTOT", zw2d(:,:,15) ) 
     417        CALL iom_put( "TZOOTOT", zw2d(:,:,16) ) 
     418         !  
     419        CALL iom_put( "FNO3PHY", zw3d(:,:,:,1) ) 
     420        CALL iom_put( "FNH4PHY", zw3d(:,:,:,2) ) 
     421        CALL iom_put( "FNH4NO3", zw3d(:,:,:,3) ) 
     422         ! 
     423       ELSE 
     424          IF( ln_diatrc ) THEN 
    424425            ! 
    425426            trc2d(:,:,jp_pcs0_2d    ) = zw2d(:,:,1)  
    … …  
    457458      IF( l_trdtrc ) THEN 
    458459         DO jl = jp_pcs0_trd, jp_pcs1_trd 
    459             CALL trd_mod_trc( trbio(:,:,:,jl), jl, kt )   ! handle the trend 
     460            CALL trd_trc( trbio(:,:,:,jl), jl, kt )   ! handle the trend 
    460461         END DO 
    461462      ENDIF 
    … …  
    467468      ENDIF 
    468469      ! 
    469       IF( ln_diatrc ) THEN 
     470      IF( ln_diatrc .OR. lk_iomput ) THEN 
    470471         CALL wrk_dealloc( jpi, jpj,     17, zw2d ) 
    471472         CALL wrk_dealloc( jpi, jpj, jpk, 3, zw3d ) 

  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zexp.F90

    r3446 r5034  
    2222   USE lbclnk 
    2323   USE prtctl_trc      ! Print control for debbuging 
    24    USE trdmod_oce 
    25    USE trdmod_trc 
     24   USE trd_oce 
     25   USE trdtrc 
    2626   USE iom 
    2727 
    … …  
    127127  
    128128      ! Oa & Ek: diagnostics depending on jpdia2d !          left as example 
    129       IF( ln_diatrc ) THEN 
    130          IF( lk_iomput ) THEN   ;   CALL iom_put( "SEDPOC" , sedpocn ) 
    131          ELSE                   ;   trc2d(:,:,jp_pcs0_2d + 18) = sedpocn(:,:) 
    132          ENDIF 
     129      IF( lk_iomput ) THEN   
     130         CALL iom_put( "SEDPOC" , sedpocn ) 
     131      ELSE 
     132         IF( ln_diatrc )           trc2d(:,:,jp_pcs0_2d + 18) = sedpocn(:,:) 
    133133      ENDIF 
    134134 
    … …  
    164164         ztrbio(:,:,:) = tra(:,:,:,jpno3) - ztrbio(:,:,:) 
    165165         jl = jp_pcs0_trd + 16 
    166          CALL trd_mod_trc( ztrbio, jl, kt )   ! handle the trend 
     166         CALL trd_trc( ztrbio, jl, kt )   ! handle the trend 
    167167         CALL wrk_dealloc( jpi, jpj, jpk, ztrbio )   ! temporary save of trends 
    168168      ENDIF 

  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zopt.F90

    r4624 r5034  
    128128              IF( etot(ji,jj,jk) >= zpar100(ji,jj) )   neln(ji,jj) = jk + 1  
    129129              !                                       ! nb. this is to ensure compatibility with 
    130               !                                       ! nmld_trc definition in trd_mld_trc_zint 
     130              !                                       ! nmld_trc definition in trd_mxl_trc_zint 
    131131           END DO 
    132132         END DO 

  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90

    r4624 r5034  
    1818   USE sms_pisces 
    1919   USE lbclnk 
    20    USE trdmod_oce 
    21    USE trdmod_trc 
     20   USE trd_oce 
     21   USE trdtrc 
    2222   USE iom 
    2323   USE prtctl_trc      ! Print control for debbuging 
    … …  
    108108      END DO 
    109109 
    110       IF( ln_diatrc ) THEN  
    111          CALL wrk_alloc( jpi, jpj, zw2d ) 
    112          zw2d(:,:) =  ztra(:,:,1) * fse3t(:,:,1) * 86400. 
    113          DO jk = 2, jpkm1 
    114             zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * fse3t(:,:,jk) * 86400. 
    115          END DO 
    116          IF( lk_iomput )  THEN 
    117            CALL iom_put( "TDETSED", zw2d ) 
    118          ELSE 
    119            trc2d(:,:,jp_pcs0_2d + 7) = zw2d(:,:) 
     110      IF( lk_iomput )  THEN 
     111         IF( iom_use( "TDETSED" ) ) THEN 
     112            CALL wrk_alloc( jpi, jpj, zw2d ) 
     113            zw2d(:,:) =  ztra(:,:,1) * fse3t(:,:,1) * 86400. 
     114            DO jk = 2, jpkm1 
     115               zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * fse3t(:,:,jk) * 86400. 
     116            END DO 
     117            CALL iom_put( "TDETSED", zw2d ) 
     118            CALL wrk_dealloc( jpi, jpj, zw2d ) 
    120119         ENDIF 
    121          CALL wrk_dealloc( jpi, jpj, zw2d ) 
     120      ELSE 
     121         IF( ln_diatrc ) THEN  
     122            CALL wrk_alloc( jpi, jpj, zw2d ) 
     123            zw2d(:,:) =  ztra(:,:,1) * fse3t(:,:,1) * 86400. 
     124            DO jk = 2, jpkm1 
     125               zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * fse3t(:,:,jk) * 86400. 
     126            END DO 
     127            trc2d(:,:,jp_pcs0_2d + 7) = zw2d(:,:) 
     128            CALL wrk_dealloc( jpi, jpj, zw2d ) 
     129         ENDIF 
    122130      ENDIF 
    123131      ! 
    … …  
    128136         ztrbio(:,:,:) = tra(:,:,:,jpdet) - ztrbio(:,:,:) 
    129137         jl = jp_pcs0_trd + 7 
    130          CALL trd_mod_trc( ztrbio, jl, kt )   ! handle the trend 
     138         CALL trd_trc( ztrbio, jl, kt )   ! handle the trend 
    131139         CALL wrk_dealloc( jpi, jpj, jpk, ztrbio ) 
    132140      ENDIF 

  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsms.F90

    r4624 r5034  
    2020   USE p2zsed 
    2121   USE p2zexp 
    22    USE trdmod_oce 
    23    USE trdmod_trc_oce 
    24    USE trdmod_trc 
    25    USE trdmld_trc 
     22   USE trd_oce 
     23   USE trdtrc_oce 
     24   USE trdtrc 
     25   USE trdmxl_trc 
    2626 
    2727   IMPLICIT NONE 
    … …  
    6161      IF( l_trdtrc ) THEN 
    6262         DO jn = jp_pcs0, jp_pcs1 
    63            CALL trd_mod_trc( tra(:,:,:,jn), jn, jptra_trd_sms, kt )   ! save trends 
     63           CALL trd_trc( tra(:,:,:,jn), jn, jptra_sms, kt )   ! save trends 
    6464         END DO 
    6565      END IF 
    6666 
    67       IF( lk_trdmld_trc )  CALL trd_mld_bio( kt )   ! trends: Mixed-layer 
     67      IF( lk_trdmxl_trc )  CALL trd_mxl_bio( kt )   ! trends: Mixed-layer 
    6868      ! 
    6969      IF ( lwm .AND. kt == nittrc000 ) CALL FLUSH    ( numonp )     ! flush output namelist PISCES 
Note: See TracChangeset for help on using the changeset viewer.