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 5038 for branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90 – NEMO

Ignore:
Timestamp:
2015-01-20T15:26:13+01:00 (9 years ago)
Author:
jamesharle
Message:

Merging branch with HEAD of the trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90

    r4609 r5038  
    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 
     
    161169      READ  ( numnatp_cfg, namlobsed, IOSTAT = ios, ERR = 902 ) 
    162170902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobsed in configuration namelist', lwp ) 
    163       WRITE ( numonp, namlobsed ) 
     171      IF(lwm) WRITE ( numonp, namlobsed ) 
    164172 
    165173      IF(lwp) THEN 
Note: See TracChangeset for help on using the changeset viewer.