Changeset 2813 for trunk/LMDZ.MARS/util/startarchive2icosa/icosa_finalize.bash

Timestamp:

Nov 8, 2022, 3:32:12 PM (2 years ago)

Author:

romain.vande

Message:

DYNAMICO MARS
Update startarchive2icosa repertory for new version of Dynamico
RV

File:

• trunk/LMDZ.MARS/util/startarchive2icosa/icosa_finalize.bash (modified) (3 diffs)

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

@@ -7 +7 @@
 # load nco module
 module purge
-source /scratch/cnt0027/lmd1167/adelavois/models/code/XIOS/arch.env
-module load nco/4.7.9-gcc-4.8.5-hdf5-1.8.18-openmpi-2.0.4
+source /scratch/cnt0027/lmd1167/adelavois/models/code/XIOS/arch.env #You have to adapt to your own path
+#module load nco/4.7.9-gcc-4.8.5-hdf5-1.8.18-openmpi-2.0.4 #Example of version of the module needed
+module load nco
 
 #conversion to netcdf 3 in order to rename and add variables, not properly done with netcdf4 files
@@ -24 +25 @@
 #Temporay soildepth array from Dynamico restartfi
 ncks -A -v soildepth \
-         ../../mars_dynamico_update_basic/startfi_0.nc startphy_nc3.nc
+         ../../mars_dynamico_update_basic/startfi_0.nc startphy_nc3.nc #You have to adapt to your own path
 
-nccopy -k 'netCDF-4' startphy_nc3.nc startfi.nc
+nccopy -k 4 startphy_nc3.nc startfi.nc
 
 #Temporary Time variable from start2archive; should be added with xios
-ncks -A -v Time \
-         start_archive_nc4.nc startfi.nc
+#ncks -A -v Time \
+#         start_archive_nc4.nc startfi.nc
+ncks -O -x -v time_counter startfi.nc startfi.nc_tmp
+ncks -A -v Time start.nc startfi.nc_tmp
+cp startfi.nc_tmp startfi.nc
 
 rearrange_startphy > rearrange_startphy.out 2>&1
@@ -47 +51 @@
      start_icosa_ref.nc start_icosa.nc
 
+# In the current version of Dynamico, the variable q for tracer is needed but unused.
+# To create this unused variable simply run this short python code
+
+#import numpy as np
+#import netCDF4 as nc
+
+#fn = 'start_icosa.nc'
+#ds = nc.Dataset(fn, 'a', format='NETCDF4')
+#nq = ds.createDimension('nq', 7)
+
+#q = ds.createVariable('q', 'f8', ('nq','lev', 'cell_mesh',))
+#nqq = ds.createVariable('nq', 'f8', ('nq',))
+#q[:,:,:]=0.
+#q.online_operation = "once"
+#q.coordinates = "lat lon"
+#ds.close()
+
 # cleanup
 rm -rf startphy_icosa_nc3.nc startphy_nc3.nc