Changeset 2500

Timestamp:

Apr 20, 2021, 12:40:10 PM (4 years ago)

Author:

abierjon

Message:

Mars GCM:
Definition of mesh longitude boundaries in inigeomphy

• at poles : west bound is set as -PI, east bound is set as PI (to ensure

that westbound < eastbound and that longitude values are
within [-PI;PI] and not [0;2PI]everywhere on the globe)

• elsewhere : rectify the index for west bounds, and make a special case for

the west bound of the first longitude (-PI).
NB: thanks to the redundant point in longitude, there is no need
of such a special case for the east bound of the last longitude

AB

File:

• trunk/LMDZ.COMMON/libf/dynphy_lonlat/inigeomphy_mod.F90 (modified) (3 diffs)

trunk/LMDZ.COMMON/libf/dynphy_lonlat/inigeomphy_mod.F90

@@ -120 +120 @@
     cufi_glo(1) = cu(1)
     cvfi_glo(1) = cv(1)
-    boundslonfi_glo(1,north_east)=0
-    boundslatfi_glo(1,north_east)=PI/2
-    boundslonfi_glo(1,north_west)=2*PI
-    boundslatfi_glo(1,north_west)=PI/2
-    boundslonfi_glo(1,south_west)=2*PI
-    boundslatfi_glo(1,south_west)=rlatv(1)
-    boundslonfi_glo(1,south_east)=0
-    boundslatfi_glo(1,south_east)=rlatv(1)
+    boundslonfi_glo(1,north_east)= PI
+    boundslatfi_glo(1,north_east)= PI/2
+    boundslonfi_glo(1,north_west)= -PI
+    boundslatfi_glo(1,north_west)= PI/2
+    boundslonfi_glo(1,south_west)= -PI
+    boundslatfi_glo(1,south_west)= rlatv(1)
+    boundslonfi_glo(1,south_east)= PI
+    boundslatfi_glo(1,south_east)= rlatv(1)
     DO j=2,jjm
       DO i=1,iim
@@ -137 +137 @@
         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)
@@ -150 +156 @@
     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
-    boundslatfi_glo(klon_glo,south_east)= -Pi/2
+    boundslonfi_glo(klon_glo,south_east)= PI
+    boundslatfi_glo(klon_glo,south_east)= -PI/2
 
     ! build airefi(), mesh area on physics grid