Changeset 2540

Timestamp:

Jun 29, 2021, 10:02:59 PM (4 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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (3 diffs)
• util/simu_MCS.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -3367 +3367 @@
 find cv.
 
-
 == 6/05/2021 == JN
 vdifc_mod update !!
@@ -3385 +3384 @@
 minor fix in co2condens4micro; fix co2 sedimentation rate outputs
 
+== 26/05/2021 == AB
+Make a cleaner initialization of some variables in simu_MCS.F90
+
 == 26/05/2021 == JN
 Bug correction in testphys1d.F, when calling "phys_start_var_init"
@@ -3403 +3405 @@
 == 22/06/2021 == CM
 co2clouds: fix mtotco2, icetotco2, vaptotco2 outputs
+
+== 29/06/2021 == AB
+Make simu_MCS.F90 compatible with DYNAMICO lon-lat output files, with a check on the latitude array order when doing the interpolation

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)