- 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/OCE/USR/usrdef_fmask.F90
r12377 r13899 58 58 !!---------------------------------------------------------------------- 59 59 ! 60 IF( TRIM( cd_cfg ) == "orca" ) THEN !== ORCA Configurations ==!60 IF( TRIM( cd_cfg ) == "orca" .OR. TRIM( cd_cfg ) == "ORCA" ) THEN !== ORCA Configurations ==! 61 61 ! 62 62 SELECT CASE ( kcfg ) … … 68 68 ! 69 69 IF(lwp) WRITE(numout,*) ' Gibraltar ' 70 ij0 = 101 ; ij1 = 101 ! Gibraltar strait : partial slip (pfmsk=0.5) 71 ii0 = 139 ; ii1 = 140 ; pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.5_wp 72 ij0 = 102 ; ij1 = 102 73 ii0 = 139 ; ii1 = 140 ; pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.5_wp 70 ij0 = 101 + nn_hls ; ij1 = 101 + nn_hls ! Gibraltar strait : partial slip (pfmsk=0.5) 71 ii0 = 139 + nn_hls - 1 ; ii1 = 140 + nn_hls - 1 72 pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.5_wp 73 ij0 = 102 + nn_hls ; ij1 = 102 + nn_hls 74 ii0 = 139 + nn_hls - 1 ; ii1 = 140 + nn_hls - 1 75 pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.5_wp 74 76 ! 75 77 IF(lwp) WRITE(numout,*) ' Bab el Mandeb ' 76 ij0 = 87 ; ij1 = 88 ! Bab el Mandeb : partial slip (pfmsk=1) 77 ii0 = 160 ; ii1 = 160 ; pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1._wp 78 ij0 = 88 ; ij1 = 88 79 ii0 = 159 ; ii1 = 159 ; pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1._wp 78 ij0 = 87 + nn_hls ; ij1 = 88 + nn_hls ! Bab el Mandeb : partial slip (pfmsk=1) 79 ii0 = 160 + nn_hls - 1 ; ii1 = 160 + nn_hls - 1 80 pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1._wp 81 ij0 = 88 + nn_hls ; ij1 = 88 + nn_hls 82 ii0 = 159 + nn_hls - 1 ; ii1 = 159 + nn_hls - 1 83 pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1._wp 80 84 ! 81 85 ! We keep this as an example but it is instable in this case … … 94 98 !!gm ! Currently these hard-wired indices relate to configuration with extend grid (jpjglo=332) 95 99 ! 96 isrow = 332 - jpjglo100 isrow = 332 - (Nj0glo + 1) ! was 332 - jpjglo -> jpjglo_old_version = Nj0glo + 1 97 101 ! 98 102 IF(lwp) WRITE(numout,*) 99 103 IF(lwp) WRITE(numout,*) ' orca_r1: increase friction near the following straits : ' 100 104 IF(lwp) WRITE(numout,*) ' Gibraltar ' 101 ii0 = 282 ; ii1 = 283 ! Gibraltar Strait 102 ij0 = 241 - isrow ; ij1 = 241 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 105 ii0 = 282 + nn_hls - 1 ; ii1 = 283 + nn_hls - 1 ! Gibraltar Strait 106 ij0 = 241 + nn_hls - isrow ; ij1 = 241 + nn_hls - isrow 107 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 103 108 ! 104 109 IF(lwp) WRITE(numout,*) ' Bhosporus ' 105 ii0 = 314 ; ii1 = 315 ! Bhosporus Strait 106 ij0 = 248 - isrow ; ij1 = 248 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 110 ii0 = 314 + nn_hls - 1 ; ii1 = 315 + nn_hls - 1 ! Bhosporus Strait 111 ij0 = 248 + nn_hls - isrow ; ij1 = 248 + nn_hls - isrow 112 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 107 113 ! 108 114 IF(lwp) WRITE(numout,*) ' Makassar (Top) ' 109 ii0 = 48 ; ii1 = 48 ! Makassar Strait (Top) 110 ij0 = 189 - isrow ; ij1 = 190 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 115 ii0 = 48 + nn_hls - 1 ; ii1 = 48 + nn_hls - 1 ! Makassar Strait (Top) 116 ij0 = 189 + nn_hls - isrow ; ij1 = 190 + nn_hls - isrow 117 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 111 118 ! 112 119 IF(lwp) WRITE(numout,*) ' Lombok ' 113 ii0 = 44 ; ii1 = 44 ! Lombok Strait 114 ij0 = 164 - isrow ; ij1 = 165 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 120 ii0 = 44 + nn_hls - 1 ; ii1 = 44 + nn_hls - 1 ! Lombok Strait 121 ij0 = 164 + nn_hls - isrow ; ij1 = 165 + nn_hls - isrow 122 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 115 123 ! 116 124 IF(lwp) WRITE(numout,*) ' Ombai ' 117 ii0 = 53 ; ii1 = 53 ! Ombai Strait 118 ij0 = 164 - isrow ; ij1 = 165 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 125 ii0 = 53 + nn_hls - 1 ; ii1 = 53 + nn_hls - 1 ! Ombai Strait 126 ij0 = 164 + nn_hls - isrow ; ij1 = 165 + nn_hls - isrow 127 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 119 128 ! 120 129 IF(lwp) WRITE(numout,*) ' Timor Passage ' 121 ii0 = 56 ; ii1 = 56 ! Timor Passage 122 ij0 = 164 - isrow ; ij1 = 165 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 130 ii0 = 56 + nn_hls - 1 ; ii1 = 56 + nn_hls - 1 ! Timor Passage 131 ij0 = 164 + nn_hls - isrow ; ij1 = 165 + nn_hls - isrow 132 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 123 133 ! 124 134 IF(lwp) WRITE(numout,*) ' West Halmahera ' 125 ii0 = 58 ; ii1 = 58 ! West Halmahera Strait 126 ij0 = 181 - isrow ; ij1 = 182 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 135 ii0 = 58 + nn_hls - 1 ; ii1 = 58 + nn_hls - 1 ! West Halmahera Strait 136 ij0 = 181 + nn_hls - isrow ; ij1 = 182 + nn_hls - isrow 137 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 127 138 ! 128 139 IF(lwp) WRITE(numout,*) ' East Halmahera ' 129 ii0 = 55 ; ii1 = 55 ! East Halmahera Strait 130 ij0 = 181 - isrow ; ij1 = 182 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 140 ii0 = 55 + nn_hls - 1 ; ii1 = 55 + nn_hls - 1 ! East Halmahera Strait 141 ij0 = 181 + nn_hls - isrow ; ij1 = 182 + nn_hls - isrow 142 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 131 143 ! 132 144 CASE DEFAULT
Note: See TracChangeset
for help on using the changeset viewer.