Changeset 2832 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Timestamp:

Nov 25, 2022, 4:05:58 PM (2 years ago)

Author:

emillour

Message:

Mars GCM:

• follow-up to r2826-2829: adapt newstart and lect_start_archive to handle "old" start_archive files where surface co2 ice is stored as "co2ice".

EM

Location:

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Files:

• lect_start_archive.F (modified) (5 diffs)
• newstart.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F

@@ -25 +25 @@
       implicit none
 
-#include "dimensions.h"
-#include "paramet.h"
-#include "comgeom2.h"
-#include "netcdf.inc"
+      include "dimensions.h"
+      include "paramet.h"
+      include "comgeom2.h"
+      include "netcdf.inc"
 c=======================================================================
 c   Declarations
@@ -196 +196 @@
       real :: tmpval ! to store a temporary variable/value
 
+      ! flag to check if CO2 surface ice is in "co2ice" (old start_archive)
+      ! then it is set  to .true. or else it is in "co2_surf" (newer 
+      ! start_archive) and then .false.
+      logical :: old_co2ice=.false. 
 c=======================================================================
 
@@ -702 +706 @@
       count=(/imold+1,jmold+1,1,0/)
        
-      ierr = NF_INQ_VARID (nid, "co2", nvarid)
-      IF (ierr .NE. NF_NOERR) THEN
-         PRINT*, "lect_start_archive: <co2> is missing"
-         CALL abort
+      ! look for CO2ice on the surface
+      ! first try the "old" co2ice field for retro-compatibility:
+      ierr = NF_INQ_VARID (nid, "co2ice", nvarid)
+      IF (ierr .EQ. NF_NOERR) THEN
+        WRITE(*,*)" Found co2ice => this is an 'old' start_archive"
+        WRITE(*,*)" will read in 'co2ice' instead of 'co2_surf'"
+        old_co2ice=.true.
+      ELSE
+        ! no 'co2ice', look for co2_surf instead
+        ierr = NF_INQ_VARID (nid, "co2ice", nvarid)
+        IF (ierr .NE. NF_NOERR) THEN
+          PRINT*, "lect_start_archive: <co2_surf> is missing"
+          PRINT*, NF_STRERROR(ierr)
+          CALL abort
+        ENDIF
       ENDIF
 #ifdef NC_DOUBLE
@@ -713 +728 @@
 #endif
       IF (ierr .NE. NF_NOERR) THEN
-         PRINT*, "lect_start_archive: Failed loading <co2>"
-         PRINT*, NF_STRERROR(ierr)
-         CALL abort
+        PRINT*, "lect_start_archive: Failed loading <co2ice>"
+        PRINT*, NF_STRERROR(ierr)
+        CALL abort
       ENDIF
 c
@@ -878 +893 @@
           if(txt.eq."co2_surf") then
             igcm_co2=iq
+            if (old_co2ice) then
+              ! CO2 surface ice has already been loaded from "co2ice"
+              cycle
+            endif
           endif
         ENDIF ! of IF (oldtracernames)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -1786 +1786 @@
      &              nsoilmx,ngridmx,llm,
      &              nqtot,dtphys,real(day_ini),0.0,cell_area,
-     &              albfi,ithfi,zmea,zstd,zsig,zgam,zthe,
-     &              hmons,summit,base)
+     &              albfi,ithfi)
       call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
      &              dtphys,hour_ini,