Changeset 3552 for trunk/LMDZ.GENERIC/libf/phystd/soil_settings.F

Timestamp:

Dec 13, 2024, 5:17:13 PM (9 days ago)

Author:

emillour

Message:

Generic PCM

Modify iostart routines to take netcdf index as an argument in order to
be used to read/write in other files than startfi.nc (preparing 1D
restart)

MM

File:

• trunk/LMDZ.GENERIC/libf/phystd/soil_settings.F (modified) (10 diffs)

trunk/LMDZ.GENERIC/libf/phystd/soil_settings.F

@@ -81 +81 @@
 ! -------------------
 ! 1.1 Start by reading how many layers of soil there are
-        dimlen=inquire_dimension_length("subsurface_layers")
+        dimlen=inquire_dimension_length(nid,"subsurface_layers")
 
         if (dimlen.ne.nsoil) then
@@ -97 +97 @@
 
         ! check if olmlayer distribution matches current one
-        call get_var("soildepth",oldmlayer,found)
+        call get_var(nid,"soildepth",oldmlayer,found)
         if (found) then
           malpha=oldmlayer(2)/oldmlayer(1)
@@ -117 +117 @@
         
 ! 1.2 Find out the # of dimensions <inertiedat> was defined as using
-      ndims=inquire_field_ndims("inertiedat")
+      ndims=inquire_field_ndims(nid,"inertiedat")
 ! 1.3 Read depths values or set olddepthdef flag and values
       if (ndims.eq.1) then ! we know that there is none
@@ -138 +138 @@
         ! read <depth> coordinate
         if (interpol) then !put values in oldmlayer
-          call get_var("soildepth",oldmlayer,found)
+          call get_var(nid,"soildepth",oldmlayer,found)
           if (.not.found) then
             write(*,*)'soil_settings: Problem while reading <soildepth>'
           endif
         else ! put values in mlayer
-          call get_var("soildepth",mlayer,found)
+          call get_var(nid,"soildepth",mlayer,found)
           print*,"mlayer",mlayer
           if (.not.found) then
@@ -205 +205 @@
        ! Read Surface thermal inertia
        allocate(surfinertia(ngrid))
-       call get_field("inertiedat",surfinertia,found)
+       call get_field(nid,"inertiedat",surfinertia,found)
        if (.not.found) then
          write(*,*) "soil_settings: Failed loading <inertiedat>"
@@ -228 +228 @@
            endif
          endif ! of if (.not.allocated(oldinertiedat))
-        call get_field("inertiedat",oldinertiedat,found)
+        call get_field(nid,"inertiedat",oldinertiedat,found)
         if (.not.found) then
           write(*,*) "soil_settings: Failed loading <inertiedat>"
@@ -234 +234 @@
         endif
        else ! put values in therm_i
-         call get_field("inertiedat",inertiedat,found)
+         call get_field(nid,"inertiedat",inertiedat,found)
          if (.not.found) then
            write(*,*) "soil_settings: Failed loading <inertiedat>"
@@ -246 +246 @@
 ! -------------------------
 !      ierr=nf90_inq_varid(nid,"tsoil",nvarid)
-      ok=inquire_field("tsoil")
+      ok=inquire_field(nid,"tsoil")
 !      if (ierr.ne.nf90_noerr) then
       if (.not.ok) then
@@ -264 +264 @@
            endif
          endif
-         call get_field("tsoil",oldtsoil,found)
+         call get_field(nid,"tsoil",oldtsoil,found)
          if (.not.found) then
            write(*,*) "soil_settings: Failed loading <tsoil>"
@@ -270 +270 @@
          endif
        else ! put values in tsoil
-         call get_field("tsoil",tsoil,found,timeindex=indextime)
+         call get_field(nid,"tsoil",tsoil,found,
+     &                  timeindex=indextime)
          if (.not.found) then
            write(*,*) "soil_settings: Failed loading <tsoil>"