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

Timestamp:

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

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

@@ -4 +4 @@
 
 real, save :: tab_cntrl_mod(100)
+
+integer,save :: nid_start ! NetCDF file identifier for startfi.nc file
 
 !$OMP THREADPRIVATE(tab_cntrl_mod)
@@ -109 +111 @@
 if (startphy_file) then
   ! open physics initial state file:
-  call open_startphy(fichnom)
+  call open_startphy(fichnom,nid_start)
 
   ! possibility to modify tab_cntrl in tabfi
@@ -124 +126 @@
 if (startphy_file) then
   ! Load surface geopotential:
-  call get_field("phisfi",phisfi,found)
+  call get_field(nid_start,"phisfi",phisfi,found)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <phisfi>",1)
@@ -136 +138 @@
 if (startphy_file) then
   ! Load bare ground albedo: (will be stored in surfdat_h)
-  call get_field("albedodat",albedodat,found)
+  call get_field(nid_start,"albedodat",albedodat,found)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <albedodat>",1)
@@ -152 +154 @@
 if (startphy_file) then
   ! Load surface albedo (for now assume it is spectrally homogeneous)
-  call get_field("albedo",albedo(:,1),found)
+  call get_field(nid_start,"albedo",albedo(:,1),found)
   if (.not.found) then
     write(*,*) modname//": Failed loading <albedo>"
@@ -173 +175 @@
 ! ZMEA
 if (startphy_file) then
-  call get_field("ZMEA",zmea,found)
+  call get_field(nid_start,"ZMEA",zmea,found)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <ZMEA>",1)
@@ -185 +187 @@
 ! ZSTD
 if (startphy_file) then
-  call get_field("ZSTD",zstd,found)
+  call get_field(nid_start,"ZSTD",zstd,found)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <ZSTD>",1)
@@ -197 +199 @@
 ! ZSIG
 if (startphy_file) then
-  call get_field("ZSIG",zsig,found)
+  call get_field(nid_start,"ZSIG",zsig,found)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <ZSIG>",1)
@@ -209 +211 @@
 ! ZGAM
 if (startphy_file) then
-  call get_field("ZGAM",zgam,found)
+  call get_field(nid_start,"ZGAM",zgam,found)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <ZGAM>",1)
@@ -221 +223 @@
 ! ZTHE
 if (startphy_file) then
-  call get_field("ZTHE",zthe,found)
+  call get_field(nid_start,"ZTHE",zthe,found)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <ZTHE>",1)
@@ -233 +235 @@
 ! Surface temperature :
 if (startphy_file) then
-  call get_field("tsurf",tsurf,found,indextime)
+  call get_field(nid_start,"tsurf",tsurf,found,indextime)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <tsurf>",1)
@@ -245 +247 @@
 ! Surface emissivity
 if (startphy_file) then
-  call get_field("emis",emis,found,indextime)
+  call get_field(nid_start,"emis",emis,found,indextime)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <emis>",1)
@@ -262 +264 @@
 if (CLFvarying) then
   if (startphy_file) then
-    call get_field("cloudfrac",cloudfrac,found,indextime)
+    call get_field(nid_start,"cloudfrac",cloudfrac,found,indextime)
     if (.not.found) then
       call abort_physic(modname,"Failed loading <cloudfrac>",1)
@@ -278 +280 @@
 if (CLFvarying) then
   if (startphy_file) then
-    call get_field("totcloudfrac",totcloudfrac,found,indextime)
+    call get_field(nid_start,"totcloudfrac",totcloudfrac,found,indextime)
     if (.not.found) then
       call abort_physic(modname,"Failed loading <totcloudfrac>",1)
@@ -293 +295 @@
 ! Height of oceanic ice (added by BC 2010)
 if (startphy_file) then
-  call get_field("hice",hice,found,indextime)
+  call get_field(nid_start,"hice",hice,found,indextime)
   if (.not.found) then
     write(*,*) "phyetat0: Failed loading <hice>"
@@ -308 +310 @@
 if (startphy_file) then
   ! nature of the surface
