Changeset 13899 for NEMO/branches/2020/tickets_icb_1900/src/SWE/dynkeg.F90
- 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/SWE/dynkeg.F90
r12983 r13899 103 103 !!an45 to be ADDED : que cas C2 - "wet points only" il suffit de x2 le terme quadratic a la coast (nn_dynkeg_adv = 2) 104 104 CASE ( nkeg_C2_wpo ) !-- Standard scheme --! 105 DO_3D _01_01(1, jpkm1 )105 DO_3D( 0, 1, 0, 1, 1, jpkm1 ) 106 106 zu = ( puu(ji-1,jj ,jk,Kmm) * puu(ji-1,jj ,jk,Kmm) & 107 107 & + puu(ji ,jj ,jk,Kmm) * puu(ji ,jj ,jk,Kmm) ) * ( 2._wp - umask(ji-1,jj,jk) * umask(ji,jj,jk) ) … … 113 113 ! 114 114 CASE ( nkeg_C2 ) !-- Standard scheme --! 115 DO_3D _01_01(1, jpkm1 )115 DO_3D( 0, 1, 0, 1, 1, jpkm1 ) 116 116 zu = puu(ji-1,jj ,jk,Kmm) * puu(ji-1,jj ,jk,Kmm) & 117 117 & + puu(ji ,jj ,jk,Kmm) * puu(ji ,jj ,jk,Kmm) … … 121 121 END_3D 122 122 CASE ( nkeg_HW ) !-- Hollingsworth scheme --! 123 DO_3D _00_00(1, jpkm1 )123 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 124 124 zu = 8._wp * ( puu(ji-1,jj ,jk,Kmm) * puu(ji-1,jj ,jk,Kmm) & 125 125 & + puu(ji ,jj ,jk,Kmm) * puu(ji ,jj ,jk,Kmm) ) & … … 137 137 END SELECT 138 138 ! 139 DO_3D _00_00(1, jpkm1 )139 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 140 140 puu(ji,jj,jk,Krhs) = puu(ji,jj,jk,Krhs) - ( zhke(ji+1,jj ,jk) - zhke(ji,jj,jk) ) / e1u(ji,jj) 141 141 pvv(ji,jj,jk,Krhs) = pvv(ji,jj,jk,Krhs) - ( zhke(ji ,jj+1,jk) - zhke(ji,jj,jk) ) / e2v(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.