Changeset 2432 for trunk/LMDZ.MARS/util/zrecast.F90

Timestamp:

Nov 14, 2020, 3:27:44 PM (4 years ago)

Author:

emillour

Message:

Mars GCM:

• Make a single README file (get rid of REAME.exec) and update "compile" script to (hopefully) better illustrate how to compile utilities
• Update utilities zrecast and lslin: zrecast.F90: Force the vertical interpolation of any variable starting by "rho" to be in log space (as it should for density, molecular concentration in mol.cm-3 etc..). Set missing value to 1.E20 lslin.F90: fix various problems in the output.

FF + EM

File:

• trunk/LMDZ.MARS/util/zrecast.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/util/zrecast.F90

@@ -61 +61 @@
 ! TN 01/2013 : Adapted for large output files with at least 2 variables > 2 GiB
 ! TN 10/2013 : Read and write controle field, if available
+! FF 10/2020 : Force interpolation in log space for variable starting by "rho"
 !
 implicit none
@@ -94 +95 @@
 integer,dimension(3) :: surfdatashape ! shape of 3D (surface+time) datasets
 
-real :: miss_val=-9.99e+33 ! special "missing value" to specify missing data
-real,parameter :: miss_val_def=-9.99e+33 ! default value for "missing value"
+real :: miss_val= 1.E+20 ! special "missing value" to specify missing data
+real,parameter :: miss_val_def= 1.E+20 ! default value for "missing value"
 character (len=64), dimension(:), allocatable :: var
 ! var(): names of variables that will be processed
@@ -1809 +1810 @@
   ! interpolate dataset onto new grid
   ! check if variable is "rho" (to set flag for interpolation below)
-  if (var(i).eq.'rho') then
+  if (var(i)(1:3).eq.'rho') then
     j=1
+    write(*,*) trim(var(i)), ' is interpolated in log scale' 
   else
     j=0
@@ -1864 +1866 @@
     ! interpolate dataset onto new grid
     ! check if variable is "rho" (to set flag for interpolation below)
-    if (var(i).eq.'rho') then
+    if (var(i)(1:3).eq.'rho') then
       j=1
+      write(*,*) trim(var(i)), ' is interpolated in log scale' 
     else
       j=0
@@ -1919 +1922 @@
   ! interpolate dataset onto new grid
   ! check if variable is "rho" (to set flag for interpolation below)
-  if (var(i).eq.'rho') then
+  if (var(i)(1:3).eq.'rho') then
     j=1
+    write(*,*) trim(var(i)), ' is interpolated in log scale' 
   else
     j=0