Changeset 3552 for trunk/LMDZ.GENERIC/libf/phystd/phyredem.F90

Timestamp:

Dec 13, 2024, 5:17:13 PM (10 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/phyredem.F90 (modified) (7 diffs)

trunk/LMDZ.GENERIC/libf/phystd/phyredem.F90

@@ -2 +2 @@
 
 implicit none
+
+INTEGER,SAVE :: nid_restart ! NetCDF file identifier for restartfi.nc file
 
 contains
@@ -82 +84 @@
   
   ! Create physics start file
-  call create_restartphy(filename)
+  call create_restartphy(filename,nid_restart)
 
 ! tab_cntrl() contains run parameters
@@ -135 +137 @@
 
   ! Write the controle array
-  call put_var("controle","Control parameters",tab_cntrl)
+  call put_var(nid_restart,"controle","Control parameters",tab_cntrl)
   
   ! Write the controle array descriptor
-  call put_var("controle_descriptor",&
+  call put_var(nid_restart,"controle_descriptor",&
                "Description of control parameters",dscrpt_tab_cntrl)
   
   ! 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)
   
   ! Soil thermal inertia
-  call put_field("inertiedat","Soil thermal inertia",ith)
+  call put_field(nid_restart,"inertiedat","Soil thermal inertia",ith)
   
   ! Close file
-  call close_restartphy
+  call close_restartphy(nid_restart)
   
 end subroutine physdem0
@@ -217 +219 @@
   
   ! Open file
-  call open_restartphy(filename)
+  call open_restartphy(filename,nid_restart)
 
   ! First variable to write must be "Time", in order to correctly
@@ -224 +226 @@
   
   ! Surface temperature
-  call put_field("tsurf","Surface temperature",tsurf)
+  call put_field(nid_restart,"tsurf","Surface temperature",tsurf)
   
   ! 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)
   
   ! Surface albedo (assume homegeneous spectral albedo for now)
-  call put_field("albedo","Surface albedo",albedo(:,1))
+  call put_field(nid_restart,"albedo","Surface albedo",albedo(:,1))
   
   ! 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 (NB: these should be optional... to be improved)
-  call put_field("cloudfrac","Cloud fraction",cloudfrac)
-  call put_field("totcloudfrac","Total fraction",totcloudfrac)
-  call put_field("hice","Height of oceanic ice",hice)
+  call put_field(nid_restart,"cloudfrac","Cloud fraction",cloudfrac)
+  call put_field(nid_restart,"totcloudfrac","Total fraction",totcloudfrac)
+  call put_field(nid_restart,"hice","Height of oceanic ice",hice)
 
  !Slab ocean 
  if(ok_slab_ocean) then
-   call put_field("rnat","Nature of surface",rnat)
-   call put_field("pctsrf_sic","Pourcentage sea ice",pctsrf_sic)
-   call put_field("tslab","Temperature slab ocean",tslab)
-   call put_field("tsea_ice","Temperature of top layer (seaice or snow)",tsea_ice)
-   call put_field("tice","Temperature of seaice",tice)
-   call put_field("sea_ice","Sea ice mass",sea_ice)
+   call put_field(nid_restart,"rnat","Nature of surface",rnat)
+   call put_field(nid_restart,"pctsrf_sic","Pourcentage sea ice",pctsrf_sic)
+   call put_field(nid_restart,"tslab","Temperature slab ocean",tslab)
+   call put_field(nid_restart,"tsea_ice","Temperature of top layer (seaice or snow)",tsea_ice)
+   call put_field(nid_restart,"tice","Temperature of seaice",tice)
+   call put_field(nid_restart,"sea_ice","Sea ice mass",sea_ice)
  endif!(ok_slab_ocean)
 
@@ -256 +258 @@
   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)
@@ -262 +264 @@
 ! Non-orographic gavity waves
 if (calllott_nonoro) then
-   call put_field("du_nonoro_gwd","Zonal wind tendency due to GW",du_nonoro_gwd)
-   call put_field("dv_nonoro_gwd","Meridional wind tendency due to GW",dv_nonoro_gwd)
-   call put_field("east_gwstress","Eastward stress profile due to GW",east_gwstress)
-   call put_field("west_gwstress","Westward stress profile due to GW",west_gwstress)
+   call put_field(nid_restart,"du_nonoro_gwd","Zonal wind tendency due to GW",du_nonoro_gwd)
+   call put_field(nid_restart,"dv_nonoro_gwd","Meridional wind tendency due to GW",dv_nonoro_gwd)
+   call put_field(nid_restart,"east_gwstress","Eastward stress profile due to GW",east_gwstress)
+   call put_field(nid_restart,"west_gwstress","Westward stress profile due to GW",west_gwstress)
 endif
 
 ! close file
-      CALL close_restartphy
+      CALL close_restartphy(nid_restart)
 !$OMP BARRIER