- Timestamp:
- 2020-11-27T17:26:33+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/tickets_icb_1900
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_icb_1900
- Property svn:externals
-
NEMO/branches/2020/tickets_icb_1900/src/TOP/TRP/trcdmp.F90
r13237 r13899 24 24 ! 25 25 USE iom 26 USE prtctl _trc! Print control for debbuging26 USE prtctl ! Print control for debbuging 27 27 28 28 IMPLICIT NONE … … 113 113 ! 114 114 CASE( 0 ) !== newtonian damping throughout the water column ==! 115 DO_3D _00_00(1, jpkm1 )115 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 116 116 ptr(ji,jj,jk,jn,Krhs) = ptr(ji,jj,jk,jn,Krhs) + restotr(ji,jj,jk) * ( ztrcdta(ji,jj,jk) - ptr(ji,jj,jk,jn,Kbb) ) 117 117 END_3D 118 118 ! 119 119 CASE ( 1 ) !== no damping in the turbocline (avt > 5 cm2/s) ==! 120 DO_3D _00_00(1, jpkm1 )120 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 121 121 IF( avt(ji,jj,jk) <= avt_c ) THEN 122 122 ptr(ji,jj,jk,jn,Krhs) = ptr(ji,jj,jk,jn,Krhs) + restotr(ji,jj,jk) * ( ztrcdta(ji,jj,jk) - ptr(ji,jj,jk,jn,Kbb) ) … … 125 125 ! 126 126 CASE ( 2 ) !== no damping in the mixed layer ==! 127 DO_3D _00_00(1, jpkm1 )127 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 128 128 IF( gdept(ji,jj,jk,Kmm) >= hmlp (ji,jj) ) THEN 129 129 ptr(ji,jj,jk,jn,Krhs) = ptr(ji,jj,jk,jn,Krhs) + restotr(ji,jj,jk) * ( ztrcdta(ji,jj,jk) - ptr(ji,jj,jk,jn,Kbb) ) … … 149 149 IF( sn_cfctl%l_prttrc ) THEN 150 150 WRITE(charout, FMT="('dmp ')") 151 CALL prt_ctl_ trc_info(charout)152 CALL prt_ctl _trc( tab4d=ptr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm, clinfo2='trd' )151 CALL prt_ctl_info( charout, cdcomp = 'top' ) 152 CALL prt_ctl( tab4d_1=ptr(:,:,:,:,Krhs), mask1=tmask, clinfo=ctrcnm, clinfo3='trd' ) 153 153 ENDIF 154 154 ! … … 205 205 !Read in mask from file 206 206 CALL iom_open ( cn_resto_tr, imask) 207 CALL iom_get ( imask, jpdom_auto glo, 'resto', restotr)207 CALL iom_get ( imask, jpdom_auto, 'resto', restotr) 208 208 CALL iom_close( imask ) 209 209 ! … … 246 246 ! ! ======================= 247 247 CASE ( 1 ) ! eORCA_R1 configuration 248 ! ! ======================= 249 isrow = 332 - jpjglo 250 ! 251 nctsi1(1) = 333 ; nctsj1(1) = 243 - isrow ! Caspian Sea 252 nctsi2(1) = 342 ; nctsj2(1) = 274 - isrow 253 ! 254 nctsi1(2) = 198 ; nctsj1(2) = 258 - isrow ! Lake Superior 255 nctsi2(2) = 204 ; nctsj2(2) = 262 - isrow 256 ! 257 nctsi1(3) = 201 ; nctsj1(3) = 250 - isrow ! Lake Michigan 258 nctsi2(3) = 203 ; nctsj2(3) = 256 - isrow 259 ! 260 nctsi1(4) = 204 ; nctsj1(4) = 252 - isrow ! Lake Huron 261 nctsi2(4) = 209 ; nctsj2(4) = 256 - isrow 262 ! 263 nctsi1(5) = 206 ; nctsj1(5) = 249 - isrow ! Lake Erie 264 nctsi2(5) = 209 ; nctsj2(5) = 251 - isrow 265 ! 266 nctsi1(6) = 210 ; nctsj1(6) = 252 - isrow ! Lake Ontario 267 nctsi2(6) = 212 ; nctsj2(6) = 252 - isrow 268 ! 269 nctsi1(7) = 321 ; nctsj1(7) = 180 - isrow ! Victoria Lake 270 nctsi2(7) = 322 ; nctsj2(7) = 189 - isrow 271 ! 272 nctsi1(8) = 297 ; nctsj1(8) = 270 - isrow ! Baltic Sea 273 nctsi2(8) = 308 ; nctsj2(8) = 293 - isrow 274 ! 275 ! ! ======================= 248 ! ! ======================= 249 ! 250 isrow = 332 - (Nj0glo + 1) ! was 332 - jpjglo -> jpjglo_old_version = Nj0glo + 1 251 ! 252 nctsi1(1) = 333 ; nctsj1(1) = 243 - isrow ! Caspian Sea 253 nctsi2(1) = 342 ; nctsj2(1) = 274 - isrow 254 ! 255 nctsi1(2) = 198 ; nctsj1(2) = 258 - isrow ! Lake Superior 256 nctsi2(2) = 204 ; nctsj2(2) = 262 - isrow 257 ! 258 nctsi1(3) = 201 ; nctsj1(3) = 250 - isrow ! Lake Michigan 259 nctsi2(3) = 203 ; nctsj2(3) = 256 - isrow 260 ! 261 nctsi1(4) = 204 ; nctsj1(4) = 252 - isrow ! Lake Huron 262 nctsi2(4) = 209 ; nctsj2(4) = 256 - isrow 263 ! 264 nctsi1(5) = 206 ; nctsj1(5) = 249 - isrow ! Lake Erie 265 nctsi2(5) = 209 ; nctsj2(5) = 251 - isrow 266 ! 267 nctsi1(6) = 210 ; nctsj1(6) = 252 - isrow ! Lake Ontario 268 nctsi2(6) = 212 ; nctsj2(6) = 252 - isrow 269 ! 270 nctsi1(7) = 321 ; nctsj1(7) = 180 - isrow ! Victoria Lake 271 nctsi2(7) = 322 ; nctsj2(7) = 189 - isrow 272 ! 273 nctsi1(8) = 297 ; nctsj1(8) = 270 - isrow ! Baltic Sea 274 nctsi2(8) = 308 ; nctsj2(8) = 293 - isrow 275 ! 276 ! ! ======================= 276 277 CASE ( 2 ) ! ORCA_R2 configuration 277 278 ! ! ======================= … … 286 287 nctsi2(3) = 181 ; nctsj2(3) = 112 287 288 ! 288 nctsi1(4) = 2 ; nctsj1(4) = 107 ! Black Sea 2 : est part of the Black Sea289 nctsi1(4) = 2 ; nctsj1(4) = 107 ! Black Sea 2 : est part of the Black Sea 289 290 nctsi2(4) = 6 ; nctsj2(4) = 112 290 291 ! 291 292 nctsi1(5) = 145 ; nctsj1(5) = 116 ! Baltic Sea 292 293 nctsi2(5) = 150 ; nctsj2(5) = 126 294 ! 293 295 ! ! ======================= 294 296 CASE ( 4 ) ! ORCA_R4 configuration … … 306 308 nctsi1(4) = 75 ; nctsj1(4) = 59 ! Baltic Sea 307 309 nctsi2(4) = 76 ; nctsj2(4) = 61 310 ! 308 311 ! ! ======================= 309 312 CASE ( 025 ) ! ORCA_R025 configuration … … 319 322 ! 320 323 ENDIF 324 ! 325 nctsi1(:) = nctsi1(:) + nn_hls - 1 ; nctsi2(:) = nctsi2(:) + nn_hls - 1 ! -1 as x-perio included in old input files 326 nctsj1(:) = nctsj1(:) + nn_hls ; nctsj2(:) = nctsj2(:) + nn_hls 321 327 ! 322 328 ! convert the position in local domain indices
Note: See TracChangeset
for help on using the changeset viewer.