Changeset 3895 for LMDZ6/trunk/libf/dynphy_lonlat

Timestamp:

May 12, 2021, 9:58:07 AM (4 years ago)

Author:

Ehouarn Millour

Message:

Correction on the values of the longitude boundaries which were offset wrt scalar grid point longitudes.
EM

File:

• LMDZ6/trunk/libf/dynphy_lonlat/inigeomphy_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/dynphy_lonlat/inigeomphy_mod.F90

@@ -93 +93 @@
   ALLOCATE(boundslat_reg(jjm+1,2))
   
-  DO i=1,iim
-   boundslon_reg(i,east)=rlonu(i+1) 
-   boundslon_reg(i,west)=rlonu(i) 
+  ! specific handling of the -180 longitude scalar grid point boundaries
+  boundslon_reg(1,east)=rlonu(1)
+  boundslon_reg(1,west)=rlonu(iim)-2*PI
+  DO i=2,iim
+   boundslon_reg(i,east)=rlonu(i) 
+   boundslon_reg(i,west)=rlonu(i-1) 
   ENDDO
 
@@ -124 +127 @@
     cufi_glo(1) = cu(1)
     cvfi_glo(1) = cv(1)
-    boundslonfi_glo(1,north_east)=0
+    boundslonfi_glo(1,north_east)=PI
     boundslatfi_glo(1,north_east)=PI/2
-    boundslonfi_glo(1,north_west)=2*PI
+    boundslonfi_glo(1,north_west)=-PI
     boundslatfi_glo(1,north_west)=PI/2
-    boundslonfi_glo(1,south_west)=2*PI
+    boundslonfi_glo(1,south_west)=-PI
     boundslatfi_glo(1,south_west)=rlatv(1)
-    boundslonfi_glo(1,south_east)=0
+    boundslonfi_glo(1,south_east)=PI
     boundslatfi_glo(1,south_east)=rlatv(1)
     DO j=2,jjm
@@ -141 +144 @@
         boundslonfi_glo(k,north_east)=rlonu(i)
         boundslatfi_glo(k,north_east)=rlatv(j-1)
-        boundslonfi_glo(k,north_west)=rlonu(i+1)
+        if (i.eq.1) then
+          ! special case for the first longitude's west bound
+          boundslonfi_glo(k,north_west)=rlonu(iim)-2*PI
+          boundslonfi_glo(k,south_west)=rlonu(iim)-2*PI
+        else
+          boundslonfi_glo(k,north_west)=rlonu(i-1)
+          boundslonfi_glo(k,south_west)=rlonu(i-1)
+        endif
         boundslatfi_glo(k,north_west)=rlatv(j-1)
-        boundslonfi_glo(k,south_west)=rlonu(i+1)
         boundslatfi_glo(k,south_west)=rlatv(j)
         boundslonfi_glo(k,south_east)=rlonu(i)
@@ -154 +163 @@
     cufi_glo(klon_glo) = cu((iim+1)*jjm+1)
     cvfi_glo(klon_glo) = cv((iim+1)*jjm-iim)
-    boundslonfi_glo(klon_glo,north_east)= 0
+    boundslonfi_glo(klon_glo,north_east)= PI
     boundslatfi_glo(klon_glo,north_east)= rlatv(jjm)
-    boundslonfi_glo(klon_glo,north_west)= 2*PI
+    boundslonfi_glo(klon_glo,north_west)= -PI
     boundslatfi_glo(klon_glo,north_west)= rlatv(jjm)
-    boundslonfi_glo(klon_glo,south_west)= 2*PI
+    boundslonfi_glo(klon_glo,south_west)= -PI
     boundslatfi_glo(klon_glo,south_west)= -PI/2
-    boundslonfi_glo(klon_glo,south_east)= 0
+    boundslonfi_glo(klon_glo,south_east)= PI
     boundslatfi_glo(klon_glo,south_east)= -Pi/2