Changeset 10954 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/TRA/tradmp.F90
- Timestamp:
- 2019-05-09T18:12:29+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/TRA/tradmp.F90
r10946 r10954 72 72 73 73 74 SUBROUTINE tra_dmp( kt, K mm, Krhs )74 SUBROUTINE tra_dmp( kt, Kbb, Kmm, Krhs ) 75 75 !!---------------------------------------------------------------------- 76 76 !! *** ROUTINE tra_dmp *** … … 91 91 !!---------------------------------------------------------------------- 92 92 INTEGER, INTENT(in) :: kt ! ocean time-step index 93 INTEGER, INTENT(in) :: K mm, Krhs ! time level indices93 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 94 94 ! 95 95 INTEGER :: ji, jj, jk, jn ! dummy loop indices … … 102 102 IF( l_trdtra ) THEN !* Save ta and sa trends 103 103 ALLOCATE( ztrdts(jpi,jpj,jpk,jpts) ) 104 ztrdts(:,:,:,:) = ts a(:,:,:,:)104 ztrdts(:,:,:,:) = ts(:,:,:,:,Krhs) 105 105 ENDIF 106 106 ! !== input T-S data at kt ==! … … 114 114 DO jj = 2, jpjm1 115 115 DO ji = fs_2, fs_jpim1 ! vector opt. 116 ts a(ji,jj,jk,jn) = tsa(ji,jj,jk,jn) + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jn) - tsb(ji,jj,jk,jn) )116 ts(ji,jj,jk,jn,Krhs) = ts(ji,jj,jk,jn,Krhs) + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jn) - ts(ji,jj,jk,jn,Kbb) ) 117 117 END DO 118 118 END DO … … 125 125 DO ji = fs_2, fs_jpim1 ! vector opt. 126 126 IF( avt(ji,jj,jk) <= avt_c ) THEN 127 ts a(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) &128 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_tem) - ts b(ji,jj,jk,jp_tem) )129 ts a(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) &130 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_sal) - ts b(ji,jj,jk,jp_sal) )127 ts(ji,jj,jk,jp_tem,Krhs) = ts(ji,jj,jk,jp_tem,Krhs) & 128 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_tem) - ts(ji,jj,jk,jp_tem,Kbb) ) 129 ts(ji,jj,jk,jp_sal,Krhs) = ts(ji,jj,jk,jp_sal,Krhs) & 130 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_sal) - ts(ji,jj,jk,jp_sal,Kbb) ) 131 131 ENDIF 132 132 END DO … … 138 138 DO jj = 2, jpjm1 139 139 DO ji = fs_2, fs_jpim1 ! vector opt. 140 IF( gdept _n(ji,jj,jk) >= hmlp (ji,jj) ) THEN141 ts a(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) &142 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_tem) - ts b(ji,jj,jk,jp_tem) )143 ts a(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) &144 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_sal) - ts b(ji,jj,jk,jp_sal) )140 IF( gdept(ji,jj,jk,Kmm) >= hmlp (ji,jj) ) THEN 141 ts(ji,jj,jk,jp_tem,Krhs) = ts(ji,jj,jk,jp_tem,Krhs) & 142 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_tem) - ts(ji,jj,jk,jp_tem,Kbb) ) 143 ts(ji,jj,jk,jp_sal,Krhs) = ts(ji,jj,jk,jp_sal,Krhs) & 144 & + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jp_sal) - ts(ji,jj,jk,jp_sal,Kbb) ) 145 145 ENDIF 146 146 END DO … … 151 151 ! 152 152 IF( l_trdtra ) THEN ! trend diagnostic 153 ztrdts(:,:,:,:) = ts a(:,:,:,:) - ztrdts(:,:,:,:)153 ztrdts(:,:,:,:) = ts(:,:,:,:,Krhs) - ztrdts(:,:,:,:) 154 154 CALL trd_tra( kt, Kmm, Krhs, 'TRA', jp_tem, jptra_dmp, ztrdts(:,:,:,jp_tem) ) 155 155 CALL trd_tra( kt, Kmm, Krhs, 'TRA', jp_sal, jptra_dmp, ztrdts(:,:,:,jp_sal) ) … … 157 157 ENDIF 158 158 ! ! Control print 159 IF(ln_ctl) CALL prt_ctl( tab3d_1=ts a(:,:,:,jp_tem), clinfo1=' dmp - Ta: ', mask1=tmask, &160 & tab3d_2=ts a(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' )159 IF(ln_ctl) CALL prt_ctl( tab3d_1=ts(:,:,:,jp_tem,Krhs), clinfo1=' dmp - Ta: ', mask1=tmask, & 160 & tab3d_2=ts(:,:,:,jp_sal,Krhs), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 161 161 ! 162 162 IF( ln_timing ) CALL timing_stop('tra_dmp')
Note: See TracChangeset
for help on using the changeset viewer.