Changeset 392 for trunk/NEMO/OPA_SRC/DYN
- Timestamp:
- 2006-03-10T11:55:45+01:00 (18 years ago)
- Location:
- trunk/NEMO/OPA_SRC/DYN
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/divcur.F90
r389 r392 130 130 IF( lp_obc_south ) hdivn(nis0 :nis1 ,njs0 :njs1 ,jk) = 0.e0 ! south 131 131 #endif 132 #if defined key_ AGRIF132 #if defined key_agrif 133 133 if ( .NOT. AGRIF_Root() ) then 134 134 IF ((nbondi == 1).OR.(nbondi == 2)) hdivn(nlci-1 , : ,jk) = 0.e0 ! east … … 334 334 IF( lp_obc_south ) hdivn(nis0 :nis1 ,njs0 :njs1 ,jk) = 0.e0 ! south 335 335 #endif 336 #if defined key_ AGRIF336 #if defined key_agrif 337 337 if ( .NOT. AGRIF_Root() ) then 338 338 IF ((nbondi == 1).OR.(nbondi == 2)) hdivn(nlci-1 , : ,jk) = 0.e0 ! east -
trunk/NEMO/OPA_SRC/DYN/dynnxt.F90
r389 r392 135 135 ! ! =============== 136 136 # endif 137 # if defined key_ AGRIF137 # if defined key_agrif 138 138 ! ! =============== 139 139 END DO ! End of slab … … 178 178 ENDIF 179 179 180 #if defined key_ AGRIF180 #if defined key_agrif 181 181 IF (.NOT.Agrif_Root()) CALL Agrif_Update_Dyn( kt ) 182 182 #endif -
trunk/NEMO/OPA_SRC/DYN/dynspg_flt.F90
r389 r392 171 171 CALL obc_vol( kt ) 172 172 #endif 173 #if defined key_ AGRIF173 #if defined key_agrif 174 174 ! Update velocities on each coarse/fine interfaces 175 175 … … 239 239 IF( nsolv == 4 ) CALL lbc_lnk_e( gcb, c_solver_pt, 1. ) 240 240 241 #if defined key_ AGRIF241 #if defined key_agrif 242 242 243 243 If (.NOT.AGRIF_ROOT()) THEN … … 315 315 END DO 316 316 317 #if defined key_ AGRIF317 #if defined key_agrif 318 318 IF (.NOT. Agrif_Root()) THEN 319 319 ! caution : grad D (fine) = grad D (coarse) at coarse/fine interface -
trunk/NEMO/OPA_SRC/DYN/dynspg_flt_jki.F90
r358 r392 37 37 USE prtctl ! Print control 38 38 USE in_out_manager ! I/O manager 39 USE agrif_opa_interp 39 40 40 41 IMPLICIT NONE … … 178 179 CALL obc_vol( kt ) 179 180 #endif 181 #if defined key_agrif 182 ! Update velocities on each coarse/fine interfaces 183 184 CALL Agrif_dyn( kt ) 185 186 #endif 180 187 #if defined key_orca_r2 181 188 IF( n_cla == 1 ) CALL dyn_spg_cla( kt ) ! Cross Land Advection (Update (ua,va)) … … 237 244 ! applied the lateral boundary conditions 238 245 IF( nsolv == 4) CALL lbc_lnk_e( gcb, c_solver_pt, 1. ) 246 247 #if defined key_agrif 248 249 If (.NOT.AGRIF_ROOT()) THEN 250 251 ! add contribution of gradient of after barotropic transport divergence 252 IF ((nbondi == -1).OR.(nbondi == 2)) gcb(3,:) = gcb(3,:) & 253 -znugdt * z2dt*laplacu(2,:)*gcdprc(3,:)*hu(2,:)*e2u(2,:) 254 IF ((nbondi == 1).OR.(nbondi == 2)) gcb(nlci-2,:) = gcb(nlci-2,:) & 255 +znugdt * z2dt*laplacu(nlci-2,:)*gcdprc(nlci-2,:)*hu(nlci-2,:)*e2u(nlci-2,:) 256 IF ((nbondj == -1).OR.(nbondj == 2)) gcb(:,3) = gcb(:,3) & 257 -znugdt * z2dt*laplacv(:,2)*gcdprc(:,3)*hv(:,2)*e1v(:,2) 258 IF ((nbondj == 1).OR.(nbondj == 2)) gcb(:,nlcj-2) = gcb(:,nlcj-2) & 259 +znugdt * z2dt*laplacv(:,nlcj-2)*gcdprc(:,nlcj-2)*hv(:,nlcj-2)*e1v(:,nlcj-2) 260 261 ENDIF 262 263 #endif 239 264 240 265 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,synchro,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, … … 310 335 END DO 311 336 337 #if defined key_agrif 338 IF (.NOT. Agrif_Root()) THEN 339 ! caution : grad D (fine) = grad D (coarse) at coarse/fine interface 340 IF ((nbondi == -1).OR.(nbondi == 2)) spgu(2,:) = znugdt * z2dt * laplacu(2,:) * umask(2,:,1) 341 IF ((nbondi == 1).OR.(nbondi == 2)) spgu(nlci-2,:) = znugdt * z2dt * laplacu(nlci-2,:) * umask(nlci-2,:,1) 342 IF ((nbondj == -1).OR.(nbondj == 2)) spgv(:,2) = znugdt * z2dt * laplacv(:,2) * vmask(:,2,1) 343 IF ((nbondj == 1).OR.(nbondj == 2)) spgv(:,nlcj-2) = znugdt * z2dt * laplacv(:,nlcj-2) * vmask(:,nlcj-2,1) 344 ENDIF 345 #endif 346 312 347 ! 7. Add the trends multiplied by z2dt to the after velocity 313 348 ! -----------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.