65
66
74
75
76
77 integer, intent(in) :: ng, tile
78 integer, intent(in) :: LBi, UBi, LBj, UBj
79 integer, intent(in) :: IminS, ImaxS, JminS, JmaxS
80 integer, intent(in) :: krhs, kstp, kout
81
82# ifdef ASSUMED_SHAPE
83 real(r8), intent(in) :: ubar(LBi:,LBj:,:)
84 real(r8), intent(in) :: vbar(LBi:,LBj:,:)
85 real(r8), intent(in) :: zeta(LBi:,LBj:,:)
86
87 real(r8), intent(inout) :: ad_ubar(LBi:,LBj:,:)
88 real(r8), intent(inout) :: ad_vbar(LBi:,LBj:,:)
89 real(r8), intent(inout) :: ad_zeta(LBi:,LBj:,:)
90# else
91 real(r8), intent(in) :: ubar(LBi:UBi,LBj:UBj,:)
92 real(r8), intent(in) :: vbar(LBi:UBi,LBj:UBj,:)
93 real(r8), intent(in) :: zeta(LBi:UBi,LBj:UBj,:)
94
95 real(r8), intent(inout) :: ad_ubar(LBi:UBi,LBj:UBj,:)
96 real(r8), intent(inout) :: ad_vbar(LBi:UBi,LBj:UBj,:)
97 real(r8), intent(inout) :: ad_zeta(LBi:UBi,LBj:UBj,:)
98# endif
99
100
101
102 integer :: Imin, Imax
103 integer :: i, j, know
104
105 real(r8) :: Ce, Cx, Zx
106 real(r8) :: bry_pgr, bry_cor, bry_str, bry_val
107 real(r8) :: cff, cff1, cff2, cff3, dt2d
108 real(r8) :: obc_in, obc_out, tau
109# if defined ATM_PRESS && defined PRESS_COMPENSATE
110 real(r8) :: OneAtm, fac
111# endif
112
113 real(r8) :: ad_Ce, ad_Cx
114 real(r8) :: ad_bry_pgr, ad_bry_cor, ad_bry_str, ad_bry_val, ad_Zx
115 real(r8) :: ad_cff, ad_cff1, ad_cff2, ad_cff3
116 real(r8) :: adfac
117
118 real(r8), dimension(IminS:ImaxS,JminS:JmaxS) :: ad_grad
119
120# include "set_bounds.h"
121
122
123
124
125
126 ad_ce=0.0_r8
127 ad_cx=0.0_r8
128 ad_zx=0.0_r8
129 ad_cff=0.0_r8
130 ad_cff1=0.0_r8
131 ad_cff2=0.0_r8
132 ad_cff3=0.0_r8
133 ad_bry_pgr=0.0_r8
134 ad_bry_cor=0.0_r8
135 ad_bry_str=0.0_r8
136 ad_bry_val=0.0_r8
137
138 ad_grad(lbi:ubi,lbj:ubj)=0.0_r8
139
140
141
142
143
144 IF (first_2d_step) THEN
145 know=krhs
148 know=krhs
150 ELSE
151 know=kstp
153 END IF
154# if defined ATM_PRESS && defined PRESS_COMPENSATE
155 oneatm=1013.25_r8
156 fac=100.0_r8/(
g*
rho0)
157# endif
158
159# if defined WET_DRY_NOT_YET
160
161
162
163
164
165
166
168 IF (
domain(ng)%NorthEast_Corner(tile))
THEN
173
174
175
176
177
178
179 END IF
180 END IF
181 IF (
domain(ng)%NorthWest_Corner(tile))
THEN
186
187
188
189
190
191
192 END IF
193 END IF
194 IF (
domain(ng)%SouthEast_Corner(tile))
THEN
199
200
201
202
203
204
205 END IF
206 END IF
207 IF (
domain(ng)%SouthWest_Corner(tile))
THEN
212
213
214
215
216
217
218 END IF
219 END IF
220 END IF
221
223 IF (
domain(ng)%Northern_Edge(tile))
THEN
224 DO i=istr,iend
227
228
229
230
231
232
233 END IF
234 END DO
235 END IF
236 IF (
domain(ng)%Southern_Edge(tile))
THEN
237 DO i=istru,iend
240
241
242
243
244
245
246 END IF
247 END DO
248 END IF
249 END IF
250
252 IF (
domain(ng)%Eastern_Edge(tile))
THEN
253 DO j=jstr,jend
256
257
258
259
260
261
262 END IF
263 END DO
264 END IF
265 IF (
domain(ng)%Western_Edge(tile))
THEN
266 DO j=jstr,jend
269
270
271
272
273
274
275 END IF
276 END DO
277 END IF
278 END IF
279# endif
280
281
282
283
284
286 IF (
domain(ng)%NorthEast_Corner(tile))
THEN
289
290
291
292
293 adfac=0.5_r8*ad_ubar(iend+1,jend+1,kout)
294 ad_ubar(iend+1,jend ,kout)=ad_ubar(iend+1,jend ,kout)+ &
295 & adfac
296 ad_ubar(iend ,jend+1,kout)=ad_ubar(iend ,jend+1,kout)+ &
297 & adfac
298 ad_ubar(iend+1,jend+1,kout)=0.0_r8
299 END IF
300 END IF
301 IF (
domain(ng)%NorthWest_Corner(tile))
THEN
304
305
306
307
308 adfac=0.5_r8*ad_ubar(istr,jend+1,kout)
309 ad_ubar(istr ,jend ,kout)=ad_ubar(istr ,jend ,kout)+ &
310 & adfac
311 ad_ubar(istr+1,jend+1,kout)=ad_ubar(istr+1,jend+1,kout)+ &
312 & adfac
313 ad_ubar(istr ,jend+1,kout)=0.0_r8
314 END IF
315 END IF
316 IF (
domain(ng)%SouthEast_Corner(tile))
THEN
319
320
321
322
323 adfac=0.5_r8*ad_ubar(iend+1,jstr-1,kout)
324 ad_ubar(iend ,jstr-1,kout)=ad_ubar(iend ,jstr-1,kout)+ &
325 & adfac
326 ad_ubar(iend+1,jstr ,kout)=ad_ubar(iend+1,jstr ,kout)+ &
327 & adfac
328 ad_ubar(iend+1,jstr-1,kout)=0.0_r8
329 END IF
330 END IF
331 IF (
domain(ng)%SouthWest_Corner(tile))
THEN
334
335
336
337
338 adfac=0.5_r8*ad_ubar(istr,jstr-1,kout)
339 ad_ubar(istr+1,jstr-1,kout)=ad_ubar(istr+1,jstr-1,kout)+ &
340 & adfac
341 ad_ubar(istr ,jstr ,kout)=ad_ubar(istr ,jstr ,kout)+ &
342 & adfac
343 ad_ubar(istr ,jstr-1,kout)=0.0_r8
344 END IF
345 END IF
346 END IF
347
348
349
350
351
352 IF (
domain(ng)%Northern_Edge(tile))
THEN
353
354
355
357 IF (
iic(ng).ne.0)
THEN
358 DO i=istru,iend
360# if defined CELERITY_READ && defined FORWARD_READ
363 obc_out=0.5_r8* &
364 & (
clima(ng)%M2nudgcof(i-1,jend+1)+ &
365 &
clima(ng)%M2nudgcof(i ,jend+1))
366 obc_in =
obcfac(ng)*obc_out
367 ELSE
370 END IF
371 IF (
boundary(ng)%ubar_north_Ce(i).lt.0.0_r8)
THEN
372 tau=obc_in
373 ELSE
374 tau=obc_out
375 END IF
376 tau=tau*dt2d
377 END IF
378# ifdef RADIATION_2D
380# else
381 cx=0.0_r8
382# endif
385# endif
386# ifdef MASKING
387
388
389
390 ad_ubar(i,jend+1,kout)=ad_ubar(i,jend+1,kout)* &
391 &
grid(ng)%umask(i,jend+1)
392# endif
394
395
396
397 ad_ubar(i,jend+1 ,know)=ad_ubar(i,jend+1 ,know)- &
398 & tau*ad_ubar(i,jend+1,kout)
399 END IF
400
401
402
403
404
405
406
407
408 adfac=ad_ubar(i,jend+1,kout)/(cff+ce)
409 ad_grad(i-1,jend+1)=ad_grad(i-1,jend+1)- &
410 & max(cx,0.0_r8)*adfac
411 ad_grad(i ,jend+1)=ad_grad(i ,jend+1)- &
412 & min(cx,0.0_r8)*adfac
413 ad_ubar(i,jend ,kout)=ad_ubar(i,jend ,kout)+ce *adfac
414 ad_ubar(i,jend+1,know)=ad_ubar(i,jend+1,know)+cff*adfac
415 ad_ubar(i,jend+1,kout)=0.0_r8
416 END IF
417 END DO
418 END IF
419
420
421
425 DO i=istru,iend
427 cff=dt2d*0.5_r8*(
grid(ng)%pn(i-1,jend)+ &
428 &
grid(ng)%pn(i ,jend))
429 cff1=sqrt(
g*0.5_r8*(
grid(ng)%h(i-1,jend)+ &
430 & zeta(i-1,jend,know)+ &
431 &
grid(ng)%h(i ,jend)+ &
432 & zeta(i ,jend,know)))
433 ce=cff*cff1
434 cff2=1.0_r8/(1.0_r8+ce)
435# ifdef MASKING
436
437
438
439 ad_ubar(i,jend+1,kout)=ad_ubar(i,jend+1,kout)* &
440 &
grid(ng)%umask(i,jend+1)
441# endif
442
443
444
445
446
447
448 adfac=cff2*ad_ubar(i,jend+1,kout)
449 ad_ubar(i,jend ,kout)=ad_ubar(i,jend ,kout)+ce*adfac
450 ad_ubar(i,jend+1,know)=ad_ubar(i,jend+1,know)+adfac
451 ad_ce=ad_ce+ubar(i,jend,kout)*adfac
452 ad_cff2=ad_cff2+ &
453 & (ubar(i,jend+1,know)+ &
454 & ce*ubar(i,jend,kout))*ad_ubar(i,jend+1,kout)
455 ad_ubar(i,jend+1,kout)=0.0_r8
456
457
458 ad_ce=ad_ce-cff2*cff2*ad_cff2
459 ad_cff2=0.0_r8
460
461
462 ad_cff1=ad_cff1+cff*ad_ce
463 ad_ce=0.0_r8
464
465
466
467
468
469 adfac=0.25_r8*
g*ad_cff1/cff1
470 grid(ng)%ad_h(i-1,jend)=
grid(ng)%ad_h(i-1,jend)+adfac
471 grid(ng)%ad_h(i ,jend)=
grid(ng)%ad_h(i ,jend)+adfac
472 ad_zeta(i-1,jend,know)=ad_zeta(i-1,jend,know)+adfac
473 ad_zeta(i ,jend,know)=ad_zeta(i ,jend,know)+adfac
474 ad_cff1=0.0_r8
475 END IF
476 END DO
477
478
479
481 DO i=istru,iend
483# ifdef MASKING
484
485
486
487 ad_ubar(i,jend+1,kout)=ad_ubar(i,jend+1,kout)* &
488 &
grid(ng)%umask(i,jend+1)
489# endif
490# ifdef ADJUST_BOUNDARY
492
493
495 & ad_ubar_north(i)+ &
496 & ad_ubar(i,jend+1,kout)
497 ad_ubar(i,jend+1,kout)=0.0_r8
498 ELSE
499
500
501 ad_ubar(i,jend+1,kout)=0.0_r8
502 END IF
503# else
504
505
506 ad_ubar(i,jend+1,kout)=0.0_r8
507# endif
508 END IF
509 END DO
510
511
512
514 DO i=istru,iend
516# ifdef MASKING
517
518
519
520 ad_ubar(i,jend+1,kout)=ad_ubar(i,jend+1,kout)* &
521 &
grid(ng)%umask(i,jend+1)
522# endif
523
524
525 ad_ubar(i,jend,kout)=ad_ubar(i,jend,kout)+ &
526 & ad_ubar(i,jend+1,kout)
527 ad_ubar(i,jend+1,kout)=0.0_r8
528 END IF
529 END DO
530
531
532
533
536 imin=istru
537 imax=iend
538 ELSE
539 imin=istr
540 imax=iendr
541 END IF
542 DO i=imin,imax
544# ifdef MASKING
545
546
547
548 ad_ubar(i,jend+1,kout)=ad_ubar(i,jend+1,kout)* &
549 &
grid(ng)%umask(i,jend+1)
550# endif
551
552
553 ad_ubar(i,jend ,kout)=ad_ubar(i,jend,kout)+ &
554 &
gamma2(ng)*ad_ubar(i,jend+1,kout)
555 ad_ubar(i,jend+1,kout)=0.0_r8
556 END IF
557 END DO
558 END IF
559 END IF
560
561
562
563
564
565 IF (
domain(ng)%Southern_Edge(tile))
THEN
566
567
568
570 IF (
iic(ng).ne.0)
THEN
571 DO i=istru,iend
573# if defined CELERITY_READ && defined FORWARD_READ
576 obc_out=0.5_r8* &
577 & (
clima(ng)%M2nudgcof(i-1,jstr-1)+ &
578 &
clima(ng)%M2nudgcof(i ,jstr-1))
579 obc_in =
obcfac(ng)*obc_out
580 ELSE
583 END IF
584 IF (
boundary(ng)%ubar_south_Ce(i).lt.0.0_r8)
THEN
585 tau=obc_in
586 ELSE
587 tau=obc_out
588 END IF
589 tau=tau*dt2d
590 END IF
591# ifdef RADIATION_2D
593# else
594 cx=0.0_r8
595# endif
598# endif
599# ifdef MASKING
600
601
602
603 ad_ubar(i,jstr-1,kout)=ad_ubar(i,jstr-1,kout)* &
604 &
grid(ng)%umask(i,jstr-1)
605# endif
607
608
609
610 ad_ubar(i,jstr-1,know)=ad_ubar(i,jstr-1,know)- &
611 & tau*ad_ubar(i,jstr-1,kout)
612 END IF
613
614
615
616
617
618
619
620
621 adfac=ad_ubar(i,jstr-1,kout)/(cff+ce)
622 ad_grad(i-1,jstr-1)=ad_grad(i-1,jstr-1)- &
623 & max(cx,0.0_r8)*adfac
624 ad_grad(i ,jstr-1)=ad_grad(i ,jstr-1)- &
625 & min(cx,0.0_r8)*adfac
626 ad_ubar(i,jstr-1,know)=ad_ubar(i,jstr-1,know)+cff*adfac
627 ad_ubar(i,jstr ,kout)=ad_ubar(i,jstr ,kout)+ce *adfac
628 ad_ubar(i,jstr-1,kout)=0.0_r8
629 END IF
630 END DO
631 END IF
632
633
634
638 DO i=istru,iend
640 cff=dt2d*0.5_r8*(
grid(ng)%pn(i-1,jstr)+ &
641 &
grid(ng)%pn(i ,jstr))
642 cff1=sqrt(
g*0.5_r8*(
grid(ng)%h(i-1,jstr)+ &
643 & zeta(i-1,jstr,know)+ &
644 &
grid(ng)%h(i ,jstr)+ &
645 & zeta(i ,jstr,know)))
646 ce=cff*cff1
647 cff2=1.0_r8/(1.0_r8+ce)
648# ifdef MASKING
649
650
651
652 ad_ubar(i,jstr-1,kout)=ad_ubar(i,jstr-1,kout)* &
653 &
grid(ng)%umask(i,jstr-1)
654# endif
655
656
657
658
659
660
661 adfac=cff2*ad_ubar(i,jstr-1,kout)
662 ad_ubar(i,jstr-1,know)=ad_ubar(i,jstr-1,know)+adfac
663 ad_ubar(i,jstr ,kout)=ad_ubar(i,jstr ,kout)+ce*adfac
664 ad_ce=ad_ce+ubar(i,jstr,kout)*adfac
665 ad_cff2=ad_cff2+ &
666 & (ubar(i,jstr-1,know)+ &
667 & ce*ubar(i,jstr,kout))*ad_ubar(i,jstr-1,kout)
668 ad_ubar(i,jstr-1,kout)=0.0_r8
669
670
671 ad_ce=ad_ce-cff2*cff2*ad_cff2
672 ad_cff2=0.0_r8
673
674
675 ad_cff1=ad_cff1+cff*ad_ce
676 ad_ce=0.0_r8
677
678
679
680
681
682 adfac=0.25_r8*
g*ad_cff1/cff1
683 grid(ng)%ad_h(i-1,jstr)=
grid(ng)%ad_h(i-1,jstr)+adfac
684 grid(ng)%ad_h(i ,jstr)=
grid(ng)%ad_h(i ,jstr)+adfac
685 ad_zeta(i-1,jstr,know)=ad_zeta(i-1,jstr,know)+adfac
686 ad_zeta(i ,jstr,know)=ad_zeta(i ,jstr,know)+adfac
687 ad_cff1=0.0_r8
688 END IF
689 END DO
690
691
692
694 DO i=istru,iend
696# ifdef MASKING
697
698
699
700 ad_ubar(i,jstr-1,kout)=ad_ubar(i,jstr-1,kout)* &
701 &
grid(ng)%umask(i,jstr-1)
702# endif
703# ifdef ADJUST_BOUNDARY
705
706
708 & ad_ubar_south(i)+ &
709 & ad_ubar(i,jstr-1,kout)
710 ad_ubar(i,jstr-1,kout)=0.0_r8
711 ELSE
712
713
714 ad_ubar(i,jstr-1,kout)=0.0_r8
715 END IF
716# else
717
718
719 ad_ubar(i,jstr-1,kout)=0.0_r8
720# endif
721 END IF
722 END DO
723
724
725
727 DO i=istru,iend
729# ifdef MASKING
730
731
732
733 ad_ubar(i,jstr-1,kout)=ad_ubar(i,jstr-1,kout)* &
734 &
grid(ng)%umask(i,jstr-1)
735# endif
736
737
738 ad_ubar(i,jstr ,kout)=ad_ubar(i,jstr ,kout)+ &
739 & ad_ubar(i,jstr-1,kout)
740 ad_ubar(i,jstr-1,kout)=0.0_r8
741 END IF
742 END DO
743
744
745
746
749 imin=istru
750 imax=iend
751 ELSE
752 imin=istr
753 imax=iendr
754 END IF
755 DO i=imin,imax
757# ifdef MASKING
758
759
760
761 ad_ubar(i,jstr-1,kout)=ad_ubar(i,jstr-1,kout)* &
762 &
grid(ng)%umask(i,jstr-1)
763# endif
764
765
766 ad_ubar(i,jstr ,kout)=ad_ubar(i,jstr,kout)+ &
767 &
gamma2(ng)*ad_ubar(i,jstr-1,kout)
768 ad_ubar(i,jstr-1,kout)=0.0_r8
769 END IF
770 END DO
771 END IF
772 END IF
773
774
775
776
777
778 IF (
domain(ng)%Eastern_Edge(tile))
THEN
779
780
781
783 IF (
iic(ng).ne.0)
THEN
784 DO j=jstr,jend
786# if defined CELERITY_READ && defined FORWARD_READ
789 obc_out=0.5_r8* &
790 & (
clima(ng)%M2nudgcof(iend ,j)+ &
791 &
clima(ng)%M2nudgcof(iend+1,j))
792 obc_in =
obcfac(ng)*obc_out
793 ELSE
796 END IF
797 IF (
boundary(ng)%ubar_east_Cx(j).lt.0.0_r8)
THEN
798 tau=obc_in
799 ELSE
800 tau=obc_out
801 END IF
802 tau=tau*dt2d
803 END IF
805# ifdef RADIATION_2D
807# else
808 ce=0.0_r8
809# endif
811# endif
812# ifdef MASKING
813
814
815
816 ad_ubar(iend+1,j,kout)=ad_ubar(iend+1,j,kout)* &
817 &
grid(ng)%umask(iend+1,j)
818# endif
820
821
822
823 ad_ubar(iend+1,j,know)=ad_ubar(iend+1,j,know)- &
824 & tau*ad_ubar(iend+1,j,kout)
825 END IF
826
827
828
829
830
831
832
833
834 adfac=ad_ubar(iend+1,j,kout)/(cff+cx)
835 ad_grad(iend+1,j )=ad_grad(iend+1,j )- &
836 & max(ce,0.0_r8)*adfac
837 ad_grad(iend+1,j+1)=ad_grad(iend+1,j+1)- &
838 & min(ce,0.0_r8)*adfac
839 ad_ubar(iend ,j,kout)=ad_ubar(iend ,j,kout)+cx *adfac
840 ad_ubar(iend+1,j,know)=ad_ubar(iend+1,j,know)+cff*adfac
841 ad_ubar(iend+1,j,kout)=0.0_r8
842 END IF
843 END DO
844 END IF
845
846
847
849 DO j=jstr,jend
851 cff=1.0_r8/(0.5_r8*(
grid(ng)%h(iend ,j)+ &
852 & zeta(iend ,j,know)+ &
853 &
grid(ng)%h(iend+1,j)+ &
854 & zeta(iend+1,j,know)))
856# ifdef MASKING
857
858
859
860 ad_ubar(iend+1,j,kout)=ad_ubar(iend+1,j,kout)* &
861 &
grid(ng)%umask(iend+1,j)
862# endif
863# ifdef ADJUST_BOUNDARY
865
866
867
870 & cx*ad_ubar(iend+1,j,kout)
871 END IF
872# endif
873# if defined ATM_PRESS && defined PRESS_COMPENSATE
874
875
876
877
878
879
880
881
882
883
884
885
886
887 adfac=cx*0.5_r8*ad_ubar(iend+1,j,kout)
888 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+adfac
889 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
890 ad_cx=ad_cx+ &
891 & (0.5_r8*(zeta(iend ,j,know)+ &
892 & zeta(iend+1,j,know)+ &
893 & fac*(
forces(ng)%Pair(iend ,j)+ &
894 &
forces(ng)%Pair(iend+1,j)- &
895 & 2.0_r8*oneatm))- &
896 &
boundary(ng)%zeta_east(j))*ad_ubar(iend+1,j,kout)
897 ad_bry_val=ad_bry_val+ad_ubar(iend+1,j,kout)
898 ad_ubar(iend+1,j,kout)=0.0_r8
899# else
900
901
902
903
904
905
906
907
908
909 adfac=cx*0.5_r8*ad_ubar(iend+1,j,kout)
910 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+adfac
911 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
912 ad_cx=ad_cx+ &
913 & (0.5_r8*(zeta(iend ,j,know)+ &
914 & zeta(iend+1,j,know))- &
915 &
boundary(ng)%zeta_east(j))*ad_ubar(iend+1,j,kout)
916 ad_bry_val=ad_bry_val+ad_ubar(iend+1,j,kout)
917 ad_ubar(iend+1,j,kout)=0.0_r8
918# endif
919
920
921 ad_cff=ad_cff+0.5_r8*
g*ad_cx/cx
922 ad_cx=0.0_r8
923
924
925
926
927
928 adfac=-cff*cff*0.5_r8*ad_cff
929 grid(ng)%ad_h(iend ,j)=
grid(ng)%ad_h(iend ,j)+adfac
930 grid(ng)%ad_h(iend+1,j)=
grid(ng)%ad_h(iend+1,j)+adfac
931 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+adfac
932 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
933 ad_cff=0.0_r8
934
935# if defined SSH_TIDES_NOT_YET && !defined UV_TIDES_NOT_YET
938 & zeta(iend,j,know))* &
939 & 0.5_r8*
grid(ng)%pm(iend,j)
940 ELSE
941 bry_pgr=-
g*(zeta(iend+1,j,know)- &
942 & zeta(iend ,j,know))* &
943 & 0.5_r8*(
grid(ng)%pm(iend ,j)+ &
944 &
grid(ng)%pm(iend+1,j))
945 END IF
946# ifdef UV_COR
947 bry_cor=0.125_r8*(vbar(iend ,j ,know)+ &
948 & vbar(iend ,j+1,know)+ &
949 & vbar(iend+1,j ,know)+ &
950 & vbar(iend+1,j+1,know))* &
951 & (
grid(ng)%f(iend ,j)+ &
952 &
grid(ng)%f(iend+1,j))
953# else
954 bry_cor=0.0_r8
955# endif
956 cff1=1.0_r8/(0.5_r8*(
grid(ng)%h(iend ,j)+ &
957 & zeta(iend ,j,know)+ &
958 &
grid(ng)%h(iend+1,j)+ &
959 & zeta(iend+1,j,know)))
960 bry_str=cff1*(
forces(ng)%sustr(iend+1,j)- &
961 &
forces(ng)%bustr(iend+1,j))
962 cx=1.0_r8/sqrt(
g*0.5_r8*(
grid(ng)%h(iend+1,j)+ &
963 & zeta(iend+1,j,know)+ &
964 &
grid(ng)%h(iend ,j)+ &
965 & zeta(iend ,j,know)))
966 cff2=
grid(ng)%om_u(iend+1,j)*cx
967
968
969
970
971
972
973
974
975 adfac=cff2*ad_bry_val
976 ad_bry_pgr=ad_bry_pgr+adfac
977 ad_bry_cor=ad_bry_cor+adfac
978 ad_bry_str=ad_bry_str+adfac
979 ad_cff2=ad_cff2+(bry_pgr+ &
980 & bry_cor+ &
981 & bry_str)*ad_bry_val
982 ad_ubar(iend,j,know)=ad_ubar(iend,j,know)+ad_bry_val
983 ad_bry_val=0.0_r8
984
985
986 ad_cx=ad_cx+
grid(ng)%om_u(iend+1,j)*ad_cff2
987 ad_cff2=0.0_r8
988
989
990
991
992
993 adfac=-cx*cx*cx*0.25_r8*
g*ad_cx
994 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+adfac
995 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
996 grid(ng)%ad_h(iend ,j)=
grid(ng)%ad_h(iend ,j)+adfac
997 grid(ng)%ad_h(iend+1,j)=
grid(ng)%ad_h(iend+1,j)+adfac
998 ad_cx=0.0_r8
999
1000
1001
1002
1003
1004 adfac=cff1*ad_bry_str
1005 forces(ng)%ad_sustr(iend+1,j)= &
1006 &
forces(ng)%ad_sustr(iend+1,j)+ &
1007 & adfac
1008 forces(ng)%tl_bustr(iend+1,j)= &
1009 &
forces(ng)%tl_bustr(iend+1,j)- &
1010 & adfac
1011 ad_cff1=ad_cff1+(
forces(ng)%sustr(iend+1,j)- &
1012 &
forces(ng)%bustr(iend+1,j))*ad_bry_str
1013 ad_bry_str=0.0_r8
1014
1015
1016
1017
1018
1019 adfac=-cff1*cff1*0.5_r8*ad_cff1
1020 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+adfac
1021 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
1022 grid(ng)%ad_h(iend ,j)=
grid(ng)%ad_h(iend ,j)+adfac
1023 grid(ng)%ad_h(iend+1,j)=
grid(ng)%ad_h(iend+1,j)+adfac
1024 ad_cff1=0.0_r8
1025# ifdef UV_COR
1026
1027
1028
1029
1030
1031
1032
1033 adfac=0.125_r8*(
grid(ng)%f(iend ,j)+ &
1034 &
grid(ng)%f(iend+1,j))*ad_bry_cor
1035 ad_vbar(iend ,j ,know)=ad_vbar(iend ,j ,know)+adfac
1036 ad_vbar(iend+1,j ,know)=ad_vbar(iend+1,j ,know)+adfac
1037 ad_vbar(iend ,j+1,know)=ad_vbar(iend ,j+1,know)+adfac
1038 ad_vbar(iend+1,j+1,know)=ad_vbar(iend+1,j+1,know)+adfac
1039 ad_bry_cor=0.0_r8
1040# else
1041
1042
1043 ad_bry_cor=0.0_r8
1044# endif
1046# ifdef ADJUST_BOUNDARY
1048
1049
1050
1051
1053 & ad_zeta_east(j)- &
1055 &
grid(ng)%pm(iend,j)* &
1056 & ad_bry_pgr
1057 END IF
1058# endif
1059
1060
1061
1062 ad_zeta(iend,j,know)=ad_zeta(iend,j,know)+ &
1063 &
g*0.5_r8*
grid(ng)%pm(iend,j)* &
1064 & ad_bry_pgr
1065 ad_bry_pgr=0.0_r8
1066 ELSE
1067
1068
1069
1070
1071
1072 adfac=-
g*0.5_r8*(
grid(ng)%pm(iend ,j)+ &
1073 &
grid(ng)%pm(iend+1,j))*ad_bry_pgr
1074 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)-adfac
1075 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
1076 ad_bry_pgr=0.0_r8
1077 END IF
1078# else
1079# ifdef ADJUST_BOUNDARY
1081
1082
1085 & ad_bry_val
1086 ad_bry_val=0.0_r8
1087 ELSE
1088
1089
1090 ad_bry_val=0.0_r8
1091 END IF
1092# else
1093
1094
1095 ad_bry_val=0.0_r8
1096# endif
1097# endif
1098 END IF
1099 END DO
1100
1101
1102
1104 DO j=jstr,jend
1106# ifdef WET_DRY_NOT_YET
1107 cff=0.5_r8*(
grid(ng)%h(iend ,j)+ &
1108 & zeta(iend ,j,know)+ &
1109 &
grid(ng)%h(iend+1,j)+ &
1110 & zeta(iend+1,j,know))
1111# else
1112 cff=0.5_r8*(
grid(ng)%h(iend ,j)+ &
1113 &
grid(ng)%h(iend+1,j))
1114# endif
1116 cx=dt2d*cff1*cff*0.5_r8*(
grid(ng)%pm(iend ,j)+ &
1117 &
grid(ng)%pm(iend+1,j))
1118 zx=(0.5_r8+cx)*zeta(iend ,j,know)+ &
1119 & (0.5_r8-cx)*zeta(iend+1,j,know)
1121 cff2=(1.0_r8-
co/cx)**2
1122 cff3=zeta(iend,j,kout)+ &
1123 & cx*zeta(iend+1,j,know)- &
1124 & (1.0_r8+cx)*zeta(iend,j,know)
1125 zx=zx+cff2*cff3
1126 END IF
1127# ifdef MASKING
1128
1129
1130
1131 ad_ubar(iend+1,j,kout)=ad_ubar(iend+1,j,kout)* &
1132 &
grid(ng)%umask(iend+1,j)
1133# endif
1134# ifdef ADJUST_BOUNDARY
1136
1137
1138
1139
1142 & 0.5_r8*cff1* &
1143 & ad_ubar(iend+1,j,kout)
1144 END IF
1145# endif
1146
1147
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157 adfac=0.5_r8*ad_ubar(iend+1,j,kout)
1158 ad_ubar(iend+1,j,know)=ad_ubar(iend+1,j,know)+ &
1159 & (1.0_r8-cx)*adfac
1160 ad_ubar(iend ,j,know)=ad_ubar(iend ,j,know)+ &
1161 & cx*adfac
1162 ad_cx=ad_cx+ &
1163 & (ubar(iend ,j,know)- &
1164 & ubar(iend+1,j,know))*adfac
1165 ad_bry_val=ad_bry_val+adfac
1166 ad_cff1=ad_cff1+ &
1167 & (zx-
boundary(ng)%zeta_east(j))*adfac
1168 ad_zx=ad_zx-cff1*adfac
1169 ad_ubar(iend+1,j,kout)=0.0_r8
1171
1172
1173
1174 ad_cff2=ad_cff2+cff3*ad_zx
1175 ad_cff3=ad_cff3+cff2*ad_zx
1176
1177
1178
1179
1180
1181
1182 ad_zeta(iend ,j,kout)=ad_zeta(iend ,j,kout)+ &
1183 & ad_cff3
1184 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)- &
1185 & (1.0_r8+cx)*ad_cff3
1186 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+ &
1187 & cx*ad_cff3
1188 ad_cx=ad_cx+ &
1189 & (zeta(iend ,j,know)+ &
1190 & zeta(iend+1,j,know))*ad_cff3
1191 ad_cff3=0.0_r8
1192
1193
1194 ad_cx=ad_cx+ &
1195 & 2.0_r8*cff2*
co*ad_cff2/(cx*cx)
1196 ad_cff2=0.0_r8
1197 END IF
1198
1199
1200
1201
1202
1203 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+ &
1204 & (0.5_r8+cx)*ad_zx
1205 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+ &
1206 & (0.5_r8-cx)*ad_zx
1207 ad_cx=ad_cx+ &
1208 & (zeta(iend ,j,know)- &
1209 & zeta(iend+1,j,know))*ad_zx
1210 ad_zx=0.0_r8
1211
1212
1213
1214
1215
1216 adfac=dt2d*0.5_r8*(
grid(ng)%pm(iend ,j)+ &
1217 &
grid(ng)%pm(iend+1,j))*ad_cx
1218 ad_cff=ad_cff+cff1*adfac
1219 ad_cff1=ad_cff1+cff*adfac
1220 ad_cx=0.0_r8
1221
1222
1223 ad_cff=ad_cff- &
1224 & 0.5_r8*cff1*ad_cff1/cff
1225 ad_cff1=0.0_r8
1226# ifdef WET_DRY_NOT_YET
1227
1228
1229
1230
1231
1232 adfac=0.5_r8*ad_cff
1233 grid(ng)%ad_h(iend ,j)=
grid(ng)%ad_h(iend ,j)+adfac
1234 grid(ng)%ad_h(iend+1,j)=
grid(ng)%ad_h(iend+1,j)+adfac
1235 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+adfac
1236 ad_zeta(iend+1,j,know)=tl_zeta(iend+1,j,know)+adfac
1237 ad_cff=0.0_r8
1238# else
1239
1240
1241
1242 adfac=0.5_r8*ad_cff
1243 grid(ng)%ad_h(iend ,j)=
grid(ng)%ad_h(iend ,j)+adfac
1244 grid(ng)%ad_h(iend+1,j)=
grid(ng)%ad_h(iend+1,j)+adfac
1245 ad_cff=0.0_r8
1246# endif
1247# if defined SSH_TIDES_NOT_YET && !defined UV_TIDES_NOT_YET
1249 bry_pgr=-
g*(
boundary(ng)%zeta_east(j)- &
1250 & zeta(iend,j,know))* &
1251 & 0.5_r8*
grid(ng)%pm(iend,j)
1252 ELSE
1253 bry_pgr=-
g*(zeta(iend+1,j,know)- &
1254 & zeta(iend ,j,know))* &
1255 & 0.5_r8*(
grid(ng)%pm(iend ,j)+ &
1256 &
grid(ng)%pm(iend+1,j))
1257 END IF
1258# ifdef UV_COR
1259 bry_cor=0.125_r8*(vbar(iend ,j ,know)+ &
1260 & vbar(iend ,j+1,know)+ &
1261 & vbar(iend+1,j ,know)+ &
1262 & vbar(iend+1,j+1,know))* &
1263 & (
grid(ng)%f(iend ,j)+ &
1264 &
grid(ng)%f(iend+1,j))
1265# else
1266 bry_cor=0.0_r8
1267# endif
1268 cff1=1.0_r8/(0.5_r8*(
grid(ng)%h(iend ,j)+ &
1269 & zeta(iend ,j,know)+ &
1270 &
grid(ng)%h(iend+1,j)+ &
1271 & zeta(iend+1,j,know)))
1272 bry_str=cff1*(
forces(ng)%sustr(iend+1,j)- &
1273 &
forces(ng)%bustr(iend+1,j))
1274 cx=1.0_r8/sqrt(
g*0.5_r8*(
grid(ng)%h(iend+1,j)+ &
1275 & zeta(iend+1,j,know)+ &
1276 &
grid(ng)%h(iend ,j)+ &
1277 & zeta(iend ,j,know)))
1278 cff2=
grid(ng)%om_u(iend+1,j)*cx
1279
1280
1281
1282
1283
1284
1285
1286
1287 adfac=cff2*ad_bry_val
1288 ad_bry_pgr=ad_bry_pgr+adfac
1289 ad_bry_cor=ad_bry_cor+adfac
1290 ad_bry_str=ad_bry_str+adfac
1291 ad_cff2=ad_cff2+(bry_pgr+ &
1292 & bry_cor+ &
1293 & bry_str)*ad_bry_val
1294 ad_ubar(iend,j,know)=ad_ubar(iend,j,know)+ad_bry_val
1295 ad_bry_val=0.0_r8
1296
1297
1298 ad_cx=ad_cx+
grid(ng)%om_u(iend+1,j)*ad_cff2
1299 ad_cff2=0.0_r8
1300
1301
1302
1303
1304
1305 adfac=-cx*cx*cx*0.25_r8*
g*ad_cx
1306 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+adfac
1307 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
1308 grid(ng)%ad_h(iend ,j)=
grid(ng)%ad_h(iend ,j)+adfac
1309 grid(ng)%ad_h(iend+1,j)=
grid(ng)%ad_h(iend+1,j)+adfac
1310 ad_cx=0.0_r8
1311
1312
1313
1314
1315
1316 adfac=cff1*ad_bry_str
1317 forces(ng)%ad_sustr(iend+1,j)= &
1318 &
forces(ng)%ad_sustr(iend+1,j)+ &
1319 & adfac
1320 forces(ng)%tl_bustr(iend+1,j)= &
1321 &
forces(ng)%tl_bustr(iend+1,j)- &
1322 & adfac
1323 ad_cff1=ad_cff1+(
forces(ng)%sustr(iend+1,j)- &
1324 &
forces(ng)%bustr(iend+1,j))*ad_bry_str
1325 ad_bry_str=0.0_r8
1326
1327
1328
1329
1330
1331 adfac=-cff1*cff1*0.5_r8*ad_cff1
1332 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)+adfac
1333 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
1334 grid(ng)%ad_h(iend ,j)=
grid(ng)%ad_h(iend ,j)+adfac
1335 grid(ng)%ad_h(iend+1,j)=
grid(ng)%ad_h(iend+1,j)+adfac
1336 ad_cff1=0.0_r8
1337# ifdef UV_COR
1338
1339
1340
1341
1342
1343
1344
1345 adfac=0.125_r8*(
grid(ng)%f(iend ,j)+ &
1346 &
grid(ng)%f(iend+1,j))*ad_bry_cor
1347 ad_vbar(iend ,j ,know)=ad_vbar(iend ,j ,know)+adfac
1348 ad_vbar(iend+1,j ,know)=ad_vbar(iend+1,j ,know)+adfac
1349 ad_vbar(iend ,j+1,know)=ad_vbar(iend ,j+1,know)+adfac
1350 ad_vbar(iend+1,j+1,know)=ad_vbar(iend+1,j+1,know)+adfac
1351 ad_bry_cor=0.0_r8
1352# else
1353
1354
1355 ad_bry_cor=0.0_r8
1356# endif
1358# ifdef ADJUST_BOUNDARY
1360
1361
1362
1363
1365 & ad_zeta_east(j)- &
1367 &
grid(ng)%pm(iend,j)* &
1368 & ad_bry_pgr
1369 END IF
1370# endif
1371
1372
1373
1374 ad_zeta(iend,j,know)=ad_zeta(iend,j,know)+ &
1375 &
g*0.5_r8*
grid(ng)%pm(iend,j)* &
1376 & ad_bry_pgr
1377 ad_bry_pgr=0.0_r8
1378 ELSE
1379
1380
1381
1382
1383
1384 adfac=-
g*0.5_r8*(
grid(ng)%pm(iend ,j)+ &
1385 &
grid(ng)%pm(iend+1,j))*ad_bry_pgr
1386 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)-adfac
1387 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
1388 ad_bry_pgr=0.0_r8
1389 END IF
1390# else
1391# ifdef ADJUST_BOUNDARY
1393
1394
1397 & ad_bry_val
1398 ad_bry_val=0.0_r8
1399 ELSE
1400
1401
1402 ad_bry_val=0.0_r8
1403 END IF
1404# else
1405
1406
1407 ad_bry_val=0.0_r8
1408# endif
1409# endif
1410 END IF
1411 END DO
1412
1413
1414
1416 DO j=jstr,jend
1418# ifdef MASKING
1419
1420
1421
1422 ad_ubar(iend+1,j,kout)=ad_ubar(iend+1,j,kout)* &
1423 &
grid(ng)%umask(iend+1,j)
1424# endif
1425# ifdef ADJUST_BOUNDARY
1427
1428
1431 & ad_ubar(iend+1,j,kout)
1432 ad_ubar(iend+1,j,kout)=0.0_r8
1433 ELSE
1434
1435
1436 ad_ubar(iend+1,j,kout)=0.0_r8
1437 END IF
1438# else
1439
1440
1441 ad_ubar(iend+1,j,kout)=0.0_r8
1442# endif
1443 END IF
1444 END DO
1445
1446
1447
1449 DO j=jstr,jend
1451# ifdef MASKING
1452
1453
1454
1455 ad_ubar(iend+1,j,kout)=ad_ubar(iend+1,j,kout)* &
1456 &
grid(ng)%umask(iend+1,j)
1457# endif
1458
1459
1460 ad_ubar(iend,j,kout)=ad_ubar(iend,j,kout)+ &
1461 & ad_ubar(iend+1,j,kout)
1462 ad_ubar(iend+1,j,kout)=0.0_r8
1463 END IF
1464 END DO
1465
1466
1467
1469 DO j=jstr,jend
1471 cff=1.0_r8/(0.5_r8*(
grid(ng)%h(iend ,j)+ &
1472 & zeta(iend ,j,know)+ &
1473 &
grid(ng)%h(iend+1,j)+ &
1474 & zeta(iend+1,j,know)))
1475# ifdef MASKING
1476
1477
1478
1479 ad_ubar(iend+1,j,kout)=ad_ubar(iend+1,j,kout)* &
1480 &
grid(ng)%umask(iend+1,j)
1481# endif
1482
1483
1484
1485
1486
1487 adfac=dt2d*ad_ubar(iend+1,j,kout)
1488 ad_bry_pgr=ad_bry_pgr+adfac
1489 ad_bry_cor=ad_bry_cor+adfac
1490 ad_bry_str=ad_bry_str+adfac
1491 ad_ubar(iend+1,j,know)=ad_ubar(iend+1,j,know)+ &
1492 & ad_ubar(iend+1,j,kout)
1493 ad_ubar(iend+1,j,kout)=0.0_r8
1494
1495
1496
1497
1498
1499 adfac=cff*ad_bry_str
1500 forces(ng)%ad_sustr(iend+1,j)= &
1501 &
forces(ng)%ad_sustr(iend+1,j)+ &
1502 & adfac
1503 forces(ng)%ad_bustr(iend+1,j)= &
1504 &
forces(ng)%ad_bustr(iend+1,j)- &
1505 & adfac
1506 ad_cff=ad_cff+(
forces(ng)%sustr(iend+1,j)- &
1507 &
forces(ng)%bustr(iend+1,j))*ad_bry_str
1508 ad_bry_str=0.0_r8
1509# ifdef UV_COR
1510
1511
1512
1513
1514
1515
1516
1517 adfac=0.125_r8*(
grid(ng)%f(iend,j)+ &
1518 &
grid(ng)%f(iend+1,j))*ad_bry_cor
1519 ad_vbar(iend, j ,know)=ad_vbar(iend ,j ,know)+adfac
1520 ad_vbar(iend+1,j ,know)=ad_vbar(iend+1,j ,know)+adfac
1521 ad_vbar(iend ,j+1,know)=ad_vbar(iend ,j+1,know)+adfac
1522 ad_vbar(iend+1,j+1,know)=ad_vbar(iend+1,j+1,know)+adfac
1523 ad_bry_cor=0.0_r8
1524# else
1525
1526
1527 ad_bry_cor=0.0_r8
1528# endif
1530# ifdef ADJUST_BOUNDARY
1532
1533
1534
1535
1537 & ad_zeta_east(j)- &
1539 &
grid(ng)%pm(iend,j)* &
1540 & ad_bry_pgr
1541 END IF
1542# endif
1543
1544
1545
1546 ad_zeta(iend,j,know)=ad_zeta(iend,j,know)+ &
1547 &
g*0.5_r8*
grid(ng)%pm(iend,j)* &
1548 & ad_bry_pgr
1549 ad_bry_pgr=0.0_r8
1550 ELSE
1551
1552
1553
1554
1555
1556 adfac=-
g*0.5_r8*(
grid(ng)%pm(iend ,j)+ &
1557 &
grid(ng)%pm(iend+1,j))*ad_bry_pgr
1558 ad_zeta(iend ,j,know)=ad_zeta(iend ,j,know)-adfac
1559 ad_zeta(iend+1,j,know)=ad_zeta(iend+1,j,know)+adfac
1560 ad_bry_pgr=0.0_r8
1561 END IF
1562 END IF
1563 END DO
1564
1565
1566
1568 DO j=jstr,jend
1570
1571
1572 ad_ubar(iend+1,j,kout)=0.0_r8
1573 END IF
1574 END DO
1575 END IF
1576 END IF
1577
1578
1579
1580
1581
1582 IF (
domain(ng)%Western_Edge(tile))
THEN
1583
1584
1585
1587 IF (
iic(ng).ne.0)
THEN
1588 DO j=jstr,jend
1590# if defined CELERITY_READ && defined FORWARD_READ
1593 obc_out=0.5_r8* &
1594 & (
clima(ng)%M2nudgcof(istr-1,j)+ &
1595 &
clima(ng)%M2nudgcof(istr ,j))
1596 obc_in =
obcfac(ng)*obc_out
1597 ELSE
1600 END IF
1601 IF (
boundary(ng)%ubar_west_Cx(j).lt.0.0_r8)
THEN
1602 tau=obc_in
1603 ELSE
1604 tau=obc_out
1605 END IF
1606 tau=tau*dt2d
1607 END IF
1609# ifdef RADIATION_2D
1611# else
1612 ce=0.0_r8
1613# endif
1615# endif
1616# ifdef MASKING
1617
1618
1619
1620 ad_ubar(istr,j,kout)=ad_ubar(istr,j,kout)* &
1621 &
grid(ng)%umask(istr,j)
1622# endif
1624
1625
1626
1627 ad_ubar(istr,j,know)=ad_ubar(istr,j,know)- &
1628 & tau*ad_ubar(istr,j,kout)
1629 END IF
1630
1631
1632
1633
1634
1635
1636
1637
1638 adfac=ad_ubar(istr,j,kout)/(cff+cx)
1639 ad_grad(istr,j )=ad_grad(istr,j )-max(ce,0.0_r8)*adfac
1640 ad_grad(istr,j+1)=ad_grad(istr,j+1)-min(ce,0.0_r8)*adfac
1641 ad_ubar(istr ,j,know)=ad_ubar(istr ,j,know)+cff*adfac
1642 ad_ubar(istr+1,j,kout)=ad_ubar(istr+1,j,kout)+cx *adfac
1643 ad_ubar(istr ,j,kout)=0.0_r8
1644 END IF
1645 END DO
1646 END IF
1647
1648
1649
1651 DO j=jstr,jend
1653 cff=1.0_r8/(0.5_r8*(
grid(ng)%h(istr-1,j)+ &
1654 & zeta(istr-1,j,know)+ &
1655 &
grid(ng)%h(istr ,j)+ &
1656 & zeta(istr ,j,know)))
1658# ifdef MASKING
1659
1660
1661
1662 ad_ubar(istr,j,kout)=ad_ubar(istr,j,kout)* &
1663 &
grid(ng)%umask(istr,j)
1664# endif
1665# ifdef ADJUST_BOUNDARY
1667
1668
1669
1672 & cx*ad_ubar(istr,j,kout)
1673 END IF
1674# endif
1675# if defined ATM_PRESS && defined PRESS_COMPENSATE
1676
1677
1678
1679
1680
1681
1682
1683
1684
1685
1686
1687
1688
1689
1690 adfac=cx*0.5_r8*ad_ubar(istr,j,kout)
1691 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)-adfac
1692 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)-adfac
1693 ad_cx=ad_cx- &
1694 & (0.5_r8*(zeta(istr-1,j,know)+ &
1695 & zeta(istr ,j,know)+ &
1696 & fac*(
forces(ng)%Pair(istr-1,j)+ &
1697 &
forces(ng)%Pair(istr ,j)- &
1698 & 2.0_r8*oneatm))- &
1699 &
boundary(ng)%zeta_west(j))*ad_ubar(istr,j,kout)
1700 ad_bry_val=ad_bry_val+ad_ubar(istr,j,kout)
1701 ad_ubar(istr,j,kout)=0.0_r8
1702# else
1703
1704
1705
1706
1707
1708
1709
1710
1711
1712 adfac=cx*0.5_r8*ad_ubar(istr,j,kout)
1713 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)-adfac
1714 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)-adfac
1715 ad_cx=ad_cx- &
1716 & (0.5_r8*(zeta(istr-1,j,know)+ &
1717 & zeta(istr ,j,know))- &
1718 &
boundary(ng)%zeta_west(j))*ad_ubar(istr,j,kout)
1719 ad_bry_val=ad_bry_val+ad_ubar(istr,j,kout)
1720 ad_ubar(istr,j,kout)=0.0_r8
1721# endif
1722
1723
1724 ad_cff=ad_cff+0.5_r8*
g*ad_cx/cx
1725 ad_cx=0.0_r8
1726
1727
1728
1729
1730
1731 adfac=-cff*cff*0.5_r8*ad_cff
1732 grid(ng)%ad_h(istr-1,j)=
grid(ng)%ad_h(istr-1,j)+adfac
1733 grid(ng)%ad_h(istr ,j)=
grid(ng)%ad_h(istr ,j)+adfac
1734 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+adfac
1735 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
1736 ad_cff=0.0_r8
1737
1738# if defined SSH_TIDES_NOT_YET && !defined UV_TIDES_NOT_YET
1740 bry_pgr=-
g*(zeta(istr,j,know)- &
1742 & 0.5_r8*
grid(ng)%pm(istr,j)
1743 ELSE
1744 bry_pgr=-
g*(zeta(istr ,j,know)- &
1745 & zeta(istr-1,j,know))* &
1746 & 0.5_r8*(
grid(ng)%pm(istr-1,j)+ &
1747 &
grid(ng)%pm(istr,j))
1748 END IF
1749# ifdef UV_COR
1750 bry_cor=0.125_r8*(vbar(istr-1,j ,know)+ &
1751 & vbar(istr-1,j+1,know)+ &
1752 & vbar(istr ,j ,know)+ &
1753 & vbar(istr ,j+1,know))* &
1754 & (
grid(ng)%f(istr-1,j)+ &
1755 &
grid(ng)%f(istr,j))
1756# else
1757 bry_cor=0.0_r8
1758# endif
1759 cff1=1.0_r8/(0.5_r8*(
grid(ng)%h(istr-1,j)+ &
1760 & zeta(istr-1,j,know)+ &
1761 &
grid(ng)%h(istr ,j)+ &
1762 & zeta(istr ,j,know)))
1763 bry_str=cff1*(
forces(ng)%sustr(istr,j)- &
1764 &
forces(ng)%bustr(istr,j))
1765 cx=1.0_r8/sqrt(
g*0.5_r8*(
grid(ng)%h(istr-1,j)+ &
1766 & zeta(istr-1,j,know)+ &
1767 &
grid(ng)%h(istr ,j)+ &
1768 & zeta(istr ,j,know)))
1769 cff2=
grid(ng)%om_u(istr,j)*cx
1770
1771
1772
1773
1774
1775
1776
1777
1778 adfac=cff2*ad_bry_val
1779 ad_bry_pgr=ad_bry_pgr+adfac
1780 ad_bry_cor=ad_bry_cor+adfac
1781 ad_bry_str=ad_bry_str+adfac
1782 ad_cff2=ad_cff2+(bry_pgr+ &
1783 & bry_cor+ &
1784 & bry_str)*ad_bry_val
1785 ad_ubar(istr+1,j,know)=ad_ubar(istr+1,j,know)+ad_bry_val
1786 ad_bry_val=0.0_r8
1787
1788
1789 ad_cx=ad_cx+
grid(ng)%om_u(istr,j)*ad_cff2
1790 ad_cff2=0.0_r8
1791
1792
1793
1794
1795
1796 adfac=-cx*cx*cx*0.25_r8*
g*ad_cx
1797 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+adfac
1798 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
1799 grid(ng)%ad_h(istr-1,j)=
grid(ng)%ad_h(istr-1,j)+adfac
1800 grid(ng)%ad_h(istr ,j)=
grid(ng)%ad_h(istr ,j)+adfac
1801 ad_cx=0.0_r8
1802
1803
1804
1805
1806
1807 adfac=cff1*ad_bry_str
1808 forces(ng)%ad_sustr(istr,j)=
forces(ng)%ad_sustr(istr,j)+ &
1809 & adfac
1810 forces(ng)%ad_bustr(istr,j)=
forces(ng)%ad_bustr(istr,j)- &
1811 & adfac
1812 ad_cff1=ad_cff1+(
forces(ng)%sustr(istr,j)- &
1813 &
forces(ng)%bustr(istr,j))*ad_bry_str
1814 ad_bry_str=0.0_r8
1815
1816
1817
1818
1819
1820 adfac=-cff1*cff1*0.5_r8*ad_cff1
1821 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+adfac
1822 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
1823 grid(ng)%ad_h(istr-1,j)=
grid(ng)%ad_h(istr-1,j)+adfac
1824 grid(ng)%ad_h(istr ,j)=
grid(ng)%ad_h(istr ,j)+adfac
1825 ad_cff1=0.0_r8
1826# ifdef UV_COR
1827
1828
1829
1830
1831
1832
1833
1834 adfac=0.125_r8*(
grid(ng)%f(istr-1,j)+ &
1835 &
grid(ng)%f(istr ,j))*ad_bry_cor
1836 ad_vbar(istr-1,j ,know)=ad_vbar(istr-1,j ,know)+adfac
1837 ad_vbar(istr ,j ,know)=ad_vbar(istr ,j ,know)+adfac
1838 ad_vbar(istr-1,j+1,know)=ad_vbar(istr-1,j+1,know)+adfac
1839 ad_vbar(istr ,j+1,know)=ad_vbar(istr ,j+1,know)+adfac
1840 ad_bry_cor=0.0_r8
1841# else
1842
1843
1844 ad_bry_cor=0.0_r8
1845# endif
1847# ifdef ADJUST_BOUNDARY
1849
1850
1851
1852
1854 & ad_zeta_west(j)+ &
1856 &
grid(ng)%pm(istr,j)* &
1857 & ad_bry_pgr
1858 END IF
1859# endif
1860
1861
1862
1863 ad_zeta(istr,j,know)=ad_zeta(istr,j,know)- &
1864 &
g*0.5_r8*
grid(ng)%pm(istr,j)* &
1865 & ad_bry_pgr
1866 ad_bry_pgr=0.0_r8
1867 ELSE
1868
1869
1870
1871
1872
1873 adfac=-
g*0.5_r8*(
grid(ng)%pm(istr-1,j)+ &
1874 &
grid(ng)%pm(istr ,j))*ad_bry_pgr
1875 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)-adfac
1876 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
1877 ad_bry_pgr=0.0_r8
1878 END IF
1879# else
1880# ifdef ADJUST_BOUNDARY
1882
1883
1886 & ad_bry_val
1887 ad_bry_val=0.0_r8
1888 ELSE
1889
1890
1891 ad_bry_val=0.0_r8
1892 END IF
1893# else
1894
1895
1896 ad_bry_val=0.0_r8
1897# endif
1898# endif
1899 END IF
1900 END DO
1901
1902
1903
1905 DO j=jstr,jend
1907# ifdef WET_DRY_NOT_YET
1908 cff=0.5_r8*(
grid(ng)%h(istr-1,j)+ &
1909 & zeta(istr-1,j,know)+ &
1910 &
grid(ng)%h(istr ,j)+ &
1911 & zeta(istr ,j,know))
1912# else
1913 cff=0.5_r8*(
grid(ng)%h(istr-1,j)+ &
1914 &
grid(ng)%h(istr ,j))
1915# endif
1917 cx=dt2d*cff1*cff*0.5_r8*(
grid(ng)%pm(istr-1,j)+ &
1918 &
grid(ng)%pm(istr ,j))
1919 zx=(0.5_r8+cx)*zeta(istr ,j,know)+ &
1920 & (0.5_r8-cx)*zeta(istr-1,j,know)
1922 cff2=(1.0_r8-
co/cx)**2
1923 cff3=zeta(istr,j,kout)+ &
1924 & cx*zeta(istr-1,j,know)- &
1925 & (1.0_r8+cx)*zeta(istr,j,know)
1926 zx=zx+cff2*cff3
1927 END IF
1928# ifdef MASKING
1929
1930
1931
1932 ad_ubar(istr,j,kout)=ad_ubar(istr,j,kout)* &
1933 &
grid(ng)%umask(istr,j)
1934# endif
1935# ifdef ADJUST_BOUNDARY
1937
1938
1939
1942 & cx*ad_ubar(istr,j,kout)
1943 END IF
1944# endif
1945
1946
1947
1948
1949
1950
1951
1952
1953
1954
1955
1956 adfac=0.5_r8*ad_ubar(istr,j,kout)
1957 ad_ubar(istr, j,know)=ad_ubar(istr ,j,know)+ &
1958 & (1.0_r8-cx)*adfac
1959 ad_ubar(istr+1,j,know)=ad_ubar(istr+1,j,know)+ &
1960 & cx*adfac
1961 ad_cx=ad_cx- &
1962 (ubar(istr ,j,know)- &
1963 & ubar(istr+1,j,know))*adfac
1964 ad_bry_val=ad_bry_val+adfac
1965 ad_cff1=ad_cff1- &
1966 & (zx-
boundary(ng)%zeta_west(j))*adfac
1967 ad_zx=ad_zx-cff1*adfac
1968 ad_ubar(istr,j,kout)=0.0_r8
1970
1971
1972
1973 ad_cff2=ad_cff2+cff3*ad_zx
1974 ad_cff3=ad_cff3+cff2*ad_zx
1975
1976
1977
1978
1979
1980
1981 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+ &
1982 & cx*ad_cff3
1983 ad_zeta(istr ,j,kout)=ad_zeta(istr ,j,kout)+ &
1984 & ad_cff3
1985 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)- &
1986 & (1.0_r8+cx)*ad_cff3
1987 ad_cx=ad_cx+ &
1988 & (zeta(istr-1,j,know)+ &
1989 & zeta(istr ,j,know))*ad_cff3
1990 ad_cff3=0.0_r8
1991
1992
1993 ad_cx=ad_cx+ &
1994 & 2.0_r8*cff2*
co*ad_cff2/(cx*cx)
1995 ad_cff2=0.0_r8
1996 END IF
1997
1998
1999
2000
2001
2002 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+ &
2003 & (0.5_r8-cx)*ad_zx
2004 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+ &
2005 & (0.5_r8+cx)*ad_zx
2006 ad_cx=ad_cx+ &
2007 & (zeta(istr ,j,know)- &
2008 & zeta(istr-1,j,know))*ad_zx
2009 ad_zx=0.0_r8
2010
2011
2012
2013
2014
2015 adfac=dt2d*0.5_r8*(
grid(ng)%pm(istr-1,j)+ &
2016 &
grid(ng)%pm(istr ,j))*ad_cx
2017 ad_cff=ad_cff+cff1*adfac
2018 ad_cff1=ad_cff1+cff*adfac
2019 ad_cx=0.0_r8
2020
2021
2022 ad_cff=ad_cff- &
2023 & 0.5_r8*cff1*ad_cff1/cff
2024 ad_cff1=0.0_r8
2025# ifdef WET_DRY_NOT_YET
2026
2027
2028
2029
2030
2031 adfac=0.5_r8*ad_cff
2032 grid(ng)%ad_h(istr-1,j)=
grid(ng)%ad_h(istr-1,j)+adfac
2033 grid(ng)%ad_h(istr ,j)=
grid(ng)%ad_h(istr ,j)+adfac
2034 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+adfac
2035 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
2036 ad_cff=0.0_r8
2037# else
2038
2039
2040
2041 adfac=0.5_r8*ad_cff
2042 grid(ng)%ad_h(istr-1,j)=
grid(ng)%ad_h(istr-1,j)+adfac
2043 grid(ng)%ad_h(istr ,j)=
grid(ng)%ad_h(istr ,j)+adfac
2044 ad_cff=0.0_r8
2045# endif
2046# if defined SSH_TIDES_NOT_YET && !defined UV_TIDES_NOT_YET
2048 bry_pgr=-
g*(zeta(istr,j,know)- &
2050 & 0.5_r8*
grid(ng)%pm(istr,j)
2051 ELSE
2052 bry_pgr=-
g*(zeta(istr ,j,know)- &
2053 & zeta(istr-1,j,know))* &
2054 & 0.5_r8*(
grid(ng)%pm(istr-1,j)+ &
2055 &
grid(ng)%pm(istr ,j))
2056 END IF
2057# ifdef UV_COR
2058 bry_cor=0.125_r8*(vbar(istr-1,j ,know)+ &
2059 & vbar(istr-1,j+1,know)+ &
2060 & vbar(istr ,j ,know)+ &
2061 & vbar(istr ,j+1,know))* &
2062 & (
grid(ng)%f(istr-1,j)+ &
2063 &
grid(ng)%f(istr ,j))
2064# else
2065 bry_cor=0.0_r8
2066# endif
2067 cff1=1.0_r8/(0.5_r8*(
grid(ng)%h(istr-1,j)+ &
2068 & zeta(istr-1,j,know)+ &
2069 &
grid(ng)%h(istr ,j)+ &
2070 & zeta(istr ,j,know)))
2071 bry_str=cff1*(
forces(ng)%sustr(istr,j)- &
2072 &
forces(ng)%bustr(istr,j))
2073 cx=1.0_r8/sqrt(
g*0.5_r8*(
grid(ng)%h(istr-1,j)+ &
2074 & zeta(istr-1,j,know)+ &
2075 &
grid(ng)%h(istr ,j)+ &
2076 & zeta(istr ,j,know)))
2077 cff2=
grid(ng)%om_u(istr,j)*cx
2078
2079
2080
2081
2082
2083
2084
2085
2086 adfac=cff2*ad_bry_val
2087 ad_bry_pgr=ad_bry_pgr+adfac
2088 ad_bry_cor=ad_bry_cor+adfac
2089 ad_bry_str=ad_bry_str+adfac
2090 ad_cff2=ad_cff2+ &
2091 & (bry_pgr+ &
2092 & bry_cor+ &
2093 & bry_str)*ad_bry_val
2094 ad_ubar(istr+1,j,know)=ad_ubar(istr+1,j,know)+ &
2095 & ad_bry_val
2096 ad_bry_val=0.0_r8
2097
2098
2099 ad_cx=ad_cx+
grid(ng)%om_u(istr,j)*ad_cff2
2100 ad_cff2=0.0_r8
2101
2102
2103
2104
2105
2106 adfac=-cx*cx*cx*0.25_r8*
g*ad_cx
2107 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+adfac
2108 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
2109 grid(ng)%ad_h(istr-1,j)=
grid(ng)%ad_h(istr-1,j)+adfac
2110 grid(ng)%ad_h(istr ,j)=
grid(ng)%ad_h(istr ,j)+adfac
2111 ad_cx=0.0_r8
2112
2113
2114
2115
2116
2117 adfac=cff1*ad_bry_str
2118 forces(ng)%ad_sustr(istr,j)=
forces(ng)%ad_sustr(istr,j)+ &
2119 & adfac
2120 forces(ng)%ad_bustr(istr,j)=
forces(ng)%ad_bustr(istr,j)- &
2121 & adfac
2122 ad_cff1=ad_cff1+(
forces(ng)%sustr(istr,j)- &
2123 &
forces(ng)%bustr(istr,j))*ad_bry_str
2124 ad_bry_str=0.0_r8
2125
2126
2127
2128
2129
2130 adfac=-cff1*cff1*0.5_r8*ad_cff1
2131 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+adfac
2132 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
2133 grid(ng)%ad_h(istr-1,j)=
grid(ng)%ad_h(istr-1,j)+adfac
2134 grid(ng)%ad_h(istr ,j)=
grid(ng)%ad_h(istr ,j)+adfac
2135 ad_cff1=0.0_r8
2136# ifdef UV_COR
2137
2138
2139
2140
2141
2142
2143
2144 adfac=0.125_r8*(
grid(ng)%f(istr-1,j)+ &
2145 &
grid(ng)%f(istr ,j))*ad_bry_cor
2146 ad_vbar(istr-1,j ,know)=ad_vbar(istr-1,j ,know)+adfac
2147 ad_vbar(istr ,j ,know)=ad_vbar(istr ,j ,know)+adfac
2148 ad_vbar(istr-1,j+1,know)=ad_vbar(istr-1,j+1,know)+adfac
2149 ad_vbar(istr ,j+1,know)=ad_vbar(istr ,j+1,know)+adfac
2150 ad_bry_cor=0.0_r8
2151# else
2152
2153
2154 ad_bry_cor=0.0_r8
2155# endif
2157# ifdef ADJUST_BOUNDARY
2159
2160
2161
2162
2164 & ad_zeta_west(j)+ &
2166 &
grid(ng)%pm(istr,j)* &
2167 & ad_bry_pgr
2168 END IF
2169# endif
2170
2171
2172
2173 ad_zeta(istr,j,know)=ad_zeta(istr,j,know)- &
2174 &
g*0.5_r8*
grid(ng)%pm(istr,j)* &
2175 & ad_bry_pgr
2176 ad_bry_pgr=0.0_r8
2177 ELSE
2178
2179
2180
2181
2182
2183 adfac=-
g*0.5_r8*(
grid(ng)%pm(istr-1,j)+ &
2184 &
grid(ng)%pm(istr ,j))*ad_bry_pgr
2185 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)-adfac
2186 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
2187 ad_bry_pgr=0.0_r8
2188 END IF
2189# else
2190# ifdef ADJUST_BOUNDARY
2192
2193
2195 & ad_ubar_west(j)+ &
2196 & ad_bry_val
2197 ad_bry_val=0.0_r8
2198 ELSE
2199
2200
2201 ad_bry_val=0.0_r8
2202 END IF
2203# else
2204
2205
2206 ad_bry_val=0.0_r8
2207# endif
2208# endif
2209 END IF
2210 END DO
2211
2212
2213
2215 DO j=jstr,jend
2217# ifdef MASKING
2218
2219
2220
2221 ad_ubar(istr,j,kout)=ad_ubar(istr,j,kout)* &
2222 &
grid(ng)%umask(istr,j)
2223# endif
2224# ifdef ADJUST_BOUNDARY
2226
2227
2230 & ad_ubar(istr,j,kout)
2231 ad_ubar(istr,j,kout)=0.0_r8
2232 ELSE
2233
2234
2235 ad_ubar(istr,j,kout)=0.0_r8
2236 END IF
2237# else
2238
2239
2240 ad_ubar(istr,j,kout)=0.0_r8
2241# endif
2242 END IF
2243 END DO
2244
2245
2246
2248 DO j=jstr,jend
2250# ifdef MASKING
2251
2252
2253
2254 ad_ubar(istr,j,kout)=ad_ubar(istr,j,kout)* &
2255 &
grid(ng)%umask(istr,j)
2256# endif
2257
2258
2259 ad_ubar(istr+1,j,kout)=ad_ubar(istr+1,j,kout)+ &
2260 & ad_ubar(istr,j,kout)
2261 ad_ubar(istr ,j,kout)=0.0_r8
2262 END IF
2263 END DO
2264
2265
2266
2268 DO j=jstr,jend
2270 cff=1.0_r8/(0.5_r8*(
grid(ng)%h(iend ,j)+ &
2271 & zeta(iend ,j,know)+ &
2272 &
grid(ng)%h(iend+1,j)+ &
2273 & zeta(iend+1,j,know)))
2274# ifdef MASKING
2275
2276
2277
2278 ad_ubar(istr,j,kout)=ad_ubar(istr,j,kout)* &
2279 &
grid(ng)%umask(istr,j)
2280# endif
2281
2282
2283
2284
2285
2286 adfac=dt2d*ad_ubar(istr,j,kout)
2287 ad_bry_pgr=ad_bry_pgr+adfac
2288 ad_bry_cor=ad_bry_cor+adfac
2289 ad_bry_str=ad_bry_str+adfac
2290 ad_ubar(istr,j,know)=ad_ubar(istr,j,know)+ &
2291 & ad_ubar(istr,j,kout)
2292 ad_ubar(istr,j,kout)=0.0_r8
2293
2294
2295
2296
2297
2298 adfac=cff*ad_bry_str
2299 forces(ng)%ad_sustr(istr,j)=
forces(ng)%ad_sustr(istr,j)+ &
2300 & adfac
2301 forces(ng)%ad_bustr(istr,j)=
forces(ng)%ad_bustr(istr,j)- &
2302 & adfac
2303 ad_cff=ad_cff+(
forces(ng)%sustr(istr,j)- &
2304 &
forces(ng)%bustr(istr,j))*ad_bry_str
2305 ad_bry_str=0.0_r8
2306
2307
2308
2309
2310
2311 adfac=-cff*cff*0.5_r8*ad_cff
2312 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)+adfac
2313 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
2314 grid(ng)%ad_h(istr-1,j)=
grid(ng)%ad_h(istr-1,j)+adfac
2315 grid(ng)%ad_h(istr ,j)=
grid(ng)%ad_h(istr ,j)+adfac
2316 ad_cff=0.0_r8
2317# ifdef UV_COR
2318
2319
2320
2321
2322
2323
2324
2325 adfac=0.125_r8*(
grid(ng)%f(istr-1,j)+ &
2326 &
grid(ng)%f(istr ,j))*ad_bry_cor
2327 ad_vbar(istr-1,j ,know)=ad_vbar(istr-1,j ,know)+adfac
2328 ad_vbar(istr ,j ,know)=ad_vbar(istr ,j ,know)+adfac
2329 ad_vbar(istr-1,j+1,know)=ad_vbar(istr-1,j+1,know)+adfac
2330 ad_vbar(istr ,j+1,know)=ad_vbar(istr ,j+1,know)+adfac
2331 ad_bry_cor=0.0_r8
2332# else
2333
2334
2335 ad_bry_cor=0.0_r8
2336# endif
2338# ifdef ADJUST_BOUNDARY
2340
2341
2342
2343
2345 & ad_zeta_west(j)+ &
2347 &
grid(ng)%pm(istr,j)* &
2348 & ad_bry_pgr
2349 END IF
2350# endif
2351
2352
2353
2354 ad_zeta(istr,j,know)=ad_zeta(istr,j,know)- &
2355 &
g*0.5_r8*
grid(ng)%pm(istr,j)* &
2356 & ad_bry_pgr
2357 ad_bry_pgr=0.0_r8
2358 ELSE
2359
2360
2361
2362
2363
2364 adfac=-
g*0.5_r8*(
grid(ng)%pm(istr-1,j)+ &
2365 &
grid(ng)%pm(istr ,j))*ad_bry_pgr
2366 ad_zeta(istr-1,j,know)=ad_zeta(istr-1,j,know)-adfac
2367 ad_zeta(istr ,j,know)=ad_zeta(istr ,j,know)+adfac
2368 ad_bry_pgr=0.0_r8
2369 END IF
2370 END IF
2371 END DO
2372
2373
2374
2376 DO j=jstr,jend
2378
2379
2380 ad_ubar(istr,j,kout)=0.0_r8
2381 END IF
2382 END DO
2383 END IF
2384 END IF
2385
2386 RETURN
type(t_boundary), dimension(:), allocatable boundary
type(t_apply), dimension(:), allocatable lbc_apply
type(t_clima), dimension(:), allocatable clima
type(t_forces), dimension(:), allocatable forces
type(t_grid), dimension(:), allocatable grid
type(t_lbc), dimension(:,:,:), allocatable ad_lbc
type(t_lbc), dimension(:,:,:), allocatable lbc
type(t_domain), dimension(:), allocatable domain
logical, dimension(:), allocatable lnudgem2clm
integer, dimension(:), allocatable iic
logical, dimension(:,:,:), allocatable lobc
logical, dimension(:), allocatable ewperiodic
logical, dimension(:), allocatable nsperiodic
logical, dimension(:), allocatable predictor_2d_step
real(dp), dimension(:), allocatable obcfac
real(r8), dimension(:), allocatable gamma2
integer, parameter isouth
real(dp), dimension(:), allocatable dtfast
real(dp), dimension(:,:), allocatable m2obc_out
integer, parameter inorth
real(dp), dimension(:,:), allocatable m2obc_in