Changeset 3761

Timestamp:

May 14, 2025, 5:14:10 PM (4 weeks ago)

Author:

afalco

Message:

Pluto: start2archive: do not write to file if existing start_archive.nc has different dimensions than start/startfi.nc/executable.
AF

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/start2archive.F (modified) (2 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/start2archive.F

@@ -106 +106 @@
       REAL vn(ip1jm,llm),un(ip1jmp1,llm)
 
+c Variables from pre-existing start_archive
+c -----------------
+      integer :: nsoilold, nlatold, nlonold, naltold
 c Autres  variables
 c -----------------
@@ -439 +442 @@
       endif
 
+      ierr=NF_INQ_DIMID(nid,"altitude",dimid)
+      if (ierr.eq.NF_NOERR) then
+        ierr=NF_INQ_DIMLEN(nid,dimid,naltold)
+            if (ierr.ne.NF_NOERR) then
+            call abort_gcm('start2archive: ',
+     &              'Failed reading altitude length', 1)
+        endif
+        if (naltold.ne.llm) then
+          print*,"start.nc:          n_levels:", llm
+          print*, "start_archive.nc: n_levels:",naltold
+          call abort_gcm('start2archive: ',
+     &              'altitude dimension from start_archive.nc differ'//
+     &              'from start/startfi.nc', 1)
+        endif
+      endif
+
+      ierr=NF_INQ_DIMID(nid,"longitude",dimid)
+      if (ierr.eq.NF_NOERR) then
+        ierr=NF_INQ_DIMLEN(nid,dimid,nlonold)
+            if (ierr.ne.NF_NOERR) then
+            call abort_gcm('start2archive: ',
+     &              'Failed reading longitude length', 1)
+        endif
+        if (nlonold.ne.iim) then
+          print*,"start.nc:          n_lon:", iim
+          print*, "start_archive.nc: n_lon:",nlonold
+          call abort_gcm('start2archive: ',
+     &              'longitude dimension from start_archive.nc differ'//
+     &              'from start/startfi.nc', 1)
+        endif
+      endif
+
+
+      ierr=NF_INQ_DIMID(nid,"latitude",dimid)
+      if (ierr.eq.NF_NOERR) then
+        ierr=NF_INQ_DIMLEN(nid,dimid,nlatold)
+            if (ierr.ne.NF_NOERR) then
+            call abort_gcm('start2archive: ',
+     &              'Failed reading latitude length', 1)
+        endif
+        if (nlatold.ne.jjm) then
+          print*,"start.nc:          n_lat:", jjm
+          print*, "start_archive.nc: n_lat:",nlatold
+          call abort_gcm('start2archive: ',
+     &              'latitude dimension from start_archive.nc differ'//
+     &              'from start/startfi.nc', 1)
+        endif
+      endif
+
+      ierr=NF_INQ_DIMID(nid,"subsurface_layers",dimid)
+      if (ierr.eq.NF_NOERR) then
+        ierr=NF_INQ_DIMLEN(nid,dimid,nsoilold)
+            if (ierr.ne.NF_NOERR) then
+            call abort_gcm('start2archive: ',
+     &              'Failed reading subsurface_layers length', 1)
+        endif
+        if (nsoilold.ne.nsoilmx) then
+          print*,"start.nc:          n_levels:", nsoilmx
+          print*, "start_archive.nc: n_levels:",nsoilold
+          call abort_gcm('start2archive: ',
+     &              'Soil dimension from start_archive.nc differ'//
+     &              'from start/startfi.nc', 1)
+        endif
+      endif
+
 c-----------------------------------------------------------------------
 c Ecriture des champs  (n2ice,emis,ps,Tsurf,T,u,v,q2,q,qsurf)