Changeset 3294 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90
- Timestamp:
- 2012-01-28T17:44:18+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90
r2777 r3294 51 51 USE thd_ice ! LIM thermodynamics 52 52 USE in_out_manager ! I/O manager 53 USE lib_mpp ! MPP library 53 USE lib_mpp ! MPP library 54 USE wrk_nemo ! work arrays 54 55 55 56 IMPLICIT NONE … … 297 298 !! ** References : Vancoppenolle et al., 2007 (in preparation) 298 299 !!------------------------------------------------------------------ 299 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released300 USE wrk_nemo, ONLY: wrk_3d_3, wrk_3d_4301 300 INTEGER :: ji, jj, jk, jl ! dummy loop index 302 301 REAL(wp) :: dummy_fac0, dummy_fac1, dummy_fac, zsal ! local scalar 303 302 REAL(wp) :: zind0, zind01, zindbal, zargtemp , zs_zero ! - - 304 !305 303 REAL(wp), POINTER, DIMENSION(:,:,:) :: z_slope_s, zalpha ! 3D pointer 306 304 !!------------------------------------------------------------------ 307 305 308 IF( wrk_in_use( 3, 3,4 ) ) THEN 309 CALL ctl_stop( 'lim_var_salprof: requested workspace arrays unavailable' ) ; RETURN 310 END IF 311 312 z_slope_s => wrk_3d_3(:,:,1:jpl) ! slope of the salinity profile 313 zalpha => wrk_3d_4(:,:,1:jpl) ! weight factor for s between s_i_0 and s_i_1 306 CALL wrk_alloc( jpi, jpj, jpl, z_slope_s, zalpha ) 314 307 315 308 !--------------------------------------- … … 390 383 ENDIF ! num_sal 391 384 ! 392 IF( wrk_not_released(3, 3,4) ) CALL ctl_stop('lim_var_salprof: failed to release workspace arrays.')385 CALL wrk_dealloc( jpi, jpj, jpl, z_slope_s, zalpha ) 393 386 ! 394 387 END SUBROUTINE lim_var_salprof … … 433 426 !! Works with 1d vectors and is used by thermodynamic modules 434 427 !!------------------------------------------------------------------- 435 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released436 USE wrk_nemo, ONLY: wrk_1d_4437 428 INTEGER, INTENT(in) :: kideb, kiut ! thickness category index 438 429 ! … … 445 436 !!--------------------------------------------------------------------- 446 437 447 IF( wrk_in_use(1, 4) ) THEN 448 CALL ctl_stop('lim_var_salprof1d : requestead workspace arrays unavailable.') ; RETURN 449 END IF 450 ! Set-up pointers to sub-arrays of workspace arrays 451 z_slope_s => wrk_1d_4 (1:jpij) 438 CALL wrk_alloc( jpij, z_slope_s ) 452 439 453 440 !--------------------------------------- … … 514 501 ENDIF 515 502 ! 516 IF( wrk_not_released(1, 4) ) CALL ctl_stop( 'lim_var_salprof1d : failed to release workspace arrays')503 CALL wrk_dealloc( jpij, z_slope_s ) 517 504 ! 518 505 END SUBROUTINE lim_var_salprof1d
Note: See TracChangeset
for help on using the changeset viewer.