Changeset 3670

Timestamp:

Mar 5, 2025, 2:13:18 PM (4 months ago)

Author:

afalco

Message:

Pluto: fixed bug in LMDZ to DYNAMICO conversion.
AF

Location:

trunk/LMDZ.PLUTO/util/startarchive2icosa

Files:

• icosa_finalize.bash (modified) (2 diffs)
• rearrange_startphy.f90 (modified) (2 diffs)

trunk/LMDZ.PLUTO/util/startarchive2icosa/icosa_finalize.bash

@@ -14 +14 @@
 
 # save start*nc from lmdz grid
-mkdir -p lmdz_starts
-cp start.nc startfi.nc lmdz_starts
 
 #conversion to netcdf 3 in order to rename and add variables, not properly done with netcdf4 files
@@ -68 +66 @@
 mv start_icosa.nc icosa_starts/start.nc
 mv startfi.nc icosa_starts/startfi.nc
-
-# get back start files from lmdz grid
-cp lmdz_starts/* .

trunk/LMDZ.PLUTO/util/startarchive2icosa/rearrange_startphy.f90

@@ -144 +144 @@
 !   ierr=NF90_INQUIRE_DIMENSION(ncid, dimids(9), len=nslope)
   ierr=NF90_INQUIRE(ncid,nVariables=nvar)
+  print*, "physical_points", physical_points, "dim", dimids(1)
+  print*, "nvertex", nvertex, "dim", dimids(2)
+  print*, "index", index_, "dim", dimids(4)
+  print*, "subsurface_layers", subsurface_layer, "dim", dimids(6)
+  print*, "lev_p1", lev_p1, "dim", dimids(7)
+  print*, "Time", Time, "dim", dimids(8)
 
 
@@ -235 +241 @@
 
     ! Processing 2D variables : soil
-    else if ((ndim==2).and.(varname_dimids(2)==subsurface_layer)) then
+    else if ((ndim==2).and.(varname_dimids(2)==dimids(6))) then
       write(*,*) "processing ",trim(varname)
       ! load field_soil