Changeset 2540 for trunk/LMDZ.MARS/util/simu_MCS.F90

Timestamp:

Jun 29, 2021, 10:02:59 PM (3 years ago)

Author:

abierjon

Message:

Mars GCM:
Make simu_MCS.F90 compatible with DYNAMICO lon-lat output files, with a check on the latitude array order when doing the interpolation

AB

File:

• trunk/LMDZ.MARS/util/simu_MCS.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/util/simu_MCS.F90

@@ -693 +693 @@
 endif
 
+!write(*,*) "gcm_vars retrieved : ",gcm_vars(1:nbvar)
+
 !===============================================================================
 ! 2.2 Definition of LT variables from obsfile to outfile
@@ -1624 +1626 @@
 
 if (lat.ne.prev_lat) then
-  ! recall that latitudes start from north pole
-  do i=1,latlen-1
-    if (latitude(i).ge.lat) then
-      ilat_inf=i
-    else
-      exit
-    endif
-  enddo
+  if (latitude(1).gt.latitude(2)) then
+  ! decreasing latitudes, from 90N to (-)90S (LMDZ-like)
+    do i=1,latlen-1
+      if (latitude(i).ge.lat) then
+        ilat_inf=i
+      else
+        exit
+      endif
+    enddo
+  else
+  ! increasing latitudes, from (-)90S to 90N (DYNAMICO-like)
+    do i=1,latlen-1
+      if (latitude(i).le.lat) then
+        ilat_inf=i
+      else
+        exit
+      endif
+    enddo
+  endif
   ilat_sup=ilat_inf+1
 endif ! of if (lat.ne.prev_lat)