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

Timestamp:

Apr 25, 2025, 10:57:05 AM (4 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/phyredem.F90 (modified) (9 diffs)

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

@@ -2 +2 @@
 
 implicit none
+
+INTEGER,SAVE :: nid_restart ! NetCDF file identifier for restartfi.nc file
 
 contains
@@ -15 +17 @@
                        iceradius, dtemisice, phisfi
   use iostart, only : open_restartphy, close_restartphy, &
-                      put_var, put_field, length
+                      put_var, put_field, length, ldscrpt, ndscrpt
   use mod_grid_phy_lmdz, only : klon_glo
   use planete_mod, only: year_day, periastr, apoastr, peri_day, &
@@ -81 +83 @@
 
   ! Create physics start file
-  call open_restartphy(filename)
+  call open_restartphy(filename,nid_restart)
 
 ! tab_cntrl() contains run parameters
@@ -133 +135 @@
   tab_cntrl(35) = volcapa ! soil volumetric heat capacity
 
-  call put_var("controle","Control parameters",tab_cntrl)
+  call put_var(nid_restart,"controle","Control parameters",tab_cntrl)
 
   ! Write the controle array descriptor
@@ -140 +142 @@
 
   ! Write the mid-layer depths
-  call put_var("soildepth","Soil mid-layer depth",mlayer)
+  call put_var(nid_restart,"soildepth","Soil mid-layer depth",mlayer)
 
   ! Write longitudes
-  call put_field("longitude","Longitudes of physics grid",lonfi)
+  call put_field(nid_restart,"longitude","Longitudes of physics grid",lonfi)
 
   ! Write latitudes
-  call put_field("latitude","Latitudes of physics grid",latfi)
+  call put_field(nid_restart,"latitude","Latitudes of physics grid",latfi)
 
   ! Write mesh areas
-  call put_field("area","Mesh area",cell_area)
+  call put_field(nid_restart,"area","Mesh area",cell_area)
 
   ! Write surface geopotential
-  call put_field("phisfi","Geopotential at the surface",phisfi)
+  call put_field(nid_restart,"phisfi","Geopotential at the surface",phisfi)
 
   ! Write surface albedo
-  !call put_field("albedodat","Albedo of bare ground",alb)
+  !call put_field(nid_restart,"albedodat","Albedo of bare ground",alb)
 
   ! Subgrid topogaphy variables
-  call put_field("ZMEA","Relief: mean relief",zmea)
-  call put_field("ZSTD","Relief: standard deviation",zstd)
-  call put_field("ZSIG","Relief: sigma parameter",zsig)
-  call put_field("ZGAM","Relief: gamma parameter",zgam)
-  call put_field("ZTHE","Relief: theta parameter",zthe)
+  call put_field(nid_restart,"ZMEA","Relief: mean relief",zmea)
+  call put_field(nid_restart,"ZSTD","Relief: standard deviation",zstd)
+  call put_field(nid_restart,"ZSIG","Relief: sigma parameter",zsig)
+  call put_field(nid_restart,"ZGAM","Relief: gamma parameter",zgam)
+  call put_field(nid_restart,"ZTHE","Relief: theta parameter",zthe)
 
   ! Close file
-  call close_restartphy
+  call close_restartphy(nid_restart)
 
 end subroutine physdem0
@@ -197 +199 @@
 
   ! Open file
-  call open_restartphy(filename)
+  call open_restartphy(filename, nid_restart)
 
   ! First variable to write must be "Time", in order to correctly
@@ -204 +206 @@
 
   ! Surface temperature
-  call put_field("tsurf","Surface temperature",tsurf)
+  call put_field(nid_restart,"tsurf","Surface temperature",tsurf)
 
   ! Soil inertia
-  call put_field("inertiedat","Soil thermal inertia",inertiesoil)
+  call put_field(nid_restart,"inertiedat","Soil thermal inertia",inertiesoil)
 
   ! Soil temperature
-  call put_field("tsoil","Soil temperature",tsoil)
+  call put_field(nid_restart,"tsoil","Soil temperature",tsoil)
 
   ! Emissivity of the surface
-  call put_field("emis","Surface emissivity",emis)
+  call put_field(nid_restart,"emis","Surface emissivity",emis)
 
   ! Albedo of the surface
-  call put_field("albedodat","Albedo of bare ground",alb)
+  call put_field(nid_restart,"albedodat","Albedo of bare ground",alb)
 
   ! Planetary Boundary Layer
-  call put_field("q2","pbl wind variance",q2)
+  call put_field(nid_restart,"q2","pbl wind variance",q2)
 
 ! cloud fraction and sea ice !AF24: removed
@@ -228 +230 @@
   if (nq>0) then
     do iq=1,nq
-      call put_field(noms(iq),"tracer on surface",qsurf(:,iq))
+      call put_field(nid_restart,noms(iq),"tracer on surface",qsurf(:,iq))
     enddo
   endif ! of if (nq>0)
@@ -235 +237 @@
 
 ! close file
-      CALL close_restartphy
+      CALL close_restartphy(nid_restart)
 !$OMP BARRIER