Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90
- Timestamp:
- 2012-01-28T17:44:18+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90
r2715 r3294 23 23 USE diaptr ! poleward transport diagnostics 24 24 USE trc_oce ! share passive tracers/Ocean variables 25 USE wrk_nemo ! Memory Allocation 26 USE timing ! Timing 25 27 26 28 … … 42 44 CONTAINS 43 45 44 SUBROUTINE tra_adv_muscl2( kt, cdtype, p2dt, pun, pvn, pwn, &46 SUBROUTINE tra_adv_muscl2( kt, kit000, cdtype, p2dt, pun, pvn, pwn, & 45 47 & ptb, ptn, pta, kjpt ) 46 48 !!---------------------------------------------------------------------- … … 59 61 !! IPSL, Sept. 2000 (http://www.lodyc.jussieu.fr/opa) 60 62 !!---------------------------------------------------------------------- 61 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 62 USE oce , ONLY: zwx => ua , zwy => va ! (ua,va) used as 3D workspace 63 USE wrk_nemo, ONLY: zslpx => wrk_3d_1 , zslpy => wrk_3d_2 ! 3D workspace 63 USE oce , ONLY: zwx => ua , zwy => va ! (ua,va) used as 3D workspace 64 64 !! 65 65 INTEGER , INTENT(in ) :: kt ! ocean time-step index 66 INTEGER , INTENT(in ) :: kit000 ! first time step index 66 67 CHARACTER(len=3) , INTENT(in ) :: cdtype ! =TRA or TRC (tracer indicator) 67 68 INTEGER , INTENT(in ) :: kjpt ! number of tracers … … 75 76 REAL(wp) :: zv, z0v, zzwy, z0w ! - - 76 77 REAL(wp) :: ztra, zbtr, zdt, zalpha ! - - 78 REAL(wp), POINTER, DIMENSION(:,:,:) :: zslpx, zslpy 77 79 !!---------------------------------------------------------------------- 78 79 IF( wrk_in_use(3, 1,2) ) THEN 80 CALL ctl_stop('tra_adv_muscl2: requested workspace arrays are unavailable') ; RETURN 81 ENDIF 82 83 IF( kt == nit000 ) THEN 80 ! 81 IF( nn_timing == 1 ) CALL timing_start('tra_adv_muscl2') 82 ! 83 CALL wrk_alloc( jpi, jpj, jpk, zslpx, zslpy ) 84 ! 85 86 IF( kt == kit000 ) THEN 84 87 IF(lwp) WRITE(numout,*) 85 88 IF(lwp) WRITE(numout,*) 'tra_adv_muscl2 : MUSCL2 advection scheme on ', cdtype … … 285 288 END DO 286 289 ! 287 IF( wrk_not_released(3, 1,2) ) CALL ctl_stop('tra_adv_muscl2: failed to release workspace arrays') 290 CALL wrk_dealloc( jpi, jpj, jpk, zslpx, zslpy ) 291 ! 292 IF( nn_timing == 1 ) CALL timing_stop('tra_adv_muscl2') 288 293 ! 289 294 END SUBROUTINE tra_adv_muscl2
Note: See TracChangeset
for help on using the changeset viewer.