Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
- Timestamp:
- 2012-01-28T17:44:18+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
r2715 r3294 20 20 USE lib_mpp ! MPP library 21 21 USE prtctl ! Print control 22 USE wrk_nemo ! Memory Allocation 23 USE timing ! Timing 22 24 23 25 IMPLICIT NONE … … 47 49 !! ** Action : (ua,va) updated with the now vorticity term trend 48 50 !!---------------------------------------------------------------------- 49 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released50 USE oce , ONLY: zfu => ta , zfv => sa ! (ta,sa) used as 3D workspace51 USE wrk_nemo, ONLY: zfu_t => wrk_3d_1 , zfv_t => wrk_3d_4 , zfu_uw =>wrk_3d_6 ! 3D workspaces52 USE wrk_nemo, ONLY: zfu_f => wrk_3d_2 , zfv_f => wrk_3d_5 , zfv_vw =>wrk_3d_753 USE wrk_nemo, ONLY: zfw => wrk_3d_354 !55 51 INTEGER, INTENT( in ) :: kt ! ocean time-step index 56 52 ! 57 53 INTEGER :: ji, jj, jk ! dummy loop indices 58 54 REAL(wp) :: zbu, zbv ! local scalars 55 REAL(wp), POINTER, DIMENSION(:,:,:) :: zfu_t, zfv_t, zfu_f, zfv_f, zfu_uw, zfv_vw, zfw 56 REAL(wp), POINTER, DIMENSION(:,:,:) :: zfu, zfv 59 57 !!---------------------------------------------------------------------- 60 58 ! 59 IF( nn_timing == 1 ) CALL timing_start('dyn_adv_cen2') 60 ! 61 CALL wrk_alloc( jpi, jpj, jpk, zfu_t, zfv_t, zfu_f, zfv_f, zfu_uw, zfv_vw, zfu, zfv, zfw ) 62 ! 61 63 IF( kt == nit000 .AND. lwp ) THEN 62 64 WRITE(numout,*) … … 64 66 WRITE(numout,*) '~~~~~~~~~~~~' 65 67 ENDIF 66 67 ! Check that global workspace arrays aren't already in use 68 IF( wrk_in_use(3, 1,2,3,4,5,6,7) ) THEN 69 CALL ctl_stop('dyn_adv_cen2 : requested workspace array unavailable') ; RETURN 70 ENDIF 71 68 ! 72 69 IF( l_trddyn ) THEN ! Save ua and va trends 73 70 zfu_uw(:,:,:) = ua(:,:,:) … … 162 159 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' ) 163 160 ! 164 IF( wrk_not_released(3, 1,2,3,4,5,6,7) ) CALL ctl_stop('dyn_adv_cen2: failed to release workspace array') 161 CALL wrk_dealloc( jpi, jpj, jpk, zfu_t, zfv_t, zfu_f, zfv_f, zfu_uw, zfv_vw, zfu, zfv, zfw ) 162 ! 163 IF( nn_timing == 1 ) CALL timing_stop('dyn_adv_cen2') 165 164 ! 166 165 END SUBROUTINE dyn_adv_cen2
Note: See TracChangeset
for help on using the changeset viewer.