Changeset 2439 for trunk

Timestamp:

Dec 22, 2020, 11:51:50 AM (4 years ago)

Author:

emillour

Message:

Common dynamics:
Add a check that input fields (read from start.nc file) indeed have
matching values along the longitude end points when they should,
and if not correct field at longitude index iimp1 to match value
at longitude index 1.
EM

File:

• trunk/LMDZ.COMMON/libf/dyn3d_common/dynetat0.F90 (modified) (7 diffs)

trunk/LMDZ.COMMON/libf/dyn3d_common/dynetat0.F90

@@ -329 +329 @@
   CHARACTER(LEN=*), INTENT(IN)  :: var
   REAL,             INTENT(OUT) :: v(:,:)
+  INTEGER :: j
   CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
   CALL err(NF90_GET_VAR(fID,vID,v),"get",var)
+  ! Check that if on iip1 points along the longitude
+  ! values indeed perfectly match at end points
+  if (size(v,1)==iip1) then
+    do j=1,size(v,2)
+      if (v(1,j).ne.v(iip1,j)) then
+        ! tell the world
+        write(*,*) "Warning: longitudinal modulo discrepency for ",trim(var)
+        write(*,*) "j=",j," v(1,j)=",v(1,j)," v(iip1,j)=",v(iip1,j)
+        ! copy over value from v(1,j) to v(iip1,j)
+        write(*,*) "setting v(iip1,j)=v(1,j)"
+        v(iip1,j)=v(1,j)
+      endif ! of if (v(1,j).ne.v(iip1,j))
+    enddo ! of do j=1,size(v,2)
+  endif ! of if (size(v,1)==iip1)
 END SUBROUTINE get_var2
 
@@ -337 +352 @@
   REAL,             INTENT(OUT) :: v(:,:)
   INTEGER, INTENT(IN) :: indextime
+  INTEGER :: j
   corner(1)=1
   corner(2)=1
@@ -345 +361 @@
   CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
   CALL err(NF90_GET_VAR(fID,vID,v,corner,edges),"get",var)
+  ! Check that if on iip1 points along the longitude
+  ! values indeed perfectly match at end points
+  if (size(v,1)==iip1) then
+    do j=1,size(v,2)
+      if (v(1,j).ne.v(iip1,j)) then
+        ! tell the world
+        write(*,*) "Warning: longitudinal modulo discrepency for ",trim(var)
+        write(*,*) "j=",j," v(1,j)=",v(1,j)," v(iip1,j)=",v(iip1,j)
+        ! copy over value from v(1,j) to v(iip1,j)
+        write(*,*) "setting v(iip1,j)=v(1,j)"
+        v(iip1,j)=v(1,j)
+      endif ! of if (v(1,j).ne.v(iip1,j))
+    enddo ! of do j=1,size(v,2)
+  endif ! of if (size(v,1)==iip1)
 END SUBROUTINE get_var2_t
 
@@ -359 +389 @@
   REAL,             INTENT(OUT) :: v(:,:,:)
   INTEGER, INTENT(IN) :: indextime
+  INTEGER :: j,k
   corner(1)=1
   corner(2)=1
@@ -369 +400 @@
   CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
   CALL err(NF90_GET_VAR(fID,vID,v,corner,edges),"get",var)
+  ! Check that values indeed perfectly match at end points along longitude
+  do k=1,llm
+    do j=1,jjp1
+      if (v(1,j,k).ne.v(iip1,j,k)) then
+        ! tell the world
+        write(*,*) "Warning: longitudinal modulo discrepency for ",trim(var)
+        write(*,*) "j=",j," k=",k
+        write(*,*) " v(1,j,k)=",v(1,j,k)," v(iip1,j,k)=",v(iip1,j,k)
+        ! copy over value from v(1,j,k) to v(iip1,j,k)
+        write(*,*) "setting v(iip1,j,k)=v(1,j,k)"
+        v(iip1,j,k)=v(1,j,k)
+      endif ! of if (v(1,j,k).ne.v(iip1,j,k))
+    enddo ! of do j=1,jjp1
+  enddo ! of do k=1,llm
 END SUBROUTINE get_var3u_t
 
@@ -375 +420 @@
   REAL,             INTENT(OUT) :: v(:,:,:)
   INTEGER, INTENT(IN) :: indextime
+  INTEGER :: j,k
   corner(1)=1
   corner(2)=1
@@ -385 +431 @@
   CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
   CALL err(NF90_GET_VAR(fID,vID,v,corner,edges),"get",var)
+  ! Check that values indeed perfectly match at end points along longitude
+  do k=1,llm
+    do j=1,jjm
+      if (v(1,j,k).ne.v(iip1,j,k)) then
+        ! tell the world
+        write(*,*) "Warning: longitudinal modulo discrepency for ",trim(var)
+        write(*,*) "j=",j," k=",k
+        write(*,*) " v(1,j,k)=",v(1,j,k)," v(iip1,j,k)=",v(iip1,j,k)
+        ! copy over value from v(1,j,k) to v(iip1,j,k)
+        write(*,*) "setting v(iip1,j,k)=v(1,j,k)"
+        v(iip1,j,k)=v(1,j,k)
+      endif ! of if (v(1,j,k).ne.v(iip1,j,k))
+    enddo ! of do j=1,jjm
+  enddo ! of do k=1,llm
 END SUBROUTINE get_var3v_t