Changeset 9753 for utils/tools/NESTING/src/agrif_connect_topo.f90
- Timestamp:
- 2018-06-06T15:28:47+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
utils/tools/NESTING/src/agrif_connect_topo.f90
r9749 r9753 791 791 INTEGER :: i,j 792 792 INTEGER :: k1 793 INTEGER :: compt 793 INTEGER :: compt, compt_oce 794 794 REAL*8, POINTER, DIMENSION(:) :: gdepw,gdept,e3w,e3t 795 795 ! … … 906 906 G0%Bathy_meter(ipt,jpt) = 0. 907 907 compt = 0 908 compt_oce = 0 908 909 DO jj = j-FLOOR(irafy/2.0),j+FLOOR(irafy/2.0)-diff 909 910 DO ii = i-FLOOR(irafx/2.0),i+FLOOR(irafx/2.0)-diff 910 911 G0%Bathy_meter(ipt,jpt) = G0%Bathy_meter(ipt,jpt) + G1%Bathy_meter(ii,jj) 912 IF (G1%Bathy_meter(ii,jj)>0.) compt_oce = compt_oce + 1 911 913 compt = compt + 1 912 914 END DO 913 915 END DO 916 ! IF (compt_oce > 0.5*irafx*irafy) G0%Bathy_meter(ipt,jpt) = G0%Bathy_meter(ipt,jpt) / FLOAT(compt) 914 917 G0%Bathy_meter(ipt,jpt) = G0%Bathy_meter(ipt,jpt) / FLOAT(compt) 915 918 IF(G0%Bathy_meter(ipt,jpt).GT.0.) G0%Bathy_meter(ipt,jpt) = MAX(G0%Bathy_meter(ipt,jpt),zmin)
Note: See TracChangeset
for help on using the changeset viewer.