- Timestamp:
- 2011-02-18T13:49:27+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/OBS/obs_readmdt.F90
r2287 r2590 81 81 !! * Modules used 82 82 USE iom 83 83 USE wrk_nemo, ONLY: wrk_use, wrk_release 84 USE wrk_nemo, ONLY: z_mdt => wrk_2d_1, & ! Array to store the MDT values 85 mdtmask => wrk_2d_2 ! Array to store the mask for the MDT 86 !! 84 87 !! * Arguments 85 88 INTEGER, INTENT(IN) :: kslano ! Number of SLA Products … … 107 110 INTEGER :: i_stat 108 111 109 REAL(wp), DIMENSION(jpi,jpj) :: &110 & z_mdt, & ! Array to store the MDT values111 & mdtmask ! Array to store the mask for the MDT112 112 REAL(wp), DIMENSION(1) :: & 113 113 & zext, & … … 129 129 & igrdj 130 130 INTEGER :: nummdt 131 !!---------------------------------------------------------------------- 132 133 IF(.NOT. wrk_use(2, 1,2))THEN 134 CALL ctl_stop('obs_rea_mdt : requested workspace array unavailable.') 135 RETURN 136 END IF 131 137 132 138 IF(lwp)WRITE(numout,*) … … 234 240 END DO 235 241 242 IF(.NOT. wrk_release(2, 1,2))THEN 243 CALL ctl_stop('obs_rea_mdt : failed to release workspace arrays.') 244 END IF 245 236 246 END SUBROUTINE obs_rea_mdt 237 247 … … 256 266 !!---------------------------------------------------------------------- 257 267 !! * Modules used 258 268 USE wrk_nemo, ONLY: wrk_use, wrk_release 269 USE wrk_nemo, ONLY: zpromsk => wrk_2d_3 270 !! 259 271 !! * Arguments 260 272 REAL(wp), DIMENSION(jpi,jpj), INTENT(INOUT) :: & … … 270 282 REAL(wp) :: zcorr_bcketa 271 283 REAL(wp) :: zcorr 272 REAL(wp), DIMENSION(jpi,jpj) :: zpromsk273 284 INTEGER :: jj 274 285 INTEGER :: ji 275 286 CHARACTER(LEN=14), PARAMETER :: & 276 287 & cpname = 'obs_offset_mdt' 277 288 !!---------------------------------------------------------------------- 289 290 IF(.NOT. wrk_use(2, 3))THEN 291 CALL ctl_stop('obs_offset_mdt : requested workspace array unavailable.') 292 RETURN 293 END IF 294 278 295 ! Initialize the local mask, for domain projection 279 296 ! Also exclude mdt points which are set to missing data … … 341 358 342 359 360 IF(.NOT. wrk_release(2, 3))THEN 361 CALL ctl_stop('obs_offset_mdt : failed to release workspace array.') 362 END IF 363 343 364 END SUBROUTINE obs_offset_mdt 344 365
Note: See TracChangeset
for help on using the changeset viewer.