Changeset 3736 for trunk/LMDZ.PLUTO/libf/phypluto/soil_settings.F90

Timestamp:

Apr 25, 2025, 10:57:05 AM (3 months ago)

Author:

afalco

Message:

Pluto: added missing controle descriptor parameters.
Update for restarts from generic (restarts in 1D)
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/soil_settings.F90 (modified) (10 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/soil_settings.F90

@@ -82 +82 @@
 ! -------------------
 ! 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
@@ -98 +98 @@
 
         ! 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)
@@ -118 +118 @@
 
 ! 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
@@ -139 +139 @@
         ! 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
             if (is_master) 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)
           if (is_master) print*,"mlayer",mlayer
           if (.not.found) then
@@ -206 +206 @@
        ! Read Surface thermal inertia
        allocate(surfinertia(ngrid))
-       call get_field("inertiedat",surfinertia,found)
+       call get_field(nid,"inertiedat",surfinertia,found)
        if (.not.found) then
          if (is_master) write(*,*) "soil_settings: Failed loading <inertiedat>"
@@ -229 +229 @@
            endif
          endif ! of if (.not.allocated(oldinertiedat))
-        call get_field("inertiedat",oldinertiedat,found)
+        call get_field(nid,"inertiedat",oldinertiedat,found)
         if (.not.found) then
           if (is_master) write(*,*) "soil_settings: Failed loading <inertiedat>"
@@ -235 +235 @@
         endif
        else ! put values in therm_i
-         call get_field("inertiedat",inertiedat,found)
+         call get_field(nid,"inertiedat",inertiedat,found)
          if (.not.found) then
            if (is_master) write(*,*) "soil_settings: Failed loading <inertiedat>"
@@ -247 +247 @@
 ! -------------------------
 !      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
@@ -265 +265 @@
            endif
          endif
-         call get_field("tsoil",oldtsoil,found)
+         call get_field(nid,"tsoil",oldtsoil,found)
          if (.not.found) then
            if (is_master) write(*,*) "soil_settings: Failed loading <tsoil>"
@@ -271 +271 @@
          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
            if (is_master) write(*,*) "soil_settings: Failed loading <tsoil>"