Changeset 2952 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Apr 28, 2023, 2:28:08 PM (2 years ago)

Author:

romain.vande

Message:

Mars PCM + newstart + start_archive: Correct writing of variables inertiesoil and fluxgeo following r2919 and r2942
RV

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• comsoil_h.F90 (modified) (2 diffs)
• iostart.F90 (modified) (3 diffs)

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

@@ -84 +84 @@
     allocate(alph(ngrid,nsoilmx-1,nslope))
     allocate(beta(ngrid,nsoilmx-1,nslope))
+    allocate(flux_geo(ngrid,nslope))
  
   end subroutine ini_comsoil_h_slope_var
@@ -99 +100 @@
     if (allocated(alph)) deallocate(alph)
     if (allocated(beta)) deallocate(beta)
+    if (allocated(flux_geo)) deallocate(flux_geo)
 
   end subroutine end_comsoil_h_slope_var

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

@@ -836 +836 @@
             ierr=NF90_REDEF(nid_restart)
 #ifdef NC_DOUBLE
-            if(field_name.eq. "tsoil") then
+            if(field_name.eq. "tsoil".or. field_name.eq. "inertiesoil") then
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
                               (/idim2,idim3,idim8,idim7/),nvarid)
@@ -844 +844 @@
             endif
 #else
-            if(field_name.eq. "tsoil") then
+            if(field_name.eq. "tsoil".or. field_name.eq. "inertiesoil") then
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
                               (/idim2,idim3,idim8,idim7/),nvarid)
@@ -1015 +1015 @@
           IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
           ierr=NF90_ENDDEF(nid_restart)
-          
           firsttime=.false.
           if(first_time1d) first_time1d=.false.