-  call get_field("rnat",rnat,found,indextime)
+  call get_field(nid_start,"rnat",rnat,found,indextime)
   if (.not.found) then
     write(*,*) "phyetat0: Failed loading <rnat>"
@@ -329 +331 @@
 if (startphy_file) then
   ! Pourcentage of sea ice cover
-  call get_field("pctsrf_sic",pctsrf_sic,found,indextime)
+  call get_field(nid_start,"pctsrf_sic",pctsrf_sic,found,indextime)
   if (.not.found) then
     write(*,*) "phyetat0: Failed loading <pctsrf_sic>"
@@ -342 +344 @@
 if (startphy_file) then
   ! Slab ocean temperature (2 layers)
-  call get_field("tslab",tslab,found,indextime)
+  call get_field(nid_start,"tslab",tslab,found,indextime)
   if (.not.found) then
     write(*,*) "phyetat0: Failed loading <tslab>"
@@ -359 +361 @@
 if (startphy_file) then
   ! Oceanic ice temperature
-  call get_field("tsea_ice",tsea_ice,found,indextime)
+  call get_field(nid_start,"tsea_ice",tsea_ice,found,indextime)
   if (.not.found) then
     write(*,*) "phyetat0: Failed loading <tsea_ice>"
@@ -372 +374 @@
 if (startphy_file) then
   ! Oceanic ice temperature
-  call get_field("tice",tice,found,indextime)
+  call get_field(nid_start,"tice",tice,found,indextime)
   if (.not.found) then
     write(*,*) "phyetat0: Failed loading <tice>"
@@ -385 +387 @@
 if (startphy_file) then
   !  Oceanic ice quantity (kg/m^2)
-  call get_field("sea_ice",sea_ice,found,indextime)
+  call get_field(nid_start,"sea_ice",sea_ice,found,indextime)
   if (.not.found) then
     write(*,*) "phyetat0: Failed loading <sea_ice>"
@@ -399 +401 @@
 ! pbl wind variance
 if (startphy_file) then
-  call get_field("q2",q2,found,indextime)
+  call get_field(nid_start,"q2",q2,found,indextime)
   if (.not.found) then
     call abort_physic(modname,"Failed loading <q2>",1)
@@ -414 +416 @@
     txt=noms(iq)
     if (startphy_file) then
-      call get_field(txt,qsurf(:,iq),found,indextime)
+      call get_field(nid_start,txt,qsurf(:,iq),found,indextime)
       if (.not.found) then
         write(*,*) "phyetat0: Failed loading <",trim(txt),">"
@@ -440 +442 @@
 ! Non-orographic gravity waves
 if (startphy_file) then
-   call get_field("du_nonoro_gwd",du_nonoro_gwd,found,indextime)
+   call get_field(nid_start,"du_nonoro_gwd",du_nonoro_gwd,found,indextime)
    if (.not.found) then
       write(*,*) "phyetat0: <du_nonoro_gwd> not in file"
@@ -451 +453 @@
 endif ! of if (startphy_file)
 if (startphy_file) then
-   call get_field("dv_nonoro_gwd",dv_nonoro_gwd,found,indextime)
+   call get_field(nid_start,"dv_nonoro_gwd",dv_nonoro_gwd,found,indextime)
    if (.not.found) then
       write(*,*) "phyetat0: <dv_nonoro_gwd> not in file"
@@ -462 +464 @@
 endif ! of if (startphy_file)
 if (startphy_file) then
-   call get_field("east_gwstress",east_gwstress,found,indextime)
+   call get_field(nid_start,"east_gwstress",east_gwstress,found,indextime)
    if (.not.found) then
       write(*,*) "phyetat0: <east_gwstress> not in file"
@@ -473 +475 @@
 endif ! of if (startphy_file)
 if (startphy_file) then
-   call get_field("west_gwstress",west_gwstress,found,indextime)
+   call get_field(nid_start,"west_gwstress",west_gwstress,found,indextime)
    if (.not.found) then
       write(*,*) "phyetat0: <west_gwstress> not in file"
@@ -486 +488 @@
 ! close file:
 !
-if (startphy_file) call close_startphy
+if (startphy_file) call close_startphy(nid_start)
 
 end subroutine phyetat0