Changeset 13067
- Timestamp:
- 2020-06-09T10:12:15+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/OCE/TRD/trdtra.F90
r10425 r13067 81 81 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in), OPTIONAL :: ptra ! now tracer variable 82 82 ! 83 INTEGER :: jk ! loop indices 83 INTEGER :: jk ! loop indices 84 INTEGER :: i01 ! 0 or 1 84 85 REAL(wp), DIMENSION(jpi,jpj,jpk) :: ztrds ! 3D workspace 85 86 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zwt, zws, ztrdt ! 3D workspace … … 89 90 IF( trd_tra_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'trd_tra : unable to allocate arrays' ) 90 91 ENDIF 91 92 ! 93 i01 = COUNT( (/ PRESENT(pun) .OR. ( ktrd /= jptra_xad .AND. ktrd /= jptra_yad .AND. ktrd /= jptra_zad ) /) ) 94 ! 92 95 IF( ctype == 'TRA' .AND. ktra == jp_tem ) THEN !== Temperature trend ==! 93 96 ! 94 SELECT CASE( ktrd )97 SELECT CASE( ktrd*i01 ) 95 98 ! ! advection: transform the advective flux into a trend 96 99 CASE( jptra_xad ) ; CALL trd_tra_adv( ptrd, pun, ptra, 'X', trdtx ) … … 111 114 IF( ctype == 'TRA' .AND. ktra == jp_sal ) THEN !== Salinity trends ==! 112 115 ! 113 SELECT CASE( ktrd )116 SELECT CASE( ktrd*i01 ) 114 117 ! ! advection: transform the advective flux into a trend 115 118 ! ! and send T & S trends to trd_tra_mng … … 162 165 IF( ctype == 'TRC' ) THEN !== passive tracer trend ==! 163 166 ! 164 SELECT CASE( ktrd )167 SELECT CASE( ktrd*i01 ) 165 168 ! ! advection: transform the advective flux into a masked trend 166 169 CASE( jptra_xad ) ; CALL trd_tra_adv( ptrd , pun , ptra, 'X', ztrds )
Note: See TracChangeset
for help on using the changeset viewer.