- Timestamp:
- 2011-02-27T17:36:24+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
r2590 r2625 29 29 # include "vectopt_loop_substitute.h90" 30 30 !!---------------------------------------------------------------------- 31 !! NEMO/OPA 3.2 , LODYC-IPSL (2009)31 !! NEMO/OPA 4.0 , NEMO Consortium (2011) 32 32 !! $Id$ 33 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)33 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 34 34 !!---------------------------------------------------------------------- 35 36 35 CONTAINS 37 36 … … 47 46 !! ** Action : (ua,va) updated with the now vorticity term trend 48 47 !!---------------------------------------------------------------------- 49 USE oce, ONLY: zfu => ta ! use ta as 3D workspace 50 USE oce, ONLY: zfv => sa ! use sa as 3D workspace 51 USE wrk_nemo, ONLY: zfu_t => wrk_3d_1, & ! 3D workspaces 52 zfu_f => wrk_3d_2, & 53 zfu_uw =>wrk_3d_3, & 54 zfv_t => wrk_3d_4, & 55 zfv_f => wrk_3d_5, & 56 zfv_vw =>wrk_3d_6, & 57 zfw => wrk_3d_7, & 58 wrk_use, wrk_release 59 IMPLICIT none 48 USE wrk_nemo, ONLY: wrk_use, wrk_release 49 USE oce , ONLY: zfu => ta ! use ta as 3D workspace 50 USE oce , ONLY: zfv => sa ! use sa as 3D workspace 51 USE wrk_nemo, ONLY: zfu_t => wrk_3d_1 , zfv_t => wrk_3d_4 , zfu_uw =>wrk_3d_6 ! 3D workspaces 52 USE wrk_nemo, ONLY: zfu_f => wrk_3d_2 , zfv_f => wrk_3d_5 , zfv_vw =>wrk_3d_7 53 USE wrk_nemo, ONLY: zfw => wrk_3d_3 60 54 !! 61 55 INTEGER, INTENT( in ) :: kt ! ocean time-step index … … 65 59 !!---------------------------------------------------------------------- 66 60 67 IF( kt == nit000 ) THEN68 IF(lwp)WRITE(numout,*)69 IF(lwp)WRITE(numout,*) 'dyn_adv_cen2 : 2nd order flux form momentum advection'70 IF(lwp)WRITE(numout,*) '~~~~~~~~~~~~'61 IF( kt == nit000 .AND. lwp ) THEN 62 WRITE(numout,*) 63 WRITE(numout,*) 'dyn_adv_cen2 : 2nd order flux form momentum advection' 64 WRITE(numout,*) '~~~~~~~~~~~~' 71 65 ENDIF 72 66 73 67 ! Check that global workspace arrays aren't already in use 74 IF( .not. wrk_use(3, 1, 2, 3, 4, 5, 6, 7) )THEN 75 IF(lwp) WRITE(numout, *) 'dyn_adv_cen2 : run-time error - global workspace arrays already in use.' 76 CALL ctl_stop('dyn_adv_cen2 : run-time error - global workspace arrays already in use.') 68 IF( .not. wrk_use(3, 1,2,3,4,5,6,7) ) THEN 69 CALL ctl_stop('dyn_adv_cen2 : requested workspace array unavailable') ; RETURN 77 70 END IF 78 71 … … 169 162 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' ) 170 163 ! 171 ! Flag that the global workspace arrays are no longer in use 172 IF( .not. wrk_release(3, 1, 2, 3, 4, 5, 6, 7) )THEN 173 IF(lwp) WRITE(numout, *) 'dyn_adv_cen2 : run-time error - failed to release global workspace arrays.' 174 END IF 164 IF( .not. wrk_release(3, 1,2,3,4,5,6,7) ) CALL ctl_stop('dyn_adv_cen2 : failed to release workspace array') 175 165 ! 176 166 END SUBROUTINE dyn_adv_cen2
Note: See TracChangeset
for help on using the changeset viewer.