Changeset 13710 for NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/TOP/TRP/trcdmp.F90
- Timestamp:
- 2020-11-02T10:56:42+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
src/TOP/TRP/trcdmp.F90 (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 8 9 9 # SETTE 10 ^/utils/CI/sette@ HEADsette10 ^/utils/CI/sette@13559 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/TOP/TRP/trcdmp.F90
r12377 r13710 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 … … 45 45 !! * Substitutions 46 46 # include "do_loop_substitute.h90" 47 # include "domzgr_substitute.h90" 47 48 !!---------------------------------------------------------------------- 48 49 !! NEMO/TOP 4.0 , NEMO Consortium (2018) … … 112 113 ! 113 114 CASE( 0 ) !== newtonian damping throughout the water column ==! 114 DO_3D _00_00(1, jpkm1 )115 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 115 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) ) 116 117 END_3D 117 118 ! 118 119 CASE ( 1 ) !== no damping in the turbocline (avt > 5 cm2/s) ==! 119 DO_3D _00_00(1, jpkm1 )120 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 120 121 IF( avt(ji,jj,jk) <= avt_c ) THEN 121 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) ) … … 124 125 ! 125 126 CASE ( 2 ) !== no damping in the mixed layer ==! 126 DO_3D _00_00(1, jpkm1 )127 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 127 128 IF( gdept(ji,jj,jk,Kmm) >= hmlp (ji,jj) ) THEN 128 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) ) … … 148 149 IF( sn_cfctl%l_prttrc ) THEN 149 150 WRITE(charout, FMT="('dmp ')") 150 CALL prt_ctl_ trc_info(charout)151 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' ) 152 153 ENDIF 153 154 ! … … 204 205 !Read in mask from file 205 206 CALL iom_open ( cn_resto_tr, imask) 206 CALL iom_get ( imask, jpdom_auto glo, 'resto', restotr)207 CALL iom_get ( imask, jpdom_auto, 'resto', restotr) 207 208 CALL iom_close( imask ) 208 209 ! … … 245 246 ! ! ======================= 246 247 CASE ( 1 ) ! eORCA_R1 configuration 247 ! ! ======================= 248 isrow = 332 - jpjglo 249 ! 250 nctsi1(1) = 333 ; nctsj1(1) = 243 - isrow ! Caspian Sea 251 nctsi2(1) = 342 ; nctsj2(1) = 274 - isrow 252 ! 253 nctsi1(2) = 198 ; nctsj1(2) = 258 - isrow ! Lake Superior 254 nctsi2(2) = 204 ; nctsj2(2) = 262 - isrow 255 ! 256 nctsi1(3) = 201 ; nctsj1(3) = 250 - isrow ! Lake Michigan 257 nctsi2(3) = 203 ; nctsj2(3) = 256 - isrow 258 ! 259 nctsi1(4) = 204 ; nctsj1(4) = 252 - isrow ! Lake Huron 260 nctsi2(4) = 209 ; nctsj2(4) = 256 - isrow 261 ! 262 nctsi1(5) = 206 ; nctsj1(5) = 249 - isrow ! Lake Erie 263 nctsi2(5) = 209 ; nctsj2(5) = 251 - isrow 264 ! 265 nctsi1(6) = 210 ; nctsj1(6) = 252 - isrow ! Lake Ontario 266 nctsi2(6) = 212 ; nctsj2(6) = 252 - isrow 267 ! 268 nctsi1(7) = 321 ; nctsj1(7) = 180 - isrow ! Victoria Lake 269 nctsi2(7) = 322 ; nctsj2(7) = 189 - isrow 270 ! 271 nctsi1(8) = 297 ; nctsj1(8) = 270 - isrow ! Baltic Sea 272 nctsi2(8) = 308 ; nctsj2(8) = 293 - isrow 273 ! 274 ! ! ======================= 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 ! ! ======================= 275 277 CASE ( 2 ) ! ORCA_R2 configuration 276 278 ! ! ======================= … … 285 287 nctsi2(3) = 181 ; nctsj2(3) = 112 286 288 ! 287 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 288 290 nctsi2(4) = 6 ; nctsj2(4) = 112 289 291 ! 290 292 nctsi1(5) = 145 ; nctsj1(5) = 116 ! Baltic Sea 291 293 nctsi2(5) = 150 ; nctsj2(5) = 126 294 ! 292 295 ! ! ======================= 293 296 CASE ( 4 ) ! ORCA_R4 configuration … … 305 308 nctsi1(4) = 75 ; nctsj1(4) = 59 ! Baltic Sea 306 309 nctsi2(4) = 76 ; nctsj2(4) = 61 310 ! 307 311 ! ! ======================= 308 312 CASE ( 025 ) ! ORCA_R025 configuration … … 318 322 ! 319 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 320 327 ! 321 328 ! convert the position in local domain indices
Note: See TracChangeset
for help on using the changeset viewer.
