Changeset 2913 for trunk/LMDZ.MARS/libf/phymars/iostart.F90

Timestamp:

Mar 14, 2023, 10:07:33 AM (2 years ago)

Author:

romain.vande

Message:

Mars PCM:
Adapt start2archive.F to the subslope parametrisation.
Small correction for some dimensions of variables.
RV

File:

• trunk/LMDZ.MARS/libf/phymars/iostart.F90 (modified) (6 diffs)

trunk/LMDZ.MARS/libf/phymars/iostart.F90

@@ -15 +15 @@
     INTEGER,SAVE :: idim7 ! "Time" dimension
     INTEGER,SAVE :: idim8 ! "nslope" dimension
-    INTEGER,SAVE :: idim9 ! "inter slope" dimension
+    INTEGER,SAVE :: idim9 ! "nslope_plus_1" dimension
     INTEGER,SAVE :: timeindex ! current time index (for time-dependent fields)
     INTEGER,PARAMETER :: length=100 ! size of tab_cntrl array
@@ -564 +564 @@
       ENDIF
 
-      ierr=NF90_DEF_DIM(nid_restart,"inter slope",nslope+1,idim9)
+      ierr=NF90_DEF_DIM(nid_restart,"inter_slope",nslope+1,idim9)
       IF (ierr/=NF90_NOERR) THEN
         write(*,*)'phyredem: problem defining inter slope dimension'
@@ -660 +660 @@
   
   REAL                           :: field_glo(klon_glo,field_size)
+  REAL                           :: field_glo_reshape(klon_glo,nsoilmx,nslope,timeindex)
   INTEGER                        :: ierr
   INTEGER                        :: nvarid
@@ -804 +805 @@
         endif ! of if (.not.present(timeindex))
 
-      ELSE IF (field_size==nsoilmx) THEN
+      ELSE IF (field_size==nsoilmx .or. field_size==nsoilmx*nslope) THEN
         ! input is a 2D "subsurface field" array
         if (.not.present(time)) then ! for a time-independent field
@@ -828 +829 @@
             ierr=NF90_REDEF(nid_restart)
 #ifdef NC_DOUBLE
+            if(field_name.eq. "tsoil") then
+            ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
+                              (/idim2,idim3,idim8,idim7/),nvarid)
+            else
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
                               (/idim2,idim3,idim7/),nvarid)
-#else
+            endif
+#else
+            if(field_name.eq. "tsoil") then
+            ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
+                              (/idim2,idim3,idim8,idim7/),nvarid)
+            else
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
                               (/idim2,idim3,idim7/),nvarid)
+            endif
 #endif
            if (ierr.ne.NF90_NOERR) then
@@ -842 +853 @@
           endif
           ! Write the variable
+          if(field_name.eq. "tsoil") then
+
+          field_glo_reshape=RESHAPE(field_glo, (/klon_glo, nsoilmx, nslope, timeindex/))
+
+          ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo_reshape,&
+                            start=(/1,1,1,timeindex/))
+          else
           ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo,&
                             start=(/1,1,timeindex/))
-
+          endif
         endif ! of if (.not.present(time))