- Timestamp:
- 2011-03-15T16:27:46+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90
r2678 r2690 74 74 !! ** Action : - pta becomes the after tracer 75 75 !!--------------------------------------------------------------------- 76 USE oce , ONLY : zwd => ua ! ua used as workspace 77 USE oce , ONLY : zws => va ! va - - 78 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 79 USE wrk_nemo, ONLY: zwi => wrk_3d_1, zwt => wrk_3d_2 ! workspace arrays 80 !! 76 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 77 USE oce , ONLY: zwd => ua , zws => va ! (ua,va) used as 3D workspace 78 USE wrk_nemo, ONLY: zwi => wrk_3d_1 , zwt => wrk_3d_2 ! 3D workspace 79 ! 81 80 INTEGER , INTENT(in ) :: kt ! ocean time-step index 82 81 CHARACTER(len=3) , INTENT(in ) :: cdtype ! =TRA or TRC (tracer indicator) … … 85 84 REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in ) :: ptb ! before and now tracer fields 86 85 REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) :: pta ! tracer trend 87 !! 88 INTEGER :: ji, jj, jk, jn ! dummy loop indices 89 REAL(wp) :: zrhs ! local scalars 90 REAL(wp) :: ze3tb, ze3tn, ze3ta ! variable vertical scale factors 86 ! 87 INTEGER :: ji, jj, jk, jn ! dummy loop indices 88 REAL(wp) :: zrhs, ze3tb, ze3tn, ze3ta ! local scalars 91 89 !!--------------------------------------------------------------------- 92 90 93 IF(wrk_in_use(3, 1,2))THEN 94 CALL ctl_stop('tra_zdf_imp : requested workspace arrays unavailable.') 95 RETURN 96 END IF 91 IF( wrk_in_use(3, 1,2) ) THEN 92 CALL ctl_stop('tra_zdf_imp : requested workspace arrays unavailable.') ; RETURN 93 ENDIF 97 94 98 95 IF( kt == nit000 ) THEN … … 114 111 ! Build matrix if temperature or salinity (only in double diffusion case) or first passive tracer 115 112 ! 116 IF( ( cdtype == 'TRA' .AND. ( ( jn == jp_tem ) .OR. ( jn == jp_sal .AND. lk_zdfddm ) ) ) .OR.&113 IF( ( cdtype == 'TRA' .AND. ( jn == jp_tem .OR. ( jn == jp_sal .AND. lk_zdfddm ) ) ) .OR. & 117 114 & ( cdtype == 'TRC' .AND. jn == 1 ) ) THEN 118 115 ! … … 231 228 ! ! ================= ! 232 229 ! 233 IF(wrk_not_released(3, 1,2))THEN 234 CALL ctl_stop('tra_zdf_imp : failed to release workspace arrays.') 235 END IF 230 IF( wrk_not_released(3, 1,2) ) CALL ctl_stop('tra_zdf_imp: failed to release workspace arrays') 236 231 ! 237 232 END SUBROUTINE tra_zdf_imp
Note: See TracChangeset
for help on using the changeset viewer.