Changeset 10025 for utils/tools/NESTING/src/agrif_connect_topo.f90

Timestamp:

2018-08-02T15:25:27+02:00 (6 years ago)

Author:

clem

Message:

nesting tools are partly rewritten (mostly for create_coordinates and bathy) to get better functionality. Now you can use the nesting to either define an agrif zoom or a regional domain (for bdy purposes). Also, the nesting tools output a domain_cfg.nc that can be directly used in NEMO4 without the need of DOMAINcfg tool. The option of median average for bathymetry interpolation still does not work properly but it's not new

File:

• utils/tools/NESTING/src/agrif_connect_topo.f90 (modified) (1 diff)

utils/tools/NESTING/src/agrif_connect_topo.f90

@@ -51 +51 @@
     TYPE(Coordinates) :: Grid     
     !
-    !
-    diff = 0
-    IF(MOD(rho,2) .EQ. 0) diff = 1
-    !
-    indx = 1 + nbghostcellsfine + CEILING(irafx/2.0) + diff      
-    indy = 1 + nbghostcellsfine + CEILING(irafy/2.0) + diff
-    bornex = 1+nbghostcellsfine + CEILING(irafx/2.0) + diff - irafx
-    borney = 1+nbghostcellsfine + CEILING(irafy/2.0) + diff - irafy
-    bornex2 = nxfin - (nbghostcellsfine) + irafx - CEILING(irafx/2.0) 
-    borney2 = nyfin - (nbghostcellsfine) + irafy - CEILING(irafy/2.0) 
-    !    
-    ALLOCATE(bathy_fin_constant(bornex-FLOOR(irafx/2.0):bornex2+FLOOR(irafx/2.0), &
-         borney-FLOOR(irafy/2.0):borney2+FLOOR(irafy/2.0)))
-    !
-    DO j = borney,borney2,irafy
-
-       jpt = jmin + 1 + nbghostcellsfine + (j-indy)/irafy
-       IF(j<=1) jpt = jmin + 1
-
-       DO i = bornex,bornex2,irafx
-
-          ipt = imin + 1 + nbghostcellsfine + (i-indx)/irafx
-          IF(i<=1) ipt = imin + 1
-          !       
-          DO jj = j-FLOOR(irafy/2.0),j+FLOOR(irafy/2.0)-diff
-             DO ii = i-FLOOR(irafx/2.0),i+FLOOR(irafx/2.0)-diff
-
-                bathy_fin_constant(ii,jj) = coarse_bathy(ipt,jpt)
-
+    IF( ln_agrif_domain ) THEN
+    !
+       diff = 0
+       IF(MOD(rho,2) .EQ. 0) diff = 1
+       !
+       indx = 1 + nbghostcellsfine + CEILING(irafx/2.0) + diff      
+       indy = 1 + nbghostcellsfine + CEILING(irafy/2.0) + diff
+       bornex = 1+nbghostcellsfine + CEILING(irafx/2.0) + diff - irafx
+       borney = 1+nbghostcellsfine + CEILING(irafy/2.0) + diff - irafy
+       bornex2 = nxfin - (nbghostcellsfine) + irafx - CEILING(irafx/2.0) 
+       borney2 = nyfin - (nbghostcellsfine) + irafy - CEILING(irafy/2.0) 
+       !    
+       ALLOCATE(bathy_fin_constant(bornex-FLOOR(irafx/2.0):bornex2+FLOOR(irafx/2.0), &
+          borney-FLOOR(irafy/2.0):borney2+FLOOR(irafy/2.0)))
+       !
+       DO j = borney,borney2,irafy
+          
+          jpt = jmin + 1 + nbghostcellsfine + (j-indy)/irafy
+          IF(j<=1) jpt = jmin + 1
+          
+          DO i = bornex,bornex2,irafx
+             
+             ipt = imin + 1 + nbghostcellsfine + (i-indx)/irafx
+             IF(i<=1) ipt = imin + 1
+             !       
+             DO jj = j-FLOOR(irafy/2.0),j+FLOOR(irafy/2.0)-diff
+                DO ii = i-FLOOR(irafx/2.0),i+FLOOR(irafx/2.0)-diff
+                   
+                   bathy_fin_constant(ii,jj) = coarse_bathy(ipt,jpt)
+                   
+                END DO
              END DO
+             
           END DO
-
-       END DO
-    END DO
+       END DO
+       
+    ELSE
+
+       ALLOCATE(bathy_fin_constant(1:nxfin,1:nyfin))
+
+       DO j = 1,nyfin-irafy+1,irafy
+          jpt = jmin + FLOOR( (j - 1.) / irafy )
+          !
+          DO i = 1,nxfin-irafx+1,irafx
+             ipt = imin + FLOOR( (i - 1.) / irafx )
+             !
+             bathy_fin_constant(i:i+irafx-1,j:j+irafy-1) = coarse_bathy(ipt,jpt)
+             !
+          END DO
+       END DO
+ 
+    ENDIF
     !
     !