Changeset 10025 for utils/tools/NESTING/src/agrif_interpolation.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_interpolation.f90 (modified) (5 diffs)

utils/tools/NESTING/src/agrif_interpolation.f90

@@ -54 +54 @@
     ! 
     !
-    CALL agrif_base_interp2(tabin,tabout,imin,jmin,typevar)
-
+    IF( ln_agrif_domain ) THEN
+       CALL agrif_base_interp2(tabin,tabout,imin,jmin,typevar)
+    ELSE
+       CALL agrif_base_interp3(tabin,tabout,imin,jmin,typevar)
+    ENDIF
+    
     !
   END SUBROUTINE agrif_interp
   !
-  !********************************************************
-  !   subroutine agrif_base_interp2       *
-  !********************************************************
   !      
   SUBROUTINE agrif_base_interp2(tabin,tabout,i_min,j_min,typevar)
@@ -189 +190 @@
   END SUBROUTINE agrif_base_interp2
   !
-  !********************************************************
-  !   subroutine polint             *
-  !********************************************************
+  !      
+  SUBROUTINE agrif_base_interp3(tabin,tabout,i_min,j_min,typevar)
+    !
+    IMPLICIT NONE
+    REAL*8,DIMENSION(:,:) :: tabin,tabout
+    INTEGER :: i_min,j_min
+    CHARACTER(*) :: typevar   
+
+    INTEGER :: nxf,nyf,zx,zy
+    INTEGER :: ji,jj,jif,jjf,jic,jjc,jdecx,jdecy
+    REAL*8 :: Ax, Bx, Ay, By
+
+    nxf = SIZE(tabout,1)
+    nyf = SIZE(tabout,2)
+
+    SELECT CASE(typevar)
+    CASE('T')
+       IF(MOD(irafx,2)==1) THEN ! odd
+          zx = 1 ; zy = 1 ; jdecx = FLOOR(irafx/2.) ; jdecy = FLOOR(irafy/2.)
+       ELSE                     ! even
+          zx = 2 ; zy = 2 ; jdecx = FLOOR(irafx/2.) ; jdecy = FLOOR(irafy/2.)
+       ENDIF
+    CASE('U')
+       IF(MOD(irafx,2)==1) THEN ! odd
+          zx = 1 ; zy = 1 ; jdecx = irafx - 1 ; jdecy = FLOOR(irafy/2.)
+       ELSE                     ! even
+          zx = 1 ; zy = 2 ; jdecx = irafx - 1 ; jdecy = FLOOR(irafy/2.)
+       ENDIF
+    CASE('V')
+       IF(MOD(irafx,2)==1) THEN ! odd
+          zx = 1 ; zy = 1 ; jdecx = FLOOR(irafx/2.) ; jdecy = irafy - 1
+       ELSE                     ! even
+          zx = 2 ; zy = 1 ; jdecx = FLOOR(irafx/2.) ; jdecy = irafy - 1
+       ENDIF
+    CASE('F')
+       IF(MOD(irafx,2)==1) THEN ! odd
+          zx = 1 ; zy = 1 ; jdecx = irafx - 1 ; jdecy = irafy - 1
+       ELSE                     ! even
+          zx = 1 ; zy = 1 ; jdecx = irafx - 1 ; jdecy = irafy - 1
+       ENDIF
+    END SELECT
+
+    
+    DO jj = 1, nyf
+
+       jjf = jj - jdecy
+       jjc = j_min + FLOOR((jjf-1.) / irafy) 
+       
+       DO ji = 1, nxf
+          
+          jif = ji - jdecx 
+          jic = i_min + FLOOR((jif-1.) / irafx) 
+          
+          Bx = MOD( zx*jif-1, zx*irafx ) / REAL(zx*irafx)
+          By = MOD( zy*jjf-1, zy*irafy ) / REAL(zy*irafy)
+          Ax = 1. - Bx
+          Ay = 1. - By
+
+          tabout(ji,jj) = ( Bx * tabin(jic+1,jjc  ) + Ax * tabin(jic,jjc  ) ) * Ay + &
+             &            ( Bx * tabin(jic+1,jjc+1) + Ax * tabin(jic,jjc+1) ) * By
+       END DO
+    END DO
+    
+    !
+  END SUBROUTINE agrif_base_interp3
+
   !       
   SUBROUTINE polint(xin,valin,n,x,val)
@@ -238 +302 @@
   END SUBROUTINE polint
   !
-  !********************************************************
-  !   subroutine polcoe             *
-  !********************************************************
   !      
   SUBROUTINE polcoe(xin,valin,n,cof)
@@ -709 +770 @@
           CALL Read_Ncdf_var('time_steps',filename,tabtemp1DInt) 
           !       print *,'timedeph = ',tabtemp1DInt
-          CALL Write_Ncdf_var('time_steps','time_counter',Child_file,tabtemp1DInt)
+          CALL Write_Ncdf_var('time_steps','time_counter',Child_file,tabtemp1DInt,'integer')
           CALL Copy_Ncdf_att('time_steps',filename,Child_file) 
           DEALLOCATE(tabtemp1DInt)
@@ -984 +1045 @@
   END SUBROUTINE Interp_Extrap_var
   !
-  !**************************************************************
-  !   end subroutine Interp_var
-  !**************************************************************          
   ! 
 END MODULE agrif_interpolation