- Timestamp:
- 2011-11-25T16:42:54+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/NST_SRC/agrif_opa_update.F90
r3104 r3184 9 9 USE in_out_manager ! I/O manager 10 10 USE lib_mpp 11 USE wrk_nemo_2 12 11 13 IMPLICIT NONE 12 14 PRIVATE … … 28 30 !! *** ROUTINE Agrif_Update_Tra *** 29 31 !!--------------------------------------------- 30 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released31 USE wrk_nemo, ONLY: wrk_4d_132 32 !! 33 33 INTEGER, INTENT(in) :: kt … … 37 37 IF((Agrif_NbStepint() .NE. (Agrif_irhot()-1)).AND.(kt /= 0)) RETURN 38 38 #if defined TWO_WAY 39 IF( wrk_in_use(4, 1) ) THEN 40 CALL ctl_stop('agrif_update_tra: ERROR: requested workspace arrays unavailable') 41 RETURN 42 END IF 43 ztab => wrk_4d_1 39 CALL wrk_alloc( jpi, jpj, jpk, jpts, ztab ) 44 40 45 41 Agrif_UseSpecialValueInUpdate = .TRUE. … … 54 50 Agrif_UseSpecialValueInUpdate = .FALSE. 55 51 56 IF( wrk_not_released(4, 1) ) THEN 57 CALL ctl_stop('Agrif_Update_Tra: ERROR: failed to release workspace arrays') 58 END IF 52 CALL wrk_dealloc( jpi, jpj, jpk, jpts, ztab ) 59 53 #endif 60 54 … … 65 59 !! *** ROUTINE Agrif_Update_Dyn *** 66 60 !!--------------------------------------------- 67 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released68 USE wrk_nemo, ONLY: wrk_2d_169 USE wrk_nemo, ONLY: wrk_3d_170 61 !! 71 62 INTEGER, INTENT(in) :: kt … … 76 67 IF ((Agrif_NbStepint() .NE. (Agrif_irhot()-1)).AND.(kt /= 0)) Return 77 68 #if defined TWO_WAY 78 ztab => wrk_3d_1 ; ztab2d => wrk_2d_1 79 IF( ( wrk_in_use(2, 1)) .OR. wrk_in_use(3, 1) )THEN 80 CALL ctl_stop('agrif_update_dyn: ERROR: requested workspace arrays unavailable') 81 RETURN 82 END IF 69 CALL wrk_alloc( jpi, jpj, ztab2d ) 70 CALL wrk_alloc( jpi, jpj, jpk, ztab ) 83 71 84 72 IF (mod(nbcline,nbclineupdate) == 0) THEN … … 100 88 Agrif_UseSpecialValueInUpdate = .FALSE. 101 89 102 IF( wrk_not_released(2, 1) .OR. wrk_not_released(3, 1) )THEN 103 CALL ctl_stop('agrif_update_dyn: ERROR: failed to release workspace arrays') 104 END IF 90 CALL wrk_dealloc( jpi, jpj, ztab2d ) 91 CALL wrk_dealloc( jpi, jpj, jpk, ztab ) 105 92 106 93 !Done in step
Note: See TracChangeset
for help on using the changeset viewer.