Changeset 5038 for branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90
- Timestamp:
- 2015-01-20T15:26:13+01:00 (9 years ago)
- 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 18 18 USE sms_pisces 19 19 USE lbclnk 20 USE trd mod_oce21 USE trd mod_trc20 USE trd_oce 21 USE trdtrc 22 22 USE iom 23 23 USE prtctl_trc ! Print control for debbuging … … 108 108 END DO 109 109 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 ) 120 119 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 122 130 ENDIF 123 131 ! … … 128 136 ztrbio(:,:,:) = tra(:,:,:,jpdet) - ztrbio(:,:,:) 129 137 jl = jp_pcs0_trd + 7 130 CALL trd_ mod_trc( ztrbio, jl, kt ) ! handle the trend138 CALL trd_trc( ztrbio, jl, kt ) ! handle the trend 131 139 CALL wrk_dealloc( jpi, jpj, jpk, ztrbio ) 132 140 ENDIF … … 161 169 READ ( numnatp_cfg, namlobsed, IOSTAT = ios, ERR = 902 ) 162 170 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobsed in configuration namelist', lwp ) 163 WRITE ( numonp, namlobsed )171 IF(lwm) WRITE ( numonp, namlobsed ) 164 172 165 173 IF(lwp) THEN
Note: See TracChangeset
for help on using the changeset viewer.