Changeset 4168 for LMDZ6/trunk

Timestamp:

Jun 9, 2022, 5:29:22 PM (2 years ago)

Author:

evignon

Message:

deuxieme phase de travail pour l'inclusion de make_sso
Ehouarn, David, Romain , Etienne

Location:

LMDZ6/trunk/tools/make_sso

Files:

• README (modified) (2 diffs)
• compile (modified) (1 diff)
• make_sso (added)
• make_sso.env (added)
• make_sso_SpherePack.f90 (modified) (1 diff)

LMDZ6/trunk/tools/make_sso/README

@@ -9 +9 @@
 - ifortran
 - netcdf
+- an adequate environment defined in make_sso.env
 
 Options:
--i input_relief_file
--res nlon nlat : of the output grid  
--v : orography variable name of input file
--m : option
+
+ ./make_sso  -i <file>  -v <var>  -res <nlon> <nlat>  -m <mask>  -f <ffac>  -t
+  <tfac>
+ 
+ Where:
+ 
+  <file>:   input file name
+  <var> :   height variable name
+  <nlon>:   number of distinct longitude points of output grid (iim)
+  <nlat>:   number of distinct latitude  points of output grid (jjm+1)
+  <mask>:   can be:  1. a mask file (o2a.nc)
+                     2. "noro": computed with grid_noro
+            default: 3. "spec": computed internally
+  <ffac>:   filtering  width: ffac*Do         (Do:output resol.)
+  <tfac>:   transition width: Di+tfac*(Do-Di) (Di: input resol.)
+ 
+ Note that latitudes grid contains both poles.
+
 
 Example to compile then run make_sso from a high-resolution Relief.nc file to
@@ -21 +36 @@
 ./compile
 file=Relief.nc
-./make_sso.e -i $file -res 180 91 -v z
+./make_sso -i $file -res 180 91 -v z
 
-Output file is ${file}_180x91
+Output file is make_sso_180x91_${file}
 

LMDZ6/trunk/tools/make_sso/compile

@@ -1 @@
-case ${HOSTNAME:0:5} in
-  irene)
-    module unload netcdf-fortran
-    module unload netcdf-c 
-    module load netcdf-c
-    module load netcdf-fortran
-    module load fftw3/gnu
-    netcdf="-L$NETCDF_LIBDIR"
-  ;;
-  cicla)
-    module load intel
-    module load netcdf4/4.3.3.1-ifort
-    NETCDF=$(which ncdump); NETCDF=${NETCDF%%/bin*}
-    NETCDFFORTRAN_INCDIR=$NETCDF/include
-    NETCDFFORTRAN_LIBDIR=$NETCDF/lib
-    netcdf=""
-  ;;
-esac
+source make_sso.env
+
 
 if [ ! -f Relief.nc ]; then
 echo Attention missing Relief.nc file
 echo I download it
-wget https://web.lmd.jussieu.fr/~lmdz/pub/3DInputData/Limit/Relief.nc
+wget --no-check-certificate -nv https://web.lmd.jussieu.fr/~lmdz/pub/3DInputData/Limit/Relief.nc
 fi
 
 if [ ! -d SPHEREPACK ]; then
 echo Missing SPHEREPACK library. I download it
-wget https://web.lmd.jussieu.fr/~lmdz/pub/import/spherepack3.2.tgz 
+wget --no-check-certificate -nv https://web.lmd.jussieu.fr/~lmdz/pub/import/spherepack3.2.tgz 
 tar -xvf spherepack3.2.tgz
+mv spherepack3.2 SPHEREPACK
 mv SPHEREPACK/make.inc SPHEREPACK/make.inc.old
+
+# copy more adequate make.inc to SPHEREPACK
 mv make.inc SPHEREPACK/.
+# copy corrected Makefile to SPHEREPACK/test
+mv SPHEREPACK/test/Makefile SPHEREPACK/test/Makefile.old
+mv  Makefile.test SPHEREPACK/test/Makefile 
+# clean up
 rm spherepack3.2.tgz
 fi

LMDZ6/trunk/tools/make_sso/make_sso_SpherePack.f90

@@ -453 +453 @@
   res_in=TRIM(int2str(nlon_in))//'x'//TRIM(int2str(nlat_in))
   res_ou=TRIM(int2str(nlon_ou))//'x'//TRIM(int2str(nlat_ou))
-  f_ou=TRIM(f_in)//'_'//TRIM(res_ou)
+  f_ou='make_sso_'//TRIM(res_ou)//'_'//TRIM(f_in)
   fnam=f_ou
   CALL nc(NF90_CREATE(f_ou,NF90_CLOBBER,fID))