- Timestamp:
- 2011-03-15T16:27:46+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90
r2633 r2690 32 32 !! 33 33 !!---------------------------------------------------------------------- 34 !! * Arguments35 34 LOGICAL, INTENT (in) :: ld_print ! If true, output arrays on numout 36 37 !! * Local variables 38 INTEGER :: ji, jj 35 ! 36 INTEGER :: ji, jj 39 37 REAL(wp) :: za00, zd_max, zetmax, zeumax, zefmax, zevmax 40 38 !!---------------------------------------------------------------------- … … 43 41 IF(lwp) WRITE(numout,*) 'ldf_dyn_c2d : 2d lateral eddy viscosity coefficient' 44 42 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 45 IF(lwp) WRITE(numout,*)46 43 47 44 ! harmonic operator (ahm1, ahm2) : ( T- and F- points) (used for laplacian operators … … 123 120 ENDIF 124 121 ENDIF 125 126 122 ! 127 123 END SUBROUTINE ldf_dyn_c2d 128 124 … … 143 139 !! 144 140 !!---------------------------------------------------------------------- 145 !! * Modules used 146 USE ldftra_oce, ONLY : aht0 147 USE wrk_nemo, ONLY: iwrk_in_use, iwrk_not_released 148 USE wrk_nemo, ONLY: icof => iwrk_2d_1 149 !! * Arguments 141 USE ldftra_oce, ONLY: aht0 142 USE wrk_nemo , ONLY: iwrk_in_use, iwrk_not_released 143 USE wrk_nemo , ONLY: icof => iwrk_2d_1 144 ! 150 145 LOGICAL, INTENT (in) :: ld_print ! If true, output arrays on numout 151 152 !! * Local variables153 INTEGER :: ji, jj, jn ! dummy loop indices154 INTEGER :: inum ! temporary logical unit155 INTEGER :: iim, ijm156 INTEGER :: ifreq, il1, il2, ij, ii146 ! 147 INTEGER :: ji, jj, jn ! dummy loop indices 148 INTEGER :: inum, iim, ijm ! local integers 149 INTEGER :: ifreq, il1, il2, ij, ii 150 REAL(wp) :: zahmeq, zcoft, zcoff, zmsk 151 CHARACTER (len=15) :: clexp 157 152 INTEGER, DIMENSION(jpidta,jpidta) :: idata 158 159 REAL(wp) :: zahmeq, zcoft, zcoff, zmsk160 161 CHARACTER (len=15) :: clexp162 153 !!---------------------------------------------------------------------- 163 154 164 155 IF( iwrk_in_use(2, 1) )THEN 165 CALL ctl_stop('ldf_dyn_c2d_orca: ERROR: requested workspace array is unavailable.') 166 RETURN 167 END IF 156 CALL ctl_stop('ldf_dyn_c2d_orca: requested workspace array is unavailable') ; RETURN 157 ENDIF 168 158 169 159 IF(lwp) WRITE(numout,*) 170 160 IF(lwp) WRITE(numout,*) 'inildf: 2d eddy viscosity coefficient' 171 161 IF(lwp) WRITE(numout,*) '~~~~~~ --' 172 IF(lwp) WRITE(numout,*) 173 IF(lwp) WRITE(numout,*) ' orca ocean model' 174 IF(lwp) WRITE(numout,*) 162 IF(lwp) WRITE(numout,*) ' orca ocean configuration' 175 163 176 164 #if defined key_antarctic … … 293 281 ENDIF 294 282 295 IF( iwrk_not_released(2, 1) )THEN 296 CALL ctl_stop('ldf_dyn_c2d_orca: ERROR: failed to release workspace array.') 283 IF( iwrk_not_released(2, 1) ) CALL ctl_stop('ldf_dyn_c2d_orca: failed to release workspace array') 297 284 END IF 298 285 ! 299 286 END SUBROUTINE ldf_dyn_c2d_orca 287 300 288 301 289 SUBROUTINE ldf_dyn_c2d_orca_R1( ld_print ) … … 314 302 !! 315 303 !!---------------------------------------------------------------------- 316 !! * Modules used 317 USE ldftra_oce, ONLY : aht0 318 USE wrk_nemo, ONLY: iwrk_in_use, iwrk_not_released 319 USE wrk_nemo, ONLY: icof => iwrk_2d_1 320 321 !! * Arguments 304 USE ldftra_oce, ONLY: aht0 305 USE wrk_nemo , ONLY: iwrk_in_use, iwrk_not_released 306 USE wrk_nemo , ONLY: icof => iwrk_2d_1 307 ! 322 308 LOGICAL, INTENT (in) :: ld_print ! If true, output arrays on numout 323 324 !! * Local variables 309 ! 325 310 INTEGER :: ji, jj, jn ! dummy loop indices 326 311 INTEGER :: inum ! temporary logical unit 327 312 INTEGER :: iim, ijm 328 313 INTEGER :: ifreq, il1, il2, ij, ii 314 REAL(wp) :: zahmeq, zcoft, zcoff, zmsk, zam20s 315 CHARACTER (len=15) :: clexp 329 316 INTEGER, DIMENSION(jpidta,jpidta) :: idata 330 331 REAL(wp) :: zahmeq, zcoft, zcoff, zmsk, zam20s 332 333 CHARACTER (len=15) :: clexp 334 !!---------------------------------------------------------------------- 335 336 IF( iwrk_in_use(2, 1) )THEN 337 CALL ctl_stop('ldf_dyn_c2d_orca_R1: ERROR: requested workspace array is unavailable.') 338 RETURN 339 END IF 317 !!---------------------------------------------------------------------- 318 319 IF( iwrk_in_use(2, 1) ) THEN 320 CALL ctl_stop('ldf_dyn_c2d_orca_R1: requested workspace array is unavailable') ; RETURN 321 ENDIF 340 322 341 323 IF(lwp) WRITE(numout,*) 342 324 IF(lwp) WRITE(numout,*) 'inildf: 2d eddy viscosity coefficient' 343 325 IF(lwp) WRITE(numout,*) '~~~~~~ --' 344 IF(lwp) WRITE(numout,*) 345 IF(lwp) WRITE(numout,*) ' orca_r1 ocean model' 346 IF(lwp) WRITE(numout,*) 326 IF(lwp) WRITE(numout,*) ' orca_r1 configuration' 347 327 348 328 #if defined key_antarctic … … 472 452 ENDIF 473 453 474 IF( iwrk_not_released(2, 1) )THEN 475 CALL ctl_stop('ldf_dyn_c2d_orca_R1: ERROR: failed to release workspace array.') 476 END IF 477 454 IF( iwrk_not_released(2, 1) ) CALL ctl_stop('ldf_dyn_c2d_orca_R1: failed to release workspace array') 455 ! 478 456 END SUBROUTINE ldf_dyn_c2d_orca_R1
Note: See TracChangeset
for help on using the changeset viewer